cube_flow_uniform_grid.h

#ifndef __CUBE_FLOW_UNIFORM_GRID__
#define __CUBE_FLOW_UNIFORM_GRID__

#include "flow_solver_case_base.h"
#include "dg/dg_base.hpp"

namespace PHiLiP {
namespace FlowSolver {

template <int dim, int nspecies, int nstate>
class CubeFlow_UniformGrid : public FlowSolverCaseBase<dim, nspecies, nstate>
{
public:
    explicit CubeFlow_UniformGrid(const Parameters::AllParameters *const parameters_input);

    virtual double get_adaptive_time_step(std::shared_ptr<DGBase<dim,nspecies,double>> dg) const override;

    virtual double get_adaptive_time_step_initial(std::shared_ptr<DGBase<dim,nspecies,double>> dg) override;

    virtual void update_maximum_local_wave_speed(DGBase<dim, nspecies, double> &dg);
 
protected:
    double maximum_local_wave_speed;

    std::shared_ptr< Physics::PhysicsBase<dim,nspecies,nstate,double> > pde_physics;

};

} // FlowSolver namespace
} // PHiLiP namespace

#endif