diff --git a/Numeric/polynomialBasis.h b/Numeric/polynomialBasis.h
index 6a05e2add178dc40a97d0dda57ab082c0d49299c..c1481612e0693ac42ecc342cf671016d8fd186ed 100644
--- a/Numeric/polynomialBasis.h
+++ b/Numeric/polynomialBasis.h
@@ -97,6 +97,10 @@ class polynomialBasis
   {
     return closures[id];
   }
+  inline const std::vector<int> &getFullClosure(int id) const
+  {
+    return fullClosures[id];
+  }
   inline int getClosureId(int iEl, int iSign=1, int iRot=0) const
   {
     return iEl + numFaces*(iSign == 1 ? 0 : 1) + 2*numFaces*iRot;