Dear all,

When trying to evaluate a function u on some mesh points with

mesh = Mesh('coarse/Omega_coarse.msh')

V = FunctionSpace(mesh, "CG",1)

V_Vec = VectorFunctionSpace(mesh, "CG", 1)

u = Function(V)

points = Function(V_vec).interpolate(SpatialCoordinate(mesh))

dat = u.at(points.dat.data_ro)

File "lin_dyn_unicell.py", line 117, in <module>

    dat = u.at(points.dat.data_ro)

  File "/Users/gidel/Documents/Firedrake/firedrake/src/firedrake/firedrake/function.py", line 585, in at

    raise PointNotInDomainError(self.function_space().mesh(), points[i].reshape(-1))

firedrake.function.PointNotInDomainError: domain <Mesh #1> does not contain point [200.          66.66666667]

How is that possible while both function spaces are defined on the same mesh? 

Note that in this example I am using the same domain to make sure the points are the same in the two cases, but I used to run this with different meshes and it worked on my laptop with an old firedrake version but not on my new machine on which I installed Firedrake yesterday. 

Thank you for your feedback,

Floriane