I have a spherical geometry consisting of of two concentric spheres where the smaller sphere is cut out from the larger sphere:
so = mshr.Sphere(dolfin.Point(0, 0, 0), router)
si = mshr.Sphere(dolfin.Point(0, 0, 0), rinner)
s = so - si
mesh = mshr.generate_mesh(s, resolution)
When I now generate the boundarymesh of mesh
using
dolfin.BoundaryMesh(mesh, 'exterior', True)
I get a mesh of the surface at router
and the one at rinner
. How can I remove the one at rinner
?