Salome HOME
IPAL54382: TC8.5.0: SMESH: Regression in mesh visualization (body fitting)
[modules/smesh.git] / src / SMESH_I / SMESH_2smeshpy.cxx
1 // Copyright (C) 2007-2016  CEA/DEN, EDF R&D, OPEN CASCADE
2 //
3 // Copyright (C) 2003-2007  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, or (at your option) any later version.
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 // File      : SMESH_2smeshpy.cxx
24 // Created   : Fri Nov 18 13:20:10 2005
25 // Author    : Edward AGAPOV (eap)
26 //
27 #include "SMESH_2smeshpy.hxx"
28
29 #include "SMESH_PythonDump.hxx"
30 #include "SMESH_NoteBook.hxx"
31 #include "SMESH_Filter_i.hxx"
32
33 #include <SALOMEDS_wrap.hxx>
34 #include <utilities.h>
35
36 #include <Resource_DataMapOfAsciiStringAsciiString.hxx>
37 #include <Resource_DataMapIteratorOfDataMapOfAsciiStringAsciiString.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 #include <LDOMParser.hxx>
46
47 #ifdef WIN32
48 #include <windows.h>
49 #else
50 #include <unistd.h>
51 #endif
52
53 IMPLEMENT_STANDARD_RTTIEXT(_pyObject          ,Standard_Transient);
54 IMPLEMENT_STANDARD_RTTIEXT(_pyCommand         ,Standard_Transient);
55 IMPLEMENT_STANDARD_RTTIEXT(_pyHypothesisReader,Standard_Transient);
56 IMPLEMENT_STANDARD_RTTIEXT(_pyGen             ,_pyObject);
57 IMPLEMENT_STANDARD_RTTIEXT(_pyMesh            ,_pyObject);
58 IMPLEMENT_STANDARD_RTTIEXT(_pySubMesh         ,_pyObject);
59 IMPLEMENT_STANDARD_RTTIEXT(_pyMeshEditor      ,_pyObject);
60 IMPLEMENT_STANDARD_RTTIEXT(_pyHypothesis      ,_pyObject);
61 IMPLEMENT_STANDARD_RTTIEXT(_pySelfEraser      ,_pyObject);
62 IMPLEMENT_STANDARD_RTTIEXT(_pyGroup           ,_pyObject);
63 IMPLEMENT_STANDARD_RTTIEXT(_pyFilter          ,_pyObject);
64 IMPLEMENT_STANDARD_RTTIEXT(_pyAlgorithm       ,_pyHypothesis);
65 IMPLEMENT_STANDARD_RTTIEXT(_pyComplexParamHypo,_pyHypothesis);
66 IMPLEMENT_STANDARD_RTTIEXT(_pyNumberOfSegmentsHyp,_pyHypothesis);
67 IMPLEMENT_STANDARD_RTTIEXT(_pyLayerDistributionHypo,_pyHypothesis);
68 IMPLEMENT_STANDARD_RTTIEXT(_pySegmentLengthAroundVertexHyp,_pyHypothesis);
69
70 using namespace std;
71 using SMESH::TPythonDump;
72
73 /*!
74  * \brief Container of commands into which the initial script is split.
75  *        It also contains data corresponding to SMESH_Gen contents
76  */
77 static Handle(_pyGen) theGen;
78
79 static TCollection_AsciiString theEmptyString;
80
81 //#define DUMP_CONVERSION
82
83 #if !defined(_DEBUG_) && defined(DUMP_CONVERSION)
84 #undef DUMP_CONVERSION
85 #endif
86
87
88 namespace {
89
90   //================================================================================
91   /*!
92    * \brief Set of TCollection_AsciiString initialized by C array of C strings
93    */
94   //================================================================================
95
96   struct TStringSet: public set<TCollection_AsciiString>
97   {
98     /*!
99      * \brief Filling. The last string must be ""
100      */
101     void Insert(const char* names[]) {
102       for ( int i = 0; names[i][0] ; ++i )
103         insert( (char*) names[i] );
104     }
105     /*!
106      * \brief Check if a string is in
107      */
108     bool Contains(const TCollection_AsciiString& name ) {
109       return find( name ) != end();
110     }
111   };
112
113   //================================================================================
114   /*!
115    * \brief Map of TCollection_AsciiString initialized by C array of C strings.
116    *        Odd items of the C array are map keys, and even items are values
117    */
118   //================================================================================
119
120   struct TStringMap: public map<TCollection_AsciiString,TCollection_AsciiString>
121   {
122     /*!
123      * \brief Filling. The last string must be ""
124      */
125     void Insert(const char* names_values[]) {
126       for ( int i = 0; names_values[i][0] ; i += 2 )
127         insert( make_pair( (char*) names_values[i], names_values[i+1] ));
128     }
129     /*!
130      * \brief Check if a string is in
131      */
132     TCollection_AsciiString Value(const TCollection_AsciiString& name ) {
133       map< _AString, _AString >::iterator it = find( name );
134       return it == end() ? "" : it->second;
135     }
136   };
137
138   //================================================================================
139   /*!
140    * \brief Returns a mesh by object
141    */
142   //================================================================================
143
144   Handle(_pyMesh) ObjectToMesh( const Handle( _pyObject )& obj )
145   {
146     if ( !obj.IsNull() )
147     {
148       if ( obj->IsKind( STANDARD_TYPE( _pyMesh )))
149         return Handle(_pyMesh)::DownCast( obj );
150       else if ( obj->IsKind( STANDARD_TYPE( _pySubMesh )))
151         return Handle(_pySubMesh)::DownCast( obj )->GetMesh();
152       else if ( obj->IsKind( STANDARD_TYPE( _pyGroup )))
153         return Handle(_pyGroup)::DownCast( obj )->GetMesh();
154     }
155     return Handle(_pyMesh)();
156   }
157
158   //================================================================================
159   /*!
160    * \brief Check if objects used as args have been created by previous commands
161    */
162   //================================================================================
163
164   void CheckObjectPresence( const Handle(_pyCommand)& cmd, set<_pyID> & presentObjects)
165   {
166     // either comment or erase a command including NotPublishedObjectName()
167     if ( cmd->GetString().Location( TPythonDump::NotPublishedObjectName(), 1, cmd->Length() ))
168     {
169       bool isResultPublished = false;
170       const int nbRes = cmd->GetNbResultValues();
171       for ( int i = 0; i < nbRes; i++ )
172       {
173         _pyID objID = cmd->GetResultValue( i+1 );
174         if ( cmd->IsStudyEntry( objID ))
175           isResultPublished = (! theGen->IsNotPublished( objID ));
176         theGen->ObjectCreationRemoved( objID ); // objID.SetName( name ) is not needed
177       }
178       if ( isResultPublished )
179         cmd->Comment();
180       else
181         cmd->Clear();
182       return;
183     }
184     // check if an Object was created in the script
185     _AString comment;
186
187     _pyID obj = cmd->GetObject();
188     if ( obj.Search( "print " ) == 1 )
189       return; // print statement
190
191     if ( !obj.IsEmpty() && obj.Value( obj.Length() ) == ')' )
192       // remove an accessor method
193       obj = _pyCommand( obj ).GetObject();
194
195     const bool isMethodCall = cmd->IsMethodCall();
196     if ( !obj.IsEmpty() && isMethodCall && !presentObjects.count( obj ) )
197     {
198       comment = "not created Object";
199       theGen->ObjectCreationRemoved( obj );
200     }
201     // check if a command has not created args
202     for ( int iArg = cmd->GetNbArgs(); iArg && comment.IsEmpty(); --iArg )
203     {
204       const _pyID& arg = cmd->GetArg( iArg );
205       if ( arg.IsEmpty() || arg.Value( 1 ) == '"' || arg.Value( 1 ) == '\'' )
206         continue;
207       list< _pyID > idList = cmd->GetStudyEntries( arg );
208       list< _pyID >::iterator id = idList.begin();
209       for ( ; id != idList.end(); ++id )
210         if ( !theGen->IsGeomObject( *id ) && !presentObjects.count( *id ))
211         {
212           comment += *id + " has not been yet created";
213           break;
214         }
215       // if ( idList.empty() && cmd->IsID( arg ) && !presentObjects.count( arg ))
216       //   comment += arg + " has not been yet created";
217     }
218     // treat result objects
219     const _pyID& result = cmd->GetResultValue();
220     if ( !result.IsEmpty() && result.Value( 1 ) != '"' && result.Value( 1 ) != '\'' )
221     {
222       list< _pyID > idList = cmd->GetStudyEntries( result );
223       list< _pyID >::iterator id = idList.begin();
224       for ( ; id != idList.end(); ++id )
225       {
226         if ( comment.IsEmpty() )
227           presentObjects.insert( *id );
228         else
229           theGen->ObjectCreationRemoved( *id ); // objID.SetName( name ) is not needed
230       }
231       if ( idList.empty() && cmd->IsID( result ))
232         presentObjects.insert( result );
233     }
234     // comment the command
235     if ( !comment.IsEmpty() )
236     {
237       cmd->Comment();
238       cmd->GetString() += " ### ";
239       cmd->GetString() += comment;
240     }
241   }
242
243   //================================================================================
244   /*!
245    * \brief Fix SMESH::FunctorType arguments of SMESH::Filter::Criterion()
246    */
247   //================================================================================
248
249   void fixFunctorType( TCollection_AsciiString& Type,
250                        TCollection_AsciiString& Compare,
251                        TCollection_AsciiString& UnaryOp,
252                        TCollection_AsciiString& BinaryOp )
253   {
254     // The problem is that dumps of old studies created using filters becomes invalid
255     // when new items are inserted in the enum SMESH::FunctorType since values
256     // of this enum are dumped as integer values.
257     // This function corrects enum values of old studies given as args (Type,Compare,...)
258     // We can find out how to correct them by value of BinaryOp which can have only two
259     // values: FT_Undefined or FT_LogicalNOT.
260     // Hereafter is the history of the enum SMESH::FunctorType since v3.0.0
261     // where PythonDump appeared
262     // v 3.0.0: FT_Undefined == 25
263     // v 3.1.0: FT_Undefined == 26, new items:
264     //   - FT_Volume3D              = 7
265     // v 4.1.2: FT_Undefined == 27, new items:
266     //   - FT_BelongToGenSurface    = 17
267     // v 5.1.1: FT_Undefined == 32, new items:
268     //   - FT_FreeNodes             = 10
269     //   - FT_FreeFaces             = 11
270     //   - FT_LinearOrQuadratic     = 23
271     //   - FT_GroupColor            = 24
272     //   - FT_ElemGeomType          = 25
273     // v 5.1.5: FT_Undefined == 33, new items:
274     //   - FT_CoplanarFaces         = 26
275     // v 6.2.0: FT_Undefined == 39, new items:
276     //   - FT_MaxElementLength2D    = 8
277     //   - FT_MaxElementLength3D    = 9
278     //   - FT_BareBorderVolume      = 25
279     //   - FT_BareBorderFace        = 26
280     //   - FT_OverConstrainedVolume = 27
281     //   - FT_OverConstrainedFace   = 28
282     // v 6.5.0: FT_Undefined == 43, new items:
283     //   - FT_EqualNodes            = 14
284     //   - FT_EqualEdges            = 15
285     //   - FT_EqualFaces            = 16
286     //   - FT_EqualVolumes          = 17
287     // v 6.6.0: FT_Undefined == 44, new items:
288     //   - FT_BallDiameter          = 37
289     // v 6.7.1: FT_Undefined == 45, new items:
290     //   - FT_EntityType            = 36
291     // v 7.3.0: FT_Undefined == 46, new items:
292     //   - FT_ConnectedElements     = 39
293     // v 7.6.0: FT_Undefined == 47, new items:
294     //   - FT_BelongToMeshGroup     = 22
295     // v 8.1.0: FT_Undefined == 48, new items:
296     //   - FT_NodeConnectivityNumber= 22
297     // v 8.5.0: FT_Undefined == 49, new items:
298     //   - FT_Deflection2D          = 22
299     //
300     // It's necessary to continue recording this history and to fill
301     // undef2newItems (see below) accordingly.
302
303     typedef map< int, vector< int > > TUndef2newItems;
304     static TUndef2newItems undef2newItems;
305     if ( undef2newItems.empty() )
306     {
307       undef2newItems[ 26 ].push_back( 7 );
308       undef2newItems[ 27 ].push_back( 17 );
309       { int items[] = { 10, 11, 23, 24, 25 };
310         undef2newItems[ 32 ].assign( items, items+5 ); }
311       undef2newItems[ 33 ].push_back( 26 );
312       { int items[] = { 8, 9, 25, 26, 27, 28 };
313         undef2newItems[ 39 ].assign( items, items+6 ); }
314       { int items[] = { 14, 15, 16, 17 };
315         undef2newItems[ 43 ].assign( items, items+4 ); }
316       undef2newItems[ 44 ].push_back( 37 );
317       undef2newItems[ 45 ].push_back( 36 );
318       undef2newItems[ 46 ].push_back( 39 );
319       undef2newItems[ 47 ].push_back( 22 );
320       undef2newItems[ 48 ].push_back( 22 );
321       undef2newItems[ 49 ].push_back( 22 );
322
323       ASSERT( undef2newItems.rbegin()->first == SMESH::FT_Undefined );
324     }
325
326     int iType     = Type.IntegerValue();
327     int iCompare  = Compare.IntegerValue();
328     int iUnaryOp  = UnaryOp.IntegerValue();
329     int iBinaryOp = BinaryOp.IntegerValue();
330
331     // find out integer value of FT_Undefined at the moment of dump
332     int oldUndefined = iBinaryOp;
333     if ( iBinaryOp < iUnaryOp ) // BinaryOp was FT_LogicalNOT
334       oldUndefined += 3;
335
336     // apply history to args
337     TUndef2newItems::const_iterator undef_items =
338       undef2newItems.upper_bound( oldUndefined );
339     if ( undef_items != undef2newItems.end() )
340     {
341       int* pArg[4] = { &iType, &iCompare, &iUnaryOp, &iBinaryOp };
342       for ( ; undef_items != undef2newItems.end(); ++undef_items )
343       {
344         const vector< int > & addedItems = undef_items->second;
345         for ( size_t i = 0; i < addedItems.size(); ++i )
346           for ( int iArg = 0; iArg < 4; ++iArg )
347           {
348             int& arg = *pArg[iArg];
349             if ( arg >= addedItems[i] )
350               arg++;
351           }
352       }
353       Type     = TCollection_AsciiString( iType     );
354       Compare  = TCollection_AsciiString( iCompare  );
355       UnaryOp  = TCollection_AsciiString( iUnaryOp  );
356       BinaryOp = TCollection_AsciiString( iBinaryOp );
357     }
358   }
359
360   //================================================================================
361   /*!
362    * \brief Replaces "SMESH.PointStruct(x,y,z)" and "SMESH.DirStruct( SMESH.PointStruct(x,y,z))"
363    *        arguments of a given command by a list "[x,y,z]" if the list is accesible
364    *        type of argument.
365    */
366   //================================================================================
367
368   void StructToList( Handle( _pyCommand)& theCommand, const bool checkMethod=true )
369   {
370     static TStringSet methodsAcceptingList;
371     if ( methodsAcceptingList.empty() ) {
372       const char * methodNames[] = {
373         "GetCriterion","Reorient2D","ExtrusionSweep","ExtrusionSweepMakeGroups0D",
374         "ExtrusionSweepMakeGroups","ExtrusionSweep0D",
375         "AdvancedExtrusion","AdvancedExtrusionMakeGroups",
376         "ExtrusionSweepObject","ExtrusionSweepObject0DMakeGroups",
377         "ExtrusionSweepObjectMakeGroups","ExtrusionSweepObject0D",
378         "ExtrusionSweepObject1D","ExtrusionSweepObject1DMakeGroups",
379         "ExtrusionSweepObject2D","ExtrusionSweepObject2DMakeGroups",
380         "ExtrusionSweepObjects","RotationSweepObjects","ExtrusionAlongPathObjects",
381         "Translate","TranslateMakeGroups","TranslateMakeMesh",
382         "TranslateObject","TranslateObjectMakeGroups", "TranslateObjectMakeMesh",
383         "ExtrusionAlongPathX","ExtrusionAlongPathObjX","SplitHexahedraIntoPrisms"
384         ,"" }; // <- mark of the end
385       methodsAcceptingList.Insert( methodNames );
386     }
387     if ( !checkMethod || methodsAcceptingList.Contains( theCommand->GetMethod() ))
388     {
389       for ( int i = theCommand->GetNbArgs(); i > 0; --i )
390       {
391         const _AString & arg = theCommand->GetArg( i );
392         if ( arg.Search( "SMESH.PointStruct" ) == 1 ||
393              arg.Search( "SMESH.DirStruct"   ) == 1 )
394         {
395           Handle(_pyCommand) workCmd = new _pyCommand( arg );
396           if ( workCmd->GetNbArgs() == 1 ) // SMESH.DirStruct( SMESH.PointStruct(x,y,z))
397           {
398             workCmd = new _pyCommand( workCmd->GetArg( 1 ) );
399           }
400           if ( workCmd->GetNbArgs() == 3 ) // SMESH.PointStruct(x,y,z)
401           {
402             _AString newArg = "[ ";
403             newArg += ( workCmd->GetArg( 1 ) + ", " +
404                         workCmd->GetArg( 2 ) + ", " +
405                         workCmd->GetArg( 3 ) + " ]");
406             theCommand->SetArg( i, newArg );
407           }
408         }
409       }
410     }
411   }
412   //================================================================================
413   /*!
414    * \brief Replaces "mesh.GetIDSource([id1,id2])" argument of a given command by
415    *        a list "[id1,id2]" if the list is an accesible type of argument.
416    */
417   //================================================================================
418
419   void GetIDSourceToList( Handle( _pyCommand)& theCommand )
420   {
421     static TStringSet methodsAcceptingList;
422     if ( methodsAcceptingList.empty() ) {
423       const char * methodNames[] = {
424         "ExportPartToMED","ExportPartToDAT","ExportPartToUNV","ExportPartToSTL",
425         "ExportCGNS","ExportGMF",
426         "Create0DElementsOnAllNodes","Reorient2D","QuadTo4Tri",
427         "ScaleMakeGroups","Scale","ScaleMakeMesh",
428         "FindCoincidentNodesOnPartBut","DoubleElements",
429         "ExtrusionSweepObjects","RotationSweepObjects","ExtrusionAlongPathObjects"
430         ,"" }; // <- mark of the end
431       methodsAcceptingList.Insert( methodNames );
432     }
433     if ( methodsAcceptingList.Contains( theCommand->GetMethod() ))
434     {
435       for ( int i = theCommand->GetNbArgs(); i > 0; --i )
436       {
437         _pyCommand argCmd( theCommand->GetArg( i ));
438         if ( argCmd.GetMethod() == "GetIDSource" &&
439              argCmd.GetNbArgs() == 2 )
440         {
441           theCommand->SetArg( i, argCmd.GetArg( 1 ));
442         }
443       }
444     }
445   }
446 }
447
448 //================================================================================
449 /*!
450  * \brief Convert a python script using commands of smeshBuilder.py
451  *  \param theScriptLines - Lines of the input script
452  *  \param theEntry2AccessorMethod - returns method names to access to
453  *         objects wrapped with python class
454  *  \param theObjectNames - names of objects
455  *  \param theRemovedObjIDs - entries of objects whose created commands were removed
456  *  \param theHistoricalDump - true means to keep all commands, false means
457  *         to exclude commands relating to objects removed from study
458  *  \retval TCollection_AsciiString - Conversion result
459  */
460 //================================================================================
461
462 void
463 SMESH_2smeshpy::ConvertScript(std::list< TCollection_AsciiString >&     theScriptLines,
464                               Resource_DataMapOfAsciiStringAsciiString& theEntry2AccessorMethod,
465                               Resource_DataMapOfAsciiStringAsciiString& theObjectNames,
466                               std::set< TCollection_AsciiString >&      theRemovedObjIDs,
467                               SALOMEDS::Study_ptr&                      theStudy,
468                               const bool                                theToKeepAllCommands)
469 {
470   std::list< TCollection_AsciiString >::iterator lineIt;
471   // process notebook variables
472   {
473     SMESH_NoteBook aNoteBook;
474
475     for ( lineIt = theScriptLines.begin(); lineIt != theScriptLines.end(); ++lineIt )
476       aNoteBook.AddCommand( *lineIt );
477
478     theScriptLines.clear();
479
480     aNoteBook.ReplaceVariables();
481
482     aNoteBook.GetResultLines( theScriptLines );
483   }
484
485   // convert to smeshBuilder.py API
486
487   theGen = new _pyGen( theEntry2AccessorMethod,
488                        theObjectNames,
489                        theRemovedObjIDs,
490                        theStudy,
491                        theToKeepAllCommands );
492
493   for ( lineIt = theScriptLines.begin(); lineIt != theScriptLines.end(); ++lineIt )
494     theGen->AddCommand( *lineIt );
495
496   theScriptLines.clear();
497
498   // finish conversion
499   theGen->Flush();
500 #ifdef DUMP_CONVERSION
501   MESSAGE_BEGIN ( std::endl << " ######## RESULT ######## " << std::endl<< std::endl );
502 #endif
503
504   // clean commands of removed objects depending on myIsPublished flag
505   theGen->ClearCommands();
506
507   // reorder commands after conversion
508   list< Handle(_pyCommand) >::iterator cmd;
509   bool orderChanges;
510   do {
511     orderChanges = false;
512     for ( cmd = theGen->GetCommands().begin(); cmd != theGen->GetCommands().end(); ++cmd )
513       if ( (*cmd)->SetDependentCmdsAfter() )
514         orderChanges = true;
515   } while ( orderChanges );
516
517   // concat commands back into a script
518   TCollection_AsciiString aPrevCmd;
519   set<_pyID> createdObjects;
520   createdObjects.insert( "smeshBuilder" );
521   createdObjects.insert( "smesh" );
522   createdObjects.insert( "theStudy" );
523   for ( cmd = theGen->GetCommands().begin(); cmd != theGen->GetCommands().end(); ++cmd )
524   {
525 #ifdef DUMP_CONVERSION
526     MESSAGE_ADD ( "## COM " << (*cmd)->GetOrderNb() << ": "<< (*cmd)->GetString() << std::endl );
527 #endif
528     if ( !(*cmd)->IsEmpty() && aPrevCmd != (*cmd)->GetString()) {
529       CheckObjectPresence( *cmd, createdObjects );
530       if ( !(*cmd)->IsEmpty() ) {
531         aPrevCmd = (*cmd)->GetString();
532         theScriptLines.push_back( aPrevCmd );
533       }
534     }
535   }
536
537   theGen->Free();
538   theGen.Nullify();
539 }
540
541 //================================================================================
542 /*!
543  * \brief _pyGen constructor
544  */
545 //================================================================================
546
547 _pyGen::_pyGen(Resource_DataMapOfAsciiStringAsciiString& theEntry2AccessorMethod,
548                Resource_DataMapOfAsciiStringAsciiString& theObjectNames,
549                std::set< TCollection_AsciiString >&      theRemovedObjIDs,
550                SALOMEDS::Study_ptr&                      theStudy,
551                const bool                                theToKeepAllCommands)
552   : _pyObject( new _pyCommand( "", 0 )),
553     myNbCommands( 0 ),
554     myID2AccessorMethod( theEntry2AccessorMethod ),
555     myObjectNames( theObjectNames ),
556     myRemovedObjIDs( theRemovedObjIDs ),
557     myNbFilters( 0 ),
558     myToKeepAllCommands( theToKeepAllCommands ),
559     myStudy( SALOMEDS::Study::_duplicate( theStudy )),
560     myGeomIDNb(0), myGeomIDIndex(-1)
561 {
562   // make that GetID() to return TPythonDump::SMESHGenName()
563   GetCreationCmd()->Clear();
564   GetCreationCmd()->GetString() = TPythonDump::SMESHGenName();
565   GetCreationCmd()->GetString() += "=";
566
567   // Find 1st digit of study entry by which a GEOM object differs from a SMESH object
568   if ( !theObjectNames.IsEmpty() && !CORBA::is_nil( theStudy ))
569   {
570     // find a GEOM entry
571     _pyID geomID;
572     SALOMEDS::SComponent_wrap geomComp = theStudy->FindComponent("GEOM");
573     if ( geomComp->_is_nil() ) return;
574     CORBA::String_var entry = geomComp->GetID();
575     geomID = entry.in();
576
577     // find a SMESH entry
578     _pyID smeshID;
579     Resource_DataMapIteratorOfDataMapOfAsciiStringAsciiString e2n( theObjectNames );
580     for ( ; e2n.More() && smeshID.IsEmpty(); e2n.Next() )
581       if ( _pyCommand::IsStudyEntry( e2n.Key() ))
582         smeshID = e2n.Key();
583
584     // find 1st difference between smeshID and geomID
585     if ( !geomID.IsEmpty() && !smeshID.IsEmpty() )
586       for ( int i = 1; i <= geomID.Length() && i <= smeshID.Length(); ++i )
587         if ( geomID.Value( i ) != smeshID.Value( i ))
588         {
589           myGeomIDNb = geomID.Value( i );
590           myGeomIDIndex = i;
591         }
592   }
593 }
594
595 //================================================================================
596 /*!
597  * \brief name of SMESH_Gen in smeshBuilder.py
598  */
599 //================================================================================
600
601 const char* _pyGen::AccessorMethod() const
602 {
603   return SMESH_2smeshpy::GenName();
604 }
605
606 //================================================================================
607 /*!
608  * \brief Convert a command using a specific converter
609  *  \param theCommand - the command to convert
610  */
611 //================================================================================
612
613 Handle(_pyCommand) _pyGen::AddCommand( const TCollection_AsciiString& theCommand)
614 {
615   // store theCommand in the sequence
616   myCommands.push_back( new _pyCommand( theCommand, ++myNbCommands ));
617
618   Handle(_pyCommand) aCommand = myCommands.back();
619 #ifdef DUMP_CONVERSION
620   MESSAGE ( "## COM " << myNbCommands << ": "<< aCommand->GetString() );
621 #endif
622
623   const _pyID& objID = aCommand->GetObject();
624
625   if ( objID.IsEmpty() )
626     return aCommand;
627
628   // Prevent moving a command creating a sub-mesh to the end of the script
629   // if the sub-mesh is used in theCommand as argument
630   // if ( _pySubMesh::CanBeArgOfMethod( aCommand->GetMethod() ))
631   // {
632   //   PlaceSubmeshAfterItsCreation( aCommand );
633   // }
634
635   // Method( SMESH.PointStruct(x,y,z)... -> Method( [x,y,z]...
636   StructToList( aCommand );
637
638   const TCollection_AsciiString& method = aCommand->GetMethod();
639
640   // not to erase _pySelfEraser's etc. used as args in some commands
641   {
642 #ifdef USE_STRING_FAMILY
643     std::list<_pyID>  objIDs;
644     if ( myKeepAgrCmdsIDs.IsInArgs( aCommand, objIDs ))
645     {
646       std::list<_pyID>::iterator objID = objIDs.begin();
647       for ( ; objID != objIDs.end(); ++objID )
648       {
649         Handle(_pyObject) obj = FindObject( *objID );
650         if ( !obj.IsNull() )
651         {
652           obj->AddArgCmd( aCommand );
653           //cout << objID << " found in " << theCommand << endl;
654         }
655       }
656     }
657 #else
658     std::list< _pyID >::const_iterator id = myKeepAgrCmdsIDs.begin();
659     for ( ; id != myKeepAgrCmdsIDs.end(); ++id )
660       if ( *id != objID && theCommand.Search( *id ) > id->Length() )
661       {
662         Handle(_pyObject) obj = FindObject( *id );
663         if ( !obj.IsNull() )
664           obj->AddArgCmd( aCommand );
665       }
666 #endif
667   }
668
669   // Find an object to process theCommand
670
671   // SMESH_Gen method?
672   if ( objID == this->GetID() || objID == SMESH_2smeshpy::GenName())
673   {
674     this->Process( aCommand );
675     //addFilterUser( aCommand, theGen ); // protect filters from clearing
676     return aCommand;
677   }
678
679   // SMESH_Mesh method?
680   map< _pyID, Handle(_pyMesh) >::iterator id_mesh = myMeshes.find( objID );
681   if ( id_mesh != myMeshes.end() )
682   {
683     //id_mesh->second->AddProcessedCmd( aCommand );
684
685     // Wrap Export*() into try-except
686     if ( aCommand->MethodStartsFrom("Export"))
687     {
688       _AString    tab = "\t";
689       _AString indent = aCommand->GetIndentation();
690       _AString tryStr = indent + "try:";
691       _AString newCmd = indent + tab + ( aCommand->GetString().ToCString() + indent.Length() );
692       _AString pasCmd = indent + tab + "pass"; // to keep valid if newCmd is erased
693       _AString excStr = indent + "except:";
694       _AString msgStr = indent + "\tprint '"; msgStr += method + "() failed. Invalid file name?'";
695
696       myCommands.insert( --myCommands.end(), new _pyCommand( tryStr, myNbCommands ));
697       aCommand->Clear();
698       aCommand->GetString() = newCmd;
699       aCommand->SetOrderNb( ++myNbCommands );
700       myCommands.push_back( new _pyCommand( pasCmd, ++myNbCommands ));
701       myCommands.push_back( new _pyCommand( excStr, ++myNbCommands ));
702       myCommands.push_back( new _pyCommand( msgStr, ++myNbCommands ));
703     }
704     // check for mesh editor object
705     if ( aCommand->GetMethod() == "GetMeshEditor" ) { // MeshEditor creation
706       _pyID editorID = aCommand->GetResultValue();
707       Handle(_pyMeshEditor) editor = new _pyMeshEditor( aCommand );
708       myMeshEditors.insert( make_pair( editorID, editor ));
709       return aCommand;
710     }
711     // check for SubMesh objects
712     else if ( aCommand->GetMethod() == "GetSubMesh" ) { // SubMesh creation
713       _pyID subMeshID = aCommand->GetResultValue();
714       Handle(_pySubMesh) subMesh = new _pySubMesh( aCommand );
715       AddObject( subMesh );
716     }
717
718     // Method( mesh.GetIDSource([id1,id2]) -> Method( [id1,id2]
719     GetIDSourceToList( aCommand );
720
721     //addFilterUser( aCommand, theGen ); // protect filters from clearing
722
723     id_mesh->second->Process( aCommand );
724     id_mesh->second->AddProcessedCmd( aCommand );
725     return aCommand;
726   }
727
728   // SMESH_MeshEditor method?
729   map< _pyID, Handle(_pyMeshEditor) >::iterator id_editor = myMeshEditors.find( objID );
730   if ( id_editor != myMeshEditors.end() )
731   {
732     // Method( mesh.GetIDSource([id1,id2]) -> Method( [id1,id2]
733     GetIDSourceToList( aCommand );
734
735     //addFilterUser( aCommand, theGen ); // protect filters from clearing
736
737     // some commands of SMESH_MeshEditor create meshes and groups
738     _pyID meshID, groups;
739     if ( method.Search("MakeMesh") != -1 )
740       meshID = aCommand->GetResultValue();
741     else if ( method == "MakeBoundaryMesh")
742       meshID = aCommand->GetResultValue(1);
743     else if ( method == "MakeBoundaryElements")
744       meshID = aCommand->GetResultValue(2);
745
746     if ( method.Search("MakeGroups") != -1      ||
747          method == "ExtrusionAlongPathX"        ||
748          method == "ExtrusionAlongPathObjX"     ||
749          method == "DoubleNodeGroupNew"         ||
750          method == "DoubleNodeGroupsNew"        ||
751          method == "DoubleNodeElemGroupNew"     ||
752          method == "DoubleNodeElemGroupsNew"    ||
753          method == "DoubleNodeElemGroup2New"    ||
754          method == "DoubleNodeElemGroups2New"   ||
755          method == "AffectedElemGroupsInRegion"
756          )
757       groups = aCommand->GetResultValue();
758     else if ( method == "MakeBoundaryMesh" )
759       groups = aCommand->GetResultValue(2);
760     else if ( method == "MakeBoundaryElements")
761       groups = aCommand->GetResultValue(3);
762     else if ( method == "Create0DElementsOnAllNodes" &&
763               aCommand->GetArg(2).Length() > 2 ) // group name != ''
764       groups = aCommand->GetResultValue();
765
766     id_editor->second->Process( aCommand );
767     id_editor->second->AddProcessedCmd( aCommand );
768
769     // create meshes
770     if ( !meshID.IsEmpty() &&
771          !myMeshes.count( meshID ) &&
772          aCommand->IsStudyEntry( meshID ))
773     {
774       _AString processedCommand = aCommand->GetString();
775       Handle(_pyMesh) mesh = new _pyMesh( aCommand, meshID );
776       CheckObjectIsReCreated( mesh );
777       myMeshes.insert( make_pair( meshID, mesh ));
778       aCommand->Clear();
779       aCommand->GetString() = processedCommand; // discard changes made by _pyMesh
780     }
781     // create groups
782     if ( !groups.IsEmpty() )
783     {
784       if ( !aCommand->IsStudyEntry( meshID ))
785         meshID = id_editor->second->GetMesh();
786       Handle(_pyMesh) mesh = myMeshes[ meshID ];
787
788       list< _pyID > idList = aCommand->GetStudyEntries( groups );
789       list< _pyID >::iterator grID = idList.begin();
790       for ( ; grID != idList.end(); ++grID )
791         if ( !myObjects.count( *grID ))
792         {
793           Handle(_pyGroup) group = new _pyGroup( aCommand, *grID );
794           AddObject( group );
795           if ( !mesh.IsNull() ) mesh->AddGroup( group );
796         }
797     }
798     return aCommand;
799   } // SMESH_MeshEditor methods
800
801   // SMESH_Hypothesis method?
802   Handle(_pyHypothesis) hyp = FindHyp( objID );
803   if ( !hyp.IsNull() && !hyp->IsAlgo() )
804   {
805     hyp->Process( aCommand );
806     hyp->AddProcessedCmd( aCommand );
807     return aCommand;
808   }
809
810   // aFilterManager.CreateFilter() ?
811   if ( aCommand->GetMethod() == "CreateFilter" )
812   {
813     // Set a more human readable name to a filter
814     // aFilter0x7fbf6c71cfb0 -> aFilter_nb
815     _pyID newID, filterID = aCommand->GetResultValue();
816     int pos = filterID.Search( "0x" );
817     if ( pos > 1 )
818       newID = (filterID.SubString(1,pos-1) + "_") + _pyID( ++myNbFilters );
819
820     Handle(_pyObject) filter( new _pyFilter( aCommand, newID ));
821     AddObject( filter );
822   }
823   // aFreeNodes0x5011f80 = aFilterManager.CreateFreeNodes() ## issue 0020976
824   else if ( theCommand.Search( "aFilterManager.Create" ) > 0 )
825   {
826     // create _pySelfEraser for functors
827     Handle(_pySelfEraser) functor = new _pySelfEraser( aCommand );
828     functor->IgnoreOwnCalls(); // to erase if not used as an argument
829     AddObject( functor );
830   }
831
832   // other object method?
833   map< _pyID, Handle(_pyObject) >::iterator id_obj = myObjects.find( objID );
834   if ( id_obj != myObjects.end() ) {
835     id_obj->second->Process( aCommand );
836     id_obj->second->AddProcessedCmd( aCommand );
837     return aCommand;
838   }
839
840   // Add access to a wrapped mesh
841   AddMeshAccessorMethod( aCommand );
842
843   // Add access to a wrapped algorithm
844   //  AddAlgoAccessorMethod( aCommand ); // ??? what if algo won't be wrapped at all ???
845
846   // PAL12227. PythonDump was not updated at proper time; result is
847   //     aCriteria.append(SMESH.Filter.Criterion(17,26,0,'L1',26,25,1e-07,SMESH.EDGE,-1))
848   // TypeError: __init__() takes exactly 11 arguments (10 given)
849   const char wrongCommand[] = "SMESH.Filter.Criterion(";
850   if ( int beg = theCommand.Location( wrongCommand, 1, theCommand.Length() ))
851   {
852     _pyCommand tmpCmd( theCommand.SubString( beg, theCommand.Length() ), -1);
853     // there must be 10 arguments, 5-th arg ThresholdID is missing,
854     const int wrongNbArgs = 9, missingArg = 5;
855     if ( tmpCmd.GetNbArgs() == wrongNbArgs )
856     {
857       for ( int i = wrongNbArgs; i > missingArg; --i )
858         tmpCmd.SetArg( i + 1, tmpCmd.GetArg( i ));
859       tmpCmd.SetArg(  missingArg, "''");
860       aCommand->GetString().Trunc( beg - 1 );
861       aCommand->GetString() += tmpCmd.GetString();
862     }
863     // IMP issue 0021014
864     // set GetCriterion(elementType,CritType,Compare,Threshold,UnaryOp,BinaryOp,Tolerance)
865     //                  1           2        3       4         5       6        7
866     // instead of "SMESH.Filter.Criterion(
867     // Type,Compare,Threshold,ThresholdStr,ThresholdID,UnaryOp,BinaryOp,Tolerance,TypeOfElement,Precision)
868     // 1    2       3         4            5           6       7        8         9             10
869     // in order to avoid the problem of type mismatch of long and FunctorType
870     const TCollection_AsciiString
871       SMESH("SMESH."), dfltFunctor("SMESH.FT_Undefined"), dfltTol("1e-07"), dfltPreci("-1");
872     TCollection_AsciiString
873       Type          = aCommand->GetArg(1),  // long
874       Compare       = aCommand->GetArg(2),  // long
875       Threshold     = aCommand->GetArg(3),  // double
876       ThresholdStr  = aCommand->GetArg(4),  // string
877       ThresholdID   = aCommand->GetArg(5),  // string
878       UnaryOp       = aCommand->GetArg(6),  // long
879       BinaryOp      = aCommand->GetArg(7),  // long
880       Tolerance     = aCommand->GetArg(8),  // double
881       TypeOfElement = aCommand->GetArg(9),  // ElementType
882       Precision     = aCommand->GetArg(10); // long
883     fixFunctorType( Type, Compare, UnaryOp, BinaryOp );
884     Type     = SMESH + SMESH::FunctorTypeToString( SMESH::FunctorType( Type.IntegerValue() ));
885     Compare  = SMESH + SMESH::FunctorTypeToString( SMESH::FunctorType( Compare.IntegerValue() ));
886     UnaryOp  = SMESH + SMESH::FunctorTypeToString( SMESH::FunctorType( UnaryOp.IntegerValue() ));
887     BinaryOp = SMESH + SMESH::FunctorTypeToString( SMESH::FunctorType( BinaryOp.IntegerValue() ));
888
889     if ( Compare == "SMESH.FT_EqualTo" )
890       Compare = "'='";
891
892     aCommand->RemoveArgs();
893     aCommand->SetObject( SMESH_2smeshpy::GenName() );
894     aCommand->SetMethod( "GetCriterion" );
895
896     aCommand->SetArg( 1, TypeOfElement );
897     aCommand->SetArg( 2, Type );
898     aCommand->SetArg( 3, Compare );
899
900     if ( Threshold.IsIntegerValue() )
901     {
902       int iGeom = Threshold.IntegerValue();
903       if ( Type == "SMESH.FT_ElemGeomType" )
904       {
905         // set SMESH.GeometryType instead of a numerical Threshold
906         const int nbTypes = SMESH::Geom_LAST;
907         const char* types[] = {
908           "Geom_POINT", "Geom_EDGE", "Geom_TRIANGLE", "Geom_QUADRANGLE", "Geom_POLYGON",
909           "Geom_TETRA", "Geom_PYRAMID", "Geom_HEXA", "Geom_PENTA", "Geom_HEXAGONAL_PRISM",
910           "Geom_POLYHEDRA", "Geom_BALL" };
911         if ( -1 < iGeom && iGeom < nbTypes )
912           Threshold = SMESH + types[ iGeom ];
913 #ifdef _DEBUG_
914         // is types complete? (compilation failure mains that enum GeometryType changed)
915         int _asrt[( sizeof(types) / sizeof(const char*) == nbTypes ) ? 2 : -1 ]; _asrt[0]=_asrt[1];
916 #endif
917       }
918       if (Type == "SMESH.FT_EntityType")
919       {
920         // set SMESH.EntityType instead of a numerical Threshold
921         const int nbTypes = SMESH::Entity_Last;
922         const char* types[] = {
923           "Entity_Node", "Entity_0D", "Entity_Edge", "Entity_Quad_Edge",
924           "Entity_Triangle", "Entity_Quad_Triangle", "Entity_BiQuad_Triangle",
925           "Entity_Quadrangle", "Entity_Quad_Quadrangle", "Entity_BiQuad_Quadrangle",
926           "Entity_Polygon", "Entity_Quad_Polygon", "Entity_Tetra", "Entity_Quad_Tetra",
927           "Entity_Pyramid", "Entity_Quad_Pyramid",
928           "Entity_Hexa", "Entity_Quad_Hexa", "Entity_TriQuad_Hexa",
929           "Entity_Penta", "Entity_Quad_Penta", "Entity_BiQuad_Penta", "Entity_Hexagonal_Prism",
930           "Entity_Polyhedra", "Entity_Quad_Polyhedra", "Entity_Ball" };
931         if ( -1 < iGeom && iGeom < nbTypes )
932           Threshold = SMESH + types[ iGeom ];
933 #ifdef _DEBUG_
934         // is 'types' complete? (compilation failure mains that enum EntityType changed)
935         int _asrt[( sizeof(types) / sizeof(const char*) == nbTypes ) ? 2 : -1 ]; _asrt[0]=_asrt[1];
936 #endif
937       }
938     }
939     if ( ThresholdID.Length() != 2 ) // neither '' nor ""
940       aCommand->SetArg( 4, ThresholdID.SubString( 2, ThresholdID.Length()-1 )); // shape entry
941     else if ( ThresholdStr.Length() != 2 )
942       aCommand->SetArg( 4, ThresholdStr );
943     else if ( ThresholdID.Length() != 2 )
944       aCommand->SetArg( 4, ThresholdID );
945     else
946       aCommand->SetArg( 4, Threshold );
947     // find the last not default arg
948     int lastDefault = 8;
949     if ( Tolerance == dfltTol ) {
950       lastDefault = 7;
951       if ( BinaryOp == dfltFunctor ) {
952         lastDefault = 6;
953         if ( UnaryOp == dfltFunctor )
954           lastDefault = 5;
955       }
956     }
957     if ( 5 < lastDefault ) aCommand->SetArg( 5, UnaryOp );
958     if ( 6 < lastDefault ) aCommand->SetArg( 6, BinaryOp );
959     if ( 7 < lastDefault ) aCommand->SetArg( 7, Tolerance );
960     if ( Precision != dfltPreci )
961     {
962       TCollection_AsciiString crit = aCommand->GetResultValue();
963       aCommand->GetString() += "; ";
964       aCommand->GetString() += crit + ".Precision = " + Precision;
965     }
966   }
967   return aCommand;
968 }
969
970 //================================================================================
971 /*!
972  * \brief Convert the command or remember it for later conversion
973  *  \param theCommand - The python command calling a method of SMESH_Gen
974  */
975 //================================================================================
976
977 void _pyGen::Process( const Handle(_pyCommand)& theCommand )
978 {
979   // there are methods to convert:
980   // CreateMesh( shape )
981   // Concatenate( [mesh1, ...], ... )
982   // CreateHypothesis( theHypType, theLibName )
983   // Compute( mesh, geom )
984   // Evaluate( mesh, geom )
985   // mesh creation
986   TCollection_AsciiString method = theCommand->GetMethod();
987
988   if ( method == "CreateMesh" || method == "CreateEmptyMesh")
989   {
990     Handle(_pyMesh) mesh = new _pyMesh( theCommand );
991     AddObject( mesh );
992     return;
993   }
994   if ( method == "CreateMeshesFromUNV" ||
995        method == "CreateMeshesFromSTL" ||
996        method == "CopyMesh" ) // command result is a mesh
997   {
998     Handle(_pyMesh) mesh = new _pyMesh( theCommand, theCommand->GetResultValue() );
999     AddObject( mesh );
1000     return;
1001   }
1002   if( method == "CreateMeshesFromMED" ||
1003       method == "CreateMeshesFromSAUV"||
1004       method == "CreateMeshesFromCGNS" ||
1005       method == "CreateMeshesFromGMF" ) // command result is ( [mesh1,mesh2], status )
1006   {
1007     std::list< _pyID > meshIDs = theCommand->GetStudyEntries( theCommand->GetResultValue() );
1008     std::list< _pyID >::iterator meshID = meshIDs.begin();
1009     for ( ; meshID != meshIDs.end(); ++meshID )
1010     {
1011       Handle(_pyMesh) mesh = new _pyMesh( theCommand, *meshID );
1012       AddObject( mesh );
1013     }
1014     if ( method == "CreateMeshesFromGMF" )
1015     {
1016       // CreateMeshesFromGMF( theFileName, theMakeRequiredGroups ) ->
1017       // CreateMeshesFromGMF( theFileName )
1018       _AString file = theCommand->GetArg(1);
1019       theCommand->RemoveArgs();
1020       theCommand->SetArg( 1, file );
1021     }
1022   }
1023
1024   // CreateHypothesis()
1025   if ( method == "CreateHypothesis" )
1026   {
1027     // issue 199929, remove standard library name (default parameter)
1028     const TCollection_AsciiString & aLibName = theCommand->GetArg( 2 );
1029     if ( aLibName.Search( "StdMeshersEngine" ) != -1 ) {
1030       // keep the first argument
1031       TCollection_AsciiString arg = theCommand->GetArg( 1 );
1032       theCommand->RemoveArgs();
1033       theCommand->SetArg( 1, arg );
1034     }
1035
1036     Handle(_pyHypothesis) hyp = _pyHypothesis::NewHypothesis( theCommand );
1037     CheckObjectIsReCreated( hyp );
1038     myHypos.insert( make_pair( hyp->GetID(), hyp ));
1039
1040     return;
1041   }
1042
1043   // smeshgen.Compute( mesh, geom ) --> mesh.Compute()
1044   if ( method == "Compute" )
1045   {
1046     const _pyID& meshID = theCommand->GetArg( 1 );
1047     map< _pyID, Handle(_pyMesh) >::iterator id_mesh = myMeshes.find( meshID );
1048     if ( id_mesh != myMeshes.end() ) {
1049       theCommand->SetObject( meshID );
1050       theCommand->RemoveArgs();
1051       id_mesh->second->Process( theCommand );
1052       id_mesh->second->AddProcessedCmd( theCommand );
1053       return;
1054     }
1055   }
1056
1057   // smeshgen.Evaluate( mesh, geom ) --> mesh.Evaluate(geom)
1058   if ( method == "Evaluate" )
1059   {
1060     const _pyID& meshID = theCommand->GetArg( 1 );
1061     map< _pyID, Handle(_pyMesh) >::iterator id_mesh = myMeshes.find( meshID );
1062     if ( id_mesh != myMeshes.end() ) {
1063       theCommand->SetObject( meshID );
1064       _pyID geom = theCommand->GetArg( 2 );
1065       theCommand->RemoveArgs();
1066       theCommand->SetArg( 1, geom );
1067       id_mesh->second->AddProcessedCmd( theCommand );
1068       return;
1069     }
1070   }
1071
1072   // objects erasing creation command if no more its commands invoked:
1073   // SMESH_Pattern, FilterManager
1074   if ( method == "GetPattern" ||
1075        method == "CreateFilterManager" ||
1076        method == "CreateMeasurements" )
1077   {
1078     Handle(_pyObject) obj = new _pySelfEraser( theCommand );
1079     if ( !AddObject( obj ) )
1080       theCommand->Clear(); // already created
1081   }
1082   // Concatenate( [mesh1, ...], ... )
1083   else if ( method == "Concatenate" || method == "ConcatenateWithGroups")
1084   {
1085     if ( method == "ConcatenateWithGroups" ) {
1086       theCommand->SetMethod( "Concatenate" );
1087       theCommand->SetArg( theCommand->GetNbArgs() + 1, "True" );
1088     }
1089     Handle(_pyMesh) mesh = new _pyMesh( theCommand, theCommand->GetResultValue() );
1090     AddObject( mesh );
1091     AddMeshAccessorMethod( theCommand );
1092   }
1093   else if ( method == "SetName" ) // SetName(obj,name)
1094   {
1095     // store theCommand as one of object commands to erase it along with the object
1096     const _pyID& objID = theCommand->GetArg( 1 );
1097     Handle(_pyObject) obj = FindObject( objID );
1098     if ( !obj.IsNull() )
1099       obj->AddProcessedCmd( theCommand );
1100   }
1101
1102   // Replace name of SMESH_Gen
1103
1104   // names of SMESH_Gen methods fully equal to methods defined in smeshBuilder.py
1105   static TStringSet smeshpyMethods;
1106   if ( smeshpyMethods.empty() ) {
1107     const char * names[] =
1108       { "SetEmbeddedMode","IsEmbeddedMode","SetCurrentStudy","GetCurrentStudy",
1109         "GetPattern","GetSubShapesId",
1110         "" }; // <- mark of array end
1111     smeshpyMethods.Insert( names );
1112   }
1113   if ( smeshpyMethods.Contains( theCommand->GetMethod() ))
1114     // smeshgen.Method() --> smesh.Method()
1115     theCommand->SetObject( SMESH_2smeshpy::SmeshpyName() );
1116   else
1117     // smeshgen.Method() --> smesh.Method()
1118     theCommand->SetObject( SMESH_2smeshpy::GenName() );
1119 }
1120
1121 //================================================================================
1122 /*!
1123  * \brief Convert the remembered commands
1124  */
1125 //================================================================================
1126
1127 void _pyGen::Flush()
1128 {
1129   // create an empty command
1130   myLastCommand = new _pyCommand();
1131
1132   map< _pyID, Handle(_pyMesh) >::iterator id_mesh;
1133   map< _pyID, Handle(_pyObject) >::iterator id_obj;
1134   map< _pyID, Handle(_pyHypothesis) >::iterator id_hyp;
1135
1136   if ( IsToKeepAllCommands() ) // historical dump
1137   {
1138     // set myIsPublished = true to all objects
1139     for ( id_mesh = myMeshes.begin(); id_mesh != myMeshes.end(); ++id_mesh )
1140       id_mesh->second->SetRemovedFromStudy( false );
1141     for ( id_hyp = myHypos.begin(); id_hyp != myHypos.end(); ++id_hyp )
1142       id_hyp->second->SetRemovedFromStudy( false );
1143     for ( id_obj = myObjects.begin(); id_obj != myObjects.end(); ++id_obj )
1144       id_obj->second->SetRemovedFromStudy( false );
1145   }
1146   else
1147   {
1148     // let hypotheses find referred objects in order to prevent clearing
1149     // not published referred hyps (it's needed for hyps like "LayerDistribution")
1150     list< Handle(_pyMesh) > fatherMeshes;
1151     for ( id_hyp = myHypos.begin(); id_hyp != myHypos.end(); ++id_hyp )
1152       if ( !id_hyp->second.IsNull() )
1153         id_hyp->second->GetReferredMeshesAndGeom( fatherMeshes );
1154   }
1155   // set myIsPublished = false to all objects depending on
1156   // meshes built on a removed geometry
1157   for ( id_mesh = myMeshes.begin(); id_mesh != myMeshes.end(); ++id_mesh )
1158     if ( id_mesh->second->IsNotGeomPublished() )
1159       id_mesh->second->SetRemovedFromStudy( true );
1160
1161   // Flush meshes
1162   for ( id_mesh = myMeshes.begin(); id_mesh != myMeshes.end(); ++id_mesh )
1163     if ( ! id_mesh->second.IsNull() )
1164       id_mesh->second->Flush();
1165
1166   // Flush hyps
1167   for ( id_hyp = myHypos.begin(); id_hyp != myHypos.end(); ++id_hyp )
1168     if ( !id_hyp->second.IsNull() ) {
1169       id_hyp->second->Flush();
1170       // smeshgen.CreateHypothesis() --> smesh.CreateHypothesis()
1171       if ( !id_hyp->second->IsWrapped() )
1172         id_hyp->second->GetCreationCmd()->SetObject( SMESH_2smeshpy::GenName() );
1173     }
1174
1175   // Flush other objects. 2 times, for objects depending on Flush() of later created objects
1176   std::list< Handle(_pyObject) >::reverse_iterator robj = myOrderedObjects.rbegin();
1177   for ( ; robj != myOrderedObjects.rend(); ++robj )
1178     if ( ! robj->IsNull() )
1179       (*robj)->Flush();
1180   std::list< Handle(_pyObject) >::iterator obj = myOrderedObjects.begin();
1181   for ( ; obj != myOrderedObjects.end(); ++obj )
1182     if ( ! obj->IsNull() )
1183       (*obj)->Flush();
1184
1185   myLastCommand->SetOrderNb( ++myNbCommands );
1186   myCommands.push_back( myLastCommand );
1187 }
1188
1189 //================================================================================
1190 /*!
1191  * \brief Prevent moving a command creating a sub-mesh to the end of the script
1192  *        if the sub-mesh is used in theCmdUsingSubmesh as argument
1193  */
1194 //================================================================================
1195
1196 void _pyGen::PlaceSubmeshAfterItsCreation( Handle(_pyCommand) theCmdUsingSubmesh ) const
1197 {
1198   // map< _pyID, Handle(_pyObject) >::const_iterator id_obj = myObjects.begin();
1199   // for ( ; id_obj != myObjects.end(); ++id_obj )
1200   // {
1201   //   if ( !id_obj->second->IsKind( STANDARD_TYPE( _pySubMesh ))) continue;
1202   //   for ( int iArg = theCmdUsingSubmesh->GetNbArgs(); iArg; --iArg )
1203   //   {
1204   //     const _pyID& arg = theCmdUsingSubmesh->GetArg( iArg );
1205   //     if ( arg.IsEmpty() || arg.Value( 1 ) == '"' || arg.Value( 1 ) == '\'' )
1206   //       continue;
1207   //     list< _pyID > idList = theCmdUsingSubmesh->GetStudyEntries( arg );
1208   //     list< _pyID >::iterator id = idList.begin();
1209   //     for ( ; id != idList.end(); ++id )
1210   //       if ( id_obj->first == *id )
1211   //         // _pySubMesh::Process() does what we need
1212   //         Handle(_pySubMesh)::DownCast( id_obj->second )->Process( theCmdUsingSubmesh );
1213   //   }
1214   // }
1215 }
1216
1217 //================================================================================
1218 /*!
1219  * \brief Clean commands of removed objects depending on myIsPublished flag
1220  */
1221 //================================================================================
1222
1223 void _pyGen::ClearCommands()
1224 {
1225   map< _pyID, Handle(_pyMesh) >::iterator id_mesh = myMeshes.begin();
1226   for ( ; id_mesh != myMeshes.end(); ++id_mesh )
1227     id_mesh->second->ClearCommands();
1228
1229   map< _pyID, Handle(_pyHypothesis) >::iterator id_hyp = myHypos.begin();
1230   for ( ; id_hyp != myHypos.end(); ++id_hyp )
1231     if ( !id_hyp->second.IsNull() )
1232       id_hyp->second->ClearCommands();
1233
1234   // Other objects. 2 times, for objects depending on ClearCommands() of later created objects
1235   std::list< Handle(_pyObject) >::reverse_iterator robj = myOrderedObjects.rbegin();
1236   for ( ; robj != myOrderedObjects.rend(); ++robj )
1237     if ( ! robj->IsNull() )
1238       (*robj)->ClearCommands();
1239   std::list< Handle(_pyObject) >::iterator obj = myOrderedObjects.begin();
1240   for ( ; obj != myOrderedObjects.end(); ++obj )
1241     if ( ! obj->IsNull() )
1242       (*obj)->ClearCommands();
1243 }
1244
1245 //================================================================================
1246 /*!
1247  * \brief Release mutual handles of objects
1248  */
1249 //================================================================================
1250
1251 void _pyGen::Free()
1252 {
1253   map< _pyID, Handle(_pyMesh) >::iterator id_mesh = myMeshes.begin();
1254   for ( ; id_mesh != myMeshes.end(); ++id_mesh )
1255     id_mesh->second->Free();
1256   myMeshes.clear();
1257
1258   map< _pyID, Handle(_pyMeshEditor) >::iterator id_ed = myMeshEditors.begin();
1259   for ( ; id_ed != myMeshEditors.end(); ++id_ed )
1260     id_ed->second->Free();
1261   myMeshEditors.clear();
1262
1263   map< _pyID, Handle(_pyObject) >::iterator id_obj = myObjects.begin();
1264   for ( ; id_obj != myObjects.end(); ++id_obj )
1265     id_obj->second->Free();
1266   myObjects.clear();
1267
1268   map< _pyID, Handle(_pyHypothesis) >::iterator id_hyp = myHypos.begin();
1269   for ( ; id_hyp != myHypos.end(); ++id_hyp )
1270     if ( !id_hyp->second.IsNull() )
1271       id_hyp->second->Free();
1272   myHypos.clear();
1273
1274   myFile2ExportedMesh.clear();
1275
1276   //myKeepAgrCmdsIDs.Print();
1277 }
1278
1279 //================================================================================
1280 /*!
1281  * \brief Add access method to mesh that is an argument
1282  *  \param theCmd - command to add access method
1283  * \retval bool - true if added
1284  */
1285 //================================================================================
1286
1287 bool _pyGen::AddMeshAccessorMethod( Handle(_pyCommand) theCmd ) const
1288 {
1289   bool added = false;
1290   map< _pyID, Handle(_pyMesh) >::const_iterator id_mesh = myMeshes.begin();
1291   for ( ; id_mesh != myMeshes.end(); ++id_mesh ) {
1292     if ( theCmd->AddAccessorMethod( id_mesh->first, id_mesh->second->AccessorMethod() ))
1293       added = true;
1294   }
1295   return added;
1296 }
1297
1298 //================================================================================
1299 /*!
1300  * \brief Add access method to algo that is an object or an argument
1301  *  \param theCmd - command to add access method
1302  * \retval bool - true if added
1303  */
1304 //================================================================================
1305
1306 bool _pyGen::AddAlgoAccessorMethod( Handle(_pyCommand) theCmd ) const
1307 {
1308   bool added = false;
1309   map< _pyID, Handle(_pyHypothesis) >::const_iterator id_hyp = myHypos.begin();
1310   for ( ; id_hyp != myHypos.end(); ++id_hyp )
1311     if ( !id_hyp->second.IsNull() &&
1312          id_hyp->second->IsAlgo() && /*(*hyp)->IsWrapped() &&*/
1313          theCmd->AddAccessorMethod( id_hyp->second->GetID(),
1314                                     id_hyp->second->AccessorMethod() ))
1315       added = true;
1316
1317   return added;
1318 }
1319
1320 //================================================================================
1321 /*!
1322  * \brief Find hypothesis by ID (entry)
1323  *  \param theHypID - The hypothesis ID
1324  * \retval Handle(_pyHypothesis) - The found hypothesis
1325  */
1326 //================================================================================
1327
1328 Handle(_pyHypothesis) _pyGen::FindHyp( const _pyID& theHypID )
1329 {
1330   map< _pyID, Handle(_pyHypothesis) >::iterator id_hyp = myHypos.find( theHypID );
1331   if ( id_hyp != myHypos.end() &&
1332        !id_hyp->second.IsNull() &&
1333        theHypID == id_hyp->second->GetID() )
1334     return id_hyp->second;
1335   return Handle(_pyHypothesis)();
1336 }
1337
1338 //================================================================================
1339 /*!
1340  * \brief Find algorithm able to create a hypothesis
1341  *  \param theGeom - The shape ID the algorithm was created on
1342  *  \param theMesh - The mesh ID that created the algorithm
1343  *  \param theHypothesis - The hypothesis the algorithm should be able to create
1344  * \retval Handle(_pyHypothesis) - The found algo
1345  */
1346 //================================================================================
1347
1348 Handle(_pyHypothesis) _pyGen::FindAlgo( const _pyID& theGeom, const _pyID& theMesh,
1349                                         const Handle(_pyHypothesis)& theHypothesis )
1350 {
1351   map< _pyID, Handle(_pyHypothesis) >::iterator id_hyp = myHypos.begin();
1352   for ( ; id_hyp != myHypos.end(); ++id_hyp )
1353     if ( !id_hyp->second.IsNull() &&
1354          id_hyp->second->IsAlgo() &&
1355          theHypothesis->CanBeCreatedBy( id_hyp->second->GetAlgoType() ) &&
1356          id_hyp->second->GetGeom() == theGeom &&
1357          id_hyp->second->GetMesh() == theMesh )
1358       return id_hyp->second;
1359   return Handle(_pyHypothesis)();
1360 }
1361
1362 //================================================================================
1363 /*!
1364  * \brief Find subMesh by ID (entry)
1365  *  \param theSubMeshID - The subMesh ID
1366  * \retval Handle(_pySubMesh) - The found subMesh
1367  */
1368 //================================================================================
1369
1370 Handle(_pySubMesh) _pyGen::FindSubMesh( const _pyID& theSubMeshID )
1371 {
1372   map< _pyID, Handle(_pyObject) >::iterator id_subMesh = myObjects.find(theSubMeshID);
1373   if ( id_subMesh != myObjects.end() )
1374     return Handle(_pySubMesh)::DownCast( id_subMesh->second );
1375   return Handle(_pySubMesh)();
1376 }
1377
1378
1379 //================================================================================
1380 /*!
1381  * \brief Change order of commands in the script
1382  *  \param theCmd1 - One command
1383  *  \param theCmd2 - Another command
1384  */
1385 //================================================================================
1386
1387 void _pyGen::ExchangeCommands( Handle(_pyCommand) theCmd1, Handle(_pyCommand) theCmd2 )
1388 {
1389   list< Handle(_pyCommand) >::iterator pos1, pos2;
1390   pos1 = find( myCommands.begin(), myCommands.end(), theCmd1 );
1391   pos2 = find( myCommands.begin(), myCommands.end(), theCmd2 );
1392   myCommands.insert( pos1, theCmd2 );
1393   myCommands.insert( pos2, theCmd1 );
1394   myCommands.erase( pos1 );
1395   myCommands.erase( pos2 );
1396
1397   int nb1 = theCmd1->GetOrderNb();
1398   theCmd1->SetOrderNb( theCmd2->GetOrderNb() );
1399   theCmd2->SetOrderNb( nb1 );
1400   //   cout << "BECOME " << theCmd1->GetOrderNb() << "\t" << theCmd1->GetString() << endl
1401   //        << "BECOME " << theCmd2->GetOrderNb() << "\t" << theCmd2->GetString() << endl << endl;
1402 }
1403
1404 //================================================================================
1405 /*!
1406  * \brief Set one command after the other
1407  *  \param theCmd - Command to move
1408  *  \param theAfterCmd - Command ater which to insert the first one
1409  */
1410 //================================================================================
1411
1412 void _pyGen::SetCommandAfter( Handle(_pyCommand) theCmd, Handle(_pyCommand) theAfterCmd )
1413 {
1414   setNeighbourCommand( theCmd, theAfterCmd, true );
1415 }
1416
1417 //================================================================================
1418 /*!
1419  * \brief Set one command before the other
1420  *  \param theCmd - Command to move
1421  *  \param theBeforeCmd - Command before which to insert the first one
1422  */
1423 //================================================================================
1424
1425 void _pyGen::SetCommandBefore( Handle(_pyCommand) theCmd, Handle(_pyCommand) theBeforeCmd )
1426 {
1427   setNeighbourCommand( theCmd, theBeforeCmd, false );
1428 }
1429
1430 //================================================================================
1431 /*!
1432  * \brief Set one command before or after the other
1433  *  \param theCmd - Command to move
1434  *  \param theOtherCmd - Command ater or before which to insert the first one
1435  */
1436 //================================================================================
1437
1438 void _pyGen::setNeighbourCommand( Handle(_pyCommand)& theCmd,
1439                                   Handle(_pyCommand)& theOtherCmd,
1440                                   const bool theIsAfter )
1441 {
1442   list< Handle(_pyCommand) >::iterator pos;
1443   pos = find( myCommands.begin(), myCommands.end(), theCmd );
1444   myCommands.erase( pos );
1445   pos = find( myCommands.begin(), myCommands.end(), theOtherCmd );
1446   myCommands.insert( (theIsAfter ? ++pos : pos), theCmd );
1447
1448   int i = 1;
1449   for ( pos = myCommands.begin(); pos != myCommands.end(); ++pos)
1450     (*pos)->SetOrderNb( i++ );
1451 }
1452
1453 //================================================================================
1454 /*!
1455  * \brief Call _pyFilter.AddUser() if a filter is used as a command arg
1456  */
1457 //================================================================================
1458
1459 // void _pyGen::addFilterUser( Handle(_pyCommand)& theCommand, const Handle(_pyObject)& user )
1460 // {
1461 // No more needed after adding _pyObject::myArgCommands
1462
1463 //   const char filterPrefix[] = "aFilter0x";
1464 //   if ( theCommand->GetString().Search( filterPrefix ) < 1 )
1465 //     return;
1466
1467 //   for ( int i = theCommand->GetNbArgs(); i > 0; --i )
1468 //   {
1469 //     const _AString & arg = theCommand->GetArg( i );
1470 //     // NOT TREATED CASE: arg == "[something, aFilter0x36a2f60]"
1471 //     if ( arg.Search( filterPrefix ) != 1 )
1472 //       continue;
1473
1474 //     Handle(_pyFilter) filter = Handle(_pyFilter)::DownCast( FindObject( arg ));
1475 //     if ( !filter.IsNull() )
1476 //     {
1477 //       filter->AddUser( user );
1478 //       if ( !filter->GetNewID().IsEmpty() )
1479 //         theCommand->SetArg( i, filter->GetNewID() );
1480 //     }
1481 //   }
1482 //}
1483
1484 //================================================================================
1485 /*!
1486  * \brief Set command be last in list of commands
1487  *  \param theCmd - Command to be last
1488  */
1489 //================================================================================
1490
1491 Handle(_pyCommand)& _pyGen::GetLastCommand()
1492 {
1493   return myLastCommand;
1494 }
1495
1496 //================================================================================
1497 /*!
1498  * \brief Set method to access to object wrapped with python class
1499  *  \param theID - The wrapped object entry
1500  *  \param theMethod - The accessor method
1501  */
1502 //================================================================================
1503
1504 void _pyGen::SetAccessorMethod(const _pyID& theID, const char* theMethod )
1505 {
1506   myID2AccessorMethod.Bind( theID, (char*) theMethod );
1507 }
1508
1509 //================================================================================
1510 /*!
1511  * \brief Generated new ID for object and assign with existing name
1512  *  \param theID - ID of existing object
1513  */
1514 //================================================================================
1515
1516 _pyID _pyGen::GenerateNewID( const _pyID& theID )
1517 {
1518   int index = 1;
1519   _pyID aNewID;
1520   do {
1521     aNewID = theID + _pyID( ":" ) + _pyID( index++ );
1522   }
1523   while ( myObjectNames.IsBound( aNewID ) );
1524
1525   if ( myObjectNames.IsBound( theID ) )
1526     myObjectNames.Bind( aNewID, ( myObjectNames.Find( theID ) + _pyID( "_" ) + _pyID( index-1 ) ) );
1527   else
1528     myObjectNames.Bind( aNewID, ( _pyID( "A" ) + aNewID ) );
1529   return aNewID;
1530 }
1531
1532 //================================================================================
1533 /*!
1534  * \brief Stores theObj in myObjects
1535  */
1536 //================================================================================
1537
1538 bool _pyGen::AddObject( Handle(_pyObject)& theObj )
1539 {
1540   if ( theObj.IsNull() ) return false;
1541
1542   CheckObjectIsReCreated( theObj );
1543
1544   bool add;
1545
1546   if ( theObj->IsKind( STANDARD_TYPE( _pyMesh ))) {
1547     add = myMeshes.insert( make_pair( theObj->GetID(),
1548                                       Handle(_pyMesh)::DownCast( theObj ))).second;
1549   }
1550   else if ( theObj->IsKind( STANDARD_TYPE( _pyMeshEditor ))) {
1551     add = myMeshEditors.insert( make_pair( theObj->GetID(),
1552                                            Handle(_pyMeshEditor)::DownCast( theObj ))).second;
1553   }
1554   else {
1555     add = myObjects.insert( make_pair( theObj->GetID(), theObj )).second;
1556     if ( add ) myOrderedObjects.push_back( theObj );
1557   }
1558   return add;
1559 }
1560
1561 //================================================================================
1562 /*!
1563  * \brief Erases an existing object with the same ID. This method should be called
1564  *        before storing theObj in _pyGen
1565  */
1566 //================================================================================
1567
1568 void _pyGen::CheckObjectIsReCreated( Handle(_pyObject)& theObj )
1569 {
1570   if ( theObj.IsNull() || !_pyCommand::IsStudyEntry( theObj->GetID() ))
1571     return;
1572
1573   const bool isHyp = theObj->IsKind( STANDARD_TYPE( _pyHypothesis ));
1574   Handle(_pyObject) existing;
1575   if( isHyp )
1576     existing = FindHyp( theObj->GetID() );
1577   else
1578     existing = FindObject( theObj->GetID() );
1579   if ( !existing.IsNull() && existing != theObj )
1580   {
1581     existing->SetRemovedFromStudy( true );
1582     existing->ClearCommands();
1583     if ( isHyp )
1584     {
1585       if ( myHypos.count( theObj->GetID() ))
1586         myHypos.erase( theObj->GetID() );
1587     }
1588     else if ( myMeshes.count( theObj->GetID() ))
1589     {
1590       myMeshes.erase( theObj->GetID() );
1591     }
1592     else if ( myObjects.count( theObj->GetID() ))
1593     {
1594       myObjects.erase( theObj->GetID() );
1595     }
1596   }
1597 }
1598
1599 //================================================================================
1600 /*!
1601  * \brief Re-register an object with other ID to make it Process() commands of
1602  * other object having this ID
1603  */
1604 //================================================================================
1605
1606 void _pyGen::SetProxyObject( const _pyID& theID, Handle(_pyObject)& theObj )
1607 {
1608   if ( theObj.IsNull() ) return;
1609
1610   if ( theObj->IsKind( STANDARD_TYPE( _pyMesh )))
1611     myMeshes.insert( make_pair( theID, Handle(_pyMesh)::DownCast( theObj )));
1612
1613   else if ( theObj->IsKind( STANDARD_TYPE( _pyMeshEditor )))
1614     myMeshEditors.insert( make_pair( theID, Handle(_pyMeshEditor)::DownCast( theObj )));
1615
1616   else
1617     myObjects.insert( make_pair( theID, theObj ));
1618 }
1619
1620 //================================================================================
1621 /*!
1622  * \brief Finds a _pyObject by ID
1623  */
1624 //================================================================================
1625
1626 Handle(_pyObject) _pyGen::FindObject( const _pyID& theObjID )  const
1627 {
1628   {
1629     map< _pyID, Handle(_pyObject) >::const_iterator id_obj = myObjects.find( theObjID );
1630     if ( id_obj != myObjects.end() )
1631       return id_obj->second;
1632   }
1633   {
1634     _pyGen* me = const_cast< _pyGen* >( this );
1635     map< _pyID, Handle(_pyMesh) >::iterator id_obj = me->myMeshes.find( theObjID );
1636     if ( id_obj != myMeshes.end() )
1637       return id_obj->second;
1638   }
1639   // {
1640   //   map< _pyID, Handle(_pyMeshEditor) >::const_iterator id_obj = myMeshEditors.find( theObjID );
1641   //   if ( id_obj != myMeshEditors.end() )
1642   //     return id_obj->second;
1643   // }
1644   return Handle(_pyObject)();
1645 }
1646
1647 //================================================================================
1648 /*!
1649  * \brief Check if a study entry is under GEOM component
1650  */
1651 //================================================================================
1652
1653 bool _pyGen::IsGeomObject(const _pyID& theObjID) const
1654 {
1655   if ( myGeomIDNb )
1656   {
1657     return ( myGeomIDIndex <= theObjID.Length() &&
1658              int( theObjID.Value( myGeomIDIndex )) == myGeomIDNb &&
1659              _pyCommand::IsStudyEntry( theObjID ));
1660   }
1661   return false;
1662 }
1663
1664 //================================================================================
1665 /*!
1666  * \brief Returns true if an object is not present in a study
1667  */
1668 //================================================================================
1669
1670 bool _pyGen::IsNotPublished(const _pyID& theObjID) const
1671 {
1672   if ( theObjID.IsEmpty() ) return false;
1673
1674   if ( myObjectNames.IsBound( theObjID ))
1675     return false; // SMESH object is in study
1676
1677   // either the SMESH object is not in study or it is a GEOM object
1678   if ( IsGeomObject( theObjID ))
1679   {
1680     SALOMEDS::SObject_wrap so = myStudy->FindObjectID( theObjID.ToCString() );
1681     if ( so->_is_nil() ) return true;
1682     CORBA::Object_var obj = so->GetObject();
1683     return CORBA::is_nil( obj );
1684   }
1685   return true; // SMESH object not in study
1686 }
1687
1688 //================================================================================
1689 /*!
1690  * \brief Add an object to myRemovedObjIDs that leads to that SetName() for
1691  *        this object is not dumped
1692  *  \param [in] theObjID - entry of the object whose creation command was eliminated
1693  */
1694 //================================================================================
1695
1696 void _pyGen::ObjectCreationRemoved(const _pyID& theObjID)
1697 {
1698   myRemovedObjIDs.insert( theObjID );
1699 }
1700
1701 //================================================================================
1702 /*!
1703  * \brief Return reader of  hypotheses of plugins
1704  */
1705 //================================================================================
1706
1707 Handle( _pyHypothesisReader ) _pyGen::GetHypothesisReader() const
1708 {
1709   if (myHypReader.IsNull() )
1710     ((_pyGen*) this)->myHypReader = new _pyHypothesisReader;
1711
1712   return myHypReader;
1713 }
1714
1715
1716 //================================================================================
1717 /*!
1718  * \brief Mesh created by SMESH_Gen
1719  */
1720 //================================================================================
1721
1722 _pyMesh::_pyMesh(const Handle(_pyCommand) theCreationCmd)
1723   : _pyObject( theCreationCmd ), myGeomNotInStudy( false )
1724 {
1725   if ( theCreationCmd->GetMethod() == "CreateMesh" && theGen->IsNotPublished( GetGeom() ))
1726     myGeomNotInStudy = true;
1727
1728   // convert my creation command --> smeshpy.Mesh(...)
1729   Handle(_pyCommand) creationCmd = GetCreationCmd();
1730   creationCmd->SetObject( SMESH_2smeshpy::SmeshpyName() );
1731   creationCmd->SetMethod( "Mesh" );
1732   theGen->SetAccessorMethod( GetID(), _pyMesh::AccessorMethod() );
1733 }
1734
1735 //================================================================================
1736 /*!
1737  * \brief Mesh created by SMESH_MeshEditor
1738  */
1739 //================================================================================
1740
1741 _pyMesh::_pyMesh(const Handle(_pyCommand) theCreationCmd, const _pyID& meshId):
1742   _pyObject(theCreationCmd,meshId), myGeomNotInStudy(false )
1743 {
1744   if ( theCreationCmd->MethodStartsFrom( "CreateMeshesFrom" ))
1745   {
1746     // this mesh depends on the exported mesh
1747     const TCollection_AsciiString& file = theCreationCmd->GetArg( 1 );
1748     if ( !file.IsEmpty() )
1749     {
1750       ExportedMeshData& exportData = theGen->FindExportedMesh( file );
1751       addFatherMesh( exportData.myMesh );
1752       if ( !exportData.myLastComputeCmd.IsNull() )
1753       {
1754         // restore cleared Compute() by which the exported mesh was generated
1755         exportData.myLastComputeCmd->GetString() = exportData.myLastComputeCmdString;
1756         // protect that Compute() cmd from clearing
1757         if ( exportData.myMesh->myLastComputeCmd == exportData.myLastComputeCmd )
1758           exportData.myMesh->myLastComputeCmd.Nullify();
1759       }
1760     }
1761   }
1762   else if ( theCreationCmd->MethodStartsFrom( "Concatenate" ))
1763   {
1764     // this mesh depends on concatenated meshes
1765     const TCollection_AsciiString& meshIDs = theCreationCmd->GetArg( 1 );
1766     list< _pyID > idList = theCreationCmd->GetStudyEntries( meshIDs );
1767     list< _pyID >::iterator meshID = idList.begin();
1768     for ( ; meshID != idList.end(); ++meshID )
1769       addFatherMesh( *meshID );
1770   }
1771   else if ( theCreationCmd->GetMethod() == "CopyMesh" )
1772   {
1773     // this mesh depends on a copied IdSource
1774     const _pyID& objID = theCreationCmd->GetArg( 1 );
1775     addFatherMesh( objID );
1776   }
1777   else if ( theCreationCmd->GetMethod().Search("MakeMesh") != -1 ||
1778             theCreationCmd->GetMethod() == "MakeBoundaryMesh" ||
1779             theCreationCmd->GetMethod() == "MakeBoundaryElements" )
1780   {
1781     // this mesh depends on a source mesh
1782     // (theCreationCmd is already Process()ed by _pyMeshEditor)
1783     const _pyID& meshID = theCreationCmd->GetObject();
1784     addFatherMesh( meshID );
1785   }
1786     
1787   // convert my creation command
1788   Handle(_pyCommand) creationCmd = GetCreationCmd();
1789   creationCmd->SetObject( SMESH_2smeshpy::SmeshpyName() );
1790   theGen->SetAccessorMethod( meshId, _pyMesh::AccessorMethod() );
1791 }
1792
1793 //================================================================================
1794 /*!
1795  * \brief Convert an IDL API command of SMESH::SMESH_Mesh to a method call of python Mesh
1796  *  \param theCommand - Engine method called for this mesh
1797  */
1798 //================================================================================
1799
1800 void _pyMesh::Process( const Handle(_pyCommand)& theCommand )
1801 {
1802   // some methods of SMESH_Mesh interface needs special conversion
1803   // to methods of Mesh python class
1804   //
1805   // 1. GetSubMesh(geom, name) + AddHypothesis(geom, algo)
1806   //     --> in Mesh_Algorithm.Create(mesh, geom, hypo, so)
1807   // 2. AddHypothesis(geom, hyp)
1808   //     --> in Mesh_Algorithm.Hypothesis(hyp, args, so)
1809   // 3. CreateGroupFromGEOM(type, name, grp)
1810   //     --> in Mesh.Group(grp, name="")
1811   // 4. ExportToMED(f, auto_groups, version)
1812   //     --> in Mesh.ExportMED( f, auto_groups, version )
1813   // 5. etc
1814
1815   const TCollection_AsciiString& method = theCommand->GetMethod();
1816   // ----------------------------------------------------------------------
1817   if ( method == "Compute" ) // in snapshot mode, clear the previous Compute()
1818   {
1819     if ( !theGen->IsToKeepAllCommands() ) // !historical
1820     {
1821       list< Handle(_pyHypothesis) >::iterator hyp;
1822       if ( !myLastComputeCmd.IsNull() )
1823       {
1824         // check if the previously computed mesh has been edited,
1825         // if so then we do not clear the previous Compute()
1826         bool toClear = true;
1827         if ( myLastComputeCmd->GetMethod() == "Compute" )
1828         {
1829           list< Handle(_pyMeshEditor)>::iterator e = myEditors.begin();
1830           for ( ; e != myEditors.end() && toClear; ++e )
1831           {
1832             list< Handle(_pyCommand)>& cmds = (*e)->GetProcessedCmds();
1833             list< Handle(_pyCommand) >::reverse_iterator cmd = cmds.rbegin();
1834             if ( cmd != cmds.rend() &&
1835                  (*cmd)->GetOrderNb() > myLastComputeCmd->GetOrderNb() )
1836               toClear = false;
1837           }
1838         }
1839         if ( toClear )
1840         {
1841           // clear hyp commands called before myLastComputeCmd
1842           for ( hyp = myHypos.begin(); hyp != myHypos.end(); ++hyp )
1843             (*hyp)->ComputeDiscarded( myLastComputeCmd );
1844
1845           myLastComputeCmd->Clear();
1846         }
1847       }
1848       myLastComputeCmd = theCommand;
1849
1850       for ( hyp = myHypos.begin(); hyp != myHypos.end(); ++hyp )
1851         (*hyp)->MeshComputed( myLastComputeCmd );
1852     }
1853     Flush();
1854   }
1855   // ----------------------------------------------------------------------
1856   else if ( method == "Clear" ) // in snapshot mode, clear all previous commands
1857   {
1858     if ( !theGen->IsToKeepAllCommands() ) // !historical
1859     {
1860       int untilCmdNb =
1861         myChildMeshes.empty() ? 0 : myChildMeshes.back()->GetCreationCmd()->GetOrderNb();
1862       // list< Handle(_pyCommand) >::reverse_iterator cmd = myProcessedCmds.rbegin();
1863       // for ( ; cmd != myProcessedCmds.rend() && (*cmd)->GetOrderNb() > untilCmdNb; ++cmd )
1864       //   (*cmd)->Clear();
1865       if ( !myLastComputeCmd.IsNull() )
1866       {
1867         list< Handle(_pyHypothesis) >::iterator hyp;
1868         for ( hyp = myHypos.begin(); hyp != myHypos.end(); ++hyp )
1869           (*hyp)->ComputeDiscarded( myLastComputeCmd );
1870
1871         myLastComputeCmd->Clear();
1872       }
1873
1874       list< Handle(_pyMeshEditor)>::iterator e = myEditors.begin();
1875       for ( ; e != myEditors.end(); ++e )
1876       {
1877         list< Handle(_pyCommand)>& cmds = (*e)->GetProcessedCmds();
1878         list< Handle(_pyCommand) >::reverse_iterator cmd = cmds.rbegin();
1879         for ( ; cmd != cmds.rend() && (*cmd)->GetOrderNb() > untilCmdNb; ++cmd )
1880           if ( !(*cmd)->IsEmpty() )
1881           {
1882             if ( (*cmd)->GetStudyEntries( (*cmd)->GetResultValue() ).empty() ) // no object created
1883               (*cmd)->Clear();
1884           }
1885       }
1886       myLastComputeCmd = theCommand; // to clear Clear() the same way as Compute()
1887     }
1888   }
1889   // ----------------------------------------------------------------------
1890   else if ( method == "GetSubMesh" ) { // collect sub-meshes of the mesh
1891     Handle(_pySubMesh) subMesh = theGen->FindSubMesh( theCommand->GetResultValue() );
1892     if ( !subMesh.IsNull() ) {
1893       subMesh->SetCreator( this );
1894       mySubmeshes.push_back( subMesh );
1895     }
1896   }
1897   // ----------------------------------------------------------------------
1898   else if ( method == "GetSubMeshes" ) { // clear as the command does nothing (0023156)
1899     theCommand->Clear();
1900   }
1901   // ----------------------------------------------------------------------
1902   else if ( method == "AddHypothesis" ) { // mesh.AddHypothesis(geom, HYPO )
1903     myAddHypCmds.push_back( theCommand );
1904     // set mesh to hypo
1905     const _pyID& hypID = theCommand->GetArg( 2 );
1906     Handle(_pyHypothesis) hyp = theGen->FindHyp( hypID );
1907     if ( !hyp.IsNull() ) {
1908       myHypos.push_back( hyp );
1909       if ( hyp->GetMesh().IsEmpty() )
1910         hyp->SetMesh( this->GetID() );
1911     }
1912   }
1913   // ----------------------------------------------------------------------
1914   else if ( method == "CreateGroup" ||
1915             method == "CreateGroupFromGEOM" ||
1916             method == "CreateGroupFromFilter" ||
1917             method == "CreateDimGroup" )
1918   {
1919     Handle(_pyGroup) group = new _pyGroup( theCommand );
1920     myGroups.push_back( group );
1921     theGen->AddObject( group );
1922   }
1923   // ----------------------------------------------------------------------
1924   // update list of groups
1925   else if ( method == "GetGroups" )
1926   {
1927     bool allGroupsRemoved = true;
1928     TCollection_AsciiString grIDs = theCommand->GetResultValue();
1929     list< _pyID >          idList = theCommand->GetStudyEntries( grIDs );
1930     list< _pyID >::iterator  grID = idList.begin();
1931     const size_t nbGroupsBefore = myGroups.size();
1932     Handle(_pyObject) obj;
1933     for ( ; grID != idList.end(); ++grID )
1934     {
1935       obj = theGen->FindObject( *grID );
1936       if ( obj.IsNull() )
1937       {
1938         Handle(_pyGroup) group = new _pyGroup( theCommand, *grID );
1939         theGen->AddObject( group );
1940         myGroups.push_back( group );
1941         obj = group;
1942       }
1943       if ( !obj->CanClear() )
1944         allGroupsRemoved = false;
1945     }
1946     if ( nbGroupsBefore == myGroups.size() ) // no new _pyGroup created
1947       obj->AddProcessedCmd( theCommand ); // to clear theCommand if all groups are removed
1948
1949     if ( !allGroupsRemoved && !theGen->IsToKeepAllCommands() )
1950     {
1951       // check if the preceding command is Compute();
1952       // if GetGroups() is just after Compute(), this can mean that the groups
1953       // were created by some algorithm and hence Compute() should not be discarded
1954       std::list< Handle(_pyCommand) >& cmdList = theGen->GetCommands();
1955       std::list< Handle(_pyCommand) >::iterator cmd = cmdList.begin();
1956       while ( (*cmd)->GetMethod() == "GetGroups" )
1957         ++cmd;
1958       if ( myLastComputeCmd == (*cmd))
1959         // protect last Compute() from clearing by the next Compute()
1960         myLastComputeCmd.Nullify();
1961     }
1962   }
1963   // ----------------------------------------------------------------------
1964   // notify a group about full removal
1965   else if ( method == "RemoveGroupWithContents" ||
1966             method == "RemoveGroup")
1967   {
1968     if ( !theGen->IsToKeepAllCommands() ) { // snapshot mode
1969       const _pyID groupID = theCommand->GetArg( 1 );
1970       Handle(_pyGroup) grp = Handle(_pyGroup)::DownCast( theGen->FindObject( groupID ));
1971       if ( !grp.IsNull() )
1972       {
1973         if ( method == "RemoveGroupWithContents" )
1974           grp->RemovedWithContents();
1975         // to clear RemoveGroup() if the group creation is cleared
1976         grp->AddProcessedCmd( theCommand );
1977       }
1978     }
1979   }
1980   // ----------------------------------------------------------------------
1981   else if ( theCommand->MethodStartsFrom( "Export" ))
1982   {
1983     if ( method == "ExportToMED" ||  // ExportToMED()  --> ExportMED()
1984          method == "ExportToMEDX" )  // ExportToMEDX() --> ExportMED()
1985     {
1986       theCommand->SetMethod( "ExportMED" );
1987       if ( theCommand->GetNbArgs() == 5 )
1988       {
1989         // ExportToMEDX(...,autoDimension) -> ExportToMEDX(...,meshPart=None,autoDimension)
1990         _AString autoDimension = theCommand->GetArg( 5 );
1991         theCommand->SetArg( 5, "None" );
1992         theCommand->SetArg( 6, autoDimension );
1993       }
1994     }
1995     else if ( method == "ExportCGNS" )
1996     { // ExportCGNS(part, ...) -> ExportCGNS(..., part)
1997       _pyID partID = theCommand->GetArg( 1 );
1998       int nbArgs = theCommand->GetNbArgs();
1999       for ( int i = 2; i <= nbArgs; ++i )
2000         theCommand->SetArg( i-1, theCommand->GetArg( i ));
2001       theCommand->SetArg( nbArgs, partID );
2002     }
2003     else if ( method == "ExportGMF" )
2004     { // ExportGMF(part,file,bool) -> ExportCGNS(file, part)
2005       _pyID partID  = theCommand->GetArg( 1 );
2006       _AString file = theCommand->GetArg( 2 );
2007       theCommand->RemoveArgs();
2008       theCommand->SetArg( 1, file );
2009       theCommand->SetArg( 2, partID );
2010     }
2011     else if ( theCommand->MethodStartsFrom( "ExportPartTo" ))
2012     { // ExportPartTo*(part, ...) -> Export*(..., part)
2013       //
2014       // remove "PartTo" from the method
2015       TCollection_AsciiString newMethod = method;
2016       newMethod.Remove( /*where=*/7, /*howmany=*/6 );
2017       theCommand->SetMethod( newMethod );
2018       // make the 1st arg be the last one (or last but three for ExportMED())
2019       _pyID partID = theCommand->GetArg( 1 );
2020       int nbArgs = theCommand->GetNbArgs() - 3 * (newMethod == "ExportMED");
2021       for ( int i = 2; i <= nbArgs; ++i )
2022         theCommand->SetArg( i-1, theCommand->GetArg( i ));
2023       theCommand->SetArg( nbArgs, partID );
2024     }
2025     // remember file name
2026     theGen->AddExportedMesh( theCommand->GetArg( 1 ),
2027                              ExportedMeshData( this, myLastComputeCmd ));
2028   }
2029   // ----------------------------------------------------------------------
2030   else if ( method == "RemoveHypothesis" ) // (geom, hyp)
2031   {
2032     _pyID hypID  = theCommand->GetArg( 2 );
2033     _pyID geomID = theCommand->GetArg( 1 );
2034     bool isLocal = ( geomID != GetGeom() );
2035
2036     // check if this mesh still has corresponding addition command
2037     Handle(_pyCommand) addCmd;
2038     list< Handle(_pyCommand) >::iterator cmd;
2039     list< Handle(_pyCommand) >* addCmds[2] = { &myAddHypCmds, &myNotConvertedAddHypCmds };
2040     for ( int i = 0; i < 2; ++i )
2041     {
2042       list< Handle(_pyCommand )> & addHypCmds = *(addCmds[i]);
2043       for ( cmd = addHypCmds.begin(); cmd != addHypCmds.end(); )
2044       {
2045         bool sameHyp = true;
2046         if ( hypID != (*cmd)->GetArg( 1 ) && hypID != (*cmd)->GetArg( 2 ))
2047           sameHyp = false; // other hyp
2048         if ( (*cmd)->GetNbArgs() == 2 &&
2049              geomID != (*cmd)->GetArg( 1 ) && geomID != (*cmd)->GetArg( 2 ))
2050           sameHyp = false; // other geom
2051         if ( (*cmd)->GetNbArgs() == 1 && isLocal )
2052           sameHyp = false; // other geom
2053         if ( sameHyp )
2054         {
2055           addCmd = *cmd;
2056           cmd    = addHypCmds.erase( cmd );
2057           if ( !theGen->IsToKeepAllCommands() /*&& CanClear()*/ ) {
2058             addCmd->Clear();
2059             theCommand->Clear();
2060           }
2061           else
2062           {
2063             // mesh.AddHypothesis(geom, hyp) --> mesh.AddHypothesis(hyp, geom=0)
2064             addCmd->RemoveArgs();
2065             addCmd->SetArg( 1, hypID );
2066             if ( isLocal )
2067               addCmd->SetArg( 2, geomID );
2068           }
2069         }
2070         else
2071         {
2072           ++cmd;
2073         }
2074       }
2075     }
2076     Handle(_pyHypothesis) hyp = theGen->FindHyp( hypID );
2077     if ( !theCommand->IsEmpty() && !hypID.IsEmpty() ) {
2078       // RemoveHypothesis(geom, hyp) --> RemoveHypothesis( hyp, geom=0 )
2079       _pyID geom = theCommand->GetArg( 1 );
2080       theCommand->RemoveArgs();
2081       theCommand->SetArg( 1, hypID );
2082       if ( geom != GetGeom() )
2083         theCommand->SetArg( 2, geom );
2084     }
2085     // remove hyp from myHypos
2086     myHypos.remove( hyp );
2087   }
2088   // check for SubMesh order commands
2089   else if ( method == "GetMeshOrder" || method == "SetMeshOrder" )
2090   {
2091     // make commands GetSubMesh() returning sub-meshes be before using sub-meshes
2092     // by GetMeshOrder() and SetMeshOrder(), since by defalut GetSubMesh()
2093     // commands are moved at the end of the script
2094     TCollection_AsciiString subIDs =
2095       ( method == "SetMeshOrder" ) ? theCommand->GetArg(1) : theCommand->GetResultValue();
2096     list< _pyID > idList = theCommand->GetStudyEntries( subIDs );
2097     list< _pyID >::iterator subID = idList.begin();
2098     for ( ; subID != idList.end(); ++subID )
2099     {
2100       Handle(_pySubMesh) subMesh = theGen->FindSubMesh( *subID );
2101       if ( !subMesh.IsNull() )
2102         subMesh->Process( theCommand ); // it moves GetSubMesh() before theCommand
2103     }
2104   }
2105   // add accessor method if necessary
2106   else
2107   {
2108     if ( NeedMeshAccess( theCommand ))
2109       // apply theCommand to the mesh wrapped by smeshpy mesh
2110       AddMeshAccess( theCommand );
2111   }
2112 }
2113
2114 //================================================================================
2115 /*!
2116  * \brief Return True if addition of accesor method is needed
2117  */
2118 //================================================================================
2119
2120 bool _pyMesh::NeedMeshAccess( const Handle(_pyCommand)& theCommand )
2121 {
2122   // names of SMESH_Mesh methods fully equal to methods of python class Mesh,
2123   // so no conversion is needed for them at all:
2124   static TStringSet sameMethods;
2125   if ( sameMethods.empty() ) {
2126     const char * names[] =
2127       { "ExportDAT","ExportUNV","ExportSTL","ExportSAUV", "RemoveGroup","RemoveGroupWithContents",
2128         "GetGroups","UnionGroups","IntersectGroups","CutGroups","CreateDimGroup","GetLog","GetId",
2129         "ClearLog","GetStudyId","HasDuplicatedGroupNamesMED","GetMEDMesh","NbNodes","NbElements",
2130         "NbEdges","NbEdgesOfOrder","NbFaces","NbFacesOfOrder","NbTriangles",
2131         "NbTrianglesOfOrder","NbQuadrangles","NbQuadranglesOfOrder","NbPolygons","NbVolumes",
2132         "NbVolumesOfOrder","NbTetras","NbTetrasOfOrder","NbHexas","NbHexasOfOrder",
2133         "NbPyramids","NbPyramidsOfOrder","NbPrisms","NbPrismsOfOrder","NbPolyhedrons",
2134         "NbSubMesh","GetElementsId","GetElementsByType","GetNodesId","GetElementType",
2135         "GetSubMeshElementsId","GetSubMeshNodesId","GetSubMeshElementType","Dump","GetNodeXYZ",
2136         "GetNodeInverseElements","GetShapeID","GetShapeIDForElem","GetElemNbNodes",
2137         "GetElemNode","IsMediumNode","IsMediumNodeOfAnyElem","ElemNbEdges","ElemNbFaces",
2138         "GetElemFaceNodes", "GetFaceNormal", "FindElementByNodes",
2139         "IsPoly","IsQuadratic","BaryCenter","GetHypothesisList", "SetAutoColor", "GetAutoColor",
2140         "Clear", "ConvertToStandalone", "GetMeshOrder", "SetMeshOrder"
2141         ,"" }; // <- mark of end
2142     sameMethods.Insert( names );
2143   }
2144
2145   return !sameMethods.Contains( theCommand->GetMethod() );
2146 }
2147
2148 //================================================================================
2149 /*!
2150  * \brief Convert creation and addition of all algos and hypos
2151  */
2152 //================================================================================
2153
2154 void _pyMesh::Flush()
2155 {
2156   {
2157     // get the meshes this mesh depends on via hypotheses
2158     list< Handle(_pyMesh) > fatherMeshes;
2159     list< Handle(_pyHypothesis) >::iterator hyp = myHypos.begin();
2160     for ( ; hyp != myHypos.end(); ++hyp )
2161       if ( ! (*hyp)->GetReferredMeshesAndGeom( fatherMeshes ))
2162         myGeomNotInStudy = true;
2163
2164     list< Handle(_pyMesh) >::iterator m = fatherMeshes.begin();
2165     for ( ; m != fatherMeshes.end(); ++m )
2166       addFatherMesh( *m );
2167     // if ( removedGeom )
2168     //     SetRemovedFromStudy(); // as referred geometry not in study
2169   }
2170   if ( myGeomNotInStudy )
2171     return;
2172
2173   list < Handle(_pyCommand) >::iterator cmd;
2174
2175   // try to convert algo addition like this:
2176   // mesh.AddHypothesis(geom, ALGO ) --> ALGO = mesh.Algo()
2177   for ( cmd = myAddHypCmds.begin(); cmd != myAddHypCmds.end(); ++cmd )
2178   {
2179     Handle(_pyCommand) addCmd = *cmd;
2180
2181     _pyID algoID = addCmd->GetArg( 2 );
2182     Handle(_pyHypothesis) algo = theGen->FindHyp( algoID );
2183     if ( algo.IsNull() || !algo->IsAlgo() )
2184       continue;
2185
2186     // check and create new algorithm instance if it is already wrapped
2187     if ( algo->IsWrapped() ) {
2188       _pyID localAlgoID = theGen->GenerateNewID( algoID );
2189       TCollection_AsciiString aNewCmdStr = addCmd->GetIndentation() + localAlgoID +
2190         TCollection_AsciiString( " = " ) + theGen->GetID() +
2191         TCollection_AsciiString( ".CreateHypothesis( \"" ) + algo->GetAlgoType() +
2192         TCollection_AsciiString( "\" )" );
2193
2194       Handle(_pyCommand) newCmd = theGen->AddCommand( aNewCmdStr );
2195       Handle(_pyAlgorithm) newAlgo = Handle(_pyAlgorithm)::DownCast(theGen->FindHyp( localAlgoID ));
2196       if ( !newAlgo.IsNull() ) {
2197         newAlgo->Assign( algo, this->GetID() );
2198         newAlgo->SetCreationCmd( newCmd );
2199         algo = newAlgo;
2200         // set algorithm creation
2201         theGen->SetCommandBefore( newCmd, addCmd );
2202         myHypos.push_back( newAlgo );
2203         if ( !myLastComputeCmd.IsNull() &&
2204              newCmd->GetOrderNb() == myLastComputeCmd->GetOrderNb() + 1)
2205           newAlgo->MeshComputed( myLastComputeCmd );
2206       }
2207       else
2208         newCmd->Clear();
2209     }
2210     _pyID geom = addCmd->GetArg( 1 );
2211     bool isLocalAlgo = ( geom != GetGeom() );
2212
2213     // try to convert
2214     if ( algo->Addition2Creation( addCmd, this->GetID() )) // OK
2215     {
2216       // wrapped algo is created after mesh creation
2217       GetCreationCmd()->AddDependantCmd( addCmd );
2218
2219       if ( isLocalAlgo ) {
2220         // mesh.AddHypothesis(geom, ALGO ) --> mesh.AlgoMethod(geom)
2221         addCmd->SetArg( addCmd->GetNbArgs() + 1,
2222                         TCollection_AsciiString( "geom=" ) + geom );
2223         // sm = mesh.GetSubMesh(geom, name) --> sm = ALGO.GetSubMesh()
2224         list < Handle(_pySubMesh) >::iterator smIt;
2225         for ( smIt = mySubmeshes.begin(); smIt != mySubmeshes.end(); ++smIt ) {
2226           Handle(_pySubMesh) subMesh = *smIt;
2227           Handle(_pyCommand) subCmd = subMesh->GetCreationCmd();
2228           if ( geom == subCmd->GetArg( 1 )) {
2229             subCmd->SetObject( algo->GetID() );
2230             subCmd->RemoveArgs();
2231             subMesh->SetCreator( algo );
2232           }
2233         }
2234       }
2235     }
2236     else // KO - ALGO was already created
2237     {
2238       // mesh.AddHypothesis(geom, ALGO) --> mesh.AddHypothesis(ALGO, geom=0)
2239       addCmd->RemoveArgs();
2240       addCmd->SetArg( 1, algoID );
2241       if ( isLocalAlgo )
2242         addCmd->SetArg( 2, geom );
2243       myNotConvertedAddHypCmds.push_back( addCmd );
2244     }
2245   }
2246
2247   // try to convert hypo addition like this:
2248   // mesh.AddHypothesis(geom, HYPO ) --> HYPO = algo.Hypo()
2249   for ( cmd = myAddHypCmds.begin(); cmd != myAddHypCmds.end(); ++cmd )
2250   {
2251     Handle(_pyCommand) addCmd = *cmd;
2252     _pyID hypID = addCmd->GetArg( 2 );
2253     Handle(_pyHypothesis) hyp = theGen->FindHyp( hypID );
2254     if ( hyp.IsNull() || hyp->IsAlgo() )
2255       continue;
2256     bool converted = hyp->Addition2Creation( addCmd, this->GetID() );
2257     if ( !converted ) {
2258       // mesh.AddHypothesis(geom, HYP) --> mesh.AddHypothesis(HYP, geom=0)
2259       _pyID geom = addCmd->GetArg( 1 );
2260       addCmd->RemoveArgs();
2261       addCmd->SetArg( 1, hypID );
2262       if ( geom != GetGeom() )
2263         addCmd->SetArg( 2, geom );
2264       myNotConvertedAddHypCmds.push_back( addCmd );
2265     }
2266   }
2267
2268   myAddHypCmds.clear();
2269   mySubmeshes.clear();
2270
2271   // flush hypotheses
2272   list< Handle(_pyHypothesis) >::iterator hyp = myHypos.begin();
2273   for ( hyp = myHypos.begin(); hyp != myHypos.end(); ++hyp )
2274     (*hyp)->Flush();
2275 }
2276
2277 //================================================================================
2278 /*!
2279  * \brief Sets myIsPublished of me and of all objects depending on me.
2280  */
2281 //================================================================================
2282
2283 void _pyMesh::SetRemovedFromStudy(const bool isRemoved)
2284 {
2285   _pyObject::SetRemovedFromStudy(isRemoved);
2286
2287   list< Handle(_pySubMesh) >::iterator sm = mySubmeshes.begin();
2288   for ( ; sm != mySubmeshes.end(); ++sm )
2289     (*sm)->SetRemovedFromStudy(isRemoved);
2290
2291   list< Handle(_pyGroup) >::iterator gr = myGroups.begin();
2292   for ( ; gr != myGroups.end(); ++gr )
2293     (*gr)->SetRemovedFromStudy(isRemoved);
2294
2295   list< Handle(_pyMesh) >::iterator m = myChildMeshes.begin();
2296   for ( ; m != myChildMeshes.end(); ++m )
2297     (*m)->SetRemovedFromStudy(isRemoved);
2298
2299   list< Handle(_pyMeshEditor)>::iterator e = myEditors.begin();
2300   for ( ; e != myEditors.end(); ++e )
2301     (*e)->SetRemovedFromStudy(isRemoved);
2302 }
2303
2304 //================================================================================
2305 /*!
2306  * \brief Return true if none of myChildMeshes is in study
2307  */
2308 //================================================================================
2309
2310 bool _pyMesh::CanClear()
2311 {
2312   if ( IsInStudy() )
2313     return false;
2314
2315   list< Handle(_pyMesh) >::iterator m = myChildMeshes.begin();
2316   for ( ; m != myChildMeshes.end(); ++m )
2317     if ( !(*m)->CanClear() )
2318       return false;
2319
2320   return true;
2321 }
2322
2323 //================================================================================
2324 /*!
2325  * \brief Clear my commands and commands of mesh editor
2326  */
2327 //================================================================================
2328
2329 void _pyMesh::ClearCommands()
2330 {
2331   if ( !CanClear() )
2332   {
2333     if ( !IsInStudy() )
2334     {
2335       // mark all sub-objects as not removed, except child meshes
2336       list< Handle(_pyMesh) > children;
2337       children.swap( myChildMeshes );
2338       SetRemovedFromStudy( false );
2339       children.swap( myChildMeshes );
2340     }
2341     return;
2342   }
2343   _pyObject::ClearCommands();
2344
2345   list< Handle(_pySubMesh) >::iterator sm = mySubmeshes.begin();
2346   for ( ; sm != mySubmeshes.end(); ++sm )
2347     (*sm)->ClearCommands();
2348   
2349   list< Handle(_pyGroup) >::iterator gr = myGroups.begin();
2350   for ( ; gr != myGroups.end(); ++gr )
2351     (*gr)->ClearCommands();
2352
2353   list< Handle(_pyMeshEditor)>::iterator e = myEditors.begin();
2354   for ( ; e != myEditors.end(); ++e )
2355     (*e)->ClearCommands();
2356 }
2357
2358 //================================================================================
2359 /*!
2360  * \brief Add a father mesh by ID
2361  */
2362 //================================================================================
2363
2364 void _pyMesh::addFatherMesh( const _pyID& meshID )
2365 {
2366   if ( !meshID.IsEmpty() && meshID != GetID() )
2367     addFatherMesh( Handle(_pyMesh)::DownCast( theGen->FindObject( meshID )));
2368 }
2369
2370 //================================================================================
2371 /*!
2372  * \brief Add a father mesh
2373  */
2374 //================================================================================
2375
2376 void _pyMesh::addFatherMesh( const Handle(_pyMesh)& mesh )
2377 {
2378   if ( !mesh.IsNull() && mesh->GetID() != GetID() )
2379   {
2380     //myFatherMeshes.push_back( mesh );
2381     mesh->myChildMeshes.push_back( this );
2382
2383     // protect last Compute() from clearing by the next Compute()
2384     mesh->myLastComputeCmd.Nullify();
2385   }
2386 }
2387
2388 //================================================================================
2389 /*!
2390  * \brief MeshEditor convert its commands to ones of mesh
2391  */
2392 //================================================================================
2393
2394 _pyMeshEditor::_pyMeshEditor(const Handle(_pyCommand)& theCreationCmd):
2395   _pyObject( theCreationCmd )
2396 {
2397   myMesh = theCreationCmd->GetObject();
2398   myCreationCmdStr = theCreationCmd->GetString();
2399   theCreationCmd->Clear();
2400
2401   Handle(_pyMesh) mesh = ObjectToMesh( theGen->FindObject( myMesh ));
2402   if ( !mesh.IsNull() )
2403     mesh->AddEditor( this );
2404 }
2405
2406 //================================================================================
2407 /*!
2408  * \brief convert its commands to ones of mesh
2409  */
2410 //================================================================================
2411
2412 void _pyMeshEditor::Process( const Handle(_pyCommand)& theCommand)
2413 {
2414   // Names of SMESH_MeshEditor methods fully equal to methods of the python class Mesh, so
2415   // commands calling these methods are converted to calls of Mesh methods without
2416   // additional modifs, only object is changed from MeshEditor to Mesh.
2417   static TStringSet sameMethods;
2418   if ( sameMethods.empty() ) {
2419     const char * names[] = {
2420       "RemoveElements","RemoveNodes","RemoveOrphanNodes",
2421       "AddNode","Add0DElement","AddEdge","AddFace","AddPolygonalFace","AddBall",
2422       "AddVolume","AddPolyhedralVolume","AddPolyhedralVolumeByFaces",
2423       "MoveNode", "MoveClosestNodeToPoint",
2424       "InverseDiag","DeleteDiag","Reorient","ReorientObject","Reorient2DBy3D",
2425       "TriToQuad","TriToQuadObject", "QuadTo4Tri", "SplitQuad","SplitQuadObject",
2426       "BestSplit","Smooth","SmoothObject","SmoothParametric","SmoothParametricObject",
2427       "ConvertToQuadratic","ConvertFromQuadratic","RenumberNodes","RenumberElements",
2428       "RotationSweep","RotationSweepObject","RotationSweepObject1D","RotationSweepObject2D",
2429       "ExtrusionSweep","AdvancedExtrusion","ExtrusionSweepObject","ExtrusionSweepObject1D",
2430       "ExtrusionByNormal", "ExtrusionSweepObject2D","ExtrusionAlongPath","ExtrusionAlongPathObject",
2431       "ExtrusionAlongPathX","ExtrusionAlongPathObject1D","ExtrusionAlongPathObject2D",
2432       "ExtrusionSweepObjects","RotationSweepObjects","ExtrusionAlongPathObjects",
2433       "Mirror","MirrorObject","Translate","TranslateObject","Rotate","RotateObject",
2434       "FindCoincidentNodes","MergeNodes","FindEqualElements","FillHole",
2435       "MergeElements","MergeEqualElements","SewFreeBorders","SewConformFreeBorders",
2436       "FindCoincidentFreeBorders", "SewCoincidentFreeBorders",
2437       "SewBorderToSide","SewSideElements","ChangeElemNodes","GetLastCreatedNodes",
2438       "GetLastCreatedElems",
2439       "MirrorMakeMesh","MirrorObjectMakeMesh","TranslateMakeMesh","TranslateObjectMakeMesh",
2440       "Scale","ScaleMakeMesh","RotateMakeMesh","RotateObjectMakeMesh","MakeBoundaryMesh",
2441       "MakeBoundaryElements", "SplitVolumesIntoTetra","SplitHexahedraIntoPrisms",
2442       "DoubleElements","DoubleNodes","DoubleNode","DoubleNodeGroup","DoubleNodeGroups",
2443       "DoubleNodeElem","DoubleNodeElemInRegion","DoubleNodeElemGroup","AffectedElemGroupsInRegion",
2444       "DoubleNodeElemGroupInRegion","DoubleNodeElemGroups","DoubleNodeElemGroupsInRegion",
2445       "DoubleNodesOnGroupBoundaries","CreateFlatElementsOnFacesGroups","CreateHoleSkin"
2446       ,"" }; // <- mark of the end
2447     sameMethods.Insert( names );
2448   }
2449
2450   // names of SMESH_MeshEditor commands in which only a method name must be replaced
2451   TStringMap diffMethods;
2452   if ( diffMethods.empty() ) {
2453     const char * orig2newName[] = {
2454       // original name --------------> new name
2455       "ExtrusionAlongPathObjX"      , "ExtrusionAlongPathX",
2456       "FindCoincidentNodesOnPartBut", "FindCoincidentNodesOnPart",
2457       "ConvertToQuadraticObject"    , "ConvertToQuadratic",
2458       "ConvertFromQuadraticObject"  , "ConvertFromQuadratic",
2459       "Create0DElementsOnAllNodes"  , "Add0DElementsToAllNodes",
2460       ""};// <- mark of the end
2461     diffMethods.Insert( orig2newName );
2462   }
2463
2464   // names of SMESH_MeshEditor methods which differ from methods of Mesh class
2465   // only by last two arguments
2466   static TStringSet diffLastTwoArgsMethods;
2467   if (diffLastTwoArgsMethods.empty() ) {
2468     const char * names[] = {
2469       "MirrorMakeGroups","MirrorObjectMakeGroups",
2470       "TranslateMakeGroups","TranslateObjectMakeGroups","ScaleMakeGroups",
2471       "RotateMakeGroups","RotateObjectMakeGroups",
2472       ""};// <- mark of the end
2473     diffLastTwoArgsMethods.Insert( names );
2474   }
2475
2476   // only a method name is to change?
2477   const TCollection_AsciiString & method = theCommand->GetMethod();
2478   bool isPyMeshMethod = sameMethods.Contains( method );
2479   if ( !isPyMeshMethod )
2480   {
2481     TCollection_AsciiString newMethod = diffMethods.Value( method );
2482     if (( isPyMeshMethod = ( newMethod.Length() > 0 )))
2483       theCommand->SetMethod( newMethod );
2484   }
2485   // ConvertToBiQuadratic(...) -> ConvertToQuadratic(...,True)
2486   if ( !isPyMeshMethod && (method == "ConvertToBiQuadratic" || method == "ConvertToBiQuadraticObject") )
2487   {
2488     isPyMeshMethod = true;
2489     theCommand->SetMethod( method.SubString( 1, 9) + method.SubString( 12, method.Length()));
2490     theCommand->SetArg( theCommand->GetNbArgs() + 1, "True" );
2491   }
2492
2493   if ( !isPyMeshMethod )
2494   {
2495     // Replace SMESH_MeshEditor "*MakeGroups" functions by the Mesh
2496     // functions with the flag "theMakeGroups = True" like:
2497     // SMESH_MeshEditor.CmdMakeGroups => Mesh.Cmd(...,True)
2498     int pos = method.Search("MakeGroups");
2499     if( pos != -1)
2500     {
2501       isPyMeshMethod = true;
2502       bool is0DmethId  = ( method == "ExtrusionSweepMakeGroups0D" );
2503       bool is0DmethObj = ( method == "ExtrusionSweepObject0DMakeGroups");
2504
2505       // 1. Remove "MakeGroups" from the Command
2506       TCollection_AsciiString aMethod = theCommand->GetMethod();
2507       int nbArgsToAdd = diffLastTwoArgsMethods.Contains(aMethod) ? 2 : 1;
2508       
2509       if(is0DmethObj)
2510         pos = pos-2;  //Remove "0D" from the Command too
2511       aMethod.Trunc(pos-1);
2512       theCommand->SetMethod(aMethod);
2513
2514       // 2. And add last "True" argument(s)
2515       while(nbArgsToAdd--)
2516         theCommand->SetArg(theCommand->GetNbArgs()+1,"True");
2517       if( is0DmethId || is0DmethObj )
2518         theCommand->SetArg(theCommand->GetNbArgs()+1,"True");
2519     }
2520   }
2521
2522   // ExtrusionSweep0D()       -> ExtrusionSweep()
2523   // ExtrusionSweepObject0D() -> ExtrusionSweepObject()
2524   if ( !isPyMeshMethod && ( method == "ExtrusionSweep0D"  ||
2525                             method == "ExtrusionSweepObject0D" ))
2526   {
2527     isPyMeshMethod = true;
2528     theCommand->SetMethod( method.SubString( 1, method.Length()-2));
2529     theCommand->SetArg(theCommand->GetNbArgs()+1,"False");  //sets flag "MakeGroups = False"
2530     theCommand->SetArg(theCommand->GetNbArgs()+1,"True");  //sets flag "IsNode = True"
2531   }
2532
2533   // DoubleNode...New(...) -> DoubleNode...(...,True)
2534   if ( !isPyMeshMethod && ( method == "DoubleNodeElemGroupNew"  ||
2535                             method == "DoubleNodeElemGroupsNew" ||
2536                             method == "DoubleNodeGroupNew"      ||
2537                             method == "DoubleNodeGroupsNew"     ||
2538                             method == "DoubleNodeElemGroup2New" ||
2539                             method == "DoubleNodeElemGroups2New"))
2540   {
2541     isPyMeshMethod = true;
2542     const int excessLen = 3 + int( method.Value( method.Length()-3 ) == '2' );
2543     theCommand->SetMethod( method.SubString( 1, method.Length()-excessLen));
2544     if ( excessLen == 3 )
2545     {
2546       theCommand->SetArg(theCommand->GetNbArgs()+1,"True");
2547     }
2548     else if ( theCommand->GetArg(4) == "0" ||
2549               theCommand->GetArg(5) == "0" )
2550     {
2551       // [ nothing, Group ] = DoubleNodeGroup2New(,,,False, True) ->
2552       // Group = DoubleNodeGroup2New(,,,False, True)
2553       _pyID groupID = theCommand->GetResultValue( 1 + int( theCommand->GetArg(4) == "0"));
2554       theCommand->SetResultValue( groupID );
2555     }
2556   }
2557   // FindAmongElementsByPoint(meshPart, x, y, z, elementType) ->
2558   // FindElementsByPoint(x, y, z, elementType, meshPart)
2559   if ( !isPyMeshMethod && method == "FindAmongElementsByPoint" )
2560   {
2561     isPyMeshMethod = true;
2562     theCommand->SetMethod( "FindElementsByPoint" );
2563     // make the 1st arg be the last one
2564     _pyID partID = theCommand->GetArg( 1 );
2565     int nbArgs = theCommand->GetNbArgs();
2566     for ( int i = 2; i <= nbArgs; ++i )
2567       theCommand->SetArg( i-1, theCommand->GetArg( i ));
2568     theCommand->SetArg( nbArgs, partID );
2569   }
2570   // Reorient2D( mesh, dir, face, point ) -> Reorient2D( mesh, dir, faceORpoint )
2571   if ( !isPyMeshMethod && method == "Reorient2D" )
2572   {
2573     isPyMeshMethod = true;
2574     _AString mesh  = theCommand->GetArg( 1 );
2575     _AString dir   = theCommand->GetArg( 2 );
2576     _AString face  = theCommand->GetArg( 3 );
2577     _AString point = theCommand->GetArg( 4 );
2578     theCommand->RemoveArgs();
2579     theCommand->SetArg( 1, mesh );
2580     theCommand->SetArg( 2, dir );
2581     if ( face.Value(1) == '-' || face.Value(1) == '0' ) // invalid: face <= 0
2582       theCommand->SetArg( 3, point );
2583     else
2584       theCommand->SetArg( 3, face );
2585   }
2586
2587   if ( method == "QuadToTri" || method == "QuadToTriObject" )
2588   {
2589     isPyMeshMethod = true;
2590     int crit_arg = theCommand->GetNbArgs();
2591     const _AString& crit = theCommand->GetArg(crit_arg);
2592     if (crit.Search("MaxElementLength2D") != -1)
2593       theCommand->SetArg(crit_arg, "");
2594   }
2595
2596   if ( isPyMeshMethod )
2597   {
2598     theCommand->SetObject( myMesh );
2599   }
2600   else
2601   {
2602     // editor creation command is needed only if any editor function is called
2603     theGen->AddMeshAccessorMethod( theCommand ); // for *Object() methods
2604     if ( !myCreationCmdStr.IsEmpty() ) {
2605       GetCreationCmd()->GetString() = myCreationCmdStr;
2606       myCreationCmdStr.Clear();
2607     }
2608   }
2609 }
2610
2611 //================================================================================
2612 /*!
2613  * \brief Return true if my mesh can be removed
2614  */
2615 //================================================================================
2616
2617 bool _pyMeshEditor::CanClear()
2618 {
2619   Handle(_pyMesh) mesh = ObjectToMesh( theGen->FindObject( myMesh ));
2620   return mesh.IsNull() ? true : mesh->CanClear();
2621 }
2622
2623 //================================================================================
2624 /*!
2625  * \brief _pyHypothesis constructor
2626  *  \param theCreationCmd -
2627  */
2628 //================================================================================
2629
2630 _pyHypothesis::_pyHypothesis(const Handle(_pyCommand)& theCreationCmd):
2631   _pyObject( theCreationCmd ), myCurCrMethod(0)
2632 {
2633   myIsAlgo = myIsWrapped = /*myIsConverted = myIsLocal = myDim = */false;
2634 }
2635
2636 //================================================================================
2637 /*!
2638  * \brief Creates algorithm or hypothesis
2639  *  \param theCreationCmd - The engine command creating a hypothesis
2640  * \retval Handle(_pyHypothesis) - Result _pyHypothesis
2641  */
2642 //================================================================================
2643
2644 Handle(_pyHypothesis) _pyHypothesis::NewHypothesis( const Handle(_pyCommand)& theCreationCmd)
2645 {
2646   // theCreationCmd: CreateHypothesis( "theHypType", "theLibName" )
2647   ASSERT (( theCreationCmd->GetMethod() == "CreateHypothesis"));
2648
2649   Handle(_pyHypothesis) hyp, algo;
2650
2651   // "theHypType"
2652   const TCollection_AsciiString & hypTypeQuoted = theCreationCmd->GetArg( 1 );
2653   if ( hypTypeQuoted.IsEmpty() )
2654     return hyp;
2655   // theHypType
2656   TCollection_AsciiString  hypType =
2657     hypTypeQuoted.SubString( 2, hypTypeQuoted.Length() - 1 );
2658
2659   algo = new _pyAlgorithm( theCreationCmd );
2660   hyp  = new _pyHypothesis( theCreationCmd );
2661
2662   if ( hypType == "NumberOfSegments" ) {
2663     hyp = new _pyNumberOfSegmentsHyp( theCreationCmd );
2664     hyp->SetConvMethodAndType( "NumberOfSegments", "Regular_1D");
2665     // arg of SetNumberOfSegments() will become the 1-st arg of hyp creation command
2666     hyp->AddArgMethod( "SetNumberOfSegments" );
2667     // arg of SetScaleFactor() will become the 2-nd arg of hyp creation command
2668     hyp->AddArgMethod( "SetScaleFactor" );
2669     hyp->AddArgMethod( "SetReversedEdges" );
2670     // same for ""CompositeSegment_1D:
2671     hyp->SetConvMethodAndType( "NumberOfSegments", "CompositeSegment_1D");
2672     hyp->AddArgMethod( "SetNumberOfSegments" );
2673     hyp->AddArgMethod( "SetScaleFactor" );
2674     hyp->AddArgMethod( "SetReversedEdges" );
2675   }
2676   else if ( hypType == "SegmentLengthAroundVertex" ) {
2677     hyp = new _pySegmentLengthAroundVertexHyp( theCreationCmd );
2678     hyp->SetConvMethodAndType( "LengthNearVertex", "Regular_1D" );
2679     hyp->AddArgMethod( "SetLength" );
2680     // same for ""CompositeSegment_1D:
2681     hyp->SetConvMethodAndType( "LengthNearVertex", "CompositeSegment_1D");
2682     hyp->AddArgMethod( "SetLength" );
2683   }
2684   else if ( hypType == "LayerDistribution2D" ) {
2685     hyp = new _pyLayerDistributionHypo( theCreationCmd, "Get2DHypothesis" );
2686     hyp->SetConvMethodAndType( "LayerDistribution", "RadialQuadrangle_1D2D");
2687   }
2688   else if ( hypType == "LayerDistribution" ) {
2689     hyp = new _pyLayerDistributionHypo( theCreationCmd, "Get3DHypothesis" );
2690     hyp->SetConvMethodAndType( "LayerDistribution", "RadialPrism_3D");
2691   }
2692   else if ( hypType == "CartesianParameters3D" ) {
2693     hyp = new _pyComplexParamHypo( theCreationCmd );
2694     hyp->SetConvMethodAndType( "SetGrid", "Cartesian_3D");
2695     for ( int iArg = 0; iArg < 4; ++iArg )
2696       hyp->setCreationArg( iArg+1, "[]");
2697     hyp->AddAccumulativeMethod( "SetGrid" );
2698     hyp->AddAccumulativeMethod( "SetGridSpacing" );
2699   }
2700   else
2701   {
2702     hyp = theGen->GetHypothesisReader()->GetHypothesis( hypType, theCreationCmd );
2703   }
2704
2705   return algo->IsValid() ? algo : hyp;
2706 }
2707
2708 //================================================================================
2709 /*!
2710  * \brief Returns true if addition of this hypothesis to a given mesh can be
2711  *        wrapped into hypothesis creation
2712  */
2713 //================================================================================
2714
2715 bool _pyHypothesis::IsWrappable(const _pyID& theMesh) const
2716 {
2717   if ( !myIsWrapped && myMesh == theMesh && IsInStudy() )
2718   {
2719     Handle(_pyObject) pyMesh = theGen->FindObject( myMesh );
2720     if ( !pyMesh.IsNull() && pyMesh->IsInStudy() )
2721       return true;
2722   }
2723   return false;
2724 }
2725
2726 //================================================================================
2727 /*!
2728  * \brief Convert the command adding a hypothesis to mesh into a smesh command
2729  *  \param theCmd - The command like mesh.AddHypothesis( geom, hypo )
2730  *  \param theAlgo - The algo that can create this hypo
2731  * \retval bool - false if the command can't be converted
2732  */
2733 //================================================================================
2734
2735 bool _pyHypothesis::Addition2Creation( const Handle(_pyCommand)& theCmd,
2736                                        const _pyID&              theMesh)
2737 {
2738   ASSERT(( theCmd->GetMethod() == "AddHypothesis" ));
2739
2740   if ( !IsWrappable( theMesh ))
2741     return false;
2742
2743   myGeom = theCmd->GetArg( 1 );
2744
2745   Handle(_pyHypothesis) algo;
2746   if ( !IsAlgo() ) {
2747     // find algo created on myGeom in theMesh
2748     algo = theGen->FindAlgo( myGeom, theMesh, this );
2749     if ( algo.IsNull() )
2750       return false;
2751     // attach hypothesis creation command to be after algo creation command
2752     // because it can be new created instance of algorithm
2753     algo->GetCreationCmd()->AddDependantCmd( theCmd );
2754   }
2755   myIsWrapped = true;
2756
2757   // mesh.AddHypothesis(geom,hyp) --> hyp = <theMesh or algo>.myCreationMethod(args)
2758   theCmd->SetResultValue( GetID() );
2759   theCmd->SetObject( IsAlgo() ? theMesh : algo->GetID());
2760   theCmd->SetMethod( IsAlgo() ? GetAlgoCreationMethod() : GetCreationMethod( algo->GetAlgoType() ));
2761   // set args (geom will be set by _pyMesh calling this method)
2762   theCmd->RemoveArgs();
2763   for ( size_t i = 0; i < myCurCrMethod->myArgs.size(); ++i ) {
2764     if ( !myCurCrMethod->myArgs[ i ].IsEmpty() )
2765       theCmd->SetArg( i+1, myCurCrMethod->myArgs[ i ]);
2766     else
2767       theCmd->SetArg( i+1, "[]");
2768   }
2769   // set a new creation command
2770   GetCreationCmd()->Clear();
2771   // replace creation command by wrapped instance
2772   // please note, that hypothesis attaches to algo creation command (see upper)
2773   SetCreationCmd( theCmd );
2774
2775
2776   // clear commands setting arg values
2777   list < Handle(_pyCommand) >::iterator argCmd = myArgCommands.begin();
2778   for ( ; argCmd != myArgCommands.end(); ++argCmd )
2779     (*argCmd)->Clear();
2780
2781   // set unknown arg commands after hypo creation
2782   Handle(_pyCommand) afterCmd = myIsWrapped ? theCmd : GetCreationCmd();
2783   list<Handle(_pyCommand)>::iterator cmd = myUnusedCommands.begin();
2784   for ( ; cmd != myUnusedCommands.end(); ++cmd ) {
2785     afterCmd->AddDependantCmd( *cmd );
2786   }
2787
2788   return myIsWrapped;
2789 }
2790
2791 //================================================================================
2792 /*!
2793  * \brief Remember hypothesis parameter values
2794  *  \param theCommand - The called hypothesis method
2795  */
2796 //================================================================================
2797
2798 void _pyHypothesis::Process( const Handle(_pyCommand)& theCommand)
2799 {
2800   ASSERT( !myIsAlgo );
2801   if ( !theGen->IsToKeepAllCommands() )
2802     rememberCmdOfParameter( theCommand );
2803   // set args
2804   bool usedCommand = false;
2805   TType2CrMethod::iterator type2meth = myAlgoType2CreationMethod.begin();
2806   for ( ; type2meth != myAlgoType2CreationMethod.end(); ++type2meth )
2807   {
2808     CreationMethod& crMethod = type2meth->second;
2809     for ( size_t i = 0; i < crMethod.myArgMethods.size(); ++i ) {
2810       if ( crMethod.myArgMethods[ i ] == theCommand->GetMethod() ) {
2811         if ( !usedCommand )
2812           myArgCommands.push_back( theCommand );
2813         usedCommand = true;
2814         while ( crMethod.myArgs.size() < i+1 )
2815           crMethod.myArgs.push_back( "None" );
2816         crMethod.myArgs[ i ] = theCommand->GetArg( crMethod.myArgNb[i] );
2817       }
2818     }
2819   }
2820   if ( !usedCommand )
2821     myUnusedCommands.push_back( theCommand );
2822 }
2823
2824 //================================================================================
2825 /*!
2826  * \brief Finish conversion
2827  */
2828 //================================================================================
2829
2830 void _pyHypothesis::Flush()
2831 {
2832   if ( !IsAlgo() )
2833   {
2834     list < Handle(_pyCommand) >::iterator cmd = myArgCommands.begin();
2835     for ( ; cmd != myArgCommands.end(); ++cmd ) {
2836       // Add access to a wrapped mesh
2837       theGen->AddMeshAccessorMethod( *cmd );
2838       // Add access to a wrapped algorithm
2839       theGen->AddAlgoAccessorMethod( *cmd );
2840     }
2841     cmd = myUnusedCommands.begin();
2842     for ( ; cmd != myUnusedCommands.end(); ++cmd ) {
2843       // Add access to a wrapped mesh
2844       theGen->AddMeshAccessorMethod( *cmd );
2845       // Add access to a wrapped algorithm
2846       theGen->AddAlgoAccessorMethod( *cmd );
2847     }
2848   }
2849   // forget previous hypothesis modifications
2850   myArgCommands.clear();
2851   myUnusedCommands.clear();
2852 }
2853
2854 //================================================================================
2855 /*!
2856  * \brief clear creation, arg and unknown commands
2857  */
2858 //================================================================================
2859
2860 void _pyHypothesis::ClearAllCommands()
2861 {
2862   GetCreationCmd()->Clear();
2863   list<Handle(_pyCommand)>::iterator cmd = myArgCommands.begin();
2864   for ( ; cmd != myArgCommands.end(); ++cmd )
2865     ( *cmd )->Clear();
2866   cmd = myUnusedCommands.begin();
2867   for ( ; cmd != myUnusedCommands.end(); ++cmd )
2868     ( *cmd )->Clear();
2869 }
2870
2871
2872 //================================================================================
2873 /*!
2874  * \brief Assign fields of theOther to me except myIsWrapped
2875  */
2876 //================================================================================
2877
2878 void _pyHypothesis::Assign( const Handle(_pyHypothesis)& theOther,
2879                             const _pyID&                 theMesh )
2880 {
2881   // myCreationCmd = theOther->myCreationCmd;
2882   myIsAlgo                  = theOther->myIsAlgo;
2883   myIsWrapped               = false;
2884   myGeom                    = theOther->myGeom;
2885   myMesh                    = theMesh;
2886   myAlgoType2CreationMethod = theOther->myAlgoType2CreationMethod;
2887   myAccumulativeMethods     = theOther->myAccumulativeMethods;
2888   //myUnusedCommands          = theOther->myUnusedCommands;
2889   // init myCurCrMethod
2890   GetCreationMethod( theOther->GetAlgoType() );
2891 }
2892
2893 //================================================================================
2894 /*!
2895  * \brief Analyze my erasability depending on myReferredObjs
2896  */
2897 //================================================================================
2898
2899 bool _pyHypothesis::CanClear()
2900 {
2901   if ( IsInStudy() )
2902   {
2903     list< Handle(_pyObject) >::iterator obj = myReferredObjs.begin();
2904     for ( ; obj != myReferredObjs.end(); ++obj )
2905       if ( (*obj)->CanClear() )
2906         return true;
2907     return false;
2908   }
2909   return true;
2910 }
2911
2912 //================================================================================
2913 /*!
2914  * \brief Clear my commands depending on usage by meshes
2915  */
2916 //================================================================================
2917
2918 void _pyHypothesis::ClearCommands()
2919 {
2920   // if ( !theGen->IsToKeepAllCommands() )
2921   // {
2922   //   bool isUsed = false;
2923   //   int lastComputeOrder = 0;
2924   //   list<Handle(_pyCommand) >::iterator cmd = myComputeCmds.begin();
2925   //   for ( ; cmd != myComputeCmds.end(); ++cmd )
2926   //     if ( ! (*cmd)->IsEmpty() )
2927   //     {
2928   //       isUsed = true;
2929   //       if ( (*cmd)->GetOrderNb() > lastComputeOrder )
2930   //         lastComputeOrder = (*cmd)->GetOrderNb();
2931   //     }
2932   //   if ( !isUsed )
2933   //   {
2934   //     SetRemovedFromStudy( true );
2935   //   }
2936   //   else
2937   //   {
2938   //     // clear my commands invoked after lastComputeOrder
2939   //     // map<TCollection_AsciiString, list< Handle(_pyCommand) > >::iterator m2c;
2940   //     // for ( m2c = myMeth2Commands.begin(); m2c != myMeth2Commands.end(); ++m2c )
2941   //     // {
2942   //     //   list< Handle(_pyCommand)> & cmds = m2c->second;
2943   //     //   if ( !cmds.empty() && cmds.back()->GetOrderNb() > lastComputeOrder )
2944   //     //     cmds.back()->Clear();
2945   //     // }
2946   //   }
2947   // }
2948   _pyObject::ClearCommands();
2949 }
2950
2951 //================================================================================
2952 /*!
2953  * \brief Find arguments that are objects like mesh, group, geometry
2954  *  \param meshes - referred meshes (directly or indirrectly)
2955  *  \retval bool - false if a referred geometry is not in the study
2956  */
2957 //================================================================================
2958
2959 bool _pyHypothesis::GetReferredMeshesAndGeom( list< Handle(_pyMesh) >& meshes )
2960 {
2961   if ( IsAlgo() ) return true;
2962
2963   bool geomPublished = true;
2964   vector< _AString > args;
2965   TType2CrMethod::iterator type2meth = myAlgoType2CreationMethod.begin();
2966   for ( ; type2meth != myAlgoType2CreationMethod.end(); ++type2meth )
2967   {
2968     CreationMethod& crMethod = type2meth->second;
2969     args.insert( args.end(), crMethod.myArgs.begin(), crMethod.myArgs.end());
2970   }
2971   list<Handle(_pyCommand)>::iterator cmd = myUnusedCommands.begin();
2972   for ( ; cmd != myUnusedCommands.end(); ++cmd ) {
2973     for ( int nb = (*cmd)->GetNbArgs(); nb; --nb )
2974       args.push_back( (*cmd)->GetArg( nb ));
2975   }
2976
2977   for ( size_t i = 0; i < args.size(); ++i )
2978   {
2979     list< _pyID > idList = _pyCommand::GetStudyEntries( args[ i ]);
2980     if ( idList.empty() && !args[ i ].IsEmpty() )
2981       idList.push_back( args[ i ]);
2982     list< _pyID >::iterator id = idList.begin();
2983     for ( ; id != idList.end(); ++id )
2984     {
2985       Handle(_pyObject)   obj = theGen->FindObject( *id );
2986       if ( obj.IsNull() ) obj = theGen->FindHyp( *id );
2987       if ( obj.IsNull() )
2988       {
2989         if ( theGen->IsGeomObject( *id ) && theGen->IsNotPublished( *id ))
2990           geomPublished = false;
2991       }
2992       else
2993       {
2994         myReferredObjs.push_back( obj );
2995         Handle(_pyMesh) mesh = ObjectToMesh( obj );
2996         if ( !mesh.IsNull() )
2997           meshes.push_back( mesh );
2998         // prevent clearing not published hyps referred e.g. by "LayerDistribution"
2999         else if ( obj->IsKind( STANDARD_TYPE( _pyHypothesis )) && this->IsInStudy() )
3000           obj->SetRemovedFromStudy( false );
3001       }
3002     }
3003   }
3004   return geomPublished;
3005 }
3006
3007 //================================================================================
3008 /*!
3009  * \brief Remember theCommand setting a parameter
3010  */
3011 //================================================================================
3012
3013 void _pyHypothesis::rememberCmdOfParameter( const Handle(_pyCommand) & theCommand )
3014 {
3015   // parameters are discriminated by method name
3016   _AString method = theCommand->GetMethod();
3017   if ( myAccumulativeMethods.count( method ))
3018     return; // this method adds values and not override the previus value
3019
3020   // discriminate commands setting different parameters via one method
3021   // by passing parameter names like e.g. SetOption("size", "0.2")
3022   if ( theCommand->GetString().FirstLocationInSet( "'\"", 1, theCommand->Length() ) &&
3023        theCommand->GetNbArgs() > 1 )
3024   {
3025     // mangle method by appending a 1st textual arg
3026     for ( int iArg = 1; iArg <= theCommand->GetNbArgs(); ++iArg )
3027     {
3028       const TCollection_AsciiString& arg = theCommand->GetArg( iArg );
3029       if ( arg.Value(1) != '\"' && arg.Value(1) != '\'' ) continue;
3030       if ( !isalpha( arg.Value(2))) continue;
3031       method += arg;
3032       break;
3033     }
3034   }
3035   // parameters are discriminated by method name
3036   list< Handle(_pyCommand)>& cmds = myMeth2Commands[ method /*theCommand->GetMethod()*/ ];
3037   if ( !cmds.empty() && !isCmdUsedForCompute( cmds.back() ))
3038   {
3039     cmds.back()->Clear(); // previous parameter value has not been used
3040     cmds.back() = theCommand;
3041   }
3042   else
3043   {
3044     cmds.push_back( theCommand );
3045   }
3046 }
3047
3048 //================================================================================
3049 /*!
3050  * \brief Return true if a setting parameter command ha been used to compute mesh
3051  */
3052 //================================================================================
3053
3054 bool _pyHypothesis::isCmdUsedForCompute( const Handle(_pyCommand) & cmd,
3055                                          _pyCommand::TAddr          avoidComputeAddr ) const
3056 {
3057   bool isUsed = false;
3058   map< _pyCommand::TAddr, list<Handle(_pyCommand) > >::const_iterator addr2cmds =
3059     myComputeAddr2Cmds.begin();
3060   for ( ; addr2cmds != myComputeAddr2Cmds.end() && !isUsed; ++addr2cmds )
3061   {
3062     if ( addr2cmds->first == avoidComputeAddr ) continue;
3063     const list<Handle(_pyCommand)> & cmds = addr2cmds->second;
3064     isUsed = ( std::find( cmds.begin(), cmds.end(), cmd ) != cmds.end() );
3065   }
3066   return isUsed;
3067 }
3068
3069 //================================================================================
3070 /*!
3071  * \brief Save commands setting parameters as they are used for a mesh computation
3072  */
3073 //================================================================================
3074
3075 void _pyHypothesis::MeshComputed( const Handle(_pyCommand)& theComputeCmd )
3076 {
3077   myComputeCmds.push_back( theComputeCmd );
3078   list<Handle(_pyCommand)>& savedCmds = myComputeAddr2Cmds[ theComputeCmd->GetAddress() ];
3079
3080   map<TCollection_AsciiString, list< Handle(_pyCommand) > >::iterator m2c;
3081   for ( m2c = myMeth2Commands.begin(); m2c != myMeth2Commands.end(); ++m2c )
3082     savedCmds.push_back( m2c->second.back() );
3083 }
3084
3085 //================================================================================
3086 /*!
3087  * \brief Clear commands setting parameters as a mesh computed using them is cleared
3088  */
3089 //================================================================================
3090
3091 void _pyHypothesis::ComputeDiscarded( const Handle(_pyCommand)& theComputeCmd )
3092 {
3093   list<Handle(_pyCommand)>& savedCmds = myComputeAddr2Cmds[ theComputeCmd->GetAddress() ];
3094
3095   list<Handle(_pyCommand)>::iterator cmd = savedCmds.begin();
3096   for ( ; cmd != savedCmds.end(); ++cmd )
3097   {
3098     // check if a cmd has been used to compute another mesh
3099     if ( isCmdUsedForCompute( *cmd, theComputeCmd->GetAddress() ))
3100       continue;
3101     // check if a cmd is a sole command setting its parameter;
3102     // don't use method name for search as it can change
3103     map<TCollection_AsciiString, list<Handle(_pyCommand)> >::iterator
3104       m2cmds = myMeth2Commands.begin();
3105     for ( ; m2cmds != myMeth2Commands.end(); ++m2cmds )
3106     {
3107       list< Handle(_pyCommand)>& cmds = m2cmds->second;
3108       list< Handle(_pyCommand)>::iterator cmdIt = std::find( cmds.begin(), cmds.end(), *cmd );