Method getPressureForceMatrix

tools/Assembler/ref/Assembler.cpp:11–29

#include <Mesh.h>

ZSparseMatrix Assembler::getPressureForceMatrix() {
    typedef Eigen::Triplet<double> T;
    std::vector<T> triplets;
    size_t dim = getDim();

    for (size_t i=0; i<m_mesh->getNbrBoundaryNodes(); i++) {
        size_t glob_idx = m_mesh->getBoundaryNode(i);
        VectorF normal = m_mesh->getBoundaryNodeNormal(i);

        for (size_t j=0; j<dim; j++) {
            triplets.push_back(T(glob_idx*dim+j  , i, normal[j]));
        }
    }

    Eigen::SparseMatrix<double> P = Eigen::SparseMatrix<double>(
            dim*m_mesh->getNbrNodes(), m_mesh->getNbrBoundaryNodes());
    P.setFromTriplets(triplets.begin(), triplets.end());
    return ZSparseMatrix(P);
}

ZSparseMatrix Assembler::getBdAreaMatrix() {
    typedef Eigen::Triplet<double> T;