Hm ...
what do you assume from you mesh and your space ?
with structured meshes and Nodal-P1 you can take your coefficient vector and fead it into Python - FFT
otherwise, you have to evaluate your GridFunction on a regular grid of points (using geometric search) first.
Joachim