build_system_matrix.h

#ifndef THESIS_ASSEMBLE_BUILD_SYSTEM_MATRIX_H
#define THESIS_ASSEMBLE_BUILD_SYSTEM_MATRIX_H
 
#include <lf/assemble/coomatrix.h>
#include <lf/assemble/dofhandler.h>
#include <lf/mesh/entity.h>
#include <lf/mesh/mesh.h>
#include <lf/mesh/utils/mesh_data_set.h>
#include <lf/quad/quad_rule.h>
 
#include <Eigen/Dense>
#include <functional>
#include <tuple>
#include <utility>
 
namespace projects::ipdg_stokes {
 
namespace assemble {
 
std::tuple<lf::assemble::COOMatrix<double>, Eigen::VectorXd>
buildSystemMatrixNoFlow(
    const std::shared_ptr<const lf::mesh::Mesh> &mesh,
    const lf::assemble::DofHandler &dofh,
    const std::function<Eigen::Vector2d(const Eigen::Vector2d &)> &f,
    const std::function<Eigen::Vector2d(const lf::mesh::Entity &)>
        &dirichlet_data,
    double sigma, const lf::quad::QuadRule &quadrule,
    bool modified_penalty = false);
 
std::tuple<lf::assemble::COOMatrix<double>, Eigen::VectorXd, Eigen::VectorXd>
buildSystemMatrixInOutFlow(
    const std::shared_ptr<const lf::mesh::Mesh> &mesh,
    const lf::assemble::DofHandler &dofh,
    const std::function<Eigen::Vector2d(const Eigen::Vector2d &)> &f,
    const std::function<Eigen::Vector2d(const lf::mesh::Entity &)>
        &dirichlet_data,
    double sigma, const lf::quad::QuadRule &quadrule,
    bool modified_penalty = false);
 
}  // end namespace assemble
 
}  // end namespace projects::ipdg_stokes
 
#endif  // THESIS_ASSEMBLE_BUILD_SYSTEM_MATRIX_H