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

Avoiding more numpy 2.0 deprecation warnings (#3049)

Browse source
This commit is contained in:
Ethan Peterson 2024-06-20 11:14:55 -07:00 • committed by GitHub
parent 78ee851990
commit 4bd0b09e60
No known key found for this signature in database
GPG key ID: B5690EEEBB952194
1 changed files with 8 additions and 8 deletions
Download patch file Download diff file Expand all files Collapse all files

16
openmc/model/surface_composite.py
View file

@ -972,19 +972,19 @@ class Polygon(CompositeSurface):
# Check if polygon is self-intersecting by comparing edges pairwise
n = len(points)
for i in range(n):
p0 = points[i, :]
p1 = points[(i + 1) % n, :]
p0 = np.append(points[i, :], 0)
p1 = np.append(points[(i + 1) % n, :], 0)
for j in range(i + 1, n):
p2 = points[j, :]
p3 = points[(j + 1) % n, :]
p2 = np.append(points[j, :], 0)
p3 = np.append(points[(j + 1) % n, :], 0)
# Compute orientation of p0 wrt p2->p3 line segment
cp0 = np.cross(p3-p0, p2-p0)
cp0 = np.cross(p3-p0, p2-p0)[-1]
# Compute orientation of p1 wrt p2->p3 line segment
cp1 = np.cross(p3-p1, p2-p1)
cp1 = np.cross(p3-p1, p2-p1)[-1]
# Compute orientation of p2 wrt p0->p1 line segment
cp2 = np.cross(p1-p2, p0-p2)
cp2 = np.cross(p1-p2, p0-p2)[-1]
# Compute orientation of p3 wrt p0->p1 line segment
cp3 = np.cross(p1-p3, p0-p3)
cp3 = np.cross(p1-p3, p0-p3)[-1]
# Group cross products in an array and find out how many are 0
cross_products = np.array([[cp0, cp1], [cp2, cp3]])