empty_RRK_base.h

#ifndef __EMPTY_RRK_BASE_H__
#define __EMPTY_RRK_BASE_H__

#include "dg/dg_base.hpp"
#include "ode_solver/runge_kutta_methods/rk_tableau_base.h"

namespace PHiLiP {
namespace ODE {

#if PHILIP_DIM==1
template <int dim, typename real, typename MeshType = dealii::Triangulation<dim>>
#else
template <int dim, typename real, typename MeshType = dealii::parallel::distributed::Triangulation<dim>>
#endif
class EmptyRRKBase
{
public:
    explicit EmptyRRKBase(std::shared_ptr<RKTableauBase<dim,real,MeshType>> /*rk_tableau_input*/);


    virtual real compute_FR_entropy_contribution(const real /*dt*/,
            std::shared_ptr<DGBase<dim,real,MeshType>>/* dg*/,
            const std::vector<dealii::LinearAlgebra::distributed::Vector<double>> &/*rk_stage*/,
            const bool /*compute_K_norm*/) const{
        return 0.0;
    };


    virtual void store_stage_solutions(const int /*istage*/,
            const dealii::LinearAlgebra::distributed::Vector<double> /*rk_stage_i*/) {
        // Do not store anything
    };


    virtual real update_relaxation_parameter(const real /*dt*/, 
            std::shared_ptr<DGBase<dim,real,MeshType>> /*dg*/,
            const std::vector<dealii::LinearAlgebra::distributed::Vector<double>> &/*rk_stage*/,
            const dealii::LinearAlgebra::distributed::Vector<double> &/*solution_update*/
            ) {
        // Return 1 such that the time step isn't modified
        return 1.0;
    };
    

};

} // ODE namespace
} // PHiLiP namespace

#endif