euler_gaussian_bump_adjoint.h

#ifndef __EULER_GAUSSIAN_BUMP_ADJOINT_H__
#define __EULER_GAUSSIAN_BUMP_ADJOINT_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 {

class BumpManifoldAdjoint: public dealii::ChartManifold<2,2,2> {
public:
    virtual dealii::Point<2> pull_back(const dealii::Point<2> &space_point) const override;
    virtual dealii::Point<2> push_forward(const dealii::Point<2> &chart_point) const override;
    virtual dealii::DerivativeForm<1,2,2> push_forward_gradient(const dealii::Point<2> &chart_point) const override;
    
    virtual std::unique_ptr<dealii::Manifold<2,2> > clone() const override;
};

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

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

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


    int run_test () const;

protected:

    // Not used?
    // double integrate_entropy_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 //__EULER_GAUSSIAN_BUMP_ADJOINT_H__