1 2 3
#define SPATIAL_DETERMINANT(gxx,gxy,gxz,gyy,gyz,gzz) \ (-(gxz)**2*(gyy) + 2*(gxy)*(gxz)*(gyz) - (gxx)*(gyz)**2 - (gxy)**2*(gzz) \ + (gxx)*(gyy)*(gzz))