696   double dJ = 0;  in _computeDeterminantAndRegularize()  local
700     dJ = 1.0;  in _computeDeterminantAndRegularize()
707     dJ = sqrt(SQU(jac[0][0]) + SQU(jac[0][1]) + SQU(jac[0][2]));  in _computeDeterminantAndRegularize()
737     dJ = sqrt(SQU(jac[0][0] * jac[1][1] - jac[0][1] * jac[1][0]) +  in _computeDeterminantAndRegularize()
757     dJ =  in _computeDeterminantAndRegularize()
764   return dJ;  in _computeDeterminantAndRegularize()
769   double dJ = 0;  in _computeDeterminantAndRegularize()  local
779     dJ = 1.0;  in _computeDeterminantAndRegularize()
786     dJ = sqrt(SQU(jac[0]) + SQU(jac[1]) + SQU(jac[2]));  in _computeDeterminantAndRegularize()
816     dJ = sqrt(SQU(jac[0] * jac[4] - jac[1] * jac[3]) +  in _computeDeterminantAndRegularize()
836     dJ = (jac[0] * jac[4] * jac[8] + jac[2] * jac[3] * jac[7] +  in _computeDeterminantAndRegularize()
842   return dJ;  in _computeDeterminantAndRegularize()