euler_split_inviscid_taylor_green_vortex.h

#ifndef __EULER_SPLIT_TAYLOR_GREEN_H__
#define __EULER_SPLIT_TAYLOR_GREEN_H__


#include <deal.II/base/tensor.h>
#include <deal.II/base/function.h>
#include <deal.II/numerics/data_out.h>
#include <deal.II/numerics/vector_tools.h>
#include <deal.II/numerics/solution_transfer.h>
#include <deal.II/base/numbers.h>
#include <deal.II/base/function_parser.h>
#include <deal.II/grid/grid_generator.h>
#include <deal.II/grid/grid_refinement.h>
#include <deal.II/grid/grid_tools.h>
#include <deal.II/grid/grid_out.h>
#include <deal.II/grid/grid_in.h>

#include <deal.II/fe/mapping_q.h>
#include "tests.h"


#include "parameters/all_parameters.h"
#include "parameters/parameters.h"
#include "physics/physics_factory.h"
#include "physics/physics.h"
#include "ode_solver/ode_solver_factory.h"

#include<fenv.h>

namespace PHiLiP {
namespace Tests {


template <int dim, int nstate>
class EulerTaylorGreen : public TestsBase
{
public:

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


    int run_test() const override;

private:

    double compute_kinetic_energy(const std::shared_ptr < DGBase<dim, double> > &dg, unsigned int poly_degree) const;

    double compute_entropy(const std::shared_ptr < DGBase<dim, double> > &dg, unsigned int poly_degree) const;

   // double compute_change_in_entropy(std::shared_ptr < DGBase<dim, double> > &dg, unsigned int poly_degree) const;
    std::array<double,2> compute_change_in_entropy(const std::shared_ptr < DGBase<dim, double> > &dg, unsigned int poly_degree) const;
    /* Based off unpublished work can show that the volume term vanishes for kinetic energy
    * from its corresponding shuffle condition.
    */
    double compute_volume_term(const std::shared_ptr < DGBase<dim, double> > &dg, unsigned int poly_degree) const;
    double get_timestep(const std::shared_ptr < DGBase<dim, double> > &dg, unsigned int poly_degree, const double delta_x) const;
};


} //Tests
} //PHiLiP

#endif