compute_integral.h

// SPDX-FileCopyrightText: © 2025 PRISMS Center at the University of Michigan
// SPDX-License-Identifier: GNU Lesser General Public Version 2.1
 
#ifndef compute_integral_h
#define compute_integral_h
 
#include <deal.II/matrix_free/matrix_free.h>
 
#include <prismspf/config.h>
#include <prismspf/core/element_volume.h>
 
PRISMS_PF_BEGIN_NAMESPACE
 
template <int dim, int degree, typename number>
class computeIntegral
{
public:
  explicit computeIntegral(const elementVolume<dim, degree, number> &_element_volume);
 
  void
  initialize(
    std::shared_ptr<dealii::MatrixFree<dim, number, dealii::VectorizedArray<number>>>
      _data);
 
  void
  compute_integral();
 
private:
  const elementVolume<dim, degree, double> &element_volume;
 
  std::shared_ptr<dealii::MatrixFree<dim, number, dealii::VectorizedArray<number>>> data;
};
class computeIntegral {…};
 
template <int dim, int degree, typename number>
computeIntegral<dim, degree, number>::computeIntegral(
  const elementVolume<dim, degree, number> &_element_volume)
  : element_volume(_element_volume)
{}
computeIntegral<dim, degree, number>::computeIntegral( {…}
 
template <int dim, int degree, typename number>
void
computeIntegral<dim, degree, number>::initialize(
  std::shared_ptr<dealii::MatrixFree<dim, number, dealii::VectorizedArray<number>>> _data)
{
  data = _data;
}
computeIntegral<dim, degree, number>::initialize( {…}
 
template <int dim, int degree, typename number>
void
computeIntegral<dim, degree, number>::compute_integral()
{}
computeIntegral<dim, degree, number>::compute_integral() {…}
 
PRISMS_PF_END_NAMESPACE
 
#endif