Method init_boundary_nodes

tools/Assembler/Mesh/TetrahedronMesh.cpp:186–209

}

void TetrahedronMesh::init_boundary_nodes() {
    const size_t dim = getDim();
    const size_t num_faces = m_mesh->get_num_faces();
    const size_t num_vertices = m_mesh->get_num_vertices();
    const size_t vertex_per_face = m_mesh->get_vertex_per_face();
    const VectorI& faces = m_mesh->get_faces();

    size_t count = 1;
    m_boundary_idx = VectorI::Zero(num_vertices);
    for (size_t i=0; i<num_faces * vertex_per_face; i++) {
        if (m_boundary_idx[faces[i]] == 0) {
            m_boundary_idx[faces[i]] = count;
            count ++;
        }
    }
    m_boundary_idx.array() -= 1;

    m_boundary.resize(count-1);
    for (size_t i=0; i<num_vertices; i++) {
        if (m_boundary_idx[i] >= 0) {
            m_boundary[m_boundary_idx[i]] = i;
        }
    }
}

void TetrahedronMesh::init_boundary_node_adjacencies() {
    std::vector<size_t> bv_adj;