[P]arallel [Hi]gh-order [Li]brary for [P]DEs  Latest
Parallel High-Order Library for PDEs through hp-adaptive Discontinuous Galerkin methods
ode_solver_factory.h
1 #ifndef __ODE_SOLVER_FACTORY__
2 #define __ODE_SOLVER_FACTORY__
3 
4 #include "dg/dg_base.hpp"
5 #include "ode_solver_base.h"
6 #include "parameters/all_parameters.h"
7 #include "reduced_order/pod_basis_base.h"
8 #include "runge_kutta_methods/rk_tableau_base.h"
9 #include "runge_kutta_methods/rk_tableau_butcher_base.h"
10 #include "runge_kutta_methods/low_storage_rk_tableau_base.h"
11 #include "relaxation_runge_kutta/empty_RRK_base.h"
12 
13 namespace PHiLiP {
14 namespace ODE {
16 
18 #if PHILIP_DIM==1
19 template <int dim, int nspecies, typename real, typename MeshType = dealii::Triangulation<dim>>
20 #else
21 template <int dim, int nspecies, typename real, typename MeshType = dealii::parallel::distributed::Triangulation<dim>>
22 #endif
23 class ODESolverFactory
24 {
25 public:
27  static std::shared_ptr<ODESolverBase<dim,nspecies,real,MeshType>> create_ODESolver(std::shared_ptr< DGBase<dim, nspecies, real, MeshType> > dg_input);
28 
30  static std::shared_ptr<ODESolverBase<dim,nspecies,real,MeshType>> create_ODESolver(std::shared_ptr< DGBase<dim, nspecies, real, MeshType> > dg_input, std::shared_ptr<ProperOrthogonalDecomposition::PODBase<dim,nspecies>> pod);
31 
33  static std::shared_ptr<ODESolverBase<dim,nspecies,real,MeshType>> create_ODESolver(std::shared_ptr< DGBase<dim, nspecies, real, MeshType> > dg_input, std::shared_ptr<ProperOrthogonalDecomposition::PODBase<dim,nspecies>> pod, Epetra_Vector weights);
34 
36  static std::shared_ptr<ODESolverBase<dim,nspecies,real,MeshType>> create_ODESolver_manual(Parameters::ODESolverParam::ODESolverEnum ode_solver_type, std::shared_ptr< DGBase<dim, nspecies, real, MeshType> > dg_input);
37 
39  static std::shared_ptr<ODESolverBase<dim,nspecies,real,MeshType>> create_ODESolver_manual(Parameters::ODESolverParam::ODESolverEnum ode_solver_type, std::shared_ptr< DGBase<dim, nspecies, real, MeshType> > dg_input, std::shared_ptr<ProperOrthogonalDecomposition::PODBase<dim,nspecies>> pod);
40 
42  static std::shared_ptr<ODESolverBase<dim,nspecies,real,MeshType>> create_ODESolver_manual(Parameters::ODESolverParam::ODESolverEnum ode_solver_type, std::shared_ptr< DGBase<dim, nspecies, real, MeshType> > dg_input, std::shared_ptr<ProperOrthogonalDecomposition::PODBase<dim,nspecies>> pod, Epetra_Vector weights);
43 
45  static void display_error_ode_solver_factory(Parameters::ODESolverParam::ODESolverEnum ode_solver_type, bool reduced_order);
46 
48  static std::shared_ptr<ODESolverBase<dim,nspecies,real,MeshType>> create_RungeKuttaODESolver(std::shared_ptr< DGBase<dim, nspecies, real, MeshType> > dg_input);
49 
51  static std::shared_ptr<ODESolverBase<dim,nspecies,real,MeshType>> create_RungeKuttaODESolver(std::shared_ptr< DGBase<dim, nspecies, real, MeshType> > dg_input, std::shared_ptr<ProperOrthogonalDecomposition::PODBase<dim,nspecies>> pod);
52 
54  static std::shared_ptr<RKTableauBase<dim,real,MeshType>> create_RKTableau(std::shared_ptr< DGBase<dim,nspecies,real,MeshType> > dg_input);
55 
57  static std::shared_ptr<EmptyRRKBase<dim,nspecies,real,MeshType>> create_RRKObject(std::shared_ptr< DGBase<dim,nspecies,real,MeshType> > dg_input,
58  std::shared_ptr<RKTableauBase<dim,real,MeshType>> rk_tableau);
59 
60 };
61 
62 } // ODE namespace
63 } // PHiLiP namespace
64 
65 #endif
PHiLiP::ODE::ODESolverFactory::create_RRKObject
static std::shared_ptr< EmptyRRKBase< dim, nspecies, real, MeshType > > create_RRKObject(std::shared_ptr< DGBase< dim, nspecies, real, MeshType > > dg_input, std::shared_ptr< RKTableauBase< dim, real, MeshType >> rk_tableau)
Creates an RRK object with specified RRK type; if no RRK is being used, creates an RRK object with em...
Definition: ode_solver_factory.cpp:329
PHiLiP
Files for the baseline physics.
Definition: ADTypes.hpp:10
PHiLiP::ODE::ODESolverFactory::create_ODESolver_manual
static std::shared_ptr< ODESolverBase< dim, nspecies, real, MeshType > > create_ODESolver_manual(Parameters::ODESolverParam::ODESolverEnum ode_solver_type, std::shared_ptr< DGBase< dim, nspecies, real, MeshType > > dg_input)
Creates either implicit or explicit ODE solver based on manual input (no POD basis given) ...
Definition: ode_solver_factory.cpp:73
PHiLiP::ODE::ODESolverFactory
Create specified ODE solver as ODESolverBase object.
Definition: ode_solver_factory.h:23
PHiLiP::ODE::ODESolverFactory::create_RungeKuttaODESolver
static std::shared_ptr< ODESolverBase< dim, nspecies, real, MeshType > > create_RungeKuttaODESolver(std::shared_ptr< DGBase< dim, nspecies, real, MeshType > > dg_input)
Creates an ODESolver object based on the specified RK method, including derived classes.
Definition: ode_solver_factory.cpp:168
PHiLiP::ODE::RKTableauBase
Base class for storing the RK method.
Definition: rk_tableau_base.h:19
PHiLiP::ODE::ODESolverFactory::display_error_ode_solver_factory
static void display_error_ode_solver_factory(Parameters::ODESolverParam::ODESolverEnum ode_solver_type, bool reduced_order)
Output error message for Implicit and Explicit solver.
Definition: ode_solver_factory.cpp:122
PHiLiP::ODE::ODESolverFactory::create_RKTableau
static std::shared_ptr< RKTableauBase< dim, real, MeshType > > create_RKTableau(std::shared_ptr< DGBase< dim, nspecies, real, MeshType > > dg_input)
Creates an RKTableau object based on the specified RK method.
Definition: ode_solver_factory.cpp:289
PHiLiP::ODE::ODESolverFactory::create_ODESolver
static std::shared_ptr< ODESolverBase< dim, nspecies, real, MeshType > > create_ODESolver(std::shared_ptr< DGBase< dim, nspecies, real, MeshType > > dg_input)
Creates either implicit or explicit ODE solver based on parameter value(no POD basis given) ...
Definition: ode_solver_factory.cpp:22
PHiLiP::DGBase
DGBase is independent of the number of state variables.
Definition: dg_base.hpp:82