fix order node of high order prism

Common/OpenFile.cpp

@@ -451,11 +451,11 @@ int MergeFile(const std::string &fileName, bool warnIfMissing)
         tmp->readMSH(fileName);
         status = GeomMeshMatcher::instance()->match(tmp2, tmp);
         delete tmp;
-	GModel::setCurrent(tmp2);
-	tmp2->setVisibility(1);
+        GModel::setCurrent(tmp2);
+        tmp2->setVisibility(1);
       }
       else
-	status = GModel::current()->readMSH(fileName);
+        status = GModel::current()->readMSH(fileName);
 #if defined(HAVE_POST)
       if(status > 1) status = PView::readMSH(fileName);
 #endif

Geo/MPrism.h

@@ -494,6 +494,10 @@ class MPrismN : public MPrism {
     }
     return "";
   }
+  virtual void reverse()
+  {
+    Msg::Error("Reverse not implemented yet for MPrismN");
+  }
   virtual void getNode(int num, double &u, double &v, double &w) const
   {
     const fullMatrix<double> &p = getFunctionSpace()->points;

benchmarks/3d/sph.geo

@@ -16,11 +16,11 @@ Circle(3) = {5,1,6};
 Circle(4) = {2,1,3};
 Circle(5) = {7,1,6};
 Circle(6) = {4,1,3};
-Line(7) = {1,5};
+//Line(7) = {1,5};
 Line(8) = {5,2};
-Line(9) = {1,6};
+//Line(9) = {1,6};
 Line(10) = {6,3};
-Line(11) = {1,7};
+//Line(11) = {1,7};
 Line(12) = {7,4};
 
 Line Loop(17) = {2,-8,-1,12};