offset_function.h

#ifndef THESIS_ASSEMBLE_OFFSET_FUNCTION_H
#define THESIS_ASSEMBLE_OFFSET_FUNCTION_H
 
#include <lf/assemble/dofhandler.h>
#include <lf/mesh/mesh.h>
#include <lf/mesh/utils/mesh_data_set.h>
 
#include <Eigen/Sparse>
#include <functional>
 
namespace projects::ipdg_stokes {
 
namespace assemble {
 
class PiecewiseBoundaryNormalJumpAssembler {
 public:
  PiecewiseBoundaryNormalJumpAssembler(
      std::shared_ptr<const lf::mesh::Mesh> mesh,
      const lf::mesh::utils::MeshDataSet<bool> &boundary);
 
  Eigen::MatrixXd Eval(const lf::mesh::Entity &entity) const;
 
  bool isActive(const lf::mesh::Entity &entity) const { return true; }
 
 private:
  const std::shared_ptr<const lf::mesh::Mesh> mesh_;
  const lf::mesh::utils::MeshDataSet<bool> &boundary_;
};
 
Eigen::VectorXd createOffsetFunction(
    const std::shared_ptr<const lf::mesh::Mesh> &mesh,
    const lf::mesh::utils::MeshDataSet<bool> &boundary,
    const lf::assemble::DofHandler &dofh,
    const std::function<Eigen::Vector2d(const lf::mesh::Entity &)>
        &dirichlet_data,
    const Eigen::SparseMatrix<double> &A);
 
}  // end namespace assemble
 
}  // end namespace projects::ipdg_stokes
 
#endif  // THESIS_ASSEMBLE_OFFSET_FUNCTION_H