dg_factory.hpp

#ifndef __DG_FACTORY_H__
#define __DG_FACTORY_H__

#include "dg_base.hpp"

namespace PHiLiP {


#if PHILIP_DIM==1 // dealii::parallel::distributed::Triangulation<dim> does not work for 1D
template <int dim, typename real, typename MeshType = dealii::Triangulation<dim>>
#else
template <int dim, typename real, typename MeshType = dealii::parallel::distributed::Triangulation<dim>>
#endif
class DGFactory
{
    using Triangulation = MeshType;
public:

    static std::shared_ptr< DGBase<dim,real,MeshType> >
        create_discontinuous_galerkin(
        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);

    static std::shared_ptr< DGBase<dim,real,MeshType> >
        create_discontinuous_galerkin(
        const Parameters::AllParameters *const parameters_input,
        const unsigned int degree,
        const unsigned int max_degree_input,
        const std::shared_ptr<Triangulation> triangulation_input);

    static std::shared_ptr< DGBase<dim,real,MeshType> >
        create_discontinuous_galerkin(
        const Parameters::AllParameters *const parameters_input,
        const unsigned int degree,
        const std::shared_ptr<Triangulation> triangulation_input);
};

} // PHiLiP namespace

#endif