rk_tableau_base.h

#ifndef __RK_TABLEAU_BASE__
#define __RK_TABLEAU_BASE__

#include <deal.II/base/conditional_ostream.h>

#include <deal.II/grid/tria.h>
#include <deal.II/distributed/shared_tria.h>
#include <deal.II/distributed/tria.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 RKTableauBase
{
public:
    RKTableauBase(const int n_rk_stages, const std::string rk_method_string_input); 

    virtual ~RKTableauBase() = default;

    double get_a(const int i, const int j) const;

    double get_b(const int i) const;

    double get_c(const int i) const;

    void set_tableau();

    const int n_rk_stages;
    
protected:

    dealii::ConditionalOStream pcout; 
    
    const std::string rk_method_string;

    dealii::Table<2,double> butcher_tableau_a;

    dealii::Table<1,double> butcher_tableau_b;
    
    dealii::Table<1,double> butcher_tableau_c;
    
    virtual void set_a() = 0;

    virtual void set_b() = 0;
    
    virtual void set_c() = 0;


};

} // ODE namespace
} // PHiLiP namespace

#endif