Method initialize_3D_grid

tools/Assembler/Materials/ElementWiseMaterial.cpp:77–92

}

void ElementWiseMaterial::initialize_3D_grid() {
    const size_t dim = m_material_mesh->get_dim();
    const size_t num_elem = m_material_mesh->get_num_voxels();
    const Float cell_size = compute_cell_size();

    m_grid = HashGrid::create(cell_size, 3);

    if (!m_material_mesh->has_attribute("voxel_centroid")) {
        m_material_mesh->add_attribute("voxel_centroid");
    }
    const VectorF& keys = m_material_mesh->get_attribute("voxel_centroid");

    for (size_t i=0; i<num_elem; i++) {
        m_grid->insert(i, keys.segment(i*dim, dim));
    }
}