How do I evaluate a grad of a solution on multiple points?

2

OK. That was fast and easy once I wrote the question.

I just used

gu=grad(u)
Gu=dolfin.project(gu)
Gu(Point(0,0))