diff --git a/Numeric/polynomialBasis.cpp b/Numeric/polynomialBasis.cpp
index da15f93cde44b302de0595cfd0882b05de1fac0c..b1e751a81417953a18de7b21ef30ba1c535d98d3 100644
--- a/Numeric/polynomialBasis.cpp
+++ b/Numeric/polynomialBasis.cpp
@@ -1526,7 +1526,7 @@ const polynomialBasis *polynomialBases::find(int tag)
   case MSH_PRI_1 :
     F.numFaces = 5;
     F.monomials = generatePascalPrism(0);
-    F.points =    gmshGeneratePointsPrism(1, false);
+    F.points =    gmshGeneratePointsPrism(0, false);
     F.parentType = TYPE_PRI;
     generateClosureOrder0(F.closures,48);
     break;