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

fix not working boundarylayer test

Browse Source
This commit is contained in:
Christopher Lackner 2020-10-27 16:53:07 +01:00
parent 47632a06aa
commit d40c05b1b1
1 changed files with 8 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
tests/pytest/test_boundarylayer.py
View File

@ -2,11 +2,15 @@
import pytest
from netgen.csg import *
def GetNSurfaceElements(mesh, boundaries):
def GetNSurfaceElements(mesh, boundaries, adjacent_domain=None):
nse_in_layer = 0
for el in mesh.Elements2D():
if mesh.GetBCName(el.index-1) in boundaries:
nse_in_layer += 1
if adjacent_domain is None:
nse_in_layer += 1
else:
if (mesh.FaceDescriptor(el.index).domin > 0 and mesh.GetMaterial(mesh.FaceDescriptor(el.index).domin) == adjacent_domain) or (mesh.FaceDescriptor(el.index).domout > 0 and mesh.GetMaterial(mesh.FaceDescriptor(el.index).domout) == adjacent_domain):
nse_in_layer += 1
return nse_in_layer
@pytest.mark.parametrize("outside", [True, False])
@ -29,7 +33,7 @@ def test_boundarylayer(outside, capfd):
assert not "elements are not matching" in capture.out
@pytest.mark.parametrize("outside", [True, False])
@pytest.mark.parametrize("version", [1]) # version 2 not working yet
@pytest.mark.parametrize("version", [1, 2]) # version 2 not working yet
def test_boundarylayer2(outside, version, capfd):
geo = CSGeometry()
top = Plane(Pnt(0,0,0.5), Vec(0,0,1))
@ -48,7 +52,7 @@ def test_boundarylayer2(outside, version, capfd):
geo.CloseSurfaces(top, bot, [])
mesh = geo.GenerateMesh()
should_ne = mesh.ne + 2 * GetNSurfaceElements(mesh, ["default"])
should_ne = mesh.ne + 2 * GetNSurfaceElements(mesh, ["default"], "part")
layersize = 0.05
mesh.BoundaryLayer("default", [0.5 * layersize, layersize], "layer", domains="part", outside=outside, grow_edges=True)
assert mesh.ne == should_ne