rk_tableau_butcher_base.h

#ifndef __RK_TABLEAU_BUTCHER_BASE__
#define __RK_TABLEAU_BUTCHER_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>

#include "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 RKTableauButcherBase: public RKTableauBase<dim,real,MeshType>
{
public:
    RKTableauButcherBase(const int n_rk_stages, const std::string rk_method_string_input); 

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

protected:

    dealii::Table<2,double> butcher_tableau_a;

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

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


};

} // ODE namespace
} // PHiLiP namespace

#endif