AdaptiveRefinement.h

#ifndef INCLUDE_ADAPTIVEREFINEMENT_H_
#define INCLUDE_ADAPTIVEREFINEMENT_H_
 
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/distributed/grid_refinement.h>
#include <deal.II/distributed/solution_transfer.h>
#include <deal.II/distributed/tria.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_system.h>
#include <deal.II/fe/fe_values.h>
#include <deal.II/lac/la_parallel_vector.h>
 
#include <core/fields.h>
#include <core/userInputParameters.h>
 
using namespace dealii;
 
template <int dim, int degree>
class AdaptiveRefinement
{
public:
  AdaptiveRefinement(
    const userInputParameters<dim>                                    &_userInputs,
    parallel::distributed::Triangulation<dim>                         &_triangulation,
    std::vector<Field<dim>>                                           &_fields,
    std::vector<dealii::LinearAlgebra::distributed::Vector<double> *> &_solutionSet,
    std::vector<parallel::distributed::SolutionTransfer<
      dim,
      dealii::LinearAlgebra::distributed::Vector<double>> *>          &_soltransSet,
    std::vector<FESystem<dim> *>                                      &_FESet,
    std::vector<DoFHandler<dim> *>                 &_dofHandlersSet_nonconst,
    std::vector<const AffineConstraints<double> *> &_constraintsDirichletSet,
    std::vector<const AffineConstraints<double> *> &_constraintsOtherSet);
 
  void
  do_adaptive_refinement(unsigned int _currentIncrement);
 
  void
  refine_grid();
 
protected:
  void
  adaptive_refinement_criterion();
 
private:
  userInputParameters<dim> userInputs;
 
  parallel::distributed::Triangulation<dim> &triangulation;
 
  std::vector<Field<dim>> &fields;
 
  std::vector<dealii::LinearAlgebra::distributed::Vector<double> *> &solutionSet;
 
  std::vector<parallel::distributed::SolutionTransfer<
    dim,
    dealii::LinearAlgebra::distributed::Vector<double>> *> &soltransSet;
 
  std::vector<FESystem<dim> *> &FESet;
 
  std::vector<DoFHandler<dim> *> &dofHandlersSet_nonconst;
 
  std::vector<const AffineConstraints<double> *> &constraintsDirichletSet;
 
  std::vector<const AffineConstraints<double> *> &constraintsOtherSet;
};
 
#endif