euler_bump_optimization.h

#ifndef __EULER_BUMP_OPTIMIZATION_H__
#define __EULER_BUMP_OPTIMIZATION_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 BumpManifold2: 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 EulerBumpOptimization: public TestsBase
{
public:
    EulerBumpOptimization () = delete;

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


    int run_test () const;

private:
    int optimize_target_bump (const unsigned int nx_ffd, const unsigned int poly_degree) const;

protected:

    //  // Integrate entropy over the entire domain to use as a functional.
    //  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