Dear all,
I would like to make an anisotropic mesh of spherical shells, with a well-resolved radial direction. I try to do this with the CloseSurfaces method, but I get Segmentation fault (core dumped) when I try. Below is an example code that results in this error. Any ideas?
All the best,
Magnus
[code]import netgen.csg as csg
import ngsolve as ng
origin = csg.Pnt(0,0,0)
inner = csg.Sphere(origin, 1.0)
outer = csg.Sphere(origin, 1.5)
shell = outer - inner
geo = csg.CSGeometry()
geo.Add(shell)
geo.CloseSurfaces(inner, outer, domain=shell)
mesh = geo.GenerateMesh(maxh=0.2)
csg.ZRefinement(mesh, geo)[/code]