Dirichlet BC for "geometric" search broken

[Anders Logg (Chalmers)]

The "geometric" search for BC nodes will fail to locate nodes on the boundary when the midpoint is not on the boundary. This is because it calls init_facets() is called inside compute_bc_geometric. That function in turn will call SubDomain::mark() which checks for the midpoint.

[Anders Logg (Chalmers)]

Fixed