diff --git a/Geo/MElement.cpp b/Geo/MElement.cpp
index 88ef313a568b584dbfb162e7d180335d887f7a40..a1474e1c00e4dada8f1390e7e2644c106d8929c9 100644
--- a/Geo/MElement.cpp
+++ b/Geo/MElement.cpp
@@ -1447,7 +1447,6 @@ MElement *MElementFactory::create(int type, std::vector<MVertex*> &v,
   case MSH_PRI_18:  return new MPrism18(v, num, part);
   case MSH_TET_20:  return new MTetrahedronN(v, 3, num, part);
   case MSH_TET_35:  return new MTetrahedronN(v, 4, num, part);
-  case MSH_TET_28:  return new MTetrahedronN(v, 5, num, part);
   case MSH_TET_56:  return new MTetrahedronN(v, 5, num, part);
   case MSH_TET_84:  return new MTetrahedronN(v, 6, num, part);
   case MSH_TET_120: return new MTetrahedronN(v, 7, num, part);
diff --git a/Geo/MTriangle.cpp b/Geo/MTriangle.cpp
index 121307e6aca2c049f4df1a78b6171cabed32cbfc..c8579331971e5efb72a32730ef665c678194fa88 100644
--- a/Geo/MTriangle.cpp
+++ b/Geo/MTriangle.cpp
@@ -142,7 +142,7 @@ const nodalBasis* MTriangle::getFunctionSpace(int order) const
 
 const JacobianBasis* MTriangle::getJacobianFuncSpace(int order) const
 {
-  if (order == -1) return BasisFactory::getNodalBasis(getTypeForMSH());
+  if (order == -1) return BasisFactory::getJacobianBasis(getTypeForMSH());
 
   switch (order) {
     case 1: return BasisFactory::getJacobianBasis(MSH_TRI_3);