diff --git a/Geo/GFaceCompound.cpp b/Geo/GFaceCompound.cpp
index 348cb6e6fbe94ffb9b5985c00e44e1a8b760fd60..37b4bbc87841f82e5a7f42d81529fd7fe6869831 100644
--- a/Geo/GFaceCompound.cpp
+++ b/Geo/GFaceCompound.cpp
@@ -661,8 +661,8 @@ bool GFaceCompound::parametrize() const
   }
   // Radial-Basis Function parametrization
   else if (_mapping == RBF){
-    printf("Parametrizing surface %d with 'rbf' \n", tag());
-    Msg::Debug("Parametrizing surface %d with 'rbf's'", tag());
+    printf("Parametrizing surface %d with 'RBF' \n", tag());
+    Msg::Debug("Parametrizing surface %d with 'RBF''", tag());
   }
 
   buildOct();  
@@ -1727,7 +1727,7 @@ GPoint GFaceCompound::point(double par1, double par2) const
     //     b102 = (2*lt->v3 + lt->v1-w31*n3)*0.333;
     //     b201 = (2*lt->v1 + lt->v3-w13*n1)*0.333;
 
-    //tagged PN trinagles (sigma=1)
+    //tagged PN triangles (sigma=1)
     double theta = 0.0;
     SVector3 d1 = lt->v1+.33*(1-theta)*(lt->v2-lt->v1);
     SVector3 d2 = lt->v2+.33*(1-theta)*(lt->v1-lt->v2);
diff --git a/Geo/STensor3.h b/Geo/STensor3.h
index b554a7edcb230aae324a38932d4b38587f3e1607..ef83e911f78e3731225cf2150656665cc3236e30 100644
--- a/Geo/STensor3.h
+++ b/Geo/STensor3.h
@@ -101,6 +101,13 @@ class SMetric3 {
     ithis.setMat(m);
     return ithis;
   }
+  double determinant () const
+  {
+    fullMatrix<double> m(3,3);
+    getMat(m);
+    double det = m.determinant();
+    return det;
+  }
   SMetric3 operator + (const SMetric3 &other) const
   {
     SMetric3 res(*this);
@@ -229,6 +236,7 @@ class STensor3 {
     ithis.setMat(m);
     return ithis;
   }
+
   STensor3 operator + (const STensor3 &other) const
   {
     STensor3 res(*this);