BoundaryMesh

Question

Is there any function in Fenics that given a mesh on a domain and a relative boundary mesh using BoundaryMesh(mesh), maps the verteices number of the domain mesh to the vertices number of the Boundarymesh and viceversa?

Thank you!

Answer 1

Hi, consider

from dolfin import *

mesh = UnitSquareMesh(10, 10)
bmesh = BoundaryMesh(mesh, 'exterior')
# Map vertices of bmesh to vertices of mesh
vb_2_v = bmesh.entity_map(0)
# Check
assert all(near(Vertex(bmesh, vb_index).point().distance(Vertex(mesh, v_index).point()), 0)
           for vb_index, v_index in enumerate(vb_2_v))
# The reverse map
v_2_vb = {v: vb for vb, v in enumerate(vb_2_v)} 

Comments

Thank you very much!!