ForkMaster/netgen
1
0
Fork 0
mirror of https://github.com/NGSolve/netgen.git synced 2025-01-12 14:10:34 +05:00
Code Issues Packages Projects Releases Wiki Activity

Merge remote-tracking branch 'bkapidani/autoPML_should_handle_default_spline_labels'

Browse Source
This commit is contained in:
Christopher Lackner 2019-12-11 10:46:28 +01:00
commit 260fd39da5
1 changed files with 8 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
python/geom2d.py
View File

@ -79,6 +79,8 @@ Returned is a dict with information to create the pml layer:
start_ndoms = ndoms = geo.GetNDomains() + 1 start_ndoms = ndoms = geo.GetNDomains() + 1
new_spline_domains = [] new_spline_domains = []
normals = {} normals = {}
duplicate_cnt = 0
for i, spline in enumerate(border): for i, spline in enumerate(border):
if i == 0: if i == 0:
global_start = Start(spline) + pml_size * spline.GetNormal(0) global_start = Start(spline) + pml_size * spline.GetNormal(0)
@ -86,6 +88,12 @@ Returned is a dict with information to create the pml layer:
next_spline = border[(i+1)%len(border)] next_spline = border[(i+1)%len(border)]
new_end = End(spline) + pml_size * spline.GetNormal(1) new_end = End(spline) + pml_size * spline.GetNormal(1)
spline_name = geo.GetBCName(spline.bc) spline_name = geo.GetBCName(spline.bc)
if "pml_" + spline_name in normals \
and normals["pml_" + spline_name] != spline.GetNormal(0):
duplicate_cnt += 1
spline_name = spline_name + "_duplicate_" + str(duplicate_cnt)
if (new_end - global_start).Norm() < tol: if (new_end - global_start).Norm() < tol:
new_spline_domains.append(ndoms) new_spline_domains.append(ndoms)
geo.Append(["line", current_start, global_start_pnt], bc="outer_" + spline_name, leftdomain = ndoms) geo.Append(["line", current_start, global_start_pnt], bc="outer_" + spline_name, leftdomain = ndoms)