shock_1d.h

#ifndef __BURGERS_STEADYSHOCK_H__
#define __BURGERS_STEADYSHOCK_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 Shock1D: public TestsBase
{
public:
    Shock1D () = delete;

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


    int run_test () const;

protected:
    void print_mesh_info(const dealii::Triangulation<dim> &triangulation,
                         const std::string &filename) const;

    static dealii::Point<dim> warp (const dealii::Point<dim> &p);


    void initialize_perturbed_solution(DGBase<dim,double> &dg, const Physics::PhysicsBase<dim,nstate,double> &physics) const;

    double integrate_solution_over_domain(DGBase<dim,double> &dg) const;
};


//   /// 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