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

check for duplicate points

Browse source
This commit is contained in:
Ethan Peterson 2023-01-30 19:40:00 -05:00
parent 4fade673a7
commit e20e648690
1 changed files with 4 additions and 1 deletions
Download patch file Download diff file Expand all files Collapse all files

5
openmc/model/surface_composite.py
View file

@ -773,13 +773,16 @@ class Polygon(CompositeSurface):
points = points[:-1, :]
check_length('points', points, 3)
if len(points) != len(np.unique(points, axis=0)):
raise ValueError('Duplicate points were detected in the Polygon input')
# Check if polygon is self-intersecting by comparing edges pairwise
n = len(points)
is_self_intersecting = False
for i in range(n):
p1x, p1y = points[i, :]
p2x, p2y = points[(i + 1) % n, :]
for j in range(i, n):
for j in range(i + 1, n):
p3x, p3y = points[j, :]
p4x, p4y = points[(j + 1) % n, :]