crp/OpenMC
1
0
Fork 0
mirror of https://github.com/openmc-dev/openmc.git synced 2026-07-27 13:45:36 -04:00
Code Issues Projects Releases Packages Wiki Activity Actions

moved coordinate axis to 0

Browse source
This commit is contained in:
Ethan Peterson 2022-11-18 11:56:27 -05:00
parent da7e06dc53
commit 5a716f1f3a
1 changed files with 6 additions and 5 deletions
Download patch file Download diff file Expand all files Collapse all files

11
openmc/mesh.py
View file

@ -172,10 +172,10 @@ class StructuredMesh(MeshBase):
-------
vertices : numpy.ndarray
Returns a numpy.ndarray representing the coordinates of the mesh
vertices with a shape equal to (dim1 + 1, ..., dimn + 1, ndim).
vertices with a shape equal to (ndim, dim1 + 1, ..., dimn + 1).
"""
return np.stack(np.meshgrid(*self._grids, indexing='ij'), axis=-1)
return np.stack(np.meshgrid(*self._grids, indexing='ij'), axis=0)
@property
def centroids(self):
@ -185,13 +185,14 @@ class StructuredMesh(MeshBase):
-------
centroids : numpy.ndarray
Returns a numpy.ndarray representing the mesh element centroid
coordinates with a shape equal to (dim1, ..., dimn, ndim).
coordinates with a shape equal to (ndim, dim1, ..., dimn). Can be
unpacked by the first dimension with xx, yy, zz = mesh.centroids
"""
ndim = self.n_dimension
vertices = self.vertices
s0 = (slice(0, -1),)*ndim + (slice(None),)
s1 = (slice(1, None),)*ndim + (slice(None),)
s0 = (slice(None),) + (slice(0, -1),)*ndim
s1 = (slice(None),) + (slice(1, None),)*ndim
return (vertices[s0] + vertices[s1]) / 2
@property