euler_entropy_waves.h

#ifndef __EULER_ENTROPY_WAVES_H__
#define __EULER_ENTROPY_WAVES_H__

#include "dg/dg_base.hpp"
#include "parameters/all_parameters.h"
#include "physics/euler.h"
#include "tests.h"

namespace PHiLiP {
namespace Tests {


template <int dim, typename real>
class EulerEntropyWavesFunction : public dealii::Function<dim,real>
{
public:

    EulerEntropyWavesFunction (const Physics::Euler<dim, dim+2, real> euler_physics, const real dimensional_density_inf);

    const Physics::Euler<dim, dim+2, real> euler_physics; 

    const real dimensional_density_inf; 
    real Q_inf; 


    real value (const dealii::Point<dim> &point, const unsigned int istate = 0) const;

private:
    dealii::Point<2> advected_location(const dealii::Point<2> old_location) const;

};

template <int dim, int nstate>
class EulerEntropyWaves: public TestsBase
{
public:
    EulerEntropyWaves () = delete;

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


    int run_test () const;

};


} // Tests namespace
} // PHiLiP namespace
#endif