numerical_flux_factory.hpp

#ifndef __NUMERICAL_FLUX_FACTORY__
#define __NUMERICAL_FLUX_FACTORY__

#include "physics/physics.h"
#include "dg/artificial_dissipation.h"

#include "convective_numerical_flux.hpp"
#include "viscous_numerical_flux.hpp"

namespace PHiLiP {
namespace NumericalFlux {

template <int dim, int nspecies, int nstate, typename real>
class NumericalFluxFactory
{
public:
    static std::unique_ptr < NumericalFluxConvective<dim,nspecies,nstate,real> >
        create_convective_numerical_flux(
            const AllParam::ConvectiveNumericalFlux conv_num_flux_type,
            const AllParam::PartialDifferentialEquation pde_type,
            const AllParam::ModelType model_type,
            std::shared_ptr<Physics::PhysicsBase<dim, nspecies, nstate, real>> physics_input);
    
    static std::unique_ptr < NumericalFluxDissipative<dim,nspecies,nstate,real> >
        create_dissipative_numerical_flux(
            const AllParam::DissipativeNumericalFlux diss_num_flux_type,
            std::shared_ptr<Physics::PhysicsBase<dim, nspecies, nstate, real>> physics_input,
            std::shared_ptr<ArtificialDissipationBase<dim, nspecies, nstate>>  artificial_dissipation_input);

protected:
    static std::unique_ptr< NumericalFluxConvective<dim,nspecies,nstate,real> > 
        create_euler_based_convective_numerical_flux(
            const AllParam::ConvectiveNumericalFlux conv_num_flux_type,
            const AllParam::PartialDifferentialEquation pde_type,
            const AllParam::ModelType model_type,
            std::shared_ptr<Physics::PhysicsBase<dim, nspecies, nstate, real>> physics_input);
};

} // NumericalFlux namespace
} // PHiLiP namespace

#endif