src/SMESH_I/SMESH_2smeshpy.cxx

   1 //  SMESH SMESH_I : idl implementation based on 'SMESH' unit's calsses
   2 //
   3 //  Copyright (C) 2003  OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
   4 //  CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
   5 //
   6 //  This library is free software; you can redistribute it and/or
   7 //  modify it under the terms of the GNU Lesser General Public
   8 //  License as published by the Free Software Foundation; either
   9 //  version 2.1 of the License.
  10 //
  11 //  This library is distributed in the hope that it will be useful,
  12 //  but WITHOUT ANY WARRANTY; without even the implied warranty of
  13 //  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  14 //  Lesser General Public License for more details.
  15 //
  16 //  You should have received a copy of the GNU Lesser General Public
  17 //  License along with this library; if not, write to the Free Software
  18 //  Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
  19 //
  20 //  See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
  21 //
  22 //
  23 //
  24 //  File   : SMESH_2D_Algo_i.hxx
  25 //  Author : Paul RASCLE, EDF
  26 //  Module : SMESH
  27 //  $Header$
  28
  29 // File      : SMESH_2smeshpy.cxx
  30 // Created   : Fri Nov 18 13:20:10 2005
  31 // Author    : Edward AGAPOV (eap)
  32
  33 #include "SMESH_2smeshpy.hxx"
  34
  35 #include "utilities.h"
  36 #include "SMESH_PythonDump.hxx"
  37 #include "Resource_DataMapOfAsciiStringAsciiString.hxx"
  38
  39 #include "SMESH_Gen_i.hxx"
  40 /* SALOME headers that include CORBA headers that include windows.h
  41  * that defines GetObject symbol as GetObjectA should stand before SALOME headers
  42  * that declare methods named GetObject - to apply the same rules of GetObject renaming
  43  * and thus to avoid mess with GetObject symbol on Windows */
  44
  45 IMPLEMENT_STANDARD_HANDLE (_pyObject          ,Standard_Transient);
  46 IMPLEMENT_STANDARD_HANDLE (_pyCommand         ,Standard_Transient);
  47 IMPLEMENT_STANDARD_HANDLE (_pyGen             ,_pyObject);
  48 IMPLEMENT_STANDARD_HANDLE (_pyMesh            ,_pyObject);
  49 IMPLEMENT_STANDARD_HANDLE (_pyMeshEditor      ,_pyObject);
  50 IMPLEMENT_STANDARD_HANDLE (_pyHypothesis      ,_pyObject);
  51 IMPLEMENT_STANDARD_HANDLE (_pyAlgorithm       ,_pyHypothesis);
  52 IMPLEMENT_STANDARD_HANDLE (_pyComplexParamHypo,_pyHypothesis);
  53 IMPLEMENT_STANDARD_HANDLE (_pyNumberOfSegmentsHyp,_pyHypothesis);
  54
  55 IMPLEMENT_STANDARD_RTTIEXT(_pyObject          ,Standard_Transient);
  56 IMPLEMENT_STANDARD_RTTIEXT(_pyCommand         ,Standard_Transient);
  57 IMPLEMENT_STANDARD_RTTIEXT(_pyGen             ,_pyObject);
  58 IMPLEMENT_STANDARD_RTTIEXT(_pyMesh            ,_pyObject);
  59 IMPLEMENT_STANDARD_RTTIEXT(_pyMeshEditor      ,_pyObject);
  60 IMPLEMENT_STANDARD_RTTIEXT(_pyHypothesis      ,_pyObject);
  61 IMPLEMENT_STANDARD_RTTIEXT(_pyAlgorithm       ,_pyHypothesis);
  62 IMPLEMENT_STANDARD_RTTIEXT(_pyComplexParamHypo,_pyHypothesis);
  63 IMPLEMENT_STANDARD_RTTIEXT(_pyNumberOfSegmentsHyp,_pyHypothesis);
  64 IMPLEMENT_STANDARD_RTTIEXT(_pyLayerDistributionHypo,_pyHypothesis);
  65 IMPLEMENT_STANDARD_RTTIEXT(_pySegmentLengthAroundVertexHyp,_pyHypothesis);
  66
  67 using namespace std;
  68 using SMESH::TPythonDump;
  69
  70 /*!
  71  * \brief Container of commands into which the initial script is split.
  72  *        It also contains data coresponding to SMESH_Gen contents
  73  */
  74 static Handle(_pyGen) theGen;
  75
  76 static TCollection_AsciiString theEmptyString;
  77
  78 //#define DUMP_CONVERSION
  79
  80 #if !defined(_DEBUG_) && defined(DUMP_CONVERSION)
  81 #undef DUMP_CONVERSION
  82 #endif
  83
  84
  85 namespace {
  86
  87   //================================================================================
  88   /*!
  89    * \brief Set of TCollection_AsciiString initialized by C array of C strings
  90    */
  91   //================================================================================
  92
  93   struct TStringSet: public set<TCollection_AsciiString>
  94   {
  95     /*!
  96      * \brief Filling. The last string must be ""
  97      */
  98     void Insert(const char* names[]) {
  99       for ( int i = 0; names[i][0] ; ++i )
 100         insert( (char*) names[i] );
 101     }
 102     /*!
 103      * \brief Check if a string is in
 104      */
 105     bool Contains(const TCollection_AsciiString& name ) {
 106       return find( name ) != end();
 107     }
 108   };
 109 }
 110
 111 //================================================================================
 112 /*!
 113  * \brief Convert python script using commands of smesh.py
 114   * \param theScript - Input script
 115   * \retval TCollection_AsciiString - Convertion result
 116   *
 117   * Class SMESH_2smeshpy declared in SMESH_PythonDump.hxx
 118  */
 119 //================================================================================
 120
 121 TCollection_AsciiString
 122 SMESH_2smeshpy::ConvertScript(const TCollection_AsciiString& theScript,
 123                               Resource_DataMapOfAsciiStringAsciiString& theEntry2AccessorMethod)
 124 {
 125   theGen = new _pyGen( theEntry2AccessorMethod );
 126
 127   // split theScript into separate commands
 128   int from = 1, end = theScript.Length(), to;
 129   while ( from < end && ( to = theScript.Location( "\n", from, end )))
 130   {
 131     if ( to != from )
 132       // cut out and store a command
 133       theGen->AddCommand( theScript.SubString( from, to - 1 ));
 134     from = to + 1;
 135   }
 136   // finish conversion
 137   theGen->Flush();
 138 #ifdef DUMP_CONVERSION
 139   MESSAGE_BEGIN ( std::endl << " ######## RESULT ######## " << std::endl<< std::endl );
 140 #endif
 141   // reorder commands after conversion
 142   list< Handle(_pyCommand) >::iterator cmd;
 143   bool orderChanges;
 144   do {
 145     orderChanges = false;
 146     for ( cmd = theGen->GetCommands().begin(); cmd != theGen->GetCommands().end(); ++cmd )
 147       if ( (*cmd)->SetDependentCmdsAfter() )
 148         orderChanges = true;
 149   } while ( orderChanges );
 150
 151   // concat commands back into a script
 152   TCollection_AsciiString aScript;
 153   for ( cmd = theGen->GetCommands().begin(); cmd != theGen->GetCommands().end(); ++cmd )
 154   {
 155 #ifdef DUMP_CONVERSION
 156     MESSAGE_ADD ( "## COM " << (*cmd)->GetOrderNb() << ": "<< (*cmd)->GetString() << std::endl );
 157 #endif
 158     if ( !(*cmd)->IsEmpty() ) {
 159       aScript += "\n";
 160       aScript += (*cmd)->GetString();
 161     }
 162   }
 163   aScript += "\n";
 164
 165   theGen.Nullify();
 166
 167   return aScript;
 168 }
 169
 170 //================================================================================
 171 /*!
 172  * \brief _pyGen constructor
 173  */
 174 //================================================================================
 175
 176 _pyGen::_pyGen(Resource_DataMapOfAsciiStringAsciiString& theEntry2AccessorMethod)
 177   : _pyObject( new _pyCommand( TPythonDump::SMESHGenName(), 0 )),
 178     myID2AccessorMethod( theEntry2AccessorMethod )
 179 {
 180   myNbCommands = 0;
 181   myHasPattern = false;
 182   // make that GetID() to return TPythonDump::SMESHGenName()
 183   GetCreationCmd()->GetString() += "=";
 184 }
 185
 186 //================================================================================
 187 /*!
 188  * \brief name of SMESH_Gen in smesh.py
 189  */
 190 //================================================================================
 191
 192 const char* _pyGen::AccessorMethod() const
 193 {
 194   return SMESH_2smeshpy::GenName();
 195 }
 196
 197 //================================================================================
 198 /*!
 199  * \brief Convert a command using a specific converter
 200   * \param theCommand - the command to convert
 201  */
 202 //================================================================================
 203
 204 Handle(_pyCommand) _pyGen::AddCommand( const TCollection_AsciiString& theCommand)
 205 {
 206   // store theCommand in the sequence
 207   myCommands.push_back( new _pyCommand( theCommand, ++myNbCommands ));
 208
 209   Handle(_pyCommand) aCommand = myCommands.back();
 210 #ifdef DUMP_CONVERSION
 211   MESSAGE ( "## COM " << myNbCommands << ": "<< aCommand->GetString() );
 212 #endif
 213
 214   _pyID objID = aCommand->GetObject();
 215
 216   if ( objID.IsEmpty() )
 217     return aCommand;
 218
 219   // SMESH_Gen method?
 220   if ( objID == this->GetID() ) {
 221     this->Process( aCommand );
 222     return aCommand;
 223   }
 224   // SMESH_Mesh method?
 225   map< _pyID, Handle(_pyMesh) >::iterator id_mesh = myMeshes.find( objID );
 226   if ( id_mesh != myMeshes.end() ) {
 227     if ( aCommand->GetMethod() == "GetMeshEditor" ) { // MeshEditor creation
 228       _pyID editorID = aCommand->GetResultValue();
 229       Handle(_pyMeshEditor) editor = new _pyMeshEditor( aCommand );
 230       myMeshEditors.insert( make_pair( editorID, editor ));
 231       return aCommand;
 232     }
 233     id_mesh->second->Process( aCommand );
 234     return aCommand;
 235   }
 236   // SMESH_MeshEditor method?
 237   map< _pyID, Handle(_pyMeshEditor) >::iterator id_editor = myMeshEditors.find( objID );
 238   if ( id_editor != myMeshEditors.end() ) {
 239     id_editor->second->Process( aCommand );
 240     return aCommand;
 241   }
 242   // SMESH_Hypothesis method?
 243   list< Handle(_pyHypothesis) >::iterator hyp = myHypos.begin();
 244   for ( ; hyp != myHypos.end(); ++hyp )
 245     if ( !(*hyp)->IsAlgo() && objID == (*hyp)->GetID() ) {
 246       (*hyp)->Process( aCommand );
 247       return aCommand;
 248     }
 249
 250   // Add access to a wrapped mesh
 251   AddMeshAccessorMethod( aCommand );
 252
 253   // Add access to a wrapped algorithm
 254   //  AddAlgoAccessorMethod( aCommand ); // ??? what if algo won't be wrapped at all ???
 255
 256   // PAL12227. PythonDump was not updated at proper time; result is
 257   //     aCriteria.append(SMESH.Filter.Criterion(17,26,0,'L1',26,25,1e-07,SMESH.EDGE,-1))
 258   // TypeError: __init__() takes exactly 11 arguments (10 given)
 259   char wrongCommand[] = "SMESH.Filter.Criterion(";
 260   if ( int beg = theCommand.Location( wrongCommand, 1, theCommand.Length() ))
 261   {
 262     _pyCommand tmpCmd( theCommand.SubString( beg, theCommand.Length() ), -1);
 263     // there must be 10 arguments, 5-th arg ThresholdID is missing,
 264     const int wrongNbArgs = 9, missingArg = 5;
 265     if ( tmpCmd.GetNbArgs() == wrongNbArgs )
 266     {
 267       for ( int i = wrongNbArgs; i > missingArg; --i )
 268         tmpCmd.SetArg( i + 1, tmpCmd.GetArg( i ));
 269       tmpCmd.SetArg(  missingArg, "''");
 270       aCommand->GetString().Trunc( beg - 1 );
 271       aCommand->GetString() += tmpCmd.GetString();
 272     }
 273   }
 274   return aCommand;
 275 }
 276
 277 //================================================================================
 278 /*!
 279  * \brief Convert the command or remember it for later conversion
 280   * \param theCommand - The python command calling a method of SMESH_Gen
 281  */
 282 //================================================================================
 283
 284 void _pyGen::Process( const Handle(_pyCommand)& theCommand )
 285 {
 286   // there are methods to convert:
 287   // CreateMesh( shape )
 288   // Concatenate( [mesh1, ...], ... )
 289   // CreateHypothesis( theHypType, theLibName )
 290   // Compute( mesh, geom )
 291   // mesh creation
 292   if ( theCommand->GetMethod() == "CreateMesh" ||
 293        theCommand->GetMethod() == "CreateEmptyMesh" ||
 294        theCommand->GetMethod() == "CreateMeshesFromUNV" ||
 295        theCommand->GetMethod() == "CreateMeshesFromSTL")
 296   {
 297     Handle(_pyMesh) mesh = new _pyMesh( theCommand );
 298     myMeshes.insert( make_pair( mesh->GetID(), mesh ));
 299     return;
 300   }
 301
 302   if(theCommand->GetMethod() == "CreateMeshesFromMED")
 303   {
 304     for(int ind = 0;ind<theCommand->GetNbResultValues();ind++)
 305     {
 306       Handle(_pyMesh) mesh = new _pyMesh( theCommand, theCommand->GetResultValue(ind));
 307       myMeshes.insert( make_pair( theCommand->GetResultValue(ind), mesh ));
 308     }
 309   }
 310
 311   // CreateHypothesis()
 312   if ( theCommand->GetMethod() == "CreateHypothesis" )
 313   {
 314     myHypos.push_back( _pyHypothesis::NewHypothesis( theCommand ));
 315     return;
 316   }
 317
 318   // smeshgen.Compute( mesh, geom ) --> mesh.Compute()
 319   if ( theCommand->GetMethod() == "Compute" )
 320   {
 321     const _pyID& meshID = theCommand->GetArg( 1 );
 322     map< _pyID, Handle(_pyMesh) >::iterator id_mesh = myMeshes.find( meshID );
 323     if ( id_mesh != myMeshes.end() ) {
 324       theCommand->SetObject( meshID );
 325       theCommand->RemoveArgs();
 326       id_mesh->second->Flush();
 327       return;
 328     }
 329   }
 330
 331   // leave only one smeshgen.GetPattern() in the script
 332   if ( theCommand->GetMethod() == "GetPattern" ) {
 333     if ( myHasPattern ) {
 334       theCommand->Clear();
 335       return;
 336     }
 337     myHasPattern = true;
 338   }
 339
 340   // Concatenate( [mesh1, ...], ... )
 341   if ( theCommand->GetMethod() == "Concatenate" ||
 342        theCommand->GetMethod() == "ConcatenateWithGroups")
 343   {
 344     AddMeshAccessorMethod( theCommand );
 345   }
 346
 347   // Replace name of SMESH_Gen
 348
 349   // names of SMESH_Gen methods fully equal to methods defined in smesh.py
 350   static TStringSet smeshpyMethods;
 351   if ( smeshpyMethods.empty() ) {
 352     const char * names[] =
 353       { "SetEmbeddedMode","IsEmbeddedMode","SetCurrentStudy","GetCurrentStudy",
 354         "GetPattern","GetSubShapesId",
 355         "" }; // <- mark of array end
 356     smeshpyMethods.Insert( names );
 357   }
 358   if ( smeshpyMethods.Contains( theCommand->GetMethod() ))
 359     // smeshgen.Method() --> smesh.Method()
 360     theCommand->SetObject( SMESH_2smeshpy::SmeshpyName() );
 361   else
 362     // smeshgen.Method() --> smesh.smesh.Method()
 363     theCommand->SetObject( SMESH_2smeshpy::GenName() );
 364 }
 365
 366 //================================================================================
 367 /*!
 368  * \brief Convert the remembered commands
 369  */
 370 //================================================================================
 371
 372 void _pyGen::Flush()
 373 {
 374   map< _pyID, Handle(_pyMesh) >::iterator id_mesh = myMeshes.begin();
 375   for ( ; id_mesh != myMeshes.end(); ++id_mesh )
 376     if ( ! id_mesh->second.IsNull() )
 377       id_mesh->second->Flush();
 378
 379   list< Handle(_pyHypothesis) >::iterator hyp = myHypos.begin();
 380   for ( ; hyp != myHypos.end(); ++hyp )
 381     if ( !hyp->IsNull() ) {
 382       (*hyp)->Flush();
 383       // smeshgen.CreateHypothesis() --> smesh.smesh.CreateHypothesis()
 384       if ( !(*hyp)->IsWrapped() )
 385         (*hyp)->GetCreationCmd()->SetObject( SMESH_2smeshpy::GenName() );
 386     }
 387 }
 388
 389 //================================================================================
 390 /*!
 391  * \brief Add access method to mesh that is an argument
 392   * \param theCmd - command to add access method
 393   * \retval bool - true if added
 394  */
 395 //================================================================================
 396
 397 bool _pyGen::AddMeshAccessorMethod( Handle(_pyCommand) theCmd ) const
 398 {
 399   bool added = false;
 400   map< _pyID, Handle(_pyMesh) >::const_iterator id_mesh = myMeshes.begin();
 401   for ( ; id_mesh != myMeshes.end(); ++id_mesh ) {
 402     if ( theCmd->AddAccessorMethod( id_mesh->first, id_mesh->second->AccessorMethod() ))
 403       added = true;
 404   }
 405   return added;
 406 }
 407
 408 //================================================================================
 409 /*!
 410  * \brief Add access method to algo that is an object or an argument
 411   * \param theCmd - command to add access method
 412   * \retval bool - true if added
 413  */
 414 //================================================================================
 415
 416 bool _pyGen::AddAlgoAccessorMethod( Handle(_pyCommand) theCmd ) const
 417 {
 418   bool added = false;
 419   list< Handle(_pyHypothesis) >::const_iterator hyp = myHypos.begin();
 420   for ( ; hyp != myHypos.end(); ++hyp ) {
 421     if ( (*hyp)->IsAlgo() && /*(*hyp)->IsWrapped() &&*/
 422          theCmd->AddAccessorMethod( (*hyp)->GetID(), (*hyp)->AccessorMethod() ))
 423       added = true;
 424   }
 425   return added;
 426 }
 427
 428 //================================================================================
 429 /*!
 430  * \brief Find hypothesis by ID (entry)
 431   * \param theHypID - The hypothesis ID
 432   * \retval Handle(_pyHypothesis) - The found hypothesis
 433  */
 434 //================================================================================
 435
 436 Handle(_pyHypothesis) _pyGen::FindHyp( const _pyID& theHypID )
 437 {
 438   list< Handle(_pyHypothesis) >::iterator hyp = myHypos.begin();
 439   for ( ; hyp != myHypos.end(); ++hyp )
 440     if ( !hyp->IsNull() && theHypID == (*hyp)->GetID() )
 441       return *hyp;
 442   return Handle(_pyHypothesis)();
 443 }
 444
 445 //================================================================================
 446 /*!
 447  * \brief Find algorithm the created algorithm
 448   * \param theGeom - The shape ID the algorithm was created on
 449   * \param theMesh - The mesh ID that created the algorithm
 450   * \param dim - The algo dimension
 451   * \retval Handle(_pyHypothesis) - The found algo
 452  */
 453 //================================================================================
 454
 455 Handle(_pyHypothesis) _pyGen::FindAlgo( const _pyID& theGeom, const _pyID& theMesh,
 456                                         const Handle(_pyHypothesis)& theHypothesis )
 457 {
 458   list< Handle(_pyHypothesis) >::iterator hyp = myHypos.begin();
 459   for ( ; hyp != myHypos.end(); ++hyp )
 460     if ( !hyp->IsNull() &&
 461          (*hyp)->IsAlgo() &&
 462          theHypothesis->CanBeCreatedBy( (*hyp)->GetAlgoType() ) &&
 463          (*hyp)->GetGeom() == theGeom &&
 464          (*hyp)->GetMesh() == theMesh )
 465       return *hyp;
 466   return 0;
 467 }
 468
 469 //================================================================================
 470 /*!
 471  * \brief Change order of commands in the script
 472   * \param theCmd1 - One command
 473   * \param theCmd2 - Another command
 474  */
 475 //================================================================================
 476
 477 void _pyGen::ExchangeCommands( Handle(_pyCommand) theCmd1, Handle(_pyCommand) theCmd2 )
 478 {
 479   list< Handle(_pyCommand) >::iterator pos1, pos2;
 480   pos1 = find( myCommands.begin(), myCommands.end(), theCmd1 );
 481   pos2 = find( myCommands.begin(), myCommands.end(), theCmd2 );
 482   myCommands.insert( pos1, theCmd2 );
 483   myCommands.insert( pos2, theCmd1 );
 484   myCommands.erase( pos1 );
 485   myCommands.erase( pos2 );
 486
 487   int nb1 = theCmd1->GetOrderNb();
 488   theCmd1->SetOrderNb( theCmd2->GetOrderNb() );
 489   theCmd2->SetOrderNb( nb1 );
 490 //   cout << "BECOME " << theCmd1->GetOrderNb() << "\t" << theCmd1->GetString() << endl
 491 //        << "BECOME " << theCmd2->GetOrderNb() << "\t" << theCmd2->GetString() << endl << endl;
 492 }
 493
 494 //================================================================================
 495 /*!
 496  * \brief Set one command after the other
 497   * \param theCmd - Command to move
 498   * \param theAfterCmd - Command ater which to insert the first one
 499  */
 500 //================================================================================
 501
 502 void _pyGen::SetCommandAfter( Handle(_pyCommand) theCmd, Handle(_pyCommand) theAfterCmd )
 503 {
 504 #ifdef _DEBUG_
 505 //cout << "SET\t" << theAfterCmd->GetString() << endl << "BEFORE\t" << theCmd->GetString() << endl<<endl;
 506 #endif
 507   list< Handle(_pyCommand) >::iterator pos;
 508   pos = find( myCommands.begin(), myCommands.end(), theCmd );
 509   myCommands.erase( pos );
 510   pos = find( myCommands.begin(), myCommands.end(), theAfterCmd );
 511   myCommands.insert( ++pos, theCmd );
 512
 513   int i = 1;
 514   for ( pos = myCommands.begin(); pos != myCommands.end(); ++pos)
 515     (*pos)->SetOrderNb( i++ );
 516 }
 517
 518 //================================================================================
 519 /*!
 520  * \brief Set method to access to object wrapped with python class
 521   * \param theID - The wrapped object entry
 522   * \param theMethod - The accessor method
 523  */
 524 //================================================================================
 525
 526 void _pyGen::SetAccessorMethod(const _pyID& theID, const char* theMethod )
 527 {
 528   myID2AccessorMethod.Bind( theID, (char*) theMethod );
 529 }
 530
 531 //================================================================================
 532 /*!
 533  * \brief Find out type of geom group
 534   * \param grpID - The geom group entry
 535   * \retval int - The type
 536  */
 537 //================================================================================
 538
 539 static bool sameGroupType( const _pyID&                   grpID,
 540                            const TCollection_AsciiString& theType)
 541 {
 542   // define group type as smesh.Mesh.Group() does
 543   int type = -1;
 544   SALOMEDS::Study_var study = SMESH_Gen_i::GetSMESHGen()->GetCurrentStudy();
 545   SALOMEDS::SObject_var aSObj = study->FindObjectID( grpID.ToCString() );
 546   if ( !aSObj->_is_nil() ) {
 547     GEOM::GEOM_Object_var aGeomObj = GEOM::GEOM_Object::_narrow( aSObj->GetObject() );
 548     if ( !aGeomObj->_is_nil() ) {
 549       switch ( aGeomObj->GetShapeType() ) {
 550       case GEOM::VERTEX: type = SMESH::NODE; break;
 551       case GEOM::EDGE:   type = SMESH::EDGE; break;
 552       case GEOM::FACE:   type = SMESH::FACE; break;
 553       case GEOM::SOLID:
 554       case GEOM::SHELL:  type = SMESH::VOLUME; break;
 555       case GEOM::COMPOUND: {
 556         GEOM::GEOM_Gen_ptr aGeomGen = SMESH_Gen_i::GetSMESHGen()->GetGeomEngine();
 557         if ( !aGeomGen->_is_nil() ) {
 558           GEOM::GEOM_IGroupOperations_var aGrpOp =
 559             aGeomGen->GetIGroupOperations( study->StudyId() );
 560           if ( !aGrpOp->_is_nil() ) {
 561             switch ( aGrpOp->GetType( aGeomObj )) {
 562             case TopAbs_VERTEX: type = SMESH::NODE; break;
 563             case TopAbs_EDGE:   type = SMESH::EDGE; break;
 564             case TopAbs_FACE:   type = SMESH::FACE; break;
 565             case TopAbs_SOLID:  type = SMESH::VOLUME; break;
 566             default:;
 567             }
 568           }
 569         }
 570       }
 571       default:;
 572       }
 573     }
 574   }
 575   if ( type < 0 ) {
 576     MESSAGE("Type of the group " << grpID << " not found");
 577     return false;
 578   }
 579   if ( theType.IsIntegerValue() )
 580     return type == theType.IntegerValue();
 581
 582   switch ( type ) {
 583   case SMESH::NODE:   return theType.Location( "NODE", 1, theType.Length() );
 584   case SMESH::EDGE:   return theType.Location( "EDGE", 1, theType.Length() );
 585   case SMESH::FACE:   return theType.Location( "FACE", 1, theType.Length() );
 586   case SMESH::VOLUME: return theType.Location( "VOLUME", 1, theType.Length() );
 587   default:;
 588   }
 589   return false;
 590 }
 591
 592 //================================================================================
 593 /*!
 594  * \brief
 595   * \param theCreationCmd -
 596  */
 597 //================================================================================
 598
 599 _pyMesh::_pyMesh(const Handle(_pyCommand) theCreationCmd):
 600   _pyObject(theCreationCmd), myHasEditor(false)
 601 {
 602   // convert my creation command
 603   Handle(_pyCommand) creationCmd = GetCreationCmd();
 604   TCollection_AsciiString str = creationCmd->GetMethod();
 605
 606   creationCmd->SetObject( SMESH_2smeshpy::SmeshpyName() );
 607   if(str != "CreateMeshesFromUNV" &&
 608      str != "CreateMeshesFromMED" &&
 609      str != "CreateMeshesFromSTL")
 610     creationCmd->SetMethod( "Mesh" );
 611
 612   theGen->SetAccessorMethod( GetID(), "GetMesh()" );
 613 }
 614
 615 //================================================================================
 616 /*!
 617  * \brief
 618   * \param theCreationCmd -
 619  */
 620 //================================================================================
 621 _pyMesh::_pyMesh(const Handle(_pyCommand) theCreationCmd, const TCollection_AsciiString& id):
 622   _pyObject(theCreationCmd), myHasEditor(false)
 623 {
 624   // convert my creation command
 625   Handle(_pyCommand) creationCmd = GetCreationCmd();
 626   creationCmd->SetObject( SMESH_2smeshpy::SmeshpyName() );
 627   theGen->SetAccessorMethod( id, "GetMesh()" );
 628 }
 629
 630 //================================================================================
 631 /*!
 632  * \brief Convert a IDL API command of SMESH::Mesh to a method call of python Mesh
 633   * \param theCommand - Engine method called for this mesh
 634  */
 635 //================================================================================
 636
 637 void _pyMesh::Process( const Handle(_pyCommand)& theCommand )
 638 {
 639   // some methods of SMESH_Mesh interface needs special conversion
 640   // to methods of Mesh python class
 641   //
 642   // 1. GetSubMesh(geom, name) + AddHypothesis(geom, algo)
 643   //     --> in Mesh_Algorithm.Create(mesh, geom, hypo, so)
 644   // 2. AddHypothesis(geom, hyp)
 645   //     --> in Mesh_Algorithm.Hypothesis(hyp, args, so)
 646   // 3. CreateGroupFromGEOM(type, name, grp)
 647   //     --> in Mesh.Group(grp, name="")
 648   // 4. ExportToMED(f, auto_groups, version)
 649   //     --> in Mesh.ExportMED( f, auto_groups, version )
 650   // 5. etc
 651
 652   const TCollection_AsciiString method = theCommand->GetMethod();
 653   // ----------------------------------------------------------------------
 654   if ( method == "GetSubMesh" ) {
 655     mySubmeshes.push_back( theCommand );
 656   }
 657   // ----------------------------------------------------------------------
 658   else if ( method == "AddHypothesis" ) { // mesh.AddHypothesis(geom, HYPO )
 659     myAddHypCmds.push_back( theCommand );
 660     // set mesh to hypo
 661     const _pyID& hypID = theCommand->GetArg( 2 );
 662     Handle(_pyHypothesis) hyp = theGen->FindHyp( hypID );
 663     if ( !hyp.IsNull() ) {
 664       myHypos.push_back( hyp );
 665       if ( hyp->GetMesh().IsEmpty() )
 666         hyp->SetMesh( this->GetID() );
 667     }
 668   }
 669   // ----------------------------------------------------------------------
 670   else if ( method == "CreateGroupFromGEOM" ) {// (type, name, grp)
 671     _pyID grp = theCommand->GetArg( 3 );
 672     if ( sameGroupType( grp, theCommand->GetArg( 1 )) ) { // --> Group(grp)
 673       theCommand->SetMethod( "Group" );
 674       theCommand->RemoveArgs();
 675       theCommand->SetArg( 1, grp );
 676     }
 677     else {
 678       AddMeshAccess( theCommand );
 679     }
 680   }
 681   // ----------------------------------------------------------------------
 682   else if ( method == "ExportToMED" ) { // ExportToMED() --> ExportMED()
 683     theCommand->SetMethod( "ExportMED" );
 684   }
 685   // ----------------------------------------------------------------------
 686   else if ( method == "CreateGroup" ) { // CreateGroup() --> CreateEmptyGroup()
 687     theCommand->SetMethod( "CreateEmptyGroup" );
 688   }
 689   // ----------------------------------------------------------------------
 690   else if ( method == "RemoveHypothesis" ) // (geom, hyp)
 691   {
 692     _pyID hypID = theCommand->GetArg( 2 );
 693
 694     // check if this mesh still has corresponding addition command
 695     bool hasAddCmd = false;
 696     list< Handle(_pyCommand) >::iterator cmd = myAddHypCmds.begin();
 697     while ( cmd != myAddHypCmds.end() )
 698     {
 699       // AddHypothesis(geom, hyp)
 700       if ( hypID == (*cmd)->GetArg( 2 )) { // erase both (add and remove) commands
 701         theCommand->Clear();
 702         (*cmd)->Clear();
 703         cmd = myAddHypCmds.erase( cmd );
 704         hasAddCmd = true;
 705       }
 706       else {
 707         ++cmd;
 708       }
 709     }
 710     Handle(_pyHypothesis) hyp = theGen->FindHyp( hypID );
 711     if ( ! hasAddCmd && hypID.Length() != 0 ) { // hypo addition already wrapped
 712       // RemoveHypothesis(geom, hyp) --> RemoveHypothesis( hyp, geom=0 )
 713       _pyID geom = theCommand->GetArg( 1 );
 714       theCommand->RemoveArgs();
 715       theCommand->SetArg( 1, hypID );
 716       if ( geom != GetGeom() )
 717         theCommand->SetArg( 2, geom );
 718     }
 719     // remove hyp from myHypos
 720     myHypos.remove( hyp );
 721   }
 722   // add accessor method if necessary
 723   else
 724   {
 725     if ( NeedMeshAccess( theCommand ))
 726       // apply theCommand to the mesh wrapped by smeshpy mesh
 727       AddMeshAccess( theCommand );
 728   }
 729 }
 730
 731 //================================================================================
 732 /*!
 733  * \brief Return True if addition of accesor method is needed
 734  */
 735 //================================================================================
 736
 737 bool _pyMesh::NeedMeshAccess( const Handle(_pyCommand)& theCommand )
 738 {
 739   // names of SMESH_Mesh methods fully equal to methods of class Mesh, so
 740   // no conversion is needed for them at all:
 741   static TStringSet sameMethods;
 742   if ( sameMethods.empty() ) {
 743     const char * names[] =
 744       { "ExportDAT","ExportUNV","ExportSTL", "RemoveGroup","RemoveGroupWithContents",
 745         "GetGroups","UnionGroups","IntersectGroups","CutGroups","GetLog","GetId","ClearLog",
 746         "GetStudyId","HasDuplicatedGroupNamesMED","GetMEDMesh","NbNodes","NbElements",
 747         "NbEdges","NbEdgesOfOrder","NbFaces","NbFacesOfOrder","NbTriangles",
 748         "NbTrianglesOfOrder","NbQuadrangles","NbQuadranglesOfOrder","NbPolygons","NbVolumes",
 749         "NbVolumesOfOrder","NbTetras","NbTetrasOfOrder","NbHexas","NbHexasOfOrder",
 750         "NbPyramids","NbPyramidsOfOrder","NbPrisms","NbPrismsOfOrder","NbPolyhedrons",
 751         "NbSubMesh","GetElementsId","GetElementsByType","GetNodesId","GetElementType",
 752         "GetSubMeshElementsId","GetSubMeshNodesId","GetSubMeshElementType","Dump","GetNodeXYZ",
 753         "GetNodeInverseElements","GetShapeID","GetShapeIDForElem","GetElemNbNodes",
 754         "GetElemNode","IsMediumNode","IsMediumNodeOfAnyElem","ElemNbEdges","ElemNbFaces",
 755         "IsPoly","IsQuadratic","BaryCenter","GetHypothesisList", "SetAutoColor", "GetAutoColor",
 756         "Clear"
 757         ,"" }; // <- mark of end
 758     sameMethods.Insert( names );
 759   }
 760
 761   return !sameMethods.Contains( theCommand->GetMethod() );
 762 }
 763
 764 //================================================================================
 765 /*!
 766  * \brief Convert creation and addition of all algos and hypos
 767  */
 768 //================================================================================
 769
 770 void _pyMesh::Flush()
 771 {
 772   list < Handle(_pyCommand) >::iterator cmd, cmd2;
 773
 774   // try to convert algo addition like this:
 775   // mesh.AddHypothesis(geom, ALGO ) --> ALGO = mesh.Algo()
 776   for ( cmd = myAddHypCmds.begin(); cmd != myAddHypCmds.end(); ++cmd )
 777   {
 778     Handle(_pyCommand) addCmd = *cmd;
 779     _pyID algoID = addCmd->GetArg( 2 );
 780     Handle(_pyHypothesis) algo = theGen->FindHyp( algoID );
 781     if ( algo.IsNull() || !algo->IsAlgo() )
 782       continue;
 783     // try to convert
 784     _pyID geom = addCmd->GetArg( 1 );
 785     bool isLocalAlgo = ( geom != GetGeom() );
 786     if ( algo->Addition2Creation( addCmd, this->GetID() )) // OK
 787     {
 788       // wrapped algo is created atfer mesh creation
 789       GetCreationCmd()->AddDependantCmd( addCmd );
 790
 791       if ( isLocalAlgo ) {
 792         // mesh.AddHypothesis(geom, ALGO ) --> mesh.AlgoMethod(geom)
 793         addCmd->SetArg( addCmd->GetNbArgs() + 1,
 794                         TCollection_AsciiString( "geom=" ) + geom );
 795         // sm = mesh.GetSubMesh(geom, name) --> sm = ALGO.GetSubMesh()
 796         for ( cmd2 = mySubmeshes.begin(); cmd2 != mySubmeshes.end(); ++cmd2 ) {
 797           Handle(_pyCommand) subCmd = *cmd2;
 798           if ( geom == subCmd->GetArg( 1 )) {
 799             subCmd->SetObject( algo->GetID() );
 800             subCmd->RemoveArgs();
 801             addCmd->AddDependantCmd( subCmd );
 802           }
 803         }
 804       }
 805     }
 806     else // KO - ALGO was already created
 807     {
 808       // mesh.AddHypothesis(geom, ALGO) --> mesh.AddHypothesis(ALGO, geom=0)
 809       addCmd->RemoveArgs();
 810       addCmd->SetArg( 1, algoID );
 811       if ( isLocalAlgo )
 812         addCmd->SetArg( 2, geom );
 813     }
 814   }
 815
 816   // try to convert hypo addition like this:
 817   // mesh.AddHypothesis(geom, HYPO ) --> HYPO = algo.Hypo()
 818   for ( cmd = myAddHypCmds.begin(); cmd != myAddHypCmds.end(); ++cmd )
 819   {
 820     Handle(_pyCommand) addCmd = *cmd;
 821     _pyID hypID = addCmd->GetArg( 2 );
 822     Handle(_pyHypothesis) hyp = theGen->FindHyp( hypID );
 823     if ( hyp.IsNull() || hyp->IsAlgo() )
 824       continue;
 825     bool converted = hyp->Addition2Creation( addCmd, this->GetID() );
 826     if ( !converted ) {
 827       // mesh.AddHypothesis(geom, HYP) --> mesh.AddHypothesis(HYP, geom=0)
 828       _pyID geom = addCmd->GetArg( 1 );
 829       addCmd->RemoveArgs();
 830       addCmd->SetArg( 1, hypID );
 831       if ( geom != GetGeom() )
 832         addCmd->SetArg( 2, geom );
 833     }
 834   }
 835
 836   // sm = mesh.GetSubMesh(geom, name) --> sm = mesh.GetMesh().GetSubMesh(geom, name)
 837 //   for ( cmd = mySubmeshes.begin(); cmd != mySubmeshes.end(); ++cmd ) {
 838 //     Handle(_pyCommand) subCmd = *cmd;
 839 //     if ( subCmd->GetNbArgs() > 0 )
 840 //       AddMeshAccess( subCmd );
 841 //   }
 842   myAddHypCmds.clear();
 843   mySubmeshes.clear();
 844
 845   // flush hypotheses
 846   list< Handle(_pyHypothesis) >::iterator hyp = myHypos.begin();
 847   for ( ; hyp != myHypos.end(); ++hyp )
 848     (*hyp)->Flush();
 849 }
 850
 851 //================================================================================
 852 /*!
 853  * \brief MeshEditor convert its commands to ones of mesh
 854  */
 855 //================================================================================
 856
 857 _pyMeshEditor::_pyMeshEditor(const Handle(_pyCommand)& theCreationCmd):
 858   _pyObject( theCreationCmd )
 859 {
 860   myMesh = theCreationCmd->GetObject();
 861   myCreationCmdStr = theCreationCmd->GetString();
 862   theCreationCmd->Clear();
 863 }
 864
 865 //================================================================================
 866 /*!
 867  * \brief convert its commands to ones of mesh
 868  */
 869 //================================================================================
 870
 871 void _pyMeshEditor::Process( const Handle(_pyCommand)& theCommand)
 872 {
 873   // names of SMESH_MeshEditor methods fully equal to methods of class Mesh, so
 874   // commands calling this methods are converted to calls of methods of Mesh
 875   static TStringSet sameMethods;
 876   if ( sameMethods.empty() ) {
 877     const char * names[] = {
 878       "RemoveElements","RemoveNodes","AddNode","AddEdge","AddFace","AddPolygonalFace",
 879       "AddVolume","AddPolyhedralVolume","AddPolyhedralVolumeByFaces","MoveNode",
 880       "InverseDiag","DeleteDiag","Reorient","ReorientObject","SplitQuad","SplitQuadObject",
 881       "BestSplit","Smooth","SmoothObject","SmoothParametric","SmoothParametricObject",
 882       "ConvertToQuadratic","ConvertFromQuadratic","RenumberNodes","RenumberElements",
 883       "RotationSweep","RotationSweepObject","ExtrusionSweep","AdvancedExtrusion",
 884       "ExtrusionSweepObject","ExtrusionSweepObject1D","ExtrusionSweepObject2D","Mirror",
 885       "MirrorObject","Translate","TranslateObject","Rotate","RotateObject",
 886       "FindCoincidentNodes","FindCoincidentNodesOnPart","MergeNodes","FindEqualElements",
 887       "MergeElements","MergeEqualElements","SewFreeBorders","SewConformFreeBorders",
 888       "SewBorderToSide","SewSideElements","ChangeElemNodes","GetLastCreatedNodes",
 889       "GetLastCreatedElems",
 890       "MirrorMakeMesh","MirrorObjectMakeMesh","TranslateMakeMesh",
 891       "TranslateObjectMakeMesh","RotateMakeMesh","RotateObjectMakeMesh"
 892       ,"" }; // <- mark of the end
 893     sameMethods.Insert( names );
 894   }
 895
 896   if ( sameMethods.Contains( theCommand->GetMethod() )) {
 897     theCommand->SetObject( myMesh );
 898
 899     // meshes made by *MakeMesh() methods are not wrapped by _pyMesh,
 900     // so let _pyMesh care of it (TMP?)
 901     if ( theCommand->GetMethod().Search("MakeMesh") != -1 )
 902       _pyMesh( new _pyCommand( theCommand->GetString(), 0 )); // for theGen->SetAccessorMethod()
 903   }
 904   else {
 905     // editor creation command is needed only if any editor function is called
 906     theGen->AddMeshAccessorMethod( theCommand ); // for *Object()
 907     if ( !myCreationCmdStr.IsEmpty() ) {
 908       GetCreationCmd()->GetString() = myCreationCmdStr;
 909       myCreationCmdStr.Clear();
 910     }
 911   }
 912 }
 913
 914 //================================================================================
 915 /*!
 916  * \brief _pyHypothesis constructor
 917   * \param theCreationCmd -
 918  */
 919 //================================================================================
 920
 921 _pyHypothesis::_pyHypothesis(const Handle(_pyCommand)& theCreationCmd):
 922   _pyObject( theCreationCmd )
 923 {
 924   myIsAlgo = myIsWrapped = /*myIsConverted = myIsLocal = myDim = */false;
 925 }
 926
 927 //================================================================================
 928 /*!
 929  * \brief Creates algorithm or hypothesis
 930   * \param theCreationCmd - The engine command creating a hypothesis
 931   * \retval Handle(_pyHypothesis) - Result _pyHypothesis
 932  */
 933 //================================================================================
 934
 935 Handle(_pyHypothesis) _pyHypothesis::NewHypothesis( const Handle(_pyCommand)& theCreationCmd)
 936 {
 937   // theCreationCmd: CreateHypothesis( "theHypType", "theLibName" )
 938   ASSERT (( theCreationCmd->GetMethod() == "CreateHypothesis"));
 939
 940   Handle(_pyHypothesis) hyp, algo;
 941
 942   // "theHypType"
 943   const TCollection_AsciiString & hypTypeQuoted = theCreationCmd->GetArg( 1 );
 944   if ( hypTypeQuoted.IsEmpty() )
 945     return hyp;
 946   // theHypType
 947   TCollection_AsciiString  hypType =
 948     hypTypeQuoted.SubString( 2, hypTypeQuoted.Length() - 1 );
 949
 950   algo = new _pyAlgorithm( theCreationCmd );
 951   hyp  = new _pyHypothesis( theCreationCmd );
 952
 953   // 1D Regular_1D ----------
 954   if ( hypType == "Regular_1D" ) {
 955     // set mesh's method creating algo,
 956     // i.e. convertion result will be "regular1d = Mesh.Segment()",
 957     // and set hypType by which algo creating a hypothesis is searched for
 958     algo->SetConvMethodAndType("Segment", hypType.ToCString());
 959   }
 960   else if ( hypType == "CompositeSegment_1D" ) {
 961     algo->SetConvMethodAndType("Segment", "Regular_1D");
 962     algo->myArgs.Append( "algo=smesh.COMPOSITE");
 963   }
 964   else if ( hypType == "LocalLength" ) {
 965     // set algo's method creating hyp, and algo type
 966     hyp->SetConvMethodAndType( "LocalLength", "Regular_1D");
 967     // set method whose 1 arg will become the 1-st arg of hyp creation command
 968     // i.e. convertion result will be "locallength = regular1d.LocalLength(<arg of SetLength()>)"
 969     hyp->AddArgMethod( "SetLength" );
 970   }
 971   else if ( hypType == "NumberOfSegments" ) {
 972     hyp = new _pyNumberOfSegmentsHyp( theCreationCmd );
 973     hyp->SetConvMethodAndType( "NumberOfSegments", "Regular_1D");
 974     // arg of SetNumberOfSegments() will become the 1-st arg of hyp creation command
 975     hyp->AddArgMethod( "SetNumberOfSegments" );
 976     // arg of SetScaleFactor() will become the 2-nd arg of hyp creation command
 977     hyp->AddArgMethod( "SetScaleFactor" );
 978   }
 979   else if ( hypType == "Arithmetic1D" ) {
 980     hyp = new _pyComplexParamHypo( theCreationCmd );
 981     hyp->SetConvMethodAndType( "Arithmetic1D", "Regular_1D");
 982   }
 983   else if ( hypType == "StartEndLength" ) {
 984     hyp = new _pyComplexParamHypo( theCreationCmd );
 985     hyp->SetConvMethodAndType( "StartEndLength", "Regular_1D");
 986   }
 987   else if ( hypType == "Deflection1D" ) {
 988     hyp->SetConvMethodAndType( "Deflection1D", "Regular_1D");
 989     hyp->AddArgMethod( "SetDeflection" );
 990   }
 991   else if ( hypType == "Propagation" ) {
 992     hyp->SetConvMethodAndType( "Propagation", "Regular_1D");
 993   }
 994   else if ( hypType == "QuadraticMesh" ) {
 995     hyp->SetConvMethodAndType( "QuadraticMesh", "Regular_1D");
 996   }
 997   else if ( hypType == "AutomaticLength" ) {
 998     hyp->SetConvMethodAndType( "AutomaticLength", "Regular_1D");
 999     hyp->AddArgMethod( "SetFineness");
1000   }
1001   else if ( hypType == "SegmentLengthAroundVertex" ) {
1002     hyp = new _pySegmentLengthAroundVertexHyp( theCreationCmd );
1003     hyp->SetConvMethodAndType( "LengthNearVertex", "Regular_1D" );
1004     hyp->AddArgMethod( "SetLength" );
1005   }
1006   // 1D Python_1D ----------
1007   else if ( hypType == "Python_1D" ) {
1008     algo->SetConvMethodAndType( "Segment", hypType.ToCString());
1009     algo->myArgs.Append( "algo=smesh.PYTHON");
1010   }
1011   else if ( hypType == "PythonSplit1D" ) {
1012     hyp->SetConvMethodAndType( "PythonSplit1D", "Python_1D");
1013     hyp->AddArgMethod( "SetNumberOfSegments");
1014     hyp->AddArgMethod( "SetPythonLog10RatioFunction");
1015   }
1016   // MEFISTO_2D ----------
1017   else if ( hypType == "MEFISTO_2D" ) { // MEFISTO_2D
1018     algo->SetConvMethodAndType( "Triangle", hypType.ToCString());
1019   }
1020   else if ( hypType == "MaxElementArea" ) {
1021     hyp->SetConvMethodAndType( "MaxElementArea", "MEFISTO_2D");
1022     hyp->SetConvMethodAndType( "MaxElementArea", "NETGEN_2D_ONLY");
1023     hyp->AddArgMethod( "SetMaxElementArea");
1024   }
1025   else if ( hypType == "LengthFromEdges" ) {
1026     hyp->SetConvMethodAndType( "LengthFromEdges", "MEFISTO_2D");
1027     hyp->SetConvMethodAndType( "LengthFromEdges", "NETGEN_2D_ONLY");
1028   }
1029   // Quadrangle_2D ----------
1030   else if ( hypType == "Quadrangle_2D" ) {
1031     algo->SetConvMethodAndType( "Quadrangle" , hypType.ToCString());
1032   }
1033   else if ( hypType == "QuadranglePreference" ) {
1034     hyp->SetConvMethodAndType( "QuadranglePreference", "Quadrangle_2D");
1035     hyp->SetConvMethodAndType( "QuadranglePreference", "NETGEN_2D_ONLY");
1036   }
1037   // NETGEN ----------
1038 //   else if ( hypType == "NETGEN_2D") { // 1D-2D
1039 //     algo->SetConvMethodAndType( "Triangle" , hypType.ToCString());
1040 //     algo->myArgs.Append( "algo=smesh.NETGEN" );
1041 //   }
1042   else if ( hypType == "NETGEN_2D_ONLY") { // 2D
1043     algo->SetConvMethodAndType( "Triangle" , hypType.ToCString());
1044     algo->myArgs.Append( "algo=smesh.NETGEN_2D" );
1045   }
1046   else if ( hypType == "NETGEN_3D") { // 3D
1047     algo->SetConvMethodAndType( "Tetrahedron" , hypType.ToCString());
1048     algo->myArgs.Append( "algo=smesh.NETGEN" );
1049   }
1050   else if ( hypType == "MaxElementVolume") {
1051     hyp->SetConvMethodAndType( "MaxElementVolume", "NETGEN_3D");
1052     hyp->AddArgMethod( "SetMaxElementVolume" );
1053   }
1054   // GHS3D_3D ----------
1055   else if ( hypType == "GHS3D_3D" ) {
1056     algo->SetConvMethodAndType( "Tetrahedron", hypType.ToCString());
1057     algo->myArgs.Append( "algo=smesh.GHS3D" );
1058   }
1059   // Hexa_3D ---------
1060   else if ( hypType == "Hexa_3D" ) {
1061     algo->SetConvMethodAndType( "Hexahedron", hypType.ToCString());
1062   }
1063   // Repetitive Projection_1D ---------
1064   else if ( hypType == "Projection_1D" ) {
1065     algo->SetConvMethodAndType( "Projection1D", hypType.ToCString());
1066   }
1067   else if ( hypType == "ProjectionSource1D" ) {
1068     hyp->SetConvMethodAndType( "SourceEdge", "Projection_1D");
1069     hyp->AddArgMethod( "SetSourceEdge");
1070     hyp->AddArgMethod( "SetSourceMesh");
1071     // 2 args of SetVertexAssociation() will become the 3-th and 4-th args of hyp creation command
1072     hyp->AddArgMethod( "SetVertexAssociation", 2 );
1073   }
1074   // Projection_2D ---------
1075   else if ( hypType == "Projection_2D" ) {
1076     algo->SetConvMethodAndType( "Projection2D", hypType.ToCString());
1077   }
1078   else if ( hypType == "ProjectionSource2D" ) {
1079     hyp->SetConvMethodAndType( "SourceFace", "Projection_2D");
1080     hyp->AddArgMethod( "SetSourceFace");
1081     hyp->AddArgMethod( "SetSourceMesh");
1082     hyp->AddArgMethod( "SetVertexAssociation", 4 );
1083   }
1084   // Projection_3D ---------
1085   else if ( hypType == "Projection_3D" ) {
1086     algo->SetConvMethodAndType( "Projection3D", hypType.ToCString());
1087   }
1088   else if ( hypType == "ProjectionSource3D" ) {
1089     hyp->SetConvMethodAndType( "SourceShape3D", "Projection_3D");
1090     hyp->AddArgMethod( "SetSource3DShape");
1091     hyp->AddArgMethod( "SetSourceMesh");
1092     hyp->AddArgMethod( "SetVertexAssociation", 4 );
1093   }
1094   // Prism_3D ---------
1095   else if ( hypType == "Prism_3D" ) {
1096     algo->SetConvMethodAndType( "Prism", hypType.ToCString());
1097   }
1098   // RadialPrism_3D ---------
1099   else if ( hypType == "RadialPrism_3D" ) {
1100     algo->SetConvMethodAndType( "Prism", hypType.ToCString());
1101   }
1102   else if ( hypType == "NumberOfLayers" ) {
1103     hyp->SetConvMethodAndType( "NumberOfLayers", "RadialPrism_3D");
1104     hyp->AddArgMethod( "SetNumberOfLayers" );
1105   }
1106   else if ( hypType == "LayerDistribution" ) {
1107     hyp = new _pyLayerDistributionHypo( theCreationCmd );
1108     hyp->SetConvMethodAndType( "LayerDistribution", "RadialPrism_3D");
1109   }
1110
1111   if ( algo->IsValid() ) {
1112     return algo;
1113   }
1114   return hyp;
1115 }
1116
1117 //================================================================================
1118 /*!
1119  * \brief Convert the command adding a hypothesis to mesh into a smesh command
1120   * \param theCmd - The command like mesh.AddHypothesis( geom, hypo )
1121   * \param theAlgo - The algo that can create this hypo
1122   * \retval bool - false if the command cant be converted
1123  */
1124 //================================================================================
1125
1126 bool _pyHypothesis::Addition2Creation( const Handle(_pyCommand)& theCmd,
1127                                        const _pyID&              theMesh)
1128 {
1129   ASSERT(( theCmd->GetMethod() == "AddHypothesis" ));
1130
1131   if ( !IsWrappable( theMesh ))
1132     return false;
1133
1134   myGeom = theCmd->GetArg( 1 );
1135
1136   Handle(_pyHypothesis) algo;
1137   if ( !IsAlgo() ) {
1138     // find algo created on myGeom in theMesh
1139     algo = theGen->FindAlgo( myGeom, theMesh, this );
1140     if ( algo.IsNull() )
1141       return false;
1142     algo->GetCreationCmd()->AddDependantCmd( theCmd );
1143   }
1144   myIsWrapped = true;
1145
1146   // mesh.AddHypothesis(geom,hyp) --> hyp = <theMesh or algo>.myCreationMethod(args)
1147   theCmd->SetResultValue( GetID() );
1148   theCmd->SetObject( IsAlgo() ? theMesh : algo->GetID());
1149   theCmd->SetMethod( IsAlgo() ? GetAlgoCreationMethod() : GetCreationMethod( algo->GetAlgoType() ));
1150   // set args
1151   theCmd->RemoveArgs();
1152   for ( int i = 1; i <= myArgs.Length(); ++i ) {
1153     if ( !myArgs( i ).IsEmpty() )
1154       theCmd->SetArg( i, myArgs( i ));
1155     else
1156       theCmd->SetArg( i, "[]");
1157   }
1158   // set a new creation command
1159   GetCreationCmd()->Clear();
1160   SetCreationCmd( theCmd );
1161
1162   // clear commands setting arg values
1163   list < Handle(_pyCommand) >::iterator argCmd = myArgCommands.begin();
1164   for ( ; argCmd != myArgCommands.end(); ++argCmd )
1165     (*argCmd)->Clear();
1166
1167   // set unknown arg commands after hypo creation
1168   Handle(_pyCommand) afterCmd = myIsWrapped ? theCmd : GetCreationCmd();
1169   list<Handle(_pyCommand)>::iterator cmd = myUnknownCommands.begin();
1170   for ( ; cmd != myUnknownCommands.end(); ++cmd ) {
1171     afterCmd->AddDependantCmd( *cmd );
1172   }
1173
1174   return myIsWrapped;
1175 }
1176
1177 //================================================================================
1178 /*!
1179  * \brief Remember hypothesis parameter values
1180  * \param theCommand - The called hypothesis method
1181  */
1182 //================================================================================
1183
1184 void _pyHypothesis::Process( const Handle(_pyCommand)& theCommand)
1185 {
1186   ASSERT( !myIsAlgo );
1187   // set args
1188   int nbArgs = 0;
1189   for ( int i = 1; i <= myArgMethods.Length(); ++i ) {
1190     if ( myArgMethods( i ) == theCommand->GetMethod() ) {
1191       while ( myArgs.Length() < nbArgs + myNbArgsByMethod( i ))
1192         myArgs.Append( "[]" );
1193       for ( int iArg = 1; iArg <= myNbArgsByMethod( i ); ++iArg )
1194         myArgs( nbArgs + iArg ) = theCommand->GetArg( iArg ); // arg value
1195       myArgCommands.push_back( theCommand );
1196       return;
1197     }
1198     nbArgs += myNbArgsByMethod( i );
1199   }
1200   myUnknownCommands.push_back( theCommand );
1201 }
1202
1203 //================================================================================
1204 /*!
1205  * \brief Finish conversion
1206  */
1207 //================================================================================
1208
1209 void _pyHypothesis::Flush()
1210 {
1211   if ( IsWrapped() ) {
1212   }
1213   else {
1214     list < Handle(_pyCommand) >::iterator cmd = myArgCommands.begin();
1215     for ( ; cmd != myArgCommands.end(); ++cmd ) {
1216       // Add access to a wrapped mesh
1217       theGen->AddMeshAccessorMethod( *cmd );
1218       // Add access to a wrapped algorithm
1219       theGen->AddAlgoAccessorMethod( *cmd );
1220     }
1221     cmd = myUnknownCommands.begin();
1222     for ( ; cmd != myUnknownCommands.end(); ++cmd ) {
1223       // Add access to a wrapped mesh
1224       theGen->AddMeshAccessorMethod( *cmd );
1225       // Add access to a wrapped algorithm
1226       theGen->AddAlgoAccessorMethod( *cmd );
1227     }
1228   }
1229   // forget previous hypothesis modifications
1230   myArgCommands.clear();
1231   myUnknownCommands.clear();
1232 }
1233
1234 //================================================================================
1235 /*!
1236  * \brief clear creation, arg and unkown commands
1237  */
1238 //================================================================================
1239
1240 void _pyHypothesis::ClearAllCommands()
1241 {
1242   GetCreationCmd()->Clear();
1243   list<Handle(_pyCommand)>::iterator cmd = myArgCommands.begin();
1244   for ( ; cmd != myArgCommands.end(); ++cmd )
1245     ( *cmd )->Clear();
1246   cmd = myUnknownCommands.begin();
1247   for ( ; cmd != myUnknownCommands.end(); ++cmd )
1248     ( *cmd )->Clear();
1249 }
1250
1251 //================================================================================
1252 /*!
1253  * \brief Remember hypothesis parameter values
1254   * \param theCommand - The called hypothesis method
1255  */
1256 //================================================================================
1257
1258 void _pyComplexParamHypo::Process( const Handle(_pyCommand)& theCommand)
1259 {
1260   // ex: hyp.SetLength(start, 1)
1261   //     hyp.SetLength(end,   0)
1262   ASSERT(( theCommand->GetMethod() == "SetLength" ));
1263   ASSERT(( theCommand->GetArg( 2 ).IsIntegerValue() ));
1264   int i = 2 - theCommand->GetArg( 2 ).IntegerValue();
1265   while ( myArgs.Length() < i )
1266     myArgs.Append( "[]" );
1267   myArgs( i ) = theCommand->GetArg( 1 ); // arg value
1268   myArgCommands.push_back( theCommand );
1269 }
1270
1271 //================================================================================
1272 /*!
1273  * \brief Convert methods of 1D hypotheses to my own methods
1274   * \param theCommand - The called hypothesis method
1275  */
1276 //================================================================================
1277
1278 void _pyLayerDistributionHypo::Process( const Handle(_pyCommand)& theCommand)
1279 {
1280   if ( theCommand->GetMethod() != "SetLayerDistribution" )
1281     return;
1282
1283   _pyID newName; // name for 1D hyp = "HypType" + "_Distribution"
1284
1285   const _pyID& hyp1dID = theCommand->GetArg( 1 );
1286   Handle(_pyHypothesis) hyp1d = theGen->FindHyp( hyp1dID );
1287   if ( hyp1d.IsNull() ) // apparently hypId changed at study restoration
1288     hyp1d = my1dHyp;
1289   else if ( !my1dHyp.IsNull() && hyp1dID != my1dHyp->GetID() ) {
1290     // 1D hypo is already set, so distribution changes and the old
1291     // 1D hypo is thrown away
1292     my1dHyp->ClearAllCommands();
1293   }
1294   my1dHyp = hyp1d;
1295   if ( my1dHyp.IsNull() )
1296     return; // something wrong :(
1297
1298   // make a new name for 1D hyp = "HypType" + "_Distribution"
1299   if ( my1dHyp->GetCreationCmd()->GetMethod() == "CreateHypothesis" ) {
1300     // not yet converted creation cmd
1301     TCollection_AsciiString hypTypeQuoted = my1dHyp->GetCreationCmd()->GetArg(1);
1302     TCollection_AsciiString hypType = hypTypeQuoted.SubString( 2, hypTypeQuoted.Length() - 1 );
1303     newName = hypType + "_Distribution";
1304     my1dHyp->GetCreationCmd()->SetResultValue( newName );
1305   }
1306   else {
1307     // already converted creation cmd
1308     newName = my1dHyp->GetCreationCmd()->GetResultValue();
1309   }
1310
1311   // as creation of 1D hyp was written later then it's edition,
1312   // we need to find all it's edition calls and process them
1313   list< Handle(_pyCommand) >& cmds = theGen->GetCommands();
1314   list< Handle(_pyCommand) >::iterator cmdIt = cmds.begin();
1315   for ( ; cmdIt != cmds.end(); ++cmdIt ) {
1316     const _pyID& objID = (*cmdIt)->GetObject();
1317     if ( objID == hyp1dID ) {
1318       my1dHyp->Process( *cmdIt );
1319       my1dHyp->GetCreationCmd()->AddDependantCmd( *cmdIt );
1320       ( *cmdIt )->SetObject( newName );
1321     }
1322   }
1323   if ( !myArgCommands.empty() )
1324     myArgCommands.front()->Clear();
1325   theCommand->SetArg( 1, newName );
1326   myArgCommands.push_back( theCommand );
1327   // copy hyp1d's creation method and args
1328 //   myCreationMethod = hyp1d->GetCreationMethod();
1329 //   myArgs           = hyp1d->GetArgs();
1330 //   // make them cleared at conversion
1331 //   myArgCommands = hyp1d->GetArgCommands();
1332
1333 //   // to be cleared at convertion only
1334 //   myArgCommands.push_back( theCommand );
1335 }
1336
1337 //================================================================================
1338 /*!
1339  * \brief
1340   * \param theAdditionCmd -
1341   * \param theMesh -
1342   * \retval bool -
1343  */
1344 //================================================================================
1345
1346 bool _pyLayerDistributionHypo::Addition2Creation( const Handle(_pyCommand)& theAdditionCmd,
1347                                                   const _pyID&              theMesh)
1348 {
1349   myIsWrapped = false;
1350
1351   if ( my1dHyp.IsNull() )
1352     return false;
1353
1354   // set "SetLayerDistribution()" after addition cmd
1355   theAdditionCmd->AddDependantCmd( myArgCommands.front() );
1356
1357   _pyID geom = theAdditionCmd->GetArg( 1 );
1358
1359   my1dHyp->SetMesh( theMesh );
1360   if ( !my1dHyp->Addition2Creation( theAdditionCmd, theMesh ))
1361     return false;
1362
1363   // clear "SetLayerDistribution()" cmd
1364   myArgCommands.front()->Clear();
1365
1366   // Convert my creation => me = RadialPrismAlgo.Get3DHypothesis()
1367
1368   // find RadialPrism algo created on <geom> for theMesh
1369   Handle(_pyHypothesis) algo = theGen->FindAlgo( geom, theMesh, this );
1370   if ( !algo.IsNull() ) {
1371     GetCreationCmd()->SetObject( algo->GetID() );
1372     GetCreationCmd()->SetMethod( "Get3DHypothesis" );
1373     GetCreationCmd()->RemoveArgs();
1374     theAdditionCmd->AddDependantCmd( GetCreationCmd() );
1375     myIsWrapped = true;
1376   }
1377   return myIsWrapped;
1378 }
1379
1380 //================================================================================
1381 /*!
1382  * \brief
1383  */
1384 //================================================================================
1385
1386 void _pyLayerDistributionHypo::Flush()
1387 {
1388   //my1dHyp.Nullify();
1389   //_pyHypothesis::Flush();
1390 }
1391
1392 //================================================================================
1393 /*!
1394  * \brief additionally to Addition2Creation, clears SetDistrType() command
1395   * \param theCmd - AddHypothesis() command
1396   * \param theMesh - mesh to which a hypothesis is added
1397   * \retval bool - convertion result
1398  */
1399 //================================================================================
1400
1401 bool _pyNumberOfSegmentsHyp::Addition2Creation( const Handle(_pyCommand)& theCmd,
1402                                                 const _pyID&              theMesh)
1403 {
1404   if ( IsWrappable( theMesh ) && myArgs.Length() > 1 ) {
1405     // scale factor (2-nd arg) is provided: clear SetDistrType(1) command
1406     bool scaleDistrType = false;
1407     list<Handle(_pyCommand)>::reverse_iterator cmd = myUnknownCommands.rbegin();
1408     for ( ; cmd != myUnknownCommands.rend(); ++cmd ) {
1409       if ( (*cmd)->GetMethod() == "SetDistrType" ) {
1410         if ( (*cmd)->GetArg( 1 ) == "1" ) {
1411           scaleDistrType = true;
1412           (*cmd)->Clear();
1413         }
1414         else if ( !scaleDistrType ) {
1415           // distribution type changed: remove scale factor from args
1416           myArgs.Remove( 2, myArgs.Length() );
1417           break;
1418         }
1419       }
1420     }
1421   }
1422   return _pyHypothesis::Addition2Creation( theCmd, theMesh );
1423 }
1424
1425 //================================================================================
1426 /*!
1427  * \brief remove repeated commands defining distribution
1428  */
1429 //================================================================================
1430
1431 void _pyNumberOfSegmentsHyp::Flush()
1432 {
1433   // find number of the last SetDistrType() command
1434   list<Handle(_pyCommand)>::reverse_iterator cmd = myUnknownCommands.rbegin();
1435   int distrTypeNb = 0;
1436   for ( ; !distrTypeNb && cmd != myUnknownCommands.rend(); ++cmd )
1437     if ( (*cmd)->GetMethod() == "SetDistrType" )
1438       distrTypeNb = (*cmd)->GetOrderNb();
1439
1440   // clear commands before the last SetDistrType()
1441   list<Handle(_pyCommand)> * cmds[2] = { &myArgCommands, &myUnknownCommands };
1442   for ( int i = 0; i < 2; ++i ) {
1443     set<TCollection_AsciiString> uniqueMethods;
1444     list<Handle(_pyCommand)> & cmdList = *cmds[i];
1445     for ( cmd = cmdList.rbegin(); cmd != cmdList.rend(); ++cmd )
1446     {
1447       bool clear = ( (*cmd)->GetOrderNb() < distrTypeNb );
1448       const TCollection_AsciiString& method = (*cmd)->GetMethod();
1449       if ( !clear || method == "SetNumberOfSegments" ) {
1450         bool isNewInSet = uniqueMethods.insert( method ).second;
1451         clear = !isNewInSet;
1452       }
1453       if ( clear )
1454         (*cmd)->Clear();
1455     }
1456     cmdList.clear();
1457   }
1458 }
1459
1460 //================================================================================
1461 /*!
1462  * \brief Convert the command adding "SegmentLengthAroundVertex" to mesh
1463  * into regular1D.LengthNearVertex( length, vertex )
1464   * \param theCmd - The command like mesh.AddHypothesis( vertex, SegmentLengthAroundVertex )
1465   * \param theMesh - The mesh needing this hypo
1466   * \retval bool - false if the command cant be converted
1467  */
1468 //================================================================================
1469
1470 bool _pySegmentLengthAroundVertexHyp::Addition2Creation( const Handle(_pyCommand)& theCmd,
1471                                                          const _pyID&              theMeshID)
1472 {
1473   if ( IsWrappable( theMeshID )) {
1474
1475     _pyID vertex = theCmd->GetArg( 1 );
1476
1477     // the problem here is that segment algo will not be found
1478     // by pyHypothesis::Addition2Creation() for <vertex>, so we try to find
1479     // geometry where segment algorithm is assigned
1480     Handle(_pyHypothesis) algo;
1481     _pyID geom = vertex;
1482     while ( algo.IsNull() && !geom.IsEmpty()) {
1483       // try to find geom as a father of <vertex>
1484       geom = FatherID( geom );
1485       algo = theGen->FindAlgo( geom, theMeshID, this );
1486     }
1487     if ( algo.IsNull() )
1488       return false; // also possible to find geom as brother of veretex...
1489     // set geom instead of vertex
1490     theCmd->SetArg( 1, geom );
1491
1492     // set vertex as a second arg
1493     if ( myArgs.Length() < 1) myArgs.Append( "1" ); // :(
1494     myArgs.Append( vertex );
1495
1496     // mesh.AddHypothesis(vertex, SegmentLengthAroundVertex) -->
1497     // theMeshID.LengthNearVertex( length, vertex )
1498     return _pyHypothesis::Addition2Creation( theCmd, theMeshID );
1499   }
1500   return false;
1501 }
1502
1503 //================================================================================
1504 /*!
1505  * \brief _pyAlgorithm constructor
1506  * \param theCreationCmd - The command like "algo = smeshgen.CreateHypothesis(type,lib)"
1507  */
1508 //================================================================================
1509
1510 _pyAlgorithm::_pyAlgorithm(const Handle(_pyCommand)& theCreationCmd)
1511   : _pyHypothesis( theCreationCmd )
1512 {
1513   myIsAlgo = true;
1514 }
1515
1516 //================================================================================
1517 /*!
1518  * \brief Convert the command adding an algorithm to mesh
1519   * \param theCmd - The command like mesh.AddHypothesis( geom, algo )
1520   * \param theMesh - The mesh needing this algo
1521   * \retval bool - false if the command cant be converted
1522  */
1523 //================================================================================
1524
1525 bool _pyAlgorithm::Addition2Creation( const Handle(_pyCommand)& theCmd,
1526                                       const _pyID&              theMeshID)
1527 {
1528   // mesh.AddHypothesis(geom,algo) --> theMeshID.myCreationMethod()
1529   if ( _pyHypothesis::Addition2Creation( theCmd, theMeshID )) {
1530     theGen->SetAccessorMethod( GetID(), "GetAlgorithm()" );
1531     return true;
1532   }
1533   return false;
1534 }
1535
1536 //================================================================================
1537 /*!
1538  * \brief Return starting position of a part of python command
1539   * \param thePartIndex - The index of command part
1540   * \retval int - Part position
1541  */
1542 //================================================================================
1543
1544 int _pyCommand::GetBegPos( int thePartIndex )
1545 {
1546   if ( IsEmpty() )
1547     return EMPTY;
1548   if ( myBegPos.Length() < thePartIndex )
1549     return UNKNOWN;
1550   return myBegPos( thePartIndex );
1551 }
1552
1553 //================================================================================
1554 /*!
1555  * \brief Store starting position of a part of python command
1556   * \param thePartIndex - The index of command part
1557   * \param thePosition - Part position
1558  */
1559 //================================================================================
1560
1561 void _pyCommand::SetBegPos( int thePartIndex, int thePosition )
1562 {
1563   while ( myBegPos.Length() < thePartIndex )
1564     myBegPos.Append( UNKNOWN );
1565   myBegPos( thePartIndex ) = thePosition;
1566 }
1567
1568 //================================================================================
1569 /*!
1570  * \brief Returns whitespace symbols at the line beginning
1571   * \retval TCollection_AsciiString - result
1572  */
1573 //================================================================================
1574
1575 TCollection_AsciiString _pyCommand::GetIndentation()
1576 {
1577   int end = 1;
1578   if ( GetBegPos( RESULT_IND ) == UNKNOWN )
1579     GetWord( myString, end, true );
1580   else
1581     end = GetBegPos( RESULT_IND );
1582   return myString.SubString( 1, end - 1 );
1583 }
1584
1585 //================================================================================
1586 /*!
1587  * \brief Return substring of python command looking like ResultValue = Obj.Meth()
1588   * \retval const TCollection_AsciiString & - ResultValue substring
1589  */
1590 //================================================================================
1591
1592 const TCollection_AsciiString & _pyCommand::GetResultValue()
1593 {
1594   if ( GetBegPos( RESULT_IND ) == UNKNOWN )
1595   {
1596     int begPos = myString.Location( "=", 1, Length() );
1597     if ( begPos )
1598       myRes = GetWord( myString, begPos, false );
1599     else
1600       begPos = EMPTY;
1601     SetBegPos( RESULT_IND, begPos );
1602   }
1603   return myRes;
1604 }
1605
1606 //================================================================================
1607 /*!
1608  * \brief Return number of python command result value ResultValue = Obj.Meth()
1609   * \retval const int
1610  */
1611 //================================================================================
1612
1613 const int _pyCommand::GetNbResultValues()
1614 {
1615   int begPos = 1;
1616   int Nb=0;
1617   int endPos = myString.Location( "=", 1, Length() );
1618   TCollection_AsciiString str = "";
1619   while ( begPos < endPos) {
1620     str = GetWord( myString, begPos, true );
1621     begPos = begPos+ str.Length();
1622     Nb++;
1623   }
1624   return (Nb-1);
1625 }
1626
1627
1628 //================================================================================
1629 /*!
1630  * \brief Return substring of python command looking like
1631  *  ResultValue1 , ResultValue1,... = Obj.Meth() with res index
1632  * \retval const TCollection_AsciiString & - ResultValue with res index substring
1633  */
1634 //================================================================================
1635 const TCollection_AsciiString & _pyCommand::GetResultValue(int res)
1636 {
1637   int begPos = 1;
1638   int Nb=0;
1639   int endPos = myString.Location( "=", 1, Length() );
1640   while ( begPos < endPos) {
1641     myRes = GetWord( myString, begPos, true );
1642     begPos = begPos + myRes.Length();
1643     Nb++;
1644     if(res == Nb){
1645       myRes.RemoveAll('[');myRes.RemoveAll(']');
1646       return myRes;
1647     }
1648     if(Nb>res)
1649       break;
1650   }
1651   return theEmptyString;
1652 }
1653
1654 //================================================================================
1655 /*!
1656  * \brief Return substring of python command looking like ResVal = Object.Meth()
1657   * \retval const TCollection_AsciiString & - Object substring
1658  */
1659 //================================================================================
1660
1661 const TCollection_AsciiString & _pyCommand::GetObject()
1662 {
1663   if ( GetBegPos( OBJECT_IND ) == UNKNOWN )
1664   {
1665     // beginning
1666     int begPos = GetBegPos( RESULT_IND ) + myRes.Length();
1667     if ( begPos < 1 )
1668       begPos = myString.Location( "=", 1, Length() ) + 1;
1669     // store
1670     myObj = GetWord( myString, begPos, true );
1671     SetBegPos( OBJECT_IND, begPos );
1672   }
1673   //SCRUTE(myObj);
1674   return myObj;
1675 }
1676
1677 //================================================================================
1678 /*!
1679  * \brief Return substring of python command looking like ResVal = Obj.Method()
1680   * \retval const TCollection_AsciiString & - Method substring
1681  */
1682 //================================================================================
1683
1684 const TCollection_AsciiString & _pyCommand::GetMethod()
1685 {
1686   if ( GetBegPos( METHOD_IND ) == UNKNOWN )
1687   {
1688     // beginning
1689     int begPos = GetBegPos( OBJECT_IND ) + myObj.Length();
1690     bool forward = true;
1691     if ( begPos < 1 ) {
1692       begPos = myString.Location( "(", 1, Length() ) - 1;
1693       forward = false;
1694     }
1695     // store
1696     myMeth = GetWord( myString, begPos, forward );
1697     SetBegPos( METHOD_IND, begPos );
1698   }
1699   //SCRUTE(myMeth);
1700   return myMeth;
1701 }
1702
1703 //================================================================================
1704 /*!
1705  * \brief Return substring of python command looking like ResVal = Obj.Meth(Arg1,...)
1706   * \retval const TCollection_AsciiString & - Arg<index> substring
1707  */
1708 //================================================================================
1709
1710 const TCollection_AsciiString & _pyCommand::GetArg( int index )
1711 {
1712   if ( GetBegPos( ARG1_IND ) == UNKNOWN )
1713   {
1714     // find all args
1715     int begPos = GetBegPos( METHOD_IND ) + myMeth.Length();
1716     if ( begPos < 1 )
1717       begPos = myString.Location( "(", 1, Length() ) + 1;
1718
1719     int i = 0, prevLen = 0;
1720     while ( begPos != EMPTY ) {
1721       begPos += prevLen;
1722       // check if we are looking at the closing parenthesis
1723       while ( begPos <= Length() && isspace( myString.Value( begPos )))
1724         ++begPos;
1725       if ( begPos > Length() || myString.Value( begPos ) == ')' )
1726         break;
1727       myArgs.Append( GetWord( myString, begPos, true, true ));
1728       SetBegPos( ARG1_IND + i, begPos );
1729       prevLen = myArgs.Last().Length();
1730       if ( prevLen == 0 )
1731         myArgs.Remove( myArgs.Length() ); // no more args
1732       i++;
1733     }
1734   }
1735   if ( myArgs.Length() < index )
1736     return theEmptyString;
1737   return myArgs( index );
1738 }
1739
1740 //================================================================================
1741 /*!
1742  * \brief Check if char is a word part
1743   * \param c - The character to check
1744   * \retval bool - The check result
1745  */
1746 //================================================================================
1747
1748 static inline bool isWord(const char c, const bool dotIsWord)
1749 {
1750   return
1751     !isspace(c) && c != ',' && c != '=' && c != ')' && c != '(' && ( dotIsWord || c != '.');
1752 }
1753
1754 //================================================================================
1755 /*!
1756  * \brief Looks for a word in the string and returns word's beginning
1757   * \param theString - The input string
1758   * \param theStartPos - The position to start the search, returning word's beginning
1759   * \param theForward - The search direction
1760   * \retval TCollection_AsciiString - The found word
1761  */
1762 //================================================================================
1763
1764 TCollection_AsciiString _pyCommand::GetWord( const TCollection_AsciiString & theString,
1765                                             int &      theStartPos,
1766                                             const bool theForward,
1767                                             const bool dotIsWord )
1768 {
1769   int beg = theStartPos, end = theStartPos;
1770   theStartPos = EMPTY;
1771   if ( beg < 1 || beg > theString.Length() )
1772     return theEmptyString;
1773
1774   if ( theForward ) { // search forward
1775     // beg
1776     while ( beg <= theString.Length() && !isWord( theString.Value( beg ), dotIsWord))
1777       ++beg;
1778     if ( beg > theString.Length() )
1779       return theEmptyString; // no word found
1780     // end
1781     end = beg + 1;
1782     while ( end <= theString.Length() && isWord( theString.Value( end ), dotIsWord))
1783       ++end;
1784     --end;
1785   }
1786   else {  // search backward
1787     // end
1788     while ( end > 0 && !isWord( theString.Value( end ), dotIsWord))
1789       --end;
1790     if ( end == 0 )
1791       return theEmptyString; // no word found
1792     beg = end - 1;
1793     while ( beg > 0 && isWord( theString.Value( beg ), dotIsWord))
1794       --beg;
1795     ++beg;
1796   }
1797   theStartPos = beg;
1798   //cout << theString << " ---- " << beg << " - " << end << endl;
1799   return theString.SubString( beg, end );
1800 }
1801
1802 //================================================================================
1803 /*!
1804  * \brief Look for position where not space char is
1805   * \param theString - The string
1806   * \param thePos - The position to search from and which returns result
1807   * \retval bool - false if there are only space after thePos in theString
1808  *
1809  *
1810  */
1811 //================================================================================
1812
1813 bool _pyCommand::SkipSpaces( const TCollection_AsciiString & theString, int & thePos )
1814 {
1815   if ( thePos < 1 || thePos > theString.Length() )
1816     return false;
1817
1818   while ( thePos <= theString.Length() && isspace( theString.Value( thePos )))
1819     ++thePos;
1820
1821   return thePos <= theString.Length();
1822 }
1823
1824 //================================================================================
1825 /*!
1826  * \brief Modify a part of the command
1827   * \param thePartIndex - The index of the part
1828   * \param thePart - The new part string
1829   * \param theOldPart - The old part
1830  */
1831 //================================================================================
1832
1833 void _pyCommand::SetPart(int thePartIndex, const TCollection_AsciiString& thePart,
1834                         TCollection_AsciiString& theOldPart)
1835 {
1836   int pos = GetBegPos( thePartIndex );
1837   if ( pos <= Length() && theOldPart != thePart)
1838   {
1839     TCollection_AsciiString seperator;
1840     if ( pos < 1 ) {
1841       pos = GetBegPos( thePartIndex + 1 );
1842       if ( pos < 1 ) return;
1843       switch ( thePartIndex ) {
1844       case RESULT_IND: seperator = " = "; break;
1845       case OBJECT_IND: seperator = "."; break;
1846       case METHOD_IND: seperator = "()"; break;
1847       default:;
1848       }
1849     }
1850     myString.Remove( pos, theOldPart.Length() );
1851     if ( !seperator.IsEmpty() )
1852       myString.Insert( pos , seperator );
1853     myString.Insert( pos, thePart );
1854     // update starting positions of the following parts
1855     int posDelta = thePart.Length() + seperator.Length() - theOldPart.Length();
1856     for ( int i = thePartIndex + 1; i <= myBegPos.Length(); ++i ) {
1857       if ( myBegPos( i ) > 0 )
1858         myBegPos( i ) += posDelta;
1859     }
1860     theOldPart = thePart;
1861   }
1862 }
1863
1864 //================================================================================
1865 /*!
1866  * \brief Set agrument
1867   * \param index - The argument index, it counts from 1
1868   * \param theArg - The argument string
1869  */
1870 //================================================================================
1871
1872 void _pyCommand::SetArg( int index, const TCollection_AsciiString& theArg)
1873 {
1874   FindAllArgs();
1875   int argInd = ARG1_IND + index - 1;
1876   int pos = GetBegPos( argInd );
1877   if ( pos < 1 ) // no index-th arg exist, append inexistent args
1878   {
1879     // find a closing parenthesis
1880     if ( int lastArgInd = GetNbArgs() ) {
1881       pos = GetBegPos( ARG1_IND + lastArgInd  - 1 ) + GetArg( lastArgInd ).Length();
1882       while ( pos > 0 && pos <= Length() && myString.Value( pos ) != ')' )
1883         ++pos;
1884     }
1885     else {
1886       pos = Length();
1887       while ( pos > 0 && myString.Value( pos ) != ')' )
1888         --pos;
1889     }
1890     if ( pos < 1 || myString.Value( pos ) != ')' ) { // no parentheses at all
1891       myString += "()";
1892       pos = Length();
1893     }
1894     while ( myArgs.Length() < index ) {
1895       if ( myArgs.Length() )
1896         myString.Insert( pos++, "," );
1897       myArgs.Append("None");
1898       myString.Insert( pos, myArgs.Last() );
1899       SetBegPos( ARG1_IND + myArgs.Length() - 1, pos );
1900       pos += myArgs.Last().Length();
1901     }
1902   }
1903   SetPart( argInd, theArg, myArgs( index ));
1904 }
1905
1906 //================================================================================
1907 /*!
1908  * \brief Empty arg list
1909  */
1910 //================================================================================
1911
1912 void _pyCommand::RemoveArgs()
1913 {
1914   if ( int pos = myString.Location( '(', 1, Length() ))
1915     myString.Trunc( pos );
1916   myString += ")";
1917   myArgs.Clear();
1918   if ( myBegPos.Length() >= ARG1_IND )
1919     myBegPos.Remove( ARG1_IND, myBegPos.Length() );
1920 }
1921
1922 //================================================================================
1923 /*!
1924  * \brief Set dependent commands after this one
1925  */
1926 //================================================================================
1927
1928 bool _pyCommand::SetDependentCmdsAfter() const
1929 {
1930   bool orderChanged = false;
1931   list< Handle(_pyCommand)>::const_reverse_iterator cmd = myDependentCmds.rbegin();
1932   for ( ; cmd != myDependentCmds.rend(); ++cmd ) {
1933     if ( (*cmd)->GetOrderNb() < GetOrderNb() ) {
1934       orderChanged = true;
1935       theGen->SetCommandAfter( *cmd, this );
1936       (*cmd)->SetDependentCmdsAfter();
1937     }
1938   }
1939   return orderChanged;
1940 }
1941 //================================================================================
1942 /*!
1943  * \brief Insert accessor method after theObjectID
1944   * \param theObjectID - id of the accessed object
1945   * \param theAcsMethod - name of the method giving access to the object
1946   * \retval bool - false if theObjectID is not found in the command string
1947  */
1948 //================================================================================
1949
1950 bool _pyCommand::AddAccessorMethod( _pyID theObjectID, const char* theAcsMethod )
1951 {
1952   if ( !theAcsMethod )
1953     return false;
1954   // start object search from the object, i.e. ignore result
1955   GetObject();
1956   int beg = GetBegPos( OBJECT_IND );
1957   if ( beg < 1 || beg > Length() )
1958     return false;
1959   bool added = false;
1960   while (( beg = myString.Location( theObjectID, beg, Length() )))
1961   {
1962     // check that theObjectID is not just a part of a longer ID
1963     int afterEnd = beg + theObjectID.Length();
1964     Standard_Character c = myString.Value( afterEnd );
1965     if ( !isalnum( c ) && c != ':' ) {
1966       // check if accessor method already present
1967       if ( c != '.' ||
1968            myString.Location( (char*) theAcsMethod, afterEnd, Length() ) != afterEnd+1) {
1969         // insertion
1970         int oldLen = Length();
1971         myString.Insert( afterEnd, (char*) theAcsMethod );
1972         myString.Insert( afterEnd, "." );
1973         // update starting positions of the parts following the modified one
1974         int posDelta = Length() - oldLen;
1975         for ( int i = 1; i <= myBegPos.Length(); ++i ) {
1976           if ( myBegPos( i ) > afterEnd )
1977             myBegPos( i ) += posDelta;
1978         }
1979         added = true;
1980       }
1981     }
1982     beg = afterEnd; // is a part - next search
1983   }
1984   return added;
1985 }
1986
1987 //================================================================================
1988 /*!
1989  * \brief Return method name giving access to an interaface object wrapped by python class
1990   * \retval const char* - method name
1991  */
1992 //================================================================================
1993
1994 const char* _pyObject::AccessorMethod() const
1995 {
1996   return 0;
1997 }
1998 //================================================================================
1999 /*!
2000  * \brief Return ID of a father
2001  */
2002 //================================================================================
2003
2004 _pyID _pyObject::FatherID(const _pyID & childID)
2005 {
2006   int colPos = childID.SearchFromEnd(':');
2007   if ( colPos > 0 )
2008     return childID.SubString( 1, colPos-1 );
2009   return "";
2010 }