Method compute_from_mesh

src/Attributes/VoxelVolumeAttribute.cpp:31–61

using namespace VoxelVolumeAttributeHelper;

void VoxelVolumeAttribute::compute_from_mesh(Mesh& mesh) {
    const size_t num_voxels = mesh.get_num_voxels();
    const size_t num_vertex_per_voxel = mesh.get_vertex_per_voxel();

    VectorF& volumes = m_values;
    volumes = VectorF::Zero(num_voxels);

    if (num_voxels > 0) {
        if (num_vertex_per_voxel == 4) {
            tbb::parallel_for(tbb::blocked_range<size_t>(0, num_voxels),
                [this, &mesh, &volumes](const tbb::blocked_range<size_t>& r) {
                    for (size_t i=r.begin(); i<r.end(); i++) {
                        volumes[i] = compute_signed_tet_volume(mesh, i);
                    }
                }
            );
        } else if (num_vertex_per_voxel == 8) {
            tbb::parallel_for(tbb::blocked_range<size_t>(0, num_voxels),
                [this, &mesh, &volumes](const tbb::blocked_range<size_t>& r) {
                    for (size_t i=r.begin(); i<r.end(); i++) {
                        volumes[i] = compute_signed_hex_volume(mesh, i);
                    }
                }
            );
        } else {
            std::cerr << "Unknown voxel type with " << num_vertex_per_voxel
                << " per voxel." << std::endl;
            return;
        }
    }
}

Float VoxelVolumeAttribute::compute_signed_tet_volume(Mesh& mesh, size_t voxel_idx) {
    VectorI voxel = mesh.get_voxel(voxel_idx);