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 void set_tableau()=0;

    const int n_rk_stages;
    
protected:

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

} // ODE namespace
} // PHiLiP namespace

#endif