fake 0D gauss rule

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);
   }
 }

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;