dof_handler.h

// SPDX-FileCopyrightText: © 2025 PRISMS Center at the University of Michigan
// SPDX-License-Identifier: GNU Lesser General Public Version 2.1
 
#ifndef dof_handler_h
#define dof_handler_h
 
#include <deal.II/base/mg_level_object.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_system.h>
 
#include <prismspf/config.h>
#include <prismspf/core/triangulation_handler.h>
#include <prismspf/user_inputs/user_input_parameters.h>
 
#include <map>
#include <vector>
 
PRISMS_PF_BEGIN_NAMESPACE
 
template <int dim>
class dofHandler
{
public:
  explicit dofHandler(const userInputParameters<dim> &_user_inputs);
 
  void
  init(const triangulationHandler<dim>                  &triangulation_handler,
       const std::map<fieldType, dealii::FESystem<dim>> &fe_system);
 
  [[nodiscard]] const std::vector<const dealii::DoFHandler<dim> *> &
  get_dof_handlers() const;
 
  [[nodiscard]] const std::map<unsigned int,
                               dealii::MGLevelObject<dealii::DoFHandler<dim>>> &
  get_mg_dof_handlers() const;
 
  [[nodiscard]] const dealii::DoFHandler<dim> &
  get_mg_dof_handler(unsigned int index, unsigned int level) const;
 
private:
  const userInputParameters<dim> *user_inputs;
 
  std::map<unsigned int, std::unique_ptr<dealii::DoFHandler<dim>>> dof_handlers;
 
  std::vector<const dealii::DoFHandler<dim> *> const_dof_handlers;
 
  bool has_multigrid = false;
 
  std::map<unsigned int, dealii::MGLevelObject<dealii::DoFHandler<dim>>> mg_dof_handlers;
};
class dofHandler {…};
 
PRISMS_PF_END_NAMESPACE
 
#endif