diff --git a/Numeric/JacobianBasis.cpp b/Numeric/JacobianBasis.cpp
index 828b8973c9971f2f5b48ae4bfece05966475cb75..a35231fd4afaa2b45fd548a994eca8011fec2656 100644
--- a/Numeric/JacobianBasis.cpp
+++ b/Numeric/JacobianBasis.cpp
@@ -376,9 +376,7 @@ void GradientBasis::mapFromIdealElement(fullMatrix<double> *dxyzdX,
 
     case TYPE_TET:
     {
-      static const double cTet[3] = {-3./2/std::sqrt(6),
-                                     -1./2/std::sqrt(2),
-                                     std::sqrt(1.5)};
+      static const double cTet[3] = {-1.5/sqrt(6), -0.75/sqrt(2), 3./sqrt(6)};
       dxyzdZ->scale(cTet[2]);
       dxyzdZ->axpy(*dxyzdX, cTet[0]);
       dxyzdZ->axpy(*dxyzdY, cTet[1]);