From 6b30ec0b7c3430e89416d5cf2f34269ac09b2953 Mon Sep 17 00:00:00 2001 From: Christopher Lackner Date: Thu, 19 Nov 2020 21:54:58 +0100 Subject: [PATCH] test case for leftdom==rightdom in geom2d --- tests/pytest/test_geom2d.py | 13 +++++++++++++ 1 file changed, 13 insertions(+) create mode 100644 tests/pytest/test_geom2d.py diff --git a/tests/pytest/test_geom2d.py b/tests/pytest/test_geom2d.py new file mode 100644 index 00000000..d40fd40a --- /dev/null +++ b/tests/pytest/test_geom2d.py @@ -0,0 +1,13 @@ +from netgen.geom2d import SplineGeometry + + +def test_leftdom_equals_rightdom(): + geo = SplineGeometry() + pnts = [(0,0), (1,0), (2,0), (2,1), (1,1), (0,1)] + gp = [geo.AppendPoint(*p) for p in pnts] + lines = [(0,1,0), (1,2,0), (2,3,0), (3,4,0), (4,5,0), (5,0,0), (1,4,1)] + for p1, p2, rd in lines: + geo.Append(["line", p1, p2], leftdomain=1, rightdomain=rd) + + mesh = geo.GenerateMesh() +