Fix: 2D Cells are now *FACES* (correct Dof for 2D - 3D interfaces) -- Nedelec...

Fix: 2D Cells are now *FACES* (correct Dof for 2D - 3D interfaces) -- Nedelec Basis are now of order 0

FunctionSpace/LineNedelecBasis.cpp

@@ -10,7 +10,7 @@ LineNedelecBasis::LineNedelecBasis(void){
   nRefSpace = refSpace->getNPermutation();
 
   // Set Basis Type //
-  order = 1;
+  order = 0;
 
   type = 1;
   dim  = 1;

FunctionSpace/QuadEdgeBasis.cpp

@@ -12,8 +12,8 @@ QuadEdgeBasis::QuadEdgeBasis(int order){
 
   nVertex = 0;
   nEdge   = 4 * (order + 1);
-  nFace   = 0;
-  nCell   = 2 * (order + 1) * order;
+  nFace   = 2 * (order + 1) * order;
+  nCell   = 0;
 
   nEdgeClosure = 2;
   nFaceClosure = 0;

FunctionSpace/QuadNodeBasis.cpp

@@ -12,8 +12,8 @@ QuadNodeBasis::QuadNodeBasis(int order){
 
   nVertex = 4;
   nEdge   = 4 * (order - 1);
-  nFace   = 0;
-  nCell   =     (order - 1) * (order - 1);
+  nFace   =     (order - 1) * (order - 1);
+  nCell   = 0;
 
   nEdgeClosure = 2;
   nFaceClosure = 0;

FunctionSpace/TetNedelecBasis.cpp

@@ -12,7 +12,7 @@ TetNedelecBasis::TetNedelecBasis(void){
     edgeV = refSpace->getAllEdge();
 
   // Set Basis Type //
-  this->order = 1;
+  this->order = 0;
 
   type = 1;
   dim  = 3;

FunctionSpace/TriEdgeBasis.cpp

@@ -20,8 +20,8 @@ TriEdgeBasis::TriEdgeBasis(unsigned int order){
 
   nVertex   = 0;
   nEdge     = 3 * (order + 1);
-  nFace     = 0;
-  nCell     = ((order - 1) * order + order - 1);
+  nFace     = ((order - 1) * order + order - 1);
+  nCell     = 0;
   nFunction = nVertex + nEdge + nFace + nCell;
 
   // Alloc Some Space //
@@ -97,7 +97,7 @@ TriEdgeBasis::TriEdgeBasis(unsigned int order){
     }
   }
 
-  // Cell Based //
+  // Face Based //
 
   // Preliminaries
   const Polynomial p = lagrange[2] * 2 - Polynomial(1, 0, 0, 0);

FunctionSpace/TriNedelecBasis.cpp

@@ -12,7 +12,7 @@ TriNedelecBasis::TriNedelecBasis(void){
     edgeV = refSpace->getAllEdge();
 
   // Set Basis Type //
-  order = 1;
+  order = 0;
 
   type = 1;
   dim  = 2;

FunctionSpace/TriNodeBasis.cpp

@@ -20,8 +20,8 @@ TriNodeBasis::TriNodeBasis(unsigned int order){
 
   nVertex   = 3;
   nEdge     = 3 * (order - 1);
-  nFace     = 0;
-  nCell     =     (order - 1) * (order - 2) / 2;
+  nFace     =     (order - 1) * (order - 2) / 2;
+  nCell     = 0;
   nFunction = nVertex + nEdge + nFace + nCell;
 
   // Alloc Some Space //
@@ -76,7 +76,7 @@ TriNodeBasis::TriNodeBasis(unsigned int order){
     }
   }
 
-  // Cell Based //
+  // Face Based //
   const Polynomial p = (lagrange[2] * 2) - Polynomial(1, 0, 0, 0);
   const int orderMinusTwo = order - 2;