positivity_preserving_tests.h

#ifndef __POSITIVITY_TESTS_H__
#define __POSITIVITY_TESTS_H__

#include "flow_solver_case_base.h"
#include "cube_flow_uniform_grid.h"
#include "dg/dg_base.hpp"

namespace PHiLiP {
namespace FlowSolver {

template <int dim, int nstate>
class PositivityPreservingTests : public CubeFlow_UniformGrid<dim, nstate>
{
#if PHILIP_DIM==1
    using Triangulation = dealii::Triangulation<PHILIP_DIM>;
 #else
    using Triangulation = dealii::parallel::distributed::Triangulation<PHILIP_DIM>;
 #endif

 public:
    explicit PositivityPreservingTests(const Parameters::AllParameters *const parameters_input);
     
    std::shared_ptr<Triangulation> generate_grid() const override;

    void display_additional_flow_case_specific_parameters() const override;

 protected:
    using CubeFlow_UniformGrid<dim, nstate>::get_adaptive_time_step;

    using CubeFlow_UniformGrid<dim, nstate>::get_adaptive_time_step_initial;

    using CubeFlow_UniformGrid<dim, nstate>::update_maximum_local_wave_speed;

    void check_positivity_density(DGBase<dim, double>& dg);

    double compute_integrated_entropy(DGBase<dim, double>& dg) const;

    const std::string unsteady_data_table_filename_with_extension;

    using FlowSolverCaseBase<dim, nstate>::compute_unsteady_data_and_write_to_table;
    void compute_unsteady_data_and_write_to_table(
        const std::shared_ptr<ODE::ODESolverBase<dim, double>> ode_solver,
        const std::shared_ptr <DGBase<dim, double>> dg,
        const std::shared_ptr<dealii::TableHandler> unsteady_data_table) override;
 
 private:
    double maximum_local_wave_speed;

    double initial_entropy;

    double previous_numerical_entropy;

    // euler physics pointer for computing physical quantities.
    std::shared_ptr < Physics::Euler<dim, nstate, double > > euler_physics;

};

} // FlowSolver namespace
} // PHiLiP namespace

#endif