parouby
/
Pamhyr2
mirror of https://gitlab.com/pamhyr/pamhyr2
1
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

debug named points in meshed profiles

master
Theophile Terraz 2025-11-28 14:12:06 +01:00
parent 0f0defb36f
commit 64c99d46d4
1 changed files with 9 additions and 18 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

27
src/Meshing/Internal.py
View File

@ -120,6 +120,7 @@ class InternalMeshing(AMeshingTool):
ptr += 1 # next profile, original ptr += 1 # next profile, original
for k in range(len(reach.profiles[ptr0].points)): for k in range(len(reach.profiles[ptr0].points)):
p.add() p.add()
p.points[k].name = reach.profiles[ptr0].points[k].name
p.points[k].x = reach.profiles[ptr0].points[k].x + \ p.points[k].x = reach.profiles[ptr0].points[k].x + \
dj*(reach.profiles[ptr].points[k].x - dj*(reach.profiles[ptr].points[k].x -
reach.profiles[ptr0].points[k].x) reach.profiles[ptr0].points[k].x)
@ -144,19 +145,11 @@ class InternalMeshing(AMeshingTool):
# len1 and 2: number of intervals # len1 and 2: number of intervals
# len1: target len # len1: target len
if tag1 == '': # left end point start1 = sect1.named_point_index(tag1)
start1 = 0 start2 = sect2.named_point_index(tag1)
start2 = 0
else:
start1 = sect1.named_point_index(tag1)
start2 = sect2.named_point_index(tag1)
if tag2 == '': # right end point end1 = sect1.named_point_index(tag2)
end1 = sect1.nb_points-1 end2 = sect2.named_point_index(tag2)
end2 = sect2.nb_points-1
else:
end1 = sect1.named_point_index(tag2)
end2 = sect2.named_point_index(tag2)
len1 = end1 - start1 len1 = end1 - start1
len2 = end2 - start2 len2 = end2 - start2
@ -192,10 +185,6 @@ class InternalMeshing(AMeshingTool):
for i in range(len1-len2): for i in range(len1-len2):
p = sect2.point(start2).copy() p = sect2.point(start2).copy()
sect2.insert_point(start2, p) sect2.insert_point(start2, p)
if tag1 == '': # left end point
sect2.point(start2).name = ''
if tag2 == '': # left end point
sect2.point(start2+1).name = ''
len2 += 1 len2 += 1
elif ltot1 < 0.0001: elif ltot1 < 0.0001:
sect2.add_npoints(len1-len2) sect2.add_npoints(len1-len2)
@ -259,8 +248,10 @@ class InternalMeshing(AMeshingTool):
p = sect2.point(start2+len2).copy() p = sect2.point(start2+len2).copy()
sect2.insert_point(start2+len2, p) sect2.insert_point(start2+len2, p)
sect2.point(start2+len2).name = '' sect2.point(start2+len2).name = ''
sect2.point(start2).name = tag1 if tag1 != "un":
sect2.point(start2+len2).name = tag2 sect2.point(start2).name = tag1
if tag2 != "np":
sect2.point(start2+len2).name = tag2
sect2.modified() sect2.modified()
def update_rk(self, reach, def update_rk(self, reach,