From b5a108eeeac323232a437fa5e0e30c4ccb270d2f Mon Sep 17 00:00:00 2001
From: eap <eap@opencascade.com>
Date: Wed, 11 May 2011 08:22:42 +0000
Subject: [PATCH] 0021238: EDF 1817 SMESH: GHS3D on quadratic meshes    fix
 crash on test_ghs3d_wo_geom_from_blqurf_quadratic_.py

---
 src/GHS3DPlugin/GHS3DPlugin_GHS3D.cxx | 28 +++++++++++++--------------
 1 file changed, 14 insertions(+), 14 deletions(-)

diff --git a/src/GHS3DPlugin/GHS3DPlugin_GHS3D.cxx b/src/GHS3DPlugin/GHS3DPlugin_GHS3D.cxx
index b64b50a..e22d8ac 100644
--- a/src/GHS3DPlugin/GHS3DPlugin_GHS3D.cxx
+++ b/src/GHS3DPlugin/GHS3DPlugin_GHS3D.cxx
@@ -1063,9 +1063,9 @@ static bool readGMFFile(const char*                     theFile,
   std::cout << "theHelper->GetMesh()->NbNodes(): " << nbMeshNodes << std::endl;
   
   const bool isQuadMesh = 
-  theHelper->GetMesh()->NbEdges( ORDER_QUADRATIC ) ||
-  theHelper->GetMesh()->NbFaces( ORDER_QUADRATIC ) ||
-  theHelper->GetMesh()->NbVolumes( ORDER_QUADRATIC );
+    theHelper->GetMesh()->NbEdges( ORDER_QUADRATIC ) ||
+    theHelper->GetMesh()->NbFaces( ORDER_QUADRATIC ) ||
+    theHelper->GetMesh()->NbVolumes( ORDER_QUADRATIC );
   std::cout << "isQuadMesh: " << isQuadMesh << std::endl;
   
   if (theHelper->GetSubShapeID() != 0)
@@ -1138,16 +1138,16 @@ static bool readGMFFile(const char*                     theFile,
       std::cout << nbInitialNodes << " from input mesh " << std::endl;
 
       // Remove orphan nodes from previous enforced mesh which was cleared
-      if ( nbElem < nbMeshNodes ) {
-        const SMDS_MeshNode* node;
-        SMDS_NodeIteratorPtr nodeIt = theMeshDS->nodesIterator();
-        while ( nodeIt->more() )
-        {
-          node = nodeIt->next();
-          if (theNodeToGhs3dIdMap.find(node) != theNodeToGhs3dIdMap.end())
-            theMeshDS->RemoveNode(node);
-        }
-      }
+//       if ( nbElem < nbMeshNodes ) {
+//         const SMDS_MeshNode* node;
+//         SMDS_NodeIteratorPtr nodeIt = theMeshDS->nodesIterator();
+//         while ( nodeIt->more() )
+//         {
+//           node = nodeIt->next();
+//           if (theNodeToGhs3dIdMap.find(node) != theNodeToGhs3dIdMap.end())
+//             theMeshDS->RemoveNode(node);
+//         }
+//       }
 
       
       int aGMFID;
@@ -1374,7 +1374,7 @@ static bool writeGMFFile(const char*                                     theMesh
                          GHS3DPlugin_Hypothesis::TGHS3DEnforcedVertexCoordsValues & theEnforcedVertices)
 {
   MESSAGE("writeGMFFile w/o geometry");
-  int idx, idxRequired = 0, idxSol;
+  int idx, idxRequired = 0, idxSol = 0;
   const int dummyint = 0;
   GHS3DPlugin_Hypothesis::TGHS3DEnforcedVertexCoordsValues::const_iterator vertexIt;
   std::vector<double> enfVertexSizes;
-- 
2.39.2