| 34 | } |
| 35 | |
| 36 | TilerEngine::FuncList MeshTilerHelper::get_3D_tiling_operators(Mesh::Ptr mesh) { |
| 37 | const size_t num_cells = mesh->get_num_voxels(); |
| 38 | const size_t num_vertex_per_cell = mesh->get_vertex_per_voxel(); |
| 39 | if (num_vertex_per_cell != 8) { |
| 40 | throw NotImplementedError("Only hex guide mesh is supported in 3D"); |
| 41 | } |
| 42 | |
| 43 | TilerEngine::FuncList operators; |
| 44 | for (size_t i=0; i<num_cells; i++) { |
| 45 | const VectorI voxel = mesh->get_voxel(i); |
| 46 | MatrixFr corners(num_vertex_per_cell, 3); |
| 47 | for (size_t j=0; j<num_vertex_per_cell; j++) { |
| 48 | corners.row(j) = mesh->get_vertex(voxel[j]).transpose(); |
| 49 | } |
| 50 | |
| 51 | operators.push_back( |
| 52 | [=](const MatrixFr& vertices) { |
| 53 | TrilinearInterpolation interpolator(corners); |
| 54 | return interpolator.interpolate_batch(vertices); |
| 55 | }); |
| 56 | } |
| 57 | return operators; |
| 58 | } |
| 59 | |
| 60 | std::vector<ParameterCommon::Variables> MeshTilerHelper::extract_face_attributes( |
| 61 | Mesh::Ptr mesh) { |
nothing calls this directly
no test coverage detected