rrk_numerical_entropy_conservation_check.h

#ifndef __RRK_NUMERICAL_ENTROPY_CONSERVATION_CHECK__
#define __RRK_NUMERICAL_ENTROPY_CONSERVATION_CHECK__

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

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

namespace PHiLiP {
namespace Tests {

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

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

    Parameters::AllParameters reinit_params(bool use_rrk, double time_step_size) const;
    
    int compare_numerical_entropy_to_initial(
            const std::unique_ptr<FlowSolver::FlowSolver<dim,nstate>> &flow_solver,
            const double initial_numerical_entropy,
            const double final_time_actual,
            bool expect_conservation
            ) const;

    int get_numerical_entropy_and_compare_to_initial(
            const Parameters::AllParameters params,
            const double numerical_entropy_initial,
            bool expect_conservation
            ) const;

};

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

#endif