Adding the storage of v_i v_j and being more careful with dx[:].

Instead of just storing v.v to compute the kinetic energy, each
element-wise v_i v_j is now stored in 1D projection for the computation
of things like the anisotropy.

grid.dx is now a vector to account for asymmetric elements.