grad(sin(pi*q(x,y,z))/sinh(pi*q(x,y,z)))