grid_study.h

#ifndef __GRID_STUDY_H__
#define __GRID_STUDY_H__

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

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

namespace PHiLiP {
namespace Tests {

template <int dim, int nstate>
class GridStudy: public TestsBase
{
public:
    GridStudy () = delete;

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


    int run_test () const;

protected:
    void print_mesh_info(const dealii::Triangulation<dim> &triangulation,
                         const std::string &filename) const;

    static dealii::Point<dim> warp (const dealii::Point<dim> &p);


    void initialize_perturbed_solution(DGBase<dim,double> &dg, const Physics::PhysicsBase<dim,nstate,double> &physics) const;

    double integrate_solution_over_domain(DGBase<dim,double> &dg) const;

    std::string get_convergence_tables_baseline_filename(const Parameters::AllParameters *const parameters_input) const;

    void write_convergence_table_to_output_file(
        const std::string error_filename_baseline, 
        const dealii::ConvergenceTable convergence_table,
        const unsigned int poly_degree) const;
};


//   /// Manufactured grid convergence
//   /** Currently the main function as all my test cases simply
//    *  check for optimal convergence of the solution
//    */
//   template<int dim>
//   int manufactured_grid_convergence (Parameters::AllParameters &parameters);

} // Tests namespace
} // PHiLiP namespace
#endif