Wrong result for TetrahedronCell::squared_distance

[Anders Logg (Chalmers)]

The following code print 0.0 while it should print 1.0:

from dolfin import *
p = Point(0.1, 0.05, -1.0)
mesh = UnitCubeMesh(8, 8, 8)
c = Cell(mesh, 0)
print c.distance(p)

[Anders Logg (Chalmers)]

Problem seems to be in TriangleCell::squared_distance which does not compute accurate results in 3D.

[Anders Logg (Chalmers)]

Fix issue 153: Wrong result for TetrahedronCell::squared_distance

→ <<cset 352ee11781ca>>

[Anders Logg (Chalmers)]

Fixed by modifying TriangleCell::squared_distance to handle 3D triangles

[Prof Garth Wells]

Removing milestone: 1.3 (automated comment)