mesh.domains().facet_domains() has gone away?

Question

In version 1.2 of DOLFIN there was such function that I could use. Now I moved to latest git revision and it seems there is not any more such function. Any idea what's wrong?

Also mesh.closest_cell() seems to be gone.

The mesh I am using has been converted by dolfin-convert from (netgen generated) diffpack to xml.

Answer 1

I'm not quite sure how mesh.domains().facet_domains() worked earlier, but the same data can be accessed from mesh.domains() and a mesh function can be created:

facet_domains = MeshFunction("size_t", mesh, D-1, mesh.domains())

The intersection code has been completely refurbished, and now use the class BoundingBoxTree:

mesh = UnitCubeMesh(4,4,4)
bbtree = BoundingBoxTree()
bbtree.build(mesh)
bbtree.compute_closest_entity(Point(2.0,0.0,0.0))