| 93 | : size(size), elementCoordinates(elementCoordinates), elementGroups(elementGroups) {} |
| 94 | |
| 95 | CellToVertexArray |
| 96 | CellToVertexArray::fromMeshReader(const seissol::geometry::MeshReader& meshReader) { |
| 97 | const auto& elements = meshReader.getElements(); |
| 98 | const auto& vertices = meshReader.getVertices(); |
| 99 | |
| 100 | return CellToVertexArray( |
| 101 | elements.size(), |
| 102 | [&](size_t index) { |
| 103 | std::array<Eigen::Vector3d, 4> verts; |
| 104 | for (size_t i = 0; i < 4; ++i) { |
| 105 | auto vindex = elements[index].vertices[i]; |
| 106 | const auto& vertex = vertices[vindex]; |
| 107 | verts[i] << vertex.coords[0], vertex.coords[1], vertex.coords[2]; |
| 108 | } |
| 109 | return verts; |
| 110 | }, |
| 111 | [&](size_t index) { return elements[index].group; }); |
| 112 | } |
| 113 | |
| 114 | #ifdef USE_HDF |
| 115 | CellToVertexArray CellToVertexArray::fromPUML(const PUML::TETPUML& mesh) { |
nothing calls this directly
no test coverage detected