euler_cylinder.h

#ifndef __EULER_CYLINDER_H__
#define __EULER_CYLINDER_H__

#include <deal.II/grid/manifold_lib.h>

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

namespace PHiLiP {
namespace Tests {

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

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

    // Warp grid into Gaussian bump
    static dealii::Point<dim> warp (const dealii::Point<dim> &p);


    int run_test () const override;

};


//   /// Manufactured grid convergence
//   /** Currently the main function as all my test cases simply
//    *  check for optimal convergence of the solution
//    */
//   template<int dim>
//   int manufactured_grid_convergence (Parameters::AllParameters &parameters);

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