Method init_boundary_face_adjacencies

tools/Assembler/Mesh/TriangleMesh.cpp:229–253

}

void TriangleMesh::init_boundary_face_adjacencies() {
    const size_t num_boundary_edges = m_boundary_edges.rows();
    const size_t num_boundary_vertices = m_boundary.size();
    std::vector<std::set<size_t> > bd_vertex_edge_adjacencies(num_boundary_vertices);
    for (size_t i=0; i<num_boundary_edges; i++) {
        const VectorI& edge = m_boundary_edges.row(i);
        bd_vertex_edge_adjacencies[edge[0]].insert(i);
        bd_vertex_edge_adjacencies[edge[1]].insert(i);
    }

    std::vector<size_t> adj;
    std::vector<size_t> adj_idx;
    adj_idx.push_back(0);
    for (size_t i=0; i<num_boundary_vertices; i++) {
        const std::set<size_t>& neighbors = bd_vertex_edge_adjacencies[i];
        adj.insert(adj.end(), neighbors.begin(), neighbors.end());
        adj_idx.push_back(adj.size());
    }

    m_boundary_node_edge_adj.resize(adj.size());
    std::copy(adj.begin(), adj.end(), m_boundary_node_edge_adj.data());

    m_boundary_node_edge_adj_idx.resize(adj_idx.size());
    std::copy(adj_idx.begin(), adj_idx.end(), m_boundary_node_edge_adj_idx.data());
}

void TriangleMesh::init_boundary_lengths() {
    const size_t num_edges = m_boundary_edges.rows();