mesh_adaptation.h

#ifndef __MESHADAPTATION_H__
#define __MESHADAPTATION_H__

#include <deal.II/distributed/grid_refinement.h>
#include <deal.II/distributed/shared_tria.h>
#include <deal.II/distributed/tria.h>
#include <deal.II/grid/grid_refinement.h>
#include <deal.II/grid/tria.h>

#include "dg/dg_base.hpp"
#include "mesh_error_estimate.h"
#include "mesh_error_factory.h"
#include "parameters/all_parameters.h"

namespace PHiLiP {


#if PHILIP_DIM==1
template <int dim, typename real, typename MeshType = dealii::Triangulation<dim>>
#else
template <int dim, typename real, typename MeshType = dealii::parallel::distributed::Triangulation<dim>>
#endif

class MeshAdaptation
{
public:

    MeshAdaptation(std::shared_ptr< DGBase<dim, real, MeshType> > dg_input, const Parameters::MeshAdaptationParam *const mesh_adaptation_param_input);

    ~MeshAdaptation() = default;

    std::unique_ptr<MeshErrorEstimateBase<dim, real, MeshType>> mesh_error;

    std::shared_ptr<DGBase<dim,real,MeshType>> dg;

    void adapt_mesh();

    int current_mesh_adaptation_cycle;

    const Parameters::MeshAdaptationParam *const mesh_adaptation_param;

protected:
    
    void fixed_fraction_isotropic_refinement_and_coarsening();
    
    void smoothness_sensor_based_hp_refinement();
    
    dealii::Vector<real> cellwise_errors;

    dealii::ConditionalOStream pcout;

};

} // namespace PHiLiP

#endif