intersect() returns invalid cell indices

[Nathan Sime]

The following code yields invalid cell indices. Should intersect be removed/updated?

from dolfin import *

mesh = UnitSquareMesh(32, 32)
for j in range(5):
    c = intersect(mesh, Point(0.0, 0.0)).intersected_cells()
    print c

output varies, for example:

[2 0]
[        0 536870912]
[0 1]
[        2 536870912]
[2818965428  307092016]

[Anders Logg (Chalmers)]

Something weird is going on. The intersect function does the following

tree = BoundingBoxTree()
tree.build(mesh)
tree.compute_entity_collisions(Point(0, 0))

from which the output is

array([1, 0], dtype=uint32)

Something is likely wrong in the SWIG layer.