parameters_manufactured_solution.h

#ifndef __PARAMETERS_MANUFACTURED_SOLUTION_H__
#define __PARAMETERS_MANUFACTURED_SOLUTION_H__

#include <deal.II/base/parameter_handler.h>
#include "parameters/parameters.h"

#include <deal.II/base/tensor.h>

namespace PHiLiP {

namespace Parameters {

class ManufacturedSolutionParam
{
public:
    bool use_manufactured_source_term;
    
    enum ManufacturedSolutionType{
        zero_solution,
        sine_solution,
        cosine_solution,
        additive_solution,
        exp_solution,
        poly_solution,
        even_poly_solution,
        atan_solution,
        boundary_layer_solution,
        s_shock_solution,
        quadratic_solution,
        example_solution,
        navah_solution_1,
        navah_solution_2,
        navah_solution_3,
        navah_solution_4,
        navah_solution_5
        };
    ManufacturedSolutionType manufactured_solution_type; 

    dealii::Tensor<2,3,double> diffusion_tensor;

    dealii::Tensor<1,3,double> advection_vector;

    double diffusion_coefficient;

    static void declare_parameters (dealii::ParameterHandler &prm);
    void parse_parameters (dealii::ParameterHandler &prm);

    static dealii::Tensor<2,3,double> get_default_diffusion_tensor();
    static dealii::Tensor<1,3,double> get_default_advection_vector();
    static double get_default_diffusion_coefficient();
};

} // Parameters namespace

} // PHiLiP namespace

#endif // __PARAMETERS_MANUFACTURED_SOLUTION_H__