trying to fix abaqus export

Geo/GModelIO_Mesh.cpp

@@ -3499,9 +3499,6 @@ int GModel::writeINP(const std::string &name, bool saveAll,
     return 0;
   }
 
-  Msg::Error("Abaqus INP export has not been tested yet -- please give us"
-             "feedback!");
-
   if(noPhysicalGroups()) saveAll = true;
 
   indexMeshVertices(saveAll);

Geo/MHexahedron.h

@@ -517,10 +517,6 @@ class MHexahedronN : public MHexahedron {
   virtual int getPolynomialOrder() const { return (int)_order; }
   virtual int getNumVertices() const { return 8 + _vs.size(); }
   virtual MVertex *getVertex(int num){ return num < 8 ? _v[num] : _vs[num - 8]; }
-  virtual MVertex *getVertexDIFF(int num)
-  {
-    throw;
-  }
   virtual int getNumEdgeVertices() const { return 12 * (_order - 1); }
   virtual int getNumFaceVertices() const { return 6 * (_order - 1)*(_order - 1); }
   virtual int getNumVolumeVertices() const { 

Geo/MTetrahedron.h

@@ -223,25 +223,19 @@ class MTetrahedron10 : public MTetrahedron {
     static const int map[10] = {0, 1, 2, 3, 4, 5, 6, 7, 9, 8};
     return getVertex(map[num]);
   }
-  virtual MVertex *getVertexDIFF(int num)
-  {
-    static const int map[10] = {0, 1, 2, 3, 4, 5, 6, 7, 9, 8};
-    return getVertex(map[num]);
-  }
+  virtual MVertex *getVertexDIFF(int num){ return getVertexBDF(num); }
+  virtual MVertex *getVertexINP(int num){ return getVertexBDF(num); }
   virtual int getNumEdgeVertices() const { return 6; }
-
   virtual void getEdgeRep(int num, double *x, double *y, double *z, SVector3 *n);
   virtual int getNumEdgesRep();
   virtual void getFaceRep(int num, double *x, double *y, double *z, SVector3 *n);
   virtual int getNumFacesRep();
-
   virtual void getEdgeVertices(const int num, std::vector<MVertex*> &v) const
   {
     v.resize(3);
     MTetrahedron::_getEdgeVertices(num, v);
     v[2] = _vs[num];
   }
-
   virtual void getFaceVertices(const int num, std::vector<MVertex*> &v) const
   {
     v.resize(6);

benchmarks/stl/mobilette.geo

@@ -26,4 +26,5 @@ EndFor
 Surface Loop(1) = {s : s + #ss[]-1};
 Volume(1) = {1};
 
-//Physical Surface(1) = {s : s + #ss[]-1};
+Physical Surface(1) = {s : s + #ss[]-1};
+Physical Volume(1) = 1;