exact_solution.h

#ifndef __EXACT_SOLUTION_H__
#define __EXACT_SOLUTION_H__

// for the exact_solution function
#include <deal.II/lac/vector.h>
#include <deal.II/base/function.h>
#include "parameters/all_parameters.h"

namespace PHiLiP {

template <int dim, int nstate, typename real>
class ExactSolutionFunction : public dealii::Function<dim,real>
{
protected:
    using dealii::Function<dim,real>::value; 

public:
    ExactSolutionFunction();

    virtual real value (const dealii::Point<dim,real> &point, const unsigned int istate = 0) const = 0;
};

template <int dim, int nstate, typename real>
class ExactSolutionFunction_Zero
        : public ExactSolutionFunction<dim,nstate,real>
{
protected:
    using dealii::Function<dim,real>::value; 

public:
    explicit ExactSolutionFunction_Zero (double time_compare);

    const double t; 

    real value (const dealii::Point<dim,real> &point, const unsigned int istate = 0) const override;
};

template <int dim, int nstate, typename real>
class ExactSolutionFunction_1DSine
        : public ExactSolutionFunction<dim,nstate,real>
{
protected:
    using dealii::Function<dim,real>::value; 

public:

    explicit ExactSolutionFunction_1DSine (double time_compare);

    const double t; 

    real value (const dealii::Point<dim,real> &point, const unsigned int istate = 0) const override;
};

template <int dim, int nstate, typename real>
class ExactSolutionFunction_IsentropicVortex
        : public ExactSolutionFunction<dim,nstate,real>
{
protected:
    using dealii::Function<dim,real>::value; 

public:

    explicit ExactSolutionFunction_IsentropicVortex (double time_compare);

    const double t; 

    real value (const dealii::Point<dim,real> &point, const unsigned int istate = 0) const override;
};

template <int dim, int nstate, typename real>
class ExactSolutionFactory
{
protected:    
    using FlowCaseEnum = Parameters::FlowSolverParam::FlowCaseType;

public:
    static std::shared_ptr<ExactSolutionFunction<dim,nstate,real>>
        create_ExactSolutionFunction(const Parameters::FlowSolverParam& flow_solver_parameters, const double time_compare);
};

} // PHiLiP namespace
#endif