| 38 | } |
| 39 | |
| 40 | Elements::Ptr Elements::adapt_boundary(Mesh::Ptr mesh) { |
| 41 | if (mesh->get_num_voxels() > 0) { |
| 42 | const size_t vertex_per_voxel = mesh->get_vertex_per_voxel(); |
| 43 | switch (vertex_per_voxel) { |
| 44 | case 4: |
| 45 | return Ptr(new TriangleElements(mesh)); |
| 46 | default: |
| 47 | std::stringstream err_msg; |
| 48 | err_msg << "Voxel with " << vertex_per_voxel |
| 49 | << " vertices is not supported yet."; |
| 50 | throw NotImplementedError(err_msg.str()); |
| 51 | } |
| 52 | } else { |
| 53 | const size_t vertex_per_face = mesh->get_vertex_per_face(); |
| 54 | switch (vertex_per_face) { |
| 55 | case 3: |
| 56 | return Ptr(new EdgeElements(mesh)); |
| 57 | default: |
| 58 | std::stringstream err_msg; |
| 59 | err_msg << "Face with " << vertex_per_face |
| 60 | << " vertices is not supported yet."; |
| 61 | throw NotImplementedError(err_msg.str()); |
| 62 | } |
| 63 | } |
| 64 | } |
nothing calls this directly
no test coverage detected