matrix_free_handler.h

// SPDX-FileCopyrightText: © 2025 PRISMS Center at the University of Michigan
// SPDX-License-Identifier: GNU Lesser General Public Version 2.1
 
#ifndef matrix_free_handler_h
#define matrix_free_handler_h
 
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/mapping.h>
#include <deal.II/lac/affine_constraints.h>
#include <deal.II/matrix_free/matrix_free.h>
 
#include <prismspf/config.h>
#include <prismspf/user_inputs/user_input_parameters.h>
 
#include <memory>
#include <vector>
 
PRISMS_PF_BEGIN_NAMESPACE
 
template <int dim, typename number = double>
class matrixfreeHandler
{
public:
  matrixfreeHandler(const userInputParameters<dim> &_user_inputs);
 
  void
  reinit(const dealii::Mapping<dim>              &mapping,
         const dealii::DoFHandler<dim>           &dof_handler,
         const dealii::AffineConstraints<number> &constraint,
         const dealii::Quadrature<1>             &quad);
 
  void
  reinit(const dealii::Mapping<dim>                                   &mapping,
         const std::vector<const dealii::DoFHandler<dim> *>           &dof_handler,
         const std::vector<const dealii::AffineConstraints<number> *> &constraint,
         const dealii::Quadrature<1>                                  &quad);
 
  void
  reinit(const dealii::Mapping<dim>                                   &mapping,
         const std::vector<const dealii::DoFHandler<dim> *>           &dof_handler,
         const std::vector<const dealii::AffineConstraints<number> *> &constraint,
         const std::vector<dealii::Quadrature<1>>                     &quad);
 
  [[nodiscard]] std::shared_ptr<dealii::MatrixFree<dim, number>>
  get_matrix_free() const;
 
private:
  const userInputParameters<dim> &user_inputs;
 
  std::shared_ptr<dealii::MatrixFree<dim, number>> matrix_free_object;
 
  typename dealii::MatrixFree<dim, number>::AdditionalData additional_data;
};
class matrixfreeHandler {…};
 
PRISMS_PF_END_NAMESPACE
 
#endif