Salome HOME
remove "using namespace std" from SMESH headers
[plugins/hexoticplugin.git] / src / HexoticPlugin / HexoticPlugin_Hexotic.cxx
index c3c54a39ce4b4cf4f0483ac1c925b1468aee80ab..73ef7d120485eab4de31bb056f7a1dad6ebfda3f 100644 (file)
@@ -1,9 +1,9 @@
-// Copyright (C) 2007-2013  CEA/DEN, EDF R&D
+// Copyright (C) 2007-2015  CEA/DEN, EDF R&D
 //
 // This library is free software; you can redistribute it and/or
 // modify it under the terms of the GNU Lesser General Public
 // License as published by the Free Software Foundation; either
-// version 2.1 of the License.
+// version 2.1 of the License, or (at your option) any later version.
 //
 // This library is distributed in the hope that it will be useful,
 // but WITHOUT ANY WARRANTY; without even the implied warranty of
@@ -31,6 +31,7 @@
 #include <sys/sysinfo.h>
 #else
 #include <errno.h>
+#include <process.h>
 #endif
 
 #ifdef _DEBUG_
@@ -98,14 +99,14 @@ HexoticPlugin_Hexotic::HexoticPlugin_Hexotic(int hypId, int studyId, SMESH_Gen*
   : SMESH_3D_Algo(hypId, studyId, gen)
 {
   MESSAGE("HexoticPlugin_Hexotic::HexoticPlugin_Hexotic");
-  _name = "Hexotic_3D";
+  _name = "MG-Hexa";
   _shapeType = (1 << TopAbs_SHELL) | (1 << TopAbs_SOLID);// 1 bit /shape type
 //   _onlyUnaryInput = false;
   _requireShape = false;
   _iShape=0;
   _nbShape=0;
   _hexoticFilesKept=false;
-  _compatibleHypothesis.push_back("Hexotic_Parameters");
+  _compatibleHypothesis.push_back( HexoticPlugin_Hypothesis::GetHypType() );
 #ifdef WITH_BLSURFPLUGIN
   _blsurfHypo = NULL;
 #endif
@@ -146,7 +147,7 @@ bool HexoticPlugin_Hexotic::CheckBLSURFHypothesis( SMESH_Mesh&         aMesh,
 
   // If a BLSURF hypothesis is applied, get it
   SMESH_HypoFilter blsurfFilter;
-  blsurfFilter.Init( blsurfFilter.HasName( "BLSURF_Parameters" ));
+  blsurfFilter.Init( blsurfFilter.HasName( BLSURFPlugin_Hypothesis::GetHypType() ));
   std::list<const SMESHDS_Hypothesis *> appliedHyps;
   aMesh.GetHypotheses( aShape, blsurfFilter, appliedHyps, false );
 
@@ -154,7 +155,7 @@ bool HexoticPlugin_Hexotic::CheckBLSURFHypothesis( SMESH_Mesh&         aMesh,
     itl = appliedHyps.begin();
     theHyp = (*itl); // use only the first hypothesis
     std::string hypName = theHyp->GetName();
-    if (hypName == "BLSURF_Parameters") {
+    if (hypName == BLSURFPlugin_Hypothesis::GetHypType()) {
       _blsurfHypo = static_cast<const BLSURFPlugin_Hypothesis*> (theHyp);
       ASSERT(_blsurfHypo);
       return true;
@@ -191,7 +192,7 @@ bool HexoticPlugin_Hexotic::CheckHypothesis( SMESH_Mesh&
   theHyp = (*itl); // use only the first hypothesis
 
   std::string hypName = theHyp->GetName();
-  if (hypName == "Hexotic_Parameters") {
+  if (hypName == HexoticPlugin_Hypothesis::GetHypType() ) {
     _hypothesis = static_cast<const HexoticPlugin_Hypothesis*> (theHyp);
     ASSERT(_hypothesis);
     aStatus = SMESH_Hypothesis::HYP_OK;
@@ -248,44 +249,6 @@ static TopoDS_Shape findShape(SMDS_MeshNode**     t_Node,
   return aShape;
 }
 
-//=======================================================================
-//function : findEdge
-//purpose  :
-//=======================================================================
-
-static int findEdge(const SMDS_MeshNode* aNode,
-                    const SMESHDS_Mesh*  aMesh,
-                    const int            nEdge,
-                    const TopoDS_Shape*  t_Edge) {
-
-  TopoDS_Shape aPntShape, foundEdge;
-  TopoDS_Vertex aVertex;
-  gp_Pnt aPnt( aNode->X(), aNode->Y(), aNode->Z() );
-
-  int foundInd, ind;
-  double nearest = RealLast(), *t_Dist;
-  double epsilon = Precision::Confusion();
-
-  t_Dist = new double[ nEdge ];
-  aPntShape = BRepBuilderAPI_MakeVertex( aPnt ).Shape();
-  aVertex   = TopoDS::Vertex( aPntShape );
-
-  for ( ind=0; ind < nEdge; ind++ ) {
-    BRepExtrema_DistShapeShape aDistance ( aVertex, t_Edge[ind] );
-    t_Dist[ind] = aDistance.Value();
-    if ( t_Dist[ind] < nearest ) {
-      nearest   = t_Dist[ind];
-      foundEdge = t_Edge[ind];
-      foundInd  = ind;
-      if ( nearest < epsilon )
-        ind = nEdge;
-    }
-  }
-
-  delete [] t_Dist;
-  return aMesh->ShapeToIndex( foundEdge );
-}
-
 //=======================================================================
 //function : getNbShape
 //purpose  :
@@ -352,7 +315,7 @@ void getShape(Mesh* mesh, Shape shape, Tab *t_Shape) {
 
 static void printWarning(const int nbExpected, std::string aString, const int nbFound) {
   cout << std::endl;
-  cout << "WARNING : " << nbExpected << " " << aString << " expected, Hexotic has found " << nbFound << std::endl;
+  cout << "WARNING : " << nbExpected << " " << aString << " expected, MG-Hexa has found " << nbFound << std::endl;
   cout << "=======" << std::endl;
   cout << std::endl;
   return;
@@ -368,43 +331,35 @@ static void removeHexoticFiles(TCollection_AsciiString file_In, TCollection_Asci
   removeFile( file_Out );
 }
 
+//=======================================================================
+//function : splitQuads
+//purpose  : splits all quadrangles into triangles
+//=======================================================================
+
+static void splitQuads(SMESH_Mesh& aMesh)
+{
+  SMESH_MeshEditor spliter( &aMesh );
+
+  TIDSortedElemSet elems;
+  SMDS_ElemIteratorPtr eIt = aMesh.GetMeshDS()->elementsIterator();
+  while( eIt->more() )
+    elems.insert( elems.end(), eIt->next() );
+  
+  spliter.QuadToTri ( elems, /*the13Diag=*/true);
+}
+
 //=======================================================================
 //function : readResult
 //purpose  : Read GMF file in case of a mesh with geometry
 //=======================================================================
 
 static bool readResult(std::string         theFile,
-#ifdef WITH_SMESH_CANCEL_COMPUTE
                        HexoticPlugin_Hexotic*  theAlgo,
-#endif
                        SMESHDS_Mesh*       theMesh,
                        const int           nbShape,
                        const TopoDS_Shape* tabShape,
                        double**            tabBox)
 {
-  // ---------------------------------
-  // Optimisation of the plugin ...
-  // Retrieve the correspondance edge --> shape
-  // (which is very costly) only when user
-  // has defined at least one group of edges
-  // which should be rare for a 3d mesh !
-  // ---------------------------------
-  
-  bool retrieve_edges = false;
-  const std::set<SMESHDS_GroupBase*>& aGroups = theMesh->GetGroups();
-  set<SMESHDS_GroupBase*>::const_iterator GrIt = aGroups.begin();
-  for (; GrIt != aGroups.end(); GrIt++)
-    {
-      SMESHDS_GroupBase* aGrp = *GrIt;
-      if ( !aGrp )
-        continue;
-      if ( aGrp->GetType() == SMDSAbs_Edge )
-        {
-          retrieve_edges = true;
-          break;
-        }
-    }
-  
   // ---------------------------------
   // Read generated elements and nodes
   // ---------------------------------
@@ -480,7 +435,7 @@ static bool readResult(std::string         theFile,
     }
 
     nbElem = 0;
-    if ( nField < (mapField.size() - 1) && nField >= 0 )
+    if ( nField < int(mapField.size() - 1) && nField >= 0 )
       fileRes >> nbElem;
 
     switch (nField) {
@@ -500,12 +455,10 @@ static bool readResult(std::string         theFile,
 
         coord = new double[nbRef];
         for ( int iElem = 0; iElem < nbElem; iElem++ ) {
-#ifdef WITH_SMESH_CANCEL_COMPUTE
           if(theAlgo->computeCanceled())
             {
               return false;
             }
-#endif
           aHexoticID = iElem + 1;
           for ( int iCoord = 0; iCoord < 3; iCoord++ )
             fileRes >> coord[ iCoord ];
@@ -530,12 +483,10 @@ static bool readResult(std::string         theFile,
         node   = new SMDS_MeshNode*[ nbRef ];
         nodeID = new int[ nbRef ];
         for ( int iElem = 0; iElem < nbElem; iElem++ ) {
-#ifdef WITH_SMESH_CANCEL_COMPUTE
           if(theAlgo->computeCanceled())
             {
               return false;
             }
-#endif
           for ( int iRef = 0; iRef < nbRef; iRef++ ) {
             fileRes >> aHexoticNodeID;                          // read nbRef aHexoticNodeID
             node[ iRef ]   = HexoticNode[ aHexoticNodeID ];
@@ -553,6 +504,10 @@ static bool readResult(std::string         theFile,
                 if ( aPnt.Distance( HexoticPnt ) < epsilon )
                   break;
               }
+              if ( nodeAssigne[ nodeID[0] ] != 0 ) { // because "Edges" go before "Corners"
+                theMesh->UnSetNodeOnShape( node[0] );
+                nodeAssigne[ nodeID[0] ] = 0;
+              }
               break;
             }
             case 4: { // "Edges"
@@ -561,10 +516,7 @@ static bool readResult(std::string         theFile,
               int iNode = 1;
               if ( nodeAssigne[ nodeID[0] ] == 0 || nodeAssigne[ nodeID[0] ] == 2 )
                 iNode = 0;
-              if(retrieve_edges)
-                shapeID = findEdge( node[iNode], theMesh, nbShapeEdge, tabEdge );
-              else
-                shapeID = 0;
+              shapeID = dummy;
               break;
             }
             case 5: { // "Ridges"
@@ -664,9 +616,7 @@ static bool readResult(std::string         theFile,
 //=======================================================================
 
 static bool readResult(std::string theFile,
-#ifdef WITH_SMESH_CANCEL_COMPUTE
                        HexoticPlugin_Hexotic*  theAlgo,
-#endif
                        SMESH_MesherHelper* theHelper)
 {
   SMESHDS_Mesh* theMesh = theHelper->GetMeshDS();
@@ -694,7 +644,15 @@ static bool readResult(std::string theFile,
   mapField["Hexahedra"]            = 7; tabRef[7] = 8; tabDummy[7] = true;
   mapField["End"]                  = 8; tabRef[8] = 0; tabDummy[8] = false;
 
-  theHelper->GetMesh()->Clear();
+  {
+    // theMesh->Clear(); -- this does not remove imported mesh
+    SMDS_ElemIteratorPtr eIt = theMesh->elementsIterator();
+    while( eIt->more() )
+      theMesh->RemoveFreeElement( eIt->next(), /*sm=*/0 );
+    SMDS_NodeIteratorPtr nIt = theMesh->nodesIterator();
+    while ( nIt->more() )
+      theMesh->RemoveFreeNode( nIt->next(), /*sm=*/0 );
+  }
 
   int nbVertices = getNbShape(theFile, "Vertices");
   HexoticNode = new SMDS_MeshNode*[ nbVertices + 1 ];
@@ -720,7 +678,7 @@ static bool readResult(std::string theFile,
     }
 
     nbElem = 0;
-    if ( nField < (mapField.size() - 1) && nField >= 0 )
+    if ( nField < int(mapField.size() - 1) && nField >= 0 )
       fileRes >> nbElem;
 
     switch (nField) {
@@ -739,12 +697,10 @@ static bool readResult(std::string theFile,
       SMDS_MeshNode * aHexoticNode;
 
       for ( int iElem = 0; iElem < nbElem; iElem++ ) {
-#ifdef WITH_SMESH_CANCEL_COMPUTE
         if(theAlgo->computeCanceled())
           {
             return false;
           }
-#endif
         aHexoticID = iElem + 1;
         for ( int iCoord = 0; iCoord < 3; iCoord++ )
           fileRes >> coord[ iCoord ];
@@ -766,12 +722,10 @@ static bool readResult(std::string theFile,
 
       for ( int iElem = 0; iElem < nbElem; iElem++ )
       {
-#ifdef WITH_SMESH_CANCEL_COMPUTE
         if(theAlgo->computeCanceled())
           {
             return false;
           }
-#endif
         for ( int iRef = 0; iRef < nbRef; iRef++ )
         {
           fileRes >> aHexoticNodeID;                          // read nbRef aHexoticNodeID
@@ -821,7 +775,7 @@ static bool readResult(std::string theFile,
 
 //=============================================================================
 /*!
- * Pass parameters to Hexotic
+ * Pass parameters to MG-Hexa
  */
 //=============================================================================
 
@@ -841,11 +795,18 @@ void HexoticPlugin_Hexotic::SetParameters(const HexoticPlugin_Hypothesis* hyp) {
     _hexoticVerbosity = hyp->GetHexoticVerbosity();
     _hexoticMaxMemory = hyp->GetHexoticMaxMemory();
     _hexoticSdMode = hyp->GetHexoticSdMode();
+    _textOptions = hyp->GetTextOptions();
     _sizeMaps = hyp->GetSizeMaps();
+    _nbLayers = hyp->GetNbLayers();
+    _firstLayerSize = hyp->GetFirstLayerSize();
+    _direction = hyp->GetDirection();
+    _growth = hyp->GetGrowth();
+    _facesWithLayers = hyp->GetFacesWithLayers();
+    _imprintedFaces = hyp->GetImprintedFaces();
   }
   else {
     cout << std::endl;
-    cout << "WARNING : The Hexotic default parameters are taken into account" << std::endl;
+    cout << "WARNING : The MG-Hexa default parameters are taken into account" << std::endl;
     cout << "=======" << std::endl;
     _hexesMinLevel = hyp->GetDefaultHexesMinLevel();
     _hexesMaxLevel = hyp->GetDefaultHexesMaxLevel();
@@ -859,7 +820,14 @@ void HexoticPlugin_Hexotic::SetParameters(const HexoticPlugin_Hypothesis* hyp) {
     _hexoticVerbosity = hyp->GetDefaultHexoticVerbosity();
     _hexoticMaxMemory = hyp->GetDefaultHexoticMaxMemory();
     _hexoticSdMode = hyp->GetDefaultHexoticSdMode();
+    _textOptions = hyp->GetDefaultTextOptions();
     _sizeMaps = hyp->GetDefaultHexoticSizeMaps();
+    _nbLayers = hyp->GetDefaultNbLayers();
+    _firstLayerSize = hyp->GetDefaultFirstLayerSize();
+    _direction = hyp->GetDefaultDirection();
+    _growth = hyp->GetDefaultGrowth();
+    _facesWithLayers = hyp->GetDefaultFacesWithLayers();
+    _imprintedFaces = hyp->GetDefaultImprintedFaces();
   }
 }
 
@@ -868,35 +836,36 @@ void HexoticPlugin_Hexotic::SetParameters(const HexoticPlugin_Hypothesis* hyp) {
 //purpose  :
 //=======================================================================
 
-static TCollection_AsciiString getTmpDir()
-{
-  TCollection_AsciiString aTmpDir;
-
-  char *Tmp_dir = getenv("SALOME_TMP_DIR");
-#ifdef WIN32
-  if(Tmp_dir == NULL) {
-    Tmp_dir = getenv("TEMP");
-    if( Tmp_dir== NULL )
-      Tmp_dir = getenv("TMP");
-#endif
-
-  if(Tmp_dir != NULL) {
-    aTmpDir = Tmp_dir;
-#ifdef WIN32
-    if(aTmpDir.Value(aTmpDir.Length()) != '\\') aTmpDir+='\\';
-#else
-    if(aTmpDir.Value(aTmpDir.Length()) != '/') aTmpDir+='/';
-#endif
-  }
-  else {
-#ifdef WIN32
-    aTmpDir = TCollection_AsciiString("C:\\");
-#else
-    aTmpDir = TCollection_AsciiString("/tmp/");
-#endif
-  }
-  return aTmpDir;
-}
+// static TCollection_AsciiString getTmpDir()
+// {
+//   TCollection_AsciiString aTmpDir;
+
+//   char *Tmp_dir = getenv("SALOME_TMP_DIR");
+// #ifdef WIN32
+//   if(Tmp_dir == NULL) {
+//     Tmp_dir = getenv("TEMP");
+//     if( Tmp_dir== NULL )
+//       Tmp_dir = getenv("TMP");
+//   }
+// #endif
+
+//   if(Tmp_dir != NULL) {
+//     aTmpDir = Tmp_dir;
+// #ifdef WIN32
+//     if(aTmpDir.Value(aTmpDir.Length()) != '\\') aTmpDir+='\\';
+// #else
+//     if(aTmpDir.Value(aTmpDir.Length()) != '/') aTmpDir+='/';
+// #endif
+//   }
+//   else {
+// #ifdef WIN32
+//     aTmpDir = TCollection_AsciiString("C:\\");
+// #else
+//     aTmpDir = TCollection_AsciiString("/tmp/");
+// #endif
+//   }
+//   return aTmpDir;
+// }
 
 //=======================================================================
 //function : getSuffix
@@ -907,18 +876,26 @@ static TCollection_AsciiString getSuffix()
 {
   TCollection_AsciiString aSuffix = "";
   aSuffix += "_";
+#ifndef WIN32
   aSuffix += getenv("USER");
+#else
+  aSuffix += getenv("USERNAME");
+#endif
   aSuffix += "_";
   aSuffix += Kernel_Utils::GetHostname().c_str();
   aSuffix += "_";
+#ifndef WIN32
   aSuffix += getpid();
+#else
+  aSuffix += _getpid();
+#endif
 
   return aSuffix;
 }
 
 //================================================================================
 /*!
- * \brief Returns a command to run Hexotic mesher
+ * \brief Returns a command to run MG-Hexa mesher
  */
 //================================================================================
 
@@ -927,7 +904,7 @@ std::string HexoticPlugin_Hexotic::getHexoticCommand(const TCollection_AsciiStri
                                                      const TCollection_AsciiString& Hexotic_SizeMap_Prefix) const
 {
   cout << std::endl;
-  cout << "Hexotic execution..." << std::endl;
+  cout << "MG-Hexa execution..." << std::endl;
   cout << _name << " parameters :" << std::endl;
   cout << "    " << _name << " Verbosity = " << _hexoticVerbosity << std::endl;
   cout << "    " << _name << " Max Memory = " << _hexoticMaxMemory << std::endl;
@@ -935,12 +912,37 @@ std::string HexoticPlugin_Hexotic::getHexoticCommand(const TCollection_AsciiStri
   cout << "    " << _name << " Segments Max Level = " << _hexesMaxLevel << std::endl;
   cout << "    " << _name << " Segments Min Size = " << _hexesMinSize << std::endl;
   cout << "    " << _name << " Segments Max Size = " << _hexesMaxSize << std::endl;
-  cout << "    " << "Hexotic can ignore ridges : " << (_hexoticIgnoreRidges ? "yes":"no") << std::endl;
-  cout << "    " << "Hexotic authorize invalide elements : " << ( _hexoticInvalidElements ? "yes":"no") << std::endl;
+  cout << "    " << "MG-Hexa can ignore ridges : " << (_hexoticIgnoreRidges ? "yes":"no") << std::endl;
+  cout << "    " << "MG-Hexa authorize invalide elements : " << ( _hexoticInvalidElements ? "yes":"no") << std::endl;
   cout << "    " << _name << " Sharp angle threshold = " << _hexoticSharpAngleThreshold << " degrees" << std::endl;
   cout << "    " << _name << " Number of threads = " << _hexoticNbProc << std::endl;
   cout << "    " << _name << " Working directory = \"" << _hexoticWorkingDirectory << "\"" << std::endl;
   cout << "    " << _name << " Sub. Dom mode = " << _hexoticSdMode << std::endl;
+  cout << "    " << _name << " Text options = \"" << _textOptions << "\"" << std::endl;
+  cout << "    " << _name << " Number of layers = " << _nbLayers << std::endl;
+  cout << "    " << _name << " Size of the first layer  = " << _firstLayerSize << std::endl;
+  cout << "    " << _name << " Direction of the layers = " << ( _direction ? "Inward" : "Outward" ) << std::endl;
+  cout << "    " << _name << " Growth = " << _growth << std::endl;
+  if (!_facesWithLayers.empty()) {
+    cout << "    " << _name << " Faces with layers = ";
+    for (size_t i = 0; i < _facesWithLayers.size(); i++)
+    {
+      cout << _facesWithLayers.at(i);
+      if ((i + 1) != _facesWithLayers.size())
+        cout << ", ";
+    }
+    cout << std::endl;
+  }
+  if (!_imprintedFaces.empty()) {
+    cout << "    " << _name << " Imprinted faces = ";
+    for (size_t i = 0; i < _imprintedFaces.size(); i++)
+    {
+      cout << _imprintedFaces.at(i);
+      if ((i + 1) != _imprintedFaces.size())
+        cout << ", ";
+    }
+    cout << std::endl;
+  }
 
   TCollection_AsciiString run_Hexotic( "mg-hexa.exe" );
 
@@ -950,11 +952,21 @@ std::string HexoticPlugin_Hexotic::getHexoticCommand(const TCollection_AsciiStri
   TCollection_AsciiString sizeMap = " --read_sizemap ";
   TCollection_AsciiString ignoreRidges = " --compute_ridges no ", invalideElements = " --allow_invalid_elements yes ";
   TCollection_AsciiString subdom = " --components ";
+#ifndef WIN32
   TCollection_AsciiString proc = " --max_number_of_threads ";
+#endif
   TCollection_AsciiString verb = " --verbose ";
   TCollection_AsciiString maxmem = " --max_memory ";
 
-  TCollection_AsciiString minLevel, maxLevel, minSize, maxSize, sharpAngle, mode, nbproc, verbosity, maxMemory;
+  TCollection_AsciiString comNbLayers = " --number_of_boundary_layers ";
+  TCollection_AsciiString comFirstLayerSize = " --height_of_the_first_layer ";
+  TCollection_AsciiString comDirection = " --boundary_layers_subdomain_direction ";
+  TCollection_AsciiString comGrowth = " --boundary_layers_geometric_progression ";
+  TCollection_AsciiString comFacesWithLayers = " --boundary_layers_surface_ids ";
+  TCollection_AsciiString comImptintedFaces = " --imprinted_surface_ids ";
+
+  TCollection_AsciiString minLevel, maxLevel, minSize, maxSize, sharpAngle, mode, nbproc, verbosity, maxMemory,
+                          textOptions, nbLayers, firstLayerSize, direction, growth, facesWithLayers, imprintedFaces;
   minLevel = _hexesMinLevel;
   maxLevel = _hexesMaxLevel;
   minSize = _hexesMinSize;
@@ -979,6 +991,23 @@ std::string HexoticPlugin_Hexotic::getHexoticCommand(const TCollection_AsciiStri
   nbproc = _hexoticNbProc;
   verbosity = _hexoticVerbosity;
   maxMemory = _hexoticMaxMemory;
+  textOptions = (" " + _textOptions + " ").c_str();
+  nbLayers = _nbLayers;
+  firstLayerSize = _firstLayerSize;
+  direction = _direction ? "1" : "-1";
+  growth = _growth;
+  for (size_t i = 0; i < _facesWithLayers.size(); i++)
+  {
+    facesWithLayers += _facesWithLayers[i];
+    if ((i + 1) != _facesWithLayers.size())
+      facesWithLayers += ",";
+  }
+  for (size_t i = 0; i < _imprintedFaces.size(); i++)
+  {
+    imprintedFaces += _imprintedFaces[i];
+    if ((i + 1) != _imprintedFaces.size())
+      imprintedFaces += ",";
+  }
 
   if (_hexoticIgnoreRidges)
     run_Hexotic +=  ignoreRidges;
@@ -1004,12 +1033,30 @@ std::string HexoticPlugin_Hexotic::getHexoticCommand(const TCollection_AsciiStri
   if (_sizeMaps.begin() != _sizeMaps.end())
     run_Hexotic += sizeMap + Hexotic_SizeMap_Prefix;
 
+  if (_nbLayers       > 0 &&
+      _firstLayerSize > 0 &&
+      _growth         > 0 &&
+      !_facesWithLayers.empty())
+  {
+    run_Hexotic += comNbLayers + nbLayers;
+    run_Hexotic += comFirstLayerSize + firstLayerSize;
+    run_Hexotic += comDirection + direction;
+    run_Hexotic += comGrowth + growth;
+    run_Hexotic += comFacesWithLayers + facesWithLayers;
+    if (!_imprintedFaces.empty())
+      run_Hexotic += comImptintedFaces + imprintedFaces;
+  }
   run_Hexotic += in + Hexotic_In + out + Hexotic_Out;
   run_Hexotic += subdom + mode;
+#ifndef WIN32
   run_Hexotic += proc + nbproc;
+#endif
   run_Hexotic += verb + verbosity;
   run_Hexotic += maxmem + maxMemory;
 
+  if (!_textOptions.empty())
+    run_Hexotic += textOptions;
+
   return run_Hexotic.ToCString();
 }
 
@@ -1147,9 +1194,13 @@ void HexoticPlugin_Hexotic::createPointsSampleFromFace( const TopoDS_Shape& aSha
   BRepMesh_IncrementalMesh M(aShape, 0.01, Standard_True);
   TopLoc_Location aLocation;
   TopoDS_Face aFace = TopoDS::Face(aShape);
-  Handle(Geom_Surface) aSurf = BRep_Tool::Surface (aFace, aLocation);
+
+  // Triangulate the face
   Handle(Poly_Triangulation) aTri = BRep_Tool::Triangulation (aFace, aLocation);
   
+  // Get the transformation associated to the face location
+  gp_Trsf aTrsf = aLocation.Transformation();
+  
   // Get triangles
   int nbTriangles = aTri->NbTriangles();
   Poly_Array1OfTriangle triangles(1,nbTriangles);
@@ -1168,6 +1219,10 @@ void HexoticPlugin_Hexotic::createPointsSampleFromFace( const TopoDS_Shape& aSha
      gp_Pnt p2 = nodes.Value(aTriangle.Value(2));
      gp_Pnt p3 = nodes.Value(aTriangle.Value(3));
      
+     p1.Transform(aTrsf);
+     p2.Transform(aTrsf);
+     p3.Transform(aTrsf);
+     
      subdivideTriangle(p1, p2, p3, theSize, thePoints);  
   }
 }
@@ -1272,7 +1327,7 @@ void HexoticPlugin_Hexotic::subdivideTriangle( const gp_Pnt& p1,
   // If the traingles share a Vertex and no Edge the distance of the mass centers 
   // to the Vertices is 2*D < S so the mass centers are distant of less than 2*S 
   
-  double threshold = sqrt( 3 ) * theSize;
+  double threshold = sqrt( 3. ) * theSize;
   
   if ( (p1.Distance(p2) > threshold ||
         p2.Distance(p3) > threshold ||
@@ -1395,16 +1450,14 @@ gp_Pnt HexoticPlugin_Hexotic::tangencyPoint(const gp_Pnt& p1,
 
 //=============================================================================
 /*!
- * Here we are going to use the Hexotic mesher
+ * Here we are going to use the MG-Hexa mesher
  */
 //=============================================================================
 
 bool HexoticPlugin_Hexotic::Compute(SMESH_Mesh&          aMesh,
-                                     const TopoDS_Shape& aShape)
+                                    const TopoDS_Shape& aShape)
 {
-#ifdef WITH_SMESH_CANCEL_COMPUTE
   _compute_canceled = false;
-#endif
   bool Ok = true;
   SMESHDS_Mesh* meshDS = aMesh.GetMeshDS();
   TCollection_AsciiString hexahedraMessage;
@@ -1415,7 +1468,7 @@ bool HexoticPlugin_Hexotic::Compute(SMESH_Mesh&          aMesh,
 
   // to prevent from displaying error message after computing,
   // SetIsAlwaysComputed( true ) to empty sub-meshes
-  vector< SMESH_subMesh* > subMeshesAlwaysComp;
+  std::vector< SMESH_subMesh* > subMeshesAlwaysComp;
   for ( int i = 0; i < _nbShape; ++i )
     if ( SMESH_subMesh* sm = aMesh.GetSubMeshContaining( aShape ))
     {
@@ -1463,15 +1516,14 @@ bool HexoticPlugin_Hexotic::Compute(SMESH_Mesh&          aMesh,
     SetParameters(_hypothesis);
 
 //     TCollection_AsciiString aTmpDir = getTmpDir();
-    TCollection_AsciiString aTmpDir = TCollection_AsciiString(_hexoticWorkingDirectory.c_str());
+    TCollection_AsciiString aTmpDir = _hexoticWorkingDirectory.c_str();
 #ifdef WIN32
     if ( aTmpDir.Value(aTmpDir.Length()) != '\\' ) aTmpDir += '\\';
 #else
     if ( aTmpDir.Value(aTmpDir.Length()) != '/' ) aTmpDir += '/';
 #endif
     TCollection_AsciiString Hexotic_In(""), Hexotic_Out, Hexotic_SizeMap_Prefix;
-    TCollection_AsciiString modeFile_In( "chmod 666 " ), modeFile_Out( "chmod 666 " );
-    TCollection_AsciiString aLogFileName = aTmpDir + "Hexotic"+getSuffix()+".log";    // log
+    TCollection_AsciiString modeFile_In( "chmod 666 " ), modeFile_Out( "chmod 666 " );    TCollection_AsciiString aLogFileName = aTmpDir + "Hexotic"+getSuffix()+".log";    // log
 
     std::map <int,int> aSmdsToHexoticIdMap;
     std::map <int,const SMDS_MeshNode*> aHexoticIdToNodeMap;
@@ -1480,16 +1532,18 @@ bool HexoticPlugin_Hexotic::Compute(SMESH_Mesh&          aMesh,
 #ifdef WITH_BLSURFPLUGIN
     bool defaultInputFile = true;
     if (_blsurfHypo && !_blsurfHypo->GetQuadAllowed()) {
-      Hexotic_In = TCollection_AsciiString(_blsurfHypo->GetGMFFile().c_str());
-      if (Hexotic_In != "")
+      Hexotic_In = _blsurfHypo->GetGMFFile().c_str();
+      if ( !Hexotic_In.IsEmpty() &&
+           SMESH_File( _blsurfHypo->GetGMFFile() ).exists() )
         defaultInputFile = false;
     }
     if (defaultInputFile) {
 #endif
       Hexotic_In  = aTmpDir + "Hexotic"+getSuffix()+"_In.mesh";
       removeHexoticFiles(Hexotic_In, Hexotic_Out);
+      splitQuads(aMesh); // quadrangles are no longer acceptable as input
       cout << std::endl;
-      cout << "Creating Hexotic input mesh file : " << Hexotic_In << std::endl;
+      cout << "Creating MG-Hexa input mesh file : " << Hexotic_In << std::endl;
       aMesh.ExportGMF(Hexotic_In.ToCString(), meshDS, true);
 #ifdef WITH_BLSURFPLUGIN
     }
@@ -1505,7 +1559,7 @@ bool HexoticPlugin_Hexotic::Compute(SMESH_Mesh&          aMesh,
     run_Hexotic += std::string(" 1> ") + aLogFileName.ToCString();  // dump into file
 
     cout << std::endl;
-    cout << "Hexotic command : " << run_Hexotic << std::endl;
+    cout << "MG-Hexa command : " << run_Hexotic << std::endl;
 
 #ifndef WIN32    
     modeFile_In += Hexotic_In;
@@ -1523,13 +1577,13 @@ bool HexoticPlugin_Hexotic::Compute(SMESH_Mesh&          aMesh,
     // --------------
 
     std::ifstream fileRes( Hexotic_Out.ToCString() );
+#ifndef WIN32  
     modeFile_Out += Hexotic_Out;
     system( modeFile_Out.ToCString() );
+#endif
     if ( ! fileRes.fail() ) {
       Ok = readResult( Hexotic_Out.ToCString(),
-#ifdef WITH_SMESH_CANCEL_COMPUTE
                        this,
-#endif
                        meshDS, _nbShape, tabShape, tabBox );
       if(Ok) {
 /*********************
@@ -1547,6 +1601,7 @@ bool HexoticPlugin_Hexotic::Compute(SMESH_Mesh&          aMesh,
       if (myError)
 */
         hexahedraMessage = "success";
+        #ifndef _DEBUG_
         removeFile(Hexotic_Out);
         removeFile(Hexotic_In);
         removeFile(aLogFileName);
@@ -1554,6 +1609,7 @@ bool HexoticPlugin_Hexotic::Compute(SMESH_Mesh&          aMesh,
         {
           removeFile( TCollection_AsciiString( sizeMapFiles[i].c_str() ) );
         }
+        #endif
       }
       else {
         hexahedraMessage = "failed"; 
@@ -1561,7 +1617,7 @@ bool HexoticPlugin_Hexotic::Compute(SMESH_Mesh&          aMesh,
     }
     else {
       hexahedraMessage = "failed";
-      cout << "Problem with Hexotic output file " << Hexotic_Out.ToCString() << std::endl;
+      cout << "Problem with MG-Hexa output file " << Hexotic_Out.ToCString() << std::endl;
       Ok = false;
       // analyse log file
       SMESH_File logFile( aLogFileName.ToCString() );
@@ -1574,11 +1630,11 @@ bool HexoticPlugin_Hexotic::Compute(SMESH_Mesh&          aMesh,
       }
 #ifndef WIN32
       if ( status > 0 && WEXITSTATUS(status) == 127 )
-        error("hexotic: command not found");
+        error("mg-hexa.exe: command not found");
 #else
       int err = errno;
       if ( status == 0 && err == ENOENT ) {
-        error("hexotic: command not found");
+        error("mg-hexa.exe: command not found");
       }
 #endif
     }
@@ -1596,10 +1652,8 @@ bool HexoticPlugin_Hexotic::Compute(SMESH_Mesh&          aMesh,
     _nbShape = 0;
     _iShape  = 0;
   }
-#ifdef WITH_SMESH_CANCEL_COMPUTE
   if(_compute_canceled)
     return error(SMESH_Comment("interruption initiated by user"));
-#endif
   return Ok;
 }
 
@@ -1616,9 +1670,7 @@ bool HexoticPlugin_Hexotic::Compute(SMESH_Mesh&          aMesh,
 
 bool HexoticPlugin_Hexotic::Compute(SMESH_Mesh & aMesh, SMESH_MesherHelper* aHelper)
 {
-#ifdef WITH_SMESH_CANCEL_COMPUTE
   _compute_canceled = false;
-#endif
 /*
   SMESH_ComputeErrorPtr myError = SMESH_ComputeError::New();
 */
@@ -1627,7 +1679,7 @@ bool HexoticPlugin_Hexotic::Compute(SMESH_Mesh & aMesh, SMESH_MesherHelper* aHel
 
   SetParameters(_hypothesis);
 
-  TCollection_AsciiString aTmpDir = getTmpDir();
+  TCollection_AsciiString aTmpDir = _hexoticWorkingDirectory.c_str();//getTmpDir();
   TCollection_AsciiString Hexotic_In, Hexotic_Out, Hexotic_SizeMap_Prefix;
   TCollection_AsciiString modeFile_In( "chmod 666 " ), modeFile_Out( "chmod 666 " );
   TCollection_AsciiString aLogFileName = aTmpDir + "Hexotic"+getSuffix()+".log";    // log
@@ -1647,8 +1699,10 @@ bool HexoticPlugin_Hexotic::Compute(SMESH_Mesh & aMesh, SMESH_MesherHelper* aHel
 
   removeHexoticFiles(Hexotic_In, Hexotic_Out);
 
+  splitQuads(aMesh); // quadrangles are no longer acceptable as input
+
   cout << std::endl;
-  cout << "Creating Hexotic input mesh file : " << Hexotic_In << std::endl;
+  cout << "Creating MG-Hexa input mesh file : " << Hexotic_In << std::endl;
   aMesh.ExportGMF(Hexotic_In.ToCString(), aHelper->GetMeshDS());
 #ifndef WIN32    
   modeFile_In += Hexotic_In;
@@ -1660,7 +1714,7 @@ bool HexoticPlugin_Hexotic::Compute(SMESH_Mesh & aMesh, SMESH_MesherHelper* aHel
   MESSAGE("HexoticPlugin_Hexotic::Compute");
 
   cout << std::endl;
-  cout << "Hexotic command : " << run_Hexotic << std::endl;
+  cout << "MG-Hexa command : " << run_Hexotic << std::endl;
   system( run_Hexotic.data() );
 
   // --------------
@@ -1672,9 +1726,7 @@ bool HexoticPlugin_Hexotic::Compute(SMESH_Mesh & aMesh, SMESH_MesherHelper* aHel
   system( modeFile_Out.ToCString() );
   if ( ! fileRes.fail() ) {
     Ok = readResult( Hexotic_Out.ToCString(),
-#ifdef WITH_SMESH_CANCEL_COMPUTE
                      this,
-#endif
                      aHelper );
     if(Ok)
 /*
@@ -1710,7 +1762,7 @@ bool HexoticPlugin_Hexotic::Compute(SMESH_Mesh & aMesh, SMESH_MesherHelper* aHel
     myError->myName = COMPERR_EXCEPTION;
 */
     hexahedraMessage = "failed";
-    cout << "Problem with Hexotic output file " << Hexotic_Out << std::endl;
+    cout << "Problem with MG-Hexa output file " << Hexotic_Out << std::endl;
     // analyse log file
     SMESH_File logFile( aLogFileName.ToCString() );
     if ( !logFile.eof() )
@@ -1720,19 +1772,17 @@ bool HexoticPlugin_Hexotic::Compute(SMESH_Mesh & aMesh, SMESH_MesherHelper* aHel
       if ( std::search( fileBeg, fileEnd, msgLic, msgLic+strlen(msgLic)) != fileEnd )
         return error("Licence problems.");
     }
-    return error(SMESH_Comment("Problem with Hexotic output file ")<<Hexotic_Out);
+    return error(SMESH_Comment("Problem with MG-Hexa output file ")<<Hexotic_Out);
   }
   cout << "Hexahedra meshing " << hexahedraMessage << std::endl;
   cout << std::endl;
 
-#ifdef WITH_SMESH_CANCEL_COMPUTE
   if(_compute_canceled)
     return error(SMESH_Comment("interruption initiated by user"));
-#endif
   removeFile(Hexotic_Out);
   removeFile(Hexotic_In);
   removeFile(aLogFileName);
-  for( int i=0; i<sizeMapFiles.size(); i++)
+  for( size_t i=0; i<sizeMapFiles.size(); i++)
   {
     removeFile( TCollection_AsciiString(sizeMapFiles[i].c_str()) );
   }
@@ -1763,17 +1813,15 @@ bool HexoticPlugin_Hexotic::Evaluate(SMESH_Mesh& aMesh,
   return true;
 }
 
-#ifdef WITH_SMESH_CANCEL_COMPUTE
 void HexoticPlugin_Hexotic::CancelCompute()
 {
   _compute_canceled = true;
-#ifdef WNT
+#ifdef WIN32
 #else
-  TCollection_AsciiString aTmpDir = getTmpDir();
+  TCollection_AsciiString aTmpDir = _hexoticWorkingDirectory.c_str(); //getTmpDir();
   TCollection_AsciiString Hexotic_In = aTmpDir + "Hexotic_In.mesh";
   TCollection_AsciiString cmd = TCollection_AsciiString("ps ux | grep ") + Hexotic_In;
   cmd += TCollection_AsciiString(" | grep -v grep | awk '{print $2}' | xargs kill -9 > /dev/null 2>&1");
   system( cmd.ToCString() );
 #endif
 }
-#endif