anisotropic_mesh_adaptation_cases.h

#ifndef __ANISOTROPICMESHADAPTATIONCASES_H__ 
#define __ANISOTROPICMESHADAPTATIONCASES_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 AnisotropicMeshAdaptationCases : public TestsBase
{
public:
    AnisotropicMeshAdaptationCases(const Parameters::AllParameters *const parameters_input,
                                       const dealii::ParameterHandler &parameter_handler_input);
    
    const dealii::ParameterHandler &parameter_handler;

    int run_test() const;

    void verify_fe_values_shape_hessian(const DGBase<dim, double> &dg) const;
    
    double evaluate_functional(std::shared_ptr<DGBase<dim,double>> dg) const;
}; 

} // Tests namespace
} // PHiLiP namespace

#endif