This is a read only copy of the old FEniCS QA forum. Please visit the new QA forum to ask questions

evaluate function fails for point inside mesh

0 votes

I know this might be me doing something wrong here so first I'd like to see if anyone can reproduce

from dolfin import *
mesh = IntervalMesh(3,0.0,1000)
f = project(Constant(1.0),mesh=mesh)
print f([1.1])

I get this error although the point is inside the mesh and I think project uses CG 1 by default.

*** Error: Unable to evaluate function at point.
*** Reason: The point is not inside the domain. Consider calling "Function::set_allow_extrapolation(true)" on this Function to allow extrapolation.
*** Where: This error was encountered inside Function.cpp.

asked Sep 6, 2015 by chaffra FEniCS User (1,830 points)
edited Sep 6, 2015 by chaffra

1 Answer

0 votes

This looks like a bug in 1D BoundingBox::compute_first_entity_collision()

answered Sep 6, 2015 by chris_richardson FEniCS Expert (31,740 points)

Confirmed bug in CollisionDetection::collides_interval_point(). I'll post a fix.

commented Sep 6, 2015 by chris_richardson FEniCS Expert (31,740 points)

Ok I see it's missing the wnorm. Thanks.

commented Sep 6, 2015 by chaffra FEniCS User (1,830 points)
...