periodic_1D_unsteady.h

#ifndef __PERIODIC_1D_UNSTEADY_H__
#define __PERIODIC_1D_UNSTEADY_H__

#include "periodic_cube_flow.h"

namespace PHiLiP {
namespace FlowSolver {

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

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

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

    double get_numerical_entropy(const std::shared_ptr <DGBase<dim, nspecies, double>> dg) const;
protected:
    void compute_unsteady_data_and_write_to_table(
            const std::shared_ptr <ODE::ODESolverBase<dim, nspecies, double>> ode_solver,
            const std::shared_ptr <DGBase<dim, nspecies, double>> dg,
            const std::shared_ptr<dealii::TableHandler> unsteady_data_table,
            const bool do_write_unsteady_data_table_file) override;
    
    std::string unsteady_data_table_filename_with_extension;
    
};


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