low_storage_rk_tableau_base.h

#ifndef __LOW_STORAGE_RK_TABLEAU_BASE__
#define __LOW_STORAGE_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>
#include <deal.II/lac/full_matrix.h>
#include <deal.II/lac/vector.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 LowStorageRKTableauBase : public RKTableauBase<dim,real,MeshType>
{
public:
    LowStorageRKTableauBase(const int n_rk_stages, const int num_delta, const std::string rk_method_string_input); 

    virtual ~LowStorageRKTableauBase() = default;

public:
    double get_gamma(const int i, const int j) const;

    double get_beta(const int i) const;

    double get_delta(const int i) const;

    double get_b_hat(const int i) const;

    void set_tableau() override;

protected:
    
    const int num_delta;

    dealii::Table<2,double> butcher_tableau_gamma;

    dealii::Table<1,double> butcher_tableau_beta;
    
    dealii::Table<1,double> butcher_tableau_delta;

    dealii::Table<1,double> butcher_tableau_b_hat;

    virtual void set_gamma() = 0;

    virtual void set_beta() = 0;
    
    virtual void set_delta() = 0;

    virtual void set_b_hat() = 0;

};

} // ODE namespace
} // PHiLiP namespace

#endif