swapped order or elements

shimwell · shimwell · commit 799098bc9a85 · 2026-04-02T11:55:39.000+02:00
diff --git a/openmc/mesh.py b/openmc/mesh.py
@@ -1839,7 +1839,7 @@ def _write_vtk_hdf5(self, filename, datasets, volume_normalization):
             )
             root.create_dataset("Dimensions", data=vertex_dims, dtype="i8")
 
-            points = vertices.reshape(-1, 3)
+            points = np.swapaxes(vertices, 0, 2).reshape(-1, 3)
             root.create_dataset("Points", data=points.astype(np.float64), dtype="f8")
 
             cell_data_group = root.create_group("CellData")
@@ -2344,7 +2344,7 @@ def _write_vtk_hdf5(self, filename, datasets, volume_normalization):
             )
             root.create_dataset("Dimensions", data=vertex_dims, dtype="i8")
 
-            points = vertices.reshape(-1, 3)
+            points = np.swapaxes(vertices, 0, 2).reshape(-1, 3)
             root.create_dataset("Points", data=points.astype(np.float64), dtype="f8")
 
             cell_data_group = root.create_group("CellData")
@@ -2785,7 +2785,7 @@ def _write_vtk_hdf5(self, filename, datasets, volume_normalization):
             )
             root.create_dataset("Dimensions", data=vertex_dims, dtype="i8")
 
-            points = vertices.reshape(-1, 3)
+            points = np.swapaxes(vertices, 0, 2).reshape(-1, 3)
             root.create_dataset("Points", data=points.astype(np.float64), dtype="f8")
 
             cell_data_group = root.create_group("CellData")

Original file line number	Diff line number	Diff line change
`@@ -1839,7 +1839,7 @@ def _write_vtk_hdf5(self, filename, datasets, volume_normalization):`
`1839`	`1839`	`)`
`1840`	`1840`	`root.create_dataset("Dimensions", data=vertex_dims, dtype="i8")`
`1841`	`1841`
`1842`		`- points = vertices.reshape(-1, 3)`
	`1842`	`+ points = np.swapaxes(vertices, 0, 2).reshape(-1, 3)`
`1843`	`1843`	`root.create_dataset("Points", data=points.astype(np.float64), dtype="f8")`
`1844`	`1844`
`1845`	`1845`	`cell_data_group = root.create_group("CellData")`
`@@ -2344,7 +2344,7 @@ def _write_vtk_hdf5(self, filename, datasets, volume_normalization):`
`2344`	`2344`	`)`
`2345`	`2345`	`root.create_dataset("Dimensions", data=vertex_dims, dtype="i8")`
`2346`	`2346`
`2347`		`- points = vertices.reshape(-1, 3)`
	`2347`	`+ points = np.swapaxes(vertices, 0, 2).reshape(-1, 3)`
`2348`	`2348`	`root.create_dataset("Points", data=points.astype(np.float64), dtype="f8")`
`2349`	`2349`
`2350`	`2350`	`cell_data_group = root.create_group("CellData")`
`@@ -2785,7 +2785,7 @@ def _write_vtk_hdf5(self, filename, datasets, volume_normalization):`
`2785`	`2785`	`)`
`2786`	`2786`	`root.create_dataset("Dimensions", data=vertex_dims, dtype="i8")`
`2787`	`2787`
`2788`		`- points = vertices.reshape(-1, 3)`
	`2788`	`+ points = np.swapaxes(vertices, 0, 2).reshape(-1, 3)`
`2789`	`2789`	`root.create_dataset("Points", data=points.astype(np.float64), dtype="f8")`
`2790`	`2790`
`2791`	`2791`	`cell_data_group = root.create_group("CellData")`