Method project_out_zero_dim

src/Geometry/MeshGeometry.cpp:26–70

}

int MeshGeometry::project_out_zero_dim() {
    if (m_dim == 2) return -1;
    assert(m_dim == 3);
    const size_t num_vertices = get_num_vertices();
    MatrixFr vertices(num_vertices, m_dim);
    std::copy(m_vertices.data(), m_vertices.data() + m_vertices.size(),
            vertices.data());
    Vector3F max_coord = vertices.colwise().maxCoeff();
    Vector3F min_coord = vertices.colwise().minCoeff();

    int zero_dim = -1;
    for (size_t i=0; i<3; i++) {
        if (max_coord[i] == min_coord[i]) {
            zero_dim = i;
            break;
        }
    }

    if (zero_dim == -1) return -1;
    m_dim = 2;
    m_vertices.resize(num_vertices * 2);
    if (zero_dim == 0) {
        // Drop X
        //std::cout << "Removing X component" << std::endl;
        for (size_t i=0; i<num_vertices; i++) {
            m_vertices[i*2  ] = vertices.row(i)[1];
            m_vertices[i*2+1] = vertices.row(i)[2];
        }
    } else if (zero_dim == 2) {
        // Drop Z
        //std::cout << "Removing Z component" << std::endl;
        for (size_t i=0; i<num_vertices; i++) {
            m_vertices[i*2  ] = vertices.row(i)[0];
            m_vertices[i*2+1] = vertices.row(i)[1];
        }
    } else {
        // Drop Y
        //std::cout << "Removing Y component" << std::endl;
        for (size_t i=0; i<num_vertices; i++) {
            m_vertices[i*2  ] = vertices.row(i)[0];
            m_vertices[i*2+1] = vertices.row(i)[2];
        }
    }
    return zero_dim;
}

void MeshGeometry::extract_faces_from_tets() {
    const VectorI& voxels = m_voxels;