diff --git a/Numeric/Gauss.cpp b/Numeric/Gauss.cpp
index 1f6e54b962b6ab780abbe15a0d5a810f09235b21..50fbb050db8d526b18d436b255943e97fd46341a 100644
--- a/Numeric/Gauss.cpp
+++ b/Numeric/Gauss.cpp
@@ -61,6 +61,11 @@ void gaussIntegration::get(int elementType, int order, fullMatrix<double> &pts,
     case TYPE_TET : pts2fullMatrix(getNGQTetPts(order), getGQTetPts(order), pts, weights); break;
     case TYPE_HEX : pts2fullMatrix(getNGQHPts(order), getGQHPts(order), pts, weights); break;
     case TYPE_PRI : pts2fullMatrix(getNGQPriPts(order), getGQPriPts(order), pts, weights); break;
+    case TYPE_PNT :
+      weights.resize(1,1);
+      weights(0,0) = 1.;
+      pts.resize(1,3);
+    break;
     default : Msg::Error("No integration rules defined for type %i", elementType);
   }
 }
diff --git a/Numeric/polynomialBasis.cpp b/Numeric/polynomialBasis.cpp
index 97d67876e45a9bf0410d1b5a27ee1c7f7432d64f..79d31aca83abeb8d822118abd92556f42adb1065 100644
--- a/Numeric/polynomialBasis.cpp
+++ b/Numeric/polynomialBasis.cpp
@@ -1032,7 +1032,7 @@ const polynomialBasis *polynomialBases::find(int tag)
     F.numFaces = 1;
     F.monomials = generate1DMonomials(0);
     F.points    = generate1DPoints(0);
-    F.parentType = MSH_PNT;
+    F.parentType = TYPE_PNT;
     break;
   case MSH_LIN_2 :
     F.numFaces = 2;