algebraic_rrk_ode_solver.h

#ifndef __ENERGY_RRK_ODESOLVER_H__
#define __ENERGY_RRK_ODESOLVER_H__

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

protected:

    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*/
            ) override;
    
    real compute_inner_product(
            const dealii::LinearAlgebra::distributed::Vector<double> &stage_i,
            const dealii::LinearAlgebra::distributed::Vector<double> &stage_j,
            std::shared_ptr<DGBase<dim,real,MeshType>> dg
            ) const;

};

} // ODE namespace
} // PHiLiP namespace

#endif