answered 2015-01-06 00:08:51 +0200
I guess
H = f.hessian().function(x,y) H(z[0,0],z[1,0])
would also work