time_refinement_study.h

#ifndef __TIME_REFINEMENT_STUDY__
#define __TIME_REFINEMENT_STUDY__

#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 TimeRefinementStudy: public TestsBase
{
public:
    TimeRefinementStudy(
            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_time_calculations;

    const double refine_ratio;

    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;

    Parameters::AllParameters reinit_params_and_refine_timestep(int refinement) const;

};

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

#endif