parameters_grid_refinement_study.h

#ifndef __PARAMETERS_GRID_REFINEMENT_STUDY_H__
#define __PARAMETERS_GRID_REFINEMENT_STUDY_H__

#include <array>

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

#include "parameters/parameters_functional.h"
#include "parameters/parameters_grid_refinement.h"
#include "parameters/parameters_manufactured_solution.h"
#include "parameters/parameters_manufactured_convergence_study.h"
#include "parameters/parameters.h"

namespace PHiLiP {

namespace Parameters {

class GridRefinementStudyParam
{
    static const unsigned int MAX_REFINEMENTS = 10;

public:
    FunctionalParam functional_param;

    ManufacturedSolutionParam manufactured_solution_param;

    std::array<GridRefinementParam, MAX_REFINEMENTS> grid_refinement_param_vector;

    unsigned int poly_degree;

    unsigned int poly_degree_max;
    unsigned int poly_degree_grid;

    unsigned int num_refinements; 

    using GridEnum = Parameters::ManufacturedConvergenceStudyParam::GridEnum;
    GridEnum grid_type;
    std::string input_grid;
    
    double grid_left; 
    double grid_right;

    unsigned int grid_size;


    bool use_interpolation;

    bool approximate_functional;
    double functional_value;

    bool output_vtk;
    bool output_adjoint_vtk;

    bool output_solution_error;
    bool output_functional_error;

    bool output_gnuplot_solution;
    bool output_gnuplot_functional;
    bool refresh_gnuplot;

    bool output_solution_time;
    bool output_adjoint_time;

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

} // Parameters namespace

} // PHiLiP namespace

#endif // __PARAMETERS_GRID_REFINEMENT_STUDY_H__