triangulation_handler.h

// SPDX-FileCopyrightText: © 2025 PRISMS Center at the University of Michigan
// SPDX-License-Identifier: GNU Lesser General Public Version 2.1
 
#ifndef triangulation_handler_h
#define triangulation_handler_h
 
#include <deal.II/distributed/tria.h>
#include <deal.II/grid/tria.h>
 
#include <prismspf/config.h>
#include <prismspf/core/solution_handler.h>
#include <prismspf/user_inputs/user_input_parameters.h>
 
#include <memory>
#include <string>
 
PRISMS_PF_BEGIN_NAMESPACE
 
template <int dim>
class triangulationHandler
{
public:
  using Triangulation =
    std::conditional_t<dim == 1,
                       dealii::Triangulation<dim>,
                       dealii::parallel::distributed::Triangulation<dim>>;
 
  explicit triangulationHandler(const userInputParameters<dim> &_user_inputs);
 
  [[nodiscard]] const Triangulation &
  get_triangulation() const;
 
  [[nodiscard]] const std::vector<std::shared_ptr<const dealii::Triangulation<dim>>> &
  get_mg_triangulation() const;
 
  [[nodiscard]] const dealii::Triangulation<dim> &
  get_mg_triangulation(unsigned int level) const;
 
  [[nodiscard]] unsigned int
  get_n_global_levels() const;
 
  [[nodiscard]] unsigned int
  get_mg_min_level() const;
 
  [[nodiscard]] unsigned int
  get_mg_max_level() const;
 
  void
  generate_mesh();
 
  void
  adaptively_refine_mesh(solutionHandler<dim> &solution_handler);
 
  void
  export_triangulation_as_vtk(const std::string &filename) const;
 
private:
  void
  mark_boundaries() const;
 
  void
  mark_periodic();
 
  const userInputParameters<dim> *user_inputs;
 
  std::shared_ptr<Triangulation> triangulation;
 
  std::vector<std::shared_ptr<const dealii::Triangulation<dim>>> coarsened_triangulations;
 
  bool has_multigrid = false;
 
  unsigned int min_level = 0;
 
  unsigned int max_level = 0;
};
class triangulationHandler {…};
 
PRISMS_PF_END_NAMESPACE
 
#endif