intersect() returns invalid cell indices
Issue #811
new
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]
Comments (2)
-
-
- removed milestone
- Log in to comment
Something weird is going on. The
intersect
function does the followingfrom which the output is
Something is likely wrong in the SWIG layer.