I’m sorry I should’ve been more specific; it does indeed work, just as
value = u(x,y)
but I need a fast way to evaluate u at many points. I really don’t want to loop over the coordinates in plain python and evaluate u pointwise, that would just take way too much time.