Something's wrong with new geometry algorithms in 1D

[Jan Blechta]

from dolfin import *
parameters['allow_extrapolation'] = True
while True:
    mesh = IntervalMesh(70, 0.06, 0.41)
    V = FunctionSpace(mesh, 'CG', 1)
    u = Function(V)
    u(0.02)

After few iterations this raises

*** Error:   Unable to complete call to function compute_closest_entity().
*** Reason:  Assertion closest_entity < std::numeric_limits<unsigned int>::max() failed.
*** Where:   This error was encountered inside /usr/users/blechta/fenics/fenics/src/dolfin/dolfin/geometry/GenericBoundingBoxTree.cpp (line 181).

but works fine with 1.2.0.

[Prof Garth Wells]

I have seen this occasionally in 1D too. Perhaps your report should make clear that your test case is1D.

[Jan Blechta]

Use IntervalMesh instead of mesh file in the snippet.

[Jan Blechta]

Fix in https://bitbucket.org/fenics-project/dolfin/pull-request/46/fix-bugs-in-boundingboxtree1d-and/diff

[Prof Garth Wells]

Fix committed.