physics_factory.h

#ifndef __PHYSICS_FACTORY__
#define __PHYSICS_FACTORY__

#include "parameters/all_parameters.h"
#include "physics.h"
#include "model.h"

namespace PHiLiP {
namespace Physics {

template <int dim, int nstate, typename real>
class PhysicsFactory
{
public:
    static std::shared_ptr< PhysicsBase<dim,nstate,real> >
        create_Physics(const Parameters::AllParameters               *const parameters_input,
                       std::shared_ptr< ModelBase<dim,nstate,real> > model_input = nullptr);

    static std::shared_ptr< PhysicsBase<dim,nstate,real> >
        create_Physics(
            const Parameters::AllParameters                              *const parameters_input,
            const Parameters::AllParameters::PartialDifferentialEquation pde_type,
            std::shared_ptr< ModelBase<dim,nstate,real> >                model_input = nullptr);

private:
    static std::shared_ptr< PhysicsBase<dim,nstate,real> >
        create_Physics_Model(
            const Parameters::AllParameters                           *const parameters_input,
            std::shared_ptr< ManufacturedSolutionFunction<dim,real> > manufactured_solution_function = nullptr,
            std::shared_ptr< ModelBase<dim,nstate,real> >             model_input = nullptr);
};


} // Physics namespace
} // PHiLiP namespace

#endif