Evaluation of function at specific point

[torsten_wendav]

When trying to evaluate a function at a specific point WITHIN the domain, I get the error:

<Point x = 0.75 y = 0.25 z = 0.125> Traceback (most recent call last): File "issue_example.py", line 7, in <module> print v(x) File "[...]/FEniCS/lib/python2.7/site-packages/dolfin/functions/function.py", line 594, in call self.eval(values,x)

from dolfin import *
mesh = UnitCubeMesh(14, 14, 14)
V = FunctionSpace(mesh, "Lagrange", 1)
v = Function(V)
x = (0.75, 0.25, 0.125)
print v(x)

[Anders Logg (Chalmers)]

The problem is a round-off error. A geometric test fails with 7.972206336921708e-16 being slightly larger than 7.142857142857143e-16. This should be fixed now.

[Anders Logg (Chalmers)]

Fix issue #168: Evaluation of function at specific point

Using both absolute and relative check in numerical tests for collisions with intervals, triangles, tetrahedrons:

const double eps = std::max(DOLFIN_EPS_LARGE, DOLFIN_EPS_LARGE*std::max(dx, std::max(dy, dz)));

→ <<cset 114c5725d7b5>>

[Prof Garth Wells]

Removing milestone: 1.3 (automated comment)