diff --git a/FunctionSpace/TetNodeBasis.cpp b/FunctionSpace/TetNodeBasis.cpp
index bd66ecdff0d5ee4e54f9d8574c7461bc343843f3..71fa5d79236cd7ac93cb787eadace310cc1c279c 100644
--- a/FunctionSpace/TetNodeBasis.cpp
+++ b/FunctionSpace/TetNodeBasis.cpp
@@ -79,8 +79,8 @@ TetNodeBasis::TetNodeBasis(size_t order){
       for(size_t l = 1; l < order; l++){
         basis[s][i] =
           new Polynomial(intLegendre[l].compose
-                         (lagrange[edgeIdx[s][e][0]] -
-                          lagrange[edgeIdx[s][e][1]]
+                         (lagrange[edgeIdx[s][e][1]] -
+                          lagrange[edgeIdx[s][e][0]]
                           ,
                           lagrange[edgeIdx[s][e][0]] +
                           lagrange[edgeIdx[s][e][1]]));
@@ -89,9 +89,6 @@ TetNodeBasis::TetNodeBasis(size_t order){
     }
   }
 
-  // Face Based //
-  // TO CHECK: Are Triangles face matching tets ?
-
   for(size_t s = 0; s < nRefSpace; s++){
     size_t i = nVertex + nEdge;
 
@@ -105,8 +102,8 @@ TetNodeBasis::TetNodeBasis(size_t order){
 
           basis[s][i] =
             new Polynomial(intLegendre[l1].compose
-                           (lagrange[faceIdx[s][f][0]] -
-                            lagrange[faceIdx[s][f][1]]
+                           (lagrange[faceIdx[s][f][1]] -
+                            lagrange[faceIdx[s][f][0]]
                             ,
                             lagrange[faceIdx[s][f][0]] +
                             lagrange[faceIdx[s][f][1]])