rrk_ode_solver_base.h

#ifndef __RRK_ODE_SOLVER_BASE_H__
#define __RRK_ODE_SOLVER_BASE_H__

#include "dg/dg_base.hpp"
#include "ode_solver/runge_kutta_ode_solver.h"
#include "ode_solver/relaxation_runge_kutta/runge_kutta_store_entropy.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 RRKODESolverBase: public RKNumEntropy<dim, real, MeshType>
{
public:
    explicit RRKODESolverBase(
            std::shared_ptr<RKTableauButcherBase<dim,real,MeshType>> rk_tableau_input);


    real relaxation_parameter = 1.0;

protected:

    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) override;
    

    virtual real compute_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*/
            ) = 0;

};

} // ODE namespace
} // PHiLiP namespace

#endif