dg_base_state.hpp

#ifndef PHILIP_DG_BASE_STATE_HPP
#define PHILIP_DG_BASE_STATE_HPP

#include <deal.II/distributed/tria.h>

#include "dg_base.hpp"
#include "parameters/all_parameters.h"
namespace PHiLiP {

/*  Contains the objects and functions that need to be templated on the number of state variables.
 */
#if PHILIP_DIM==1 // dealii::parallel::distributed::Triangulation<dim> does not work for 1D
template <int dim, int nspecies, int nstate, typename real, typename MeshType = dealii::Triangulation<dim>>
#else
template <int dim, int nspecies, int nstate, typename real, typename MeshType = dealii::parallel::distributed::Triangulation<dim>>
#endif
class DGBaseState : public DGBase<dim, nspecies, real, MeshType>
{
protected:
    using Triangulation = typename DGBase<dim,nspecies,real,MeshType>::Triangulation;

public:
    using DGBase<dim,nspecies,real,MeshType>::all_parameters; 

    DGBaseState(
        const Parameters::AllParameters *const parameters_input,
        const unsigned int degree,
        const unsigned int max_degree_input,
        const unsigned int grid_degree_input,
        const std::shared_ptr<Triangulation> triangulation_input);

    std::shared_ptr < Physics::PhysicsBase<dim, nspecies, nstate, real > > pde_physics_double;
    std::shared_ptr < Physics::ModelBase<dim, nspecies, nstate, real > > pde_model_double;
    std::unique_ptr < NumericalFlux::NumericalFluxConvective<dim, nspecies, nstate, real > > conv_num_flux_double;
    std::unique_ptr < NumericalFlux::NumericalFluxDissipative<dim, nspecies, nstate, real > > diss_num_flux_double;
    std::shared_ptr <ArtificialDissipationBase<dim,nspecies,nstate>> artificial_dissip;

    std::shared_ptr < Physics::PhysicsBase<dim, nspecies, nstate, FadType > > pde_physics_fad;
    std::shared_ptr < Physics::ModelBase<dim, nspecies, nstate, FadType > > pde_model_fad;
    std::unique_ptr < NumericalFlux::NumericalFluxConvective<dim, nspecies, nstate, FadType > > conv_num_flux_fad;
    std::unique_ptr < NumericalFlux::NumericalFluxDissipative<dim, nspecies, nstate, FadType > > diss_num_flux_fad;

    std::shared_ptr < Physics::PhysicsBase<dim, nspecies, nstate, RadType > > pde_physics_rad;
    std::shared_ptr < Physics::ModelBase<dim, nspecies, nstate, RadType > > pde_model_rad;
    std::unique_ptr < NumericalFlux::NumericalFluxConvective<dim, nspecies, nstate, RadType > > conv_num_flux_rad;
    std::unique_ptr < NumericalFlux::NumericalFluxDissipative<dim, nspecies, nstate, RadType > > diss_num_flux_rad;

    std::shared_ptr < Physics::PhysicsBase<dim, nspecies, nstate, FadFadType > > pde_physics_fad_fad;
    std::shared_ptr < Physics::ModelBase<dim, nspecies, nstate, FadFadType > > pde_model_fad_fad;
    std::unique_ptr < NumericalFlux::NumericalFluxConvective<dim, nspecies, nstate, FadFadType > > conv_num_flux_fad_fad;
    std::unique_ptr < NumericalFlux::NumericalFluxDissipative<dim, nspecies, nstate, FadFadType > > diss_num_flux_fad_fad;

    std::shared_ptr < Physics::PhysicsBase<dim, nspecies, nstate, RadFadType > > pde_physics_rad_fad;
    std::shared_ptr < Physics::ModelBase<dim, nspecies, nstate, RadFadType > > pde_model_rad_fad;
    std::unique_ptr < NumericalFlux::NumericalFluxConvective<dim, nspecies, nstate, RadFadType > > conv_num_flux_rad_fad;
    std::unique_ptr < NumericalFlux::NumericalFluxDissipative<dim, nspecies, nstate, RadFadType > > diss_num_flux_rad_fad;

    void set_physics(
        std::shared_ptr< Physics::PhysicsBase<dim, nspecies, nstate, real       > > pde_physics_double_input,
        std::shared_ptr< Physics::PhysicsBase<dim, nspecies, nstate, FadType    > > pde_physics_fad_input,
        std::shared_ptr< Physics::PhysicsBase<dim, nspecies, nstate, RadType    > > pde_physics_rad_input,
        std::shared_ptr< Physics::PhysicsBase<dim, nspecies, nstate, FadFadType > > pde_physics_fad_fad_input,
        std::shared_ptr< Physics::PhysicsBase<dim, nspecies, nstate, RadFadType > > pde_physics_rad_fad_input);

    virtual void allocate_model_variables();

    virtual void update_model_variables();

    void set_unsteady_model_time_step(const double time_step);

    void set_use_auxiliary_eq();

    void set_store_vol_flux_nodes();

    void set_store_surf_flux_nodes();

protected:

    real evaluate_CFL (std::vector< std::array<real,nstate> > soln_at_q, const real artificial_dissipation, const real cell_diameter, const unsigned int cell_degree);


    void reset_numerical_fluxes();

}; // end of DGBaseState class

}  // namespace PHiLiP

#endif  // PHILIP_DG_BASE_STATE_HPP