diff --git a/Geo/MHexahedron.cpp b/Geo/MHexahedron.cpp
index 38a392ac0f8151f6c733f50d9b50e7c65854c1d5..66cb52a2ec6a442b55001d494d7834a04fee16dc 100644
--- a/Geo/MHexahedron.cpp
+++ b/Geo/MHexahedron.cpp
@@ -124,8 +124,6 @@ const polynomialBasis* MHexahedron::getFunctionSpace(int o) const
 
   int nv = getNumVolumeVertices();
 
-  printf("nv = %d\n", nv);
-
   if ((nv == 0) && (o == -1)) {
     switch (order) {
     case 0: return polynomialBases::find(MSH_HEX_1);
@@ -146,7 +144,7 @@ const polynomialBasis* MHexahedron::getFunctionSpace(int o) const
     case 0: return polynomialBases::find(MSH_HEX_1);
     case 1: return polynomialBases::find(MSH_HEX_8);
     case 2: return polynomialBases::find(MSH_HEX_27);
-    case 3: printf("BBBBBBBBBBB\n"); return polynomialBases::find(MSH_HEX_64);
+    case 3: return polynomialBases::find(MSH_HEX_64);
     case 4: return polynomialBases::find(MSH_HEX_125);
     case 5: return polynomialBases::find(MSH_HEX_216);
     case 6: return polynomialBases::find(MSH_HEX_343);