From 06f806261a8cf12245a49f2cdc9e16028157b439 Mon Sep 17 00:00:00 2001
From: Christophe Geuzaine <cgeuzaine@ulg.ac.be>
Date: Fri, 20 May 2011 17:02:06 +0000
Subject: [PATCH] trying to fix abaqus export

---
 Geo/GModelIO_Mesh.cpp        |  3 ---
 Geo/MHexahedron.h            |  4 ----
 Geo/MTetrahedron.h           | 10 ++--------
 benchmarks/stl/mobilette.geo |  3 ++-
 4 files changed, 4 insertions(+), 16 deletions(-)

diff --git a/Geo/GModelIO_Mesh.cpp b/Geo/GModelIO_Mesh.cpp
index 3953916624..fe9ca23284 100644
--- a/Geo/GModelIO_Mesh.cpp
+++ b/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);
diff --git a/Geo/MHexahedron.h b/Geo/MHexahedron.h
index b54b3597a7..d7cb5336c3 100644
--- a/Geo/MHexahedron.h
+++ b/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 { 
diff --git a/Geo/MTetrahedron.h b/Geo/MTetrahedron.h
index f666db4a18..93e2a93e0a 100644
--- a/Geo/MTetrahedron.h
+++ b/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);
diff --git a/benchmarks/stl/mobilette.geo b/benchmarks/stl/mobilette.geo
index 8a10a2a6e5..43a248c0cd 100644
--- a/benchmarks/stl/mobilette.geo
+++ b/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;
-- 
GitLab