h_refinement_study_isentropic_vortex.h

#ifndef __H_REFINEMENT_STUDY_ISENTROPIC_VORTEX__
#define __H_REFINEMENT_STUDY_ISENTROPIC_VORTEX__

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

#include "dg/dg_base.hpp"
#include "tests.h"

namespace PHiLiP {
namespace Tests {

template <int dim, int nstate>
class HRefinementStudyIsentropicVortex: public TestsBase
{
public:
    HRefinementStudyIsentropicVortex(
            const Parameters::AllParameters *const parameters_input,
            const dealii::ParameterHandler &parameter_handler_input);

    const dealii::ParameterHandler &parameter_handler;
    
    int run_test () const override;
protected:
    const int n_calculations;

    const double refine_ratio;

    void calculate_Lp_error_at_final_time_wrt_function(double &Lp_error_density, 
            double &Lp_error_pressure,
            std::shared_ptr<DGBase<dim,double>> dg,
            const Parameters::AllParameters parameters,
            double final_time, 
            int norm_p) const;

    Parameters::AllParameters reinit_params_and_refine(int refinement) const;

};

} // End of Tests namespace
} // End of PHiLiP namespace

#endif