piecewise_const_element_matrix_provider.h

#ifndef THESIS_ASSEMBLE_PIECEWISE_CONST_ELEMENT_MATRIX_PROVIDER_H
#define THESIS_ASSEMBLE_PIECEWISE_CONST_ELEMENT_MATRIX_PROVIDER_H
 
#include <lf/mesh/utils/mesh_data_set.h>
 
#include <Eigen/Dense>
 
namespace projects::ipdg_stokes {
 
namespace assemble {
 
class PiecewiseConstElementMatrixProvider {
 public:
  PiecewiseConstElementMatrixProvider(
      double sigma, const lf::mesh::utils::MeshDataSet<bool> &boundary,
      bool modified = false);
 
  Eigen::MatrixXd Eval(const lf::mesh::Entity &entity) const;
 
  bool isActive(const lf::mesh::Entity &entity) const { return true; }
 
 private:
  const double sigma_;
  const lf::mesh::utils::MeshDataSet<bool> &boundary_;
  const bool modified_;
};
 
}  // end namespace assemble
 
}  // end namespace projects::ipdg_stokes
 
#endif  // THESIS_ASSEMBLE_PIECEWISE_CONST_ELEMENT_MATRIX_PROVIDER_H