grid_refinement_study.h

#ifndef __GRID_REFINEMENT_STUDY_H__
#define __GRID_REFINEMENT_STUDY_H__

#include "dg/dg_base.hpp"
#include "grid_refinement/gnu_out.h"
#include "parameters/all_parameters.h"
#include "physics/model.h"
#include "physics/physics.h"
#include "tests.h"

namespace PHiLiP {

namespace Tests {

#if PHILIP_DIM==1 // dealii::parallel::distributed::Triangulation<dim> does not work for 1D
template <int dim, int nstate, typename MeshType = dealii::Triangulation<dim>>
#else
template <int dim, int nstate, typename MeshType = dealii::parallel::distributed::Triangulation<dim>>
#endif
class GridRefinementStudy : public TestsBase
{
public:
    GridRefinementStudy() = delete;

    explicit GridRefinementStudy(
        const Parameters::AllParameters *const parameters_input);

    int run_test() const;

    void get_grid(
        const std::shared_ptr<MeshType>&            grid,
        const Parameters::GridRefinementStudyParam& grs_param) const;

    double approximate_exact_functional(
        const std::shared_ptr<Physics::PhysicsBase<dim,nstate,double>>& physics_double,
        const std::shared_ptr<Physics::PhysicsBase<dim,nstate,Sacado::Fad::DFad<double>>>& physics_adtype,
        const Parameters::AllParameters& param,
        const Parameters::GridRefinementStudyParam& grs_param) const;
};

// constructs the mesh for the test
template <typename MeshType>
class MeshFactory
{
public:
    static std::shared_ptr<MeshType>
    create_MeshType(const MPI_Comm mpi_communicator);
};

// function to perform the formatted output to gnuplot (of the solution error)
void output_gnufig_solution(
    PHiLiP::GridRefinement::GnuFig<double> &gf);

// function to perform the formatted output to gnuplot (of the funcitonal error)
void output_gnufig_functional(
    PHiLiP::GridRefinement::GnuFig<double> &gf);

} // Tests namespace

} // PHiLiP namespace

#endif // __GRID_REFINEMENT_STUDY_H__