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) I get the following error: 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