Method compute_from_hex_mesh

src/Attributes/VertexValanceAttribute.cpp:89–144

}

void VertexValanceAttribute::compute_from_hex_mesh(Mesh& mesh) {
    const size_t num_vertices = mesh.get_num_vertices();
    const size_t num_voxels = mesh.get_num_voxels();
    const size_t num_vertex_per_voxel = mesh.get_vertex_per_voxel();
    assert(num_vertex_per_voxel == 8);
    // Vertex ordering
    //  3 _________ 2
    //   /:       /|
    // 7/ :      / |
    // +--------+6 |
    // |  :     |  |
    // | 0:.....|..|1
    // | ;      | /
    // +--------+/
    // 4         5

    VectorF& vertex_valance = m_values;
    vertex_valance = VectorF::Zero(num_vertices);

    std::set<Duplet> edges;
    for (size_t i=0; i<num_voxels; i++) {
        VectorI voxel = mesh.get_voxel(i);
        Duplet edge[12] = {
            {voxel[0], voxel[1]},
            {voxel[1], voxel[2]},
            {voxel[2], voxel[3]},
            {voxel[3], voxel[0]},
            {voxel[4], voxel[5]},
            {voxel[5], voxel[6]},
            {voxel[6], voxel[7]},
            {voxel[7], voxel[4]},
            {voxel[0], voxel[4]},
            {voxel[1], voxel[5]},
            {voxel[2], voxel[6]},
            {voxel[3], voxel[7]} };

        edges.insert(edge[0]);
        edges.insert(edge[1]);
        edges.insert(edge[2]);
        edges.insert(edge[3]);
        edges.insert(edge[4]);
        edges.insert(edge[5]);
        edges.insert(edge[6]);
        edges.insert(edge[7]);
        edges.insert(edge[8]);
        edges.insert(edge[9]);
        edges.insert(edge[10]);
        edges.insert(edge[11]);
    }

    for (auto edge : edges) {
        const auto& data = edge.get_ori_data();
        vertex_valance[data[0]] ++;
        vertex_valance[data[1]] ++;
    }
}