piecewise_const_element_vector_provider.h

#ifndef THESIS_ASSEMBLE_PIECEWISE_CONST_ELEMENT_VECTOR_PROVIDER_H
#define THESIS_ASSEMBLE_PIECEWISE_CONST_ELEMENT_VECTOR_PROVIDER_H
 
#include <lf/mesh/entity.h>
#include <lf/mesh/utils/mesh_data_set.h>
#include <lf/quad/quad_rule.h>
 
#include <Eigen/Dense>
#include <functional>
#include <utility>
 
namespace projects::ipdg_stokes {
 
namespace assemble {
 
class PiecewiseConstElementVectorProvider {
 public:
  PiecewiseConstElementVectorProvider(
      double sigma, std::function<Eigen::Vector2d(const Eigen::Vector2d &)> f,
      lf::quad::QuadRule quadrule,
      const lf::mesh::utils::MeshDataSet<bool> &boundary,
      const lf::mesh::utils::MeshDataSet<Eigen::Vector2d> &dirichlet_data);
 
  Eigen::VectorXd Eval(const lf::mesh::Entity &entity) const;
 
  bool isActive(const lf::mesh::Entity &entity) const { return true; }
 
 private:
  const double sigma_;
  const std::function<Eigen::Vector2d(const Eigen::Vector2d &)> f_;
  const lf::quad::QuadRule quadrule_;
  const lf::mesh::utils::MeshDataSet<bool> &boundary_;
  const lf::mesh::utils::MeshDataSet<Eigen::Vector2d> &dirichlet_data_;
};
 
}  // end namespace assemble
 
}  // end namespace projects::ipdg_stokes
 
#endif  // THESIS_ASSEMBLE_PIECEWISE_CONST_ELEMENT_VECTOR_PROVIDER_H