periodic_entropy_tests.h

#ifndef __PERIODIC_ENTROPY_TESTS_H__
#define __PERIODIC_ENTROPY_TESTS_H__

#include "periodic_cube_flow.h"
#include "physics/euler.h"

namespace PHiLiP {
namespace FlowSolver {

template <int dim, int nstate>
class PeriodicEntropyTests : public PeriodicCubeFlow<dim,nstate>
{
public:

    explicit PeriodicEntropyTests(const Parameters::AllParameters *const parameters_input);

    double compute_entropy(const std::shared_ptr <DGBase<dim, double>> dg) const;


    double get_constant_time_step(std::shared_ptr<DGBase<dim,double>> dg) const override;
protected:

    enum IntegratedQuantityEnum { kinetic_energy, max_wave_speed, numerical_entropy};


    double compute_integrated_quantities(DGBase<dim, double> &dg,
            IntegratedQuantityEnum quantity, 
            const int overintegrate=10 // Overintegrate for KE, don't for num. entropy
            ) const;

    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;
    
    std::string unsteady_data_table_filename_with_extension;
    
    double initial_entropy;

    double previous_numerical_entropy;
    
    double previous_time=0;

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

};


} // FlowSolver namespace
} // PHiLiP namespace
#endif