Method compute_from_mesh

src/Attributes/VoxelCircumRadiusAttribute.cpp:9–38

using namespace PyMesh;

void VoxelCircumRadiusAttribute::compute_from_mesh(Mesh& mesh) {
    const size_t dim = mesh.get_dim();
    const size_t num_voxels = mesh.get_num_voxels();
    const size_t vertex_per_voxel = mesh.get_vertex_per_voxel();

    if (num_voxels > 0 && vertex_per_voxel != 4) {
        throw NotImplementedError(
                "Voxel circumradius is defined for 3D tet mesh only.");
    }
    if (dim != 3) {
        throw RuntimeError("Voxel circumradius is defined for 3D tet mesh only.");
    }

    if (!mesh.has_attribute("voxel_circumcenter")) {
        mesh.add_attribute("voxel_circumcenter");
    }
    const auto& circumcenter = mesh.get_attribute("voxel_circumcenter");

    VectorF& circumradius = m_values;
    circumradius.resize(num_voxels);
    const auto& voxels = mesh.get_voxels();
    const auto& vertices = mesh.get_vertices();

    for (size_t i=0; i<num_voxels; i++) {
        Vector4I voxel = voxels.segment<4>(i*4);
        Vector3F v0 = vertices.segment<3>(voxel[0]*3);
        Vector3F center = circumcenter.segment<3>(i*3);
        circumradius[i] = (v0 - center).norm();
    }
}