general_refinement_study.h

#ifndef __TIME_REFINEMENT_STUDY__
#define __TIME_REFINEMENT_STUDY__

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

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

namespace PHiLiP {
namespace Tests {

template <int dim, int nstate>
class GeneralRefinementStudy: public TestsBase
{
public:
    enum RefinementType { timestep, cell_length }; // in the future, can also add refinement by poly_degree here

    GeneralRefinementStudy(
            const Parameters::AllParameters *const parameters_input,
            const dealii::ParameterHandler &parameter_handler_input,
            const RefinementType refinement_type_input
            );

    const dealii::ParameterHandler &parameter_handler;

    int run_test () const override;
protected:
    const RefinementType refinement_type;
    
    const int n_calculations;

    const double refine_ratio;


    int run_refinement_study_and_write_result(const Parameters::AllParameters *parameters_in, const double expected_order, const bool check_only_last_refinement = false, const bool append_to_file = false) const;

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

    virtual std::tuple<double,int> process_and_write_conv_tables(std::shared_ptr<FlowSolver::FlowSolver<dim,nstate>> flow_solver, 
            const Parameters::AllParameters params, 
            double L2_error_old, 
            std::shared_ptr<dealii::ConvergenceTable> convergence_table,
            int refinement,
            const double expected_order) const;
    Parameters::AllParameters reinit_params_and_refine(const Parameters::AllParameters *parameters_in, int refinement, const RefinementType how) const;

};

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

#endif