periodic_cube_flow.h

#ifndef __PERIODIC_CUBE_FLOW_H__
#define __PERIODIC_CUBE_FLOW_H__

#include "flow_solver_case_base.h"
#include "cube_flow_uniform_grid.h"

namespace PHiLiP {
namespace FlowSolver {

#if PHILIP_DIM==1
using Triangulation = dealii::Triangulation<PHILIP_DIM>;
#else
using Triangulation = dealii::parallel::distributed::Triangulation<PHILIP_DIM>;
#endif

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

    std::shared_ptr<Triangulation> generate_grid() const override;

protected:
    const int number_of_cells_per_direction; 
    const double domain_left; 
    const double domain_right; 
    const double domain_size; 

    virtual void display_additional_flow_case_specific_parameters() const override;

    void display_grid_parameters() const;
};

} // FlowSolver namespace
} // PHiLiP namespace
#endif