timing the vortex computation, where the bottleneck is the Jacobian computation
jacobian (unstable)
backup
1. fix a bug in at_vertex() when referencing the vertex data given (i,j,k).  2. remove UnitJacobian since the way to do normalization may be wrong.  3. VTK integration now supports at_vert() and phy_coord()