diff --git a/Numeric/CondNumBasis.cpp b/Numeric/CondNumBasis.cpp
index 0b8d2063c34fc2e2d829893de59bf03253dfad71..cf46578ae835c1f96609b4cb66e8db53e113ece2 100644
--- a/Numeric/CondNumBasis.cpp
+++ b/Numeric/CondNumBasis.cpp
@@ -316,17 +316,16 @@ CondNumBasis::CondNumBasis(int tag, int cnOrder) :
     _condNumOrder(cnOrder >= 0 ? cnOrder : condNumOrder(tag))
 {
   const bool serendip = false;
-  FuncSpaceData data(true, tag, _condNumOrder, &serendip);
 
+  FuncSpaceData data = ( ElementType::ParentTypeFromTag(tag) == TYPE_PYR )? FuncSpaceData(true, tag, true, _condNumOrder+2, _condNumOrder, &serendip) : FuncSpaceData(true, tag, _condNumOrder, &serendip);
   _gradBasis = BasisFactory::getGradientBasis(data);
-
   fullMatrix<double> lagPoints;                                  // Sampling points
   gmshGeneratePoints(data, lagPoints);
   _nCondNumNodes = lagPoints.size1();
   _nMapNodes = BasisFactory::getNodalBasis(tag)->getNumShapeFunctions();
 
   // Store shape function gradients of mapping at condition number nodes
-  _gradBasis = BasisFactory::getGradientBasis(tag, _condNumOrder);
+  _gradBasis = BasisFactory::getGradientBasis(data);
 
   // Compute shape function gradients of primary mapping at barycenter,
   // in order to compute normal to straight element