Method getBdAreaMatrix

tools/Assembler/ref/Assembler.cpp:31–50

}

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

    size_t num_bd_vertices = m_mesh->getNbrBoundaryNodes();
    for (size_t i=0; i<num_bd_vertices; i++) {
        VectorI bd_faces = m_mesh->getBoundaryNodeAdjacentBoundaryFaces(i);
        double area = 0.0;
        for (size_t j=0; j<bd_faces.size(); j++) {
            area += m_mesh->getBoundaryFaceArea(bd_faces[j]);
        }
        area /= double(dim);
        triplets.push_back(T(i, i, area));
    }

    Eigen::SparseMatrix<double> Pb = Eigen::SparseMatrix<double>(num_bd_vertices, num_bd_vertices);
    Pb.setFromTriplets(triplets.begin(), triplets.end());
    return ZSparseMatrix(Pb);
}

// Item is basically a <key, value> pair.
typedef std::pair<int, double> Item;