time_refinement_study_reference.h

#ifndef __TIME_REFINEMENT_STUDY_REFERENCE__
#define __TIME_REFINEMENT_STUDY_REFERENCE__

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

    Parameters::AllParameters reinit_params_and_refine_timestep(int refinement) const;

    Parameters::AllParameters reinit_params_for_reference_solution(int number_of_timesteps, double final_time) const;

    dealii::LinearAlgebra::distributed::Vector<double> calculate_reference_solution(double final_time) const;
    
    double calculate_L2_error_at_final_time_wrt_reference(
            std::shared_ptr<DGBase<dim,double>> dg,
            const Parameters::AllParameters parameters, 
            double final_time_actual,
            dealii::LinearAlgebra::distributed::Vector<double> reference_solution) const;
};

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

#endif