invm_handler.h

// SPDX-FileCopyrightText: © 2025 PRISMS Center at the University of Michigan
// SPDX-License-Identifier: GNU Lesser General Public Version 2.1
 
#ifndef invm_handler_h
#define invm_handler_h
 
#include <deal.II/lac/la_parallel_vector.h>
#include <deal.II/matrix_free/matrix_free.h>
 
#include <prismspf/config.h>
#include <prismspf/core/variable_attributes.h>
#include <prismspf/user_inputs/user_input_parameters.h>
 
#include <map>
#include <memory>
 
PRISMS_PF_BEGIN_NAMESPACE
 
template <int dim, int degree, typename number = double>
class invmHandler
{
public:
  using VectorType = dealii::LinearAlgebra::distributed::Vector<number>;
  using size_type  = dealii::VectorizedArray<number>;
 
  explicit invmHandler(
    const std::map<unsigned int, variableAttributes> &_variable_attributes);
 
  void
  initialize(std::shared_ptr<dealii::MatrixFree<dim, number, size_type>> _data);
 
  void
  compute_invm();
 
  void
  recompute_invm();
 
  [[nodiscard]] const VectorType &
  get_invm(const unsigned int &index) const;
 
  void
  clear();
 
private:
  void
  compute_scalar_invm();
 
  void
  compute_vector_invm();
 
  const std::map<unsigned int, variableAttributes> *variable_attributes;
 
  std::shared_ptr<dealii::MatrixFree<dim, number, size_type>> data;
 
  VectorType invm_scalar;
 
  VectorType invm_vector;
 
  bool scalar_needed = false;
 
  bool vector_needed = false;
 
  unsigned int scalar_index = numbers::invalid_index;
 
  unsigned int vector_index = numbers::invalid_index;
 
  number tolerance = 1.0e-15;
};
class invmHandler {…};
 
PRISMS_PF_END_NAMESPACE
 
#endif