Method getStrainStressMatrix

tools/Assembler/ref/Assembler3D.cpp:304–321

}

ZSparseMatrix Assembler3D::getStrainStressMatrix() {
    typedef Eigen::Triplet<double> T;
    std::vector<T> triplets;

    for (size_t i=0; i<m_mesh->getNbrElements(); i++) {
        for (size_t j=0; j<6; j++) {
            for (size_t k=0; k<6; k++) {
                if (m_D(j,k) != 0) {
                    triplets.push_back(T(i*6+j, i*6+k, m_D(j,k)));
                }
            }
        }
    }

    Eigen::SparseMatrix<double> D = Eigen::SparseMatrix<double>(6*m_mesh->getNbrElements(), 6*m_mesh->getNbrElements());
    D.setFromTriplets(triplets.begin(), triplets.end());
    return ZSparseMatrix(D);
}

ZSparseMatrix Assembler3D::getBdLaplacianMatrix() {
    typedef Eigen::Triplet<double> T;