Segfault when creating 1D mesh

from netgen.meshing import Element0D, Element1D, Element2D, MeshPoint, Pnt from netgen.meshing import Mesh as NetGenMesh import ngsolve N = 10 ngmesh = NetGenMesh(dim=1) pids = [] # Add points to the mesh. for i in range(N + 1): pids.append(ngmesh.Add(MeshPoint(Pnt(i / N, 0, 0)))) # Left half of the domain is material 1, right half material 2. for i in range(N): ngmesh.Add(Element1D([pids[i], pids[i + 1]], index=1 if i <= 5 else 2)) ngmesh.SetMaterial(1, "domain_left") ngmesh.SetMaterial(2, "domain_right") # Add BC to the mesh. ngmesh.Add(Element0D(pids[0], index=1)) # Added index is 1 ngmesh.Add(Element0D(pids[N], index=2)) ngmesh.SetBCName(0, "l") # Index is retrieved as 0 ngmesh.SetBCName(1, "r") # Create NGSolve mesh and refine it once uniformly. mesh = ngsolve.Mesh(ngmesh) for el in mesh.Elements(): mesh.SetRefinementFlag(el, True) mesh.Refine()

$ netgen mwe.py NETGEN-6.2-dev Developed by Joachim Schoeberl at 2010-xxxx Vienna University of Technology 2006-2010 RWTH Aachen University 1996-2006 Johannes Kepler University Linz optfile ./ng.opt does not exist - using default values togl-version : 2 loading ngsolve library NGSolve-6.2.1907 Using Lapack Including sparse direct solver UMFPACK Running parallel using 8 thread(s) importing NGSolve-6.2.1907 (should) load python file 'mwe.py' Mesh bisection Segmentation fault: 11

idx_doml = ngmesh.AddRegion("domain_left", dim=1) idx_domr = ngmesh.AddRegion("domain_right", dim=1) for i in range(N): ngmesh.Add(Element1D([pids[i], pids[i + 1]], index=idx_doml if i < 5 else idx_domr))