Method init_voxel_adjacencies

src/Connectivity/MeshConnectivity.cpp:326–393

}

void MeshConnectivity::init_voxel_adjacencies(Mesh* mesh) {
    if (voxel_adjacencies_computed()) return;
    if (!vertex_adjacencies_computed())
        init_vertex_adjacencies(mesh);

    const size_t num_voxels = mesh->get_num_voxels();
    const size_t vertex_per_face = mesh->get_vertex_per_face();
    const size_t vertex_per_voxel = mesh->get_vertex_per_voxel();

    std::vector<IndexSet> neighbor_faces(num_voxels);
    std::vector<IndexSet> neighbor_voxels(num_voxels);

    for (size_t i=0; i<num_voxels; i++) {
        const VectorI voxel = mesh->get_voxel(i);
        IndexSet& voxel_neighbors = neighbor_faces[i];

        assert(voxel.size() == vertex_per_voxel);
        std::map<int, int> face_counter;
        for (size_t j=0; j<vertex_per_voxel; j++) {
            const VectorI neighbors = get_vertex_adjacent_faces(voxel[j]);
            const size_t num_neighbors = neighbors.size();
            for (size_t k=0; k<num_neighbors; k++) {
                if (face_counter.find(neighbors[k]) == face_counter.end()) {
                    face_counter[neighbors[k]] = 1;
                } else {
                    face_counter[neighbors[k]] += 1;
                }
            }
        }
        for (const auto& entry : face_counter) {
            if (entry.second == vertex_per_face) {
                voxel_neighbors.insert(entry.first);
            }
        }
    }

    for (size_t i=0; i<num_voxels; i++) {
        VectorI voxel = mesh->get_voxel(i);
        IndexSet& voxel_neighbors = neighbor_voxels[i];

        std::map<int, int> voxel_counter;
        for (size_t j=0; j<vertex_per_voxel; j++) {
            VectorI neighbors = get_vertex_adjacent_voxels(voxel[j]);
            size_t num_neighbors = neighbors.size();
            for (size_t k=0; k<num_neighbors; k++) {
                if (voxel_counter.find(neighbors[k]) == voxel_counter.end()) {
                    voxel_counter[neighbors[k]] = 1;
                } else {
                    voxel_counter[neighbors[k]] += 1;
                }
            }
        }
        for (const auto& entry : voxel_counter) {
            if (entry.second == vertex_per_face) {
                voxel_neighbors.insert(entry.first);
            }
        }
    }