+void _pyGen::SetAccessorMethod(const _pyID& theID, const char* theMethod )
+{
+ myID2AccessorMethod.Bind( theID, (char*) theMethod );
+}
+
+//================================================================================
+/*!
+ * \brief Generated new ID for object and assign with existing name
+ * \param theID - ID of existing object
+ */
+//================================================================================
+
+_pyID _pyGen::GenerateNewID( const _pyID& theID )
+{
+ int index = 1;
+ _pyID aNewID;
+ do {
+ aNewID = theID + _pyID( ":" ) + _pyID( index++ );
+ }
+ while ( myObjectNames.IsBound( aNewID ) );
+
+ myObjectNames.Bind( aNewID, myObjectNames.IsBound( theID )
+ ? (myObjectNames.Find( theID ) + _pyID( "_" ) + _pyID( index-1 ))
+ : _pyID( "A" ) + aNewID );
+ return aNewID;
+}
+
+//================================================================================
+/*!
+ * \brief Stores theObj in myObjects
+ */
+//================================================================================
+
+void _pyGen::AddObject( Handle(_pyObject)& theObj )
+{
+ if ( theObj.IsNull() ) return;
+
+ if ( theObj->IsKind( STANDARD_TYPE( _pyMesh )))
+ myMeshes.insert( make_pair( theObj->GetID(), Handle(_pyMesh)::DownCast( theObj )));
+
+ else if ( theObj->IsKind( STANDARD_TYPE( _pyMeshEditor )))
+ myMeshEditors.insert( make_pair( theObj->GetID(), Handle(_pyMeshEditor)::DownCast( theObj )));
+
+ else
+ myObjects.insert( make_pair( theObj->GetID(), theObj ));
+}
+
+//================================================================================
+/*!
+ * \brief Re-register an object with other ID to make it Process() commands of
+ * other object having this ID
+ */
+//================================================================================
+
+void _pyGen::SetProxyObject( const _pyID& theID, Handle(_pyObject)& theObj )
+{
+ if ( theObj.IsNull() ) return;
+
+ if ( theObj->IsKind( STANDARD_TYPE( _pyMesh )))
+ myMeshes.insert( make_pair( theID, Handle(_pyMesh)::DownCast( theObj )));
+
+ else if ( theObj->IsKind( STANDARD_TYPE( _pyMeshEditor )))
+ myMeshEditors.insert( make_pair( theID, Handle(_pyMeshEditor)::DownCast( theObj )));
+
+ else
+ myObjects.insert( make_pair( theID, theObj ));
+}
+
+//================================================================================
+/*!
+ * \brief Finds a _pyObject by ID
+ */
+//================================================================================
+
+Handle(_pyObject) _pyGen::FindObject( const _pyID& theObjID ) const
+{
+ {
+ map< _pyID, Handle(_pyObject) >::const_iterator id_obj = myObjects.find( theObjID );
+ if ( id_obj != myObjects.end() )
+ return id_obj->second;
+ }
+ {
+ map< _pyID, Handle(_pyMesh) >::const_iterator id_obj = myMeshes.find( theObjID );
+ if ( id_obj != myMeshes.end() )
+ return id_obj->second;
+ }
+ // {
+ // map< _pyID, Handle(_pyMeshEditor) >::const_iterator id_obj = myMeshEditors.find( theObjID );
+ // if ( id_obj != myMeshEditors.end() )
+ // return id_obj->second;
+ // }
+ return Handle(_pyObject)();
+}
+
+//================================================================================
+/*!
+ * \brief Check if a study entry is under GEOM component
+ */
+//================================================================================
+
+bool _pyGen::IsGeomObject(const _pyID& theObjID) const
+{
+ if ( myGeomIDNb )
+ {
+ return ( myGeomIDIndex <= theObjID.Length() &&
+ int( theObjID.Value( myGeomIDIndex )) == myGeomIDNb &&
+ _pyCommand::IsStudyEntry( theObjID ));
+ }
+ return false;
+}
+
+//================================================================================
+/*!
+ * \brief Returns true if an object is not present in a study
+ */
+//================================================================================
+
+bool _pyGen::IsNotPublished(const _pyID& theObjID) const
+{
+ if ( theObjID.IsEmpty() ) return false;
+
+ if ( myObjectNames.IsBound( theObjID ))
+ return false; // SMESH object is in study
+
+ // either the SMESH object is not in study or it is a GEOM object
+ if ( IsGeomObject( theObjID ))
+ {
+ SALOMEDS::SObject_wrap so = myStudy->FindObjectID( theObjID.ToCString() );
+ if ( so->_is_nil() ) return true;
+ CORBA::Object_var obj = so->GetObject();
+ return CORBA::is_nil( obj );
+ }
+ return true; // SMESH object not in study
+}
+
+//================================================================================
+/*!
+ * \brief Remove object name from myObjectNames that leads to that SetName() for
+ * this object is not dumped
+ * \param [in] theObjID - entry of the object whose creation command was eliminated
+ */
+//================================================================================
+
+void _pyGen::ObjectCreationRemoved(const _pyID& theObjID)
+{
+ myRemovedObjIDs.insert( theObjID );
+}
+
+//================================================================================
+/*!
+ * \brief Return reader of hypotheses of plugins
+ */
+//================================================================================
+
+Handle( _pyHypothesisReader ) _pyGen::GetHypothesisReader() const
+{
+ if (myHypReader.IsNull() )
+ ((_pyGen*) this)->myHypReader = new _pyHypothesisReader;
+
+ return myHypReader;
+}
+
+
+//================================================================================
+/*!
+ * \brief Mesh created by SMESH_Gen
+ */
+//================================================================================
+
+_pyMesh::_pyMesh(const Handle(_pyCommand) theCreationCmd)
+ : _pyObject( theCreationCmd ), myGeomNotInStudy( false )
+{
+ if ( theCreationCmd->GetMethod() == "CreateMesh" && theGen->IsNotPublished( GetGeom() ))
+ myGeomNotInStudy = true;
+
+ // convert my creation command --> smeshpy.Mesh(...)
+ Handle(_pyCommand) creationCmd = GetCreationCmd();
+ creationCmd->SetObject( SMESH_2smeshpy::SmeshpyName() );
+ creationCmd->SetMethod( "Mesh" );
+ theGen->SetAccessorMethod( GetID(), _pyMesh::AccessorMethod() );
+}
+
+//================================================================================
+/*!
+ * \brief Mesh created by SMESH_MeshEditor
+ */
+//================================================================================
+
+_pyMesh::_pyMesh(const Handle(_pyCommand) theCreationCmd, const _pyID& meshId):
+ _pyObject(theCreationCmd,meshId), myGeomNotInStudy(false )
+{
+ if ( theCreationCmd->MethodStartsFrom( "CreateMeshesFrom" ))
+ {
+ // this mesh depends on the exported mesh
+ const TCollection_AsciiString& file = theCreationCmd->GetArg( 1 );
+ if ( !file.IsEmpty() )
+ {
+ ExportedMeshData& exportData = theGen->FindExportedMesh( file );
+ addFatherMesh( exportData.myMesh );
+ if ( !exportData.myLastComputeCmd.IsNull() )
+ {
+ // restore cleared Compute() by which the exported mesh was generated
+ exportData.myLastComputeCmd->GetString() = exportData.myLastComputeCmdString;
+ // protect that Compute() cmd from clearing
+ if ( exportData.myMesh->myLastComputeCmd == exportData.myLastComputeCmd )
+ exportData.myMesh->myLastComputeCmd.Nullify();
+ }
+ }
+ }
+ else if ( theCreationCmd->MethodStartsFrom( "Concatenate" ))
+ {
+ // this mesh depends on concatenated meshes
+ const TCollection_AsciiString& meshIDs = theCreationCmd->GetArg( 1 );
+ list< _pyID > idList = theCreationCmd->GetStudyEntries( meshIDs );
+ list< _pyID >::iterator meshID = idList.begin();
+ for ( ; meshID != idList.end(); ++meshID )
+ addFatherMesh( *meshID );
+ }
+ else if ( theCreationCmd->GetMethod() == "CopyMesh" )
+ {
+ // this mesh depends on a copied IdSource
+ const _pyID& objID = theCreationCmd->GetArg( 1 );
+ addFatherMesh( objID );
+ }
+ else if ( theCreationCmd->GetMethod().Search("MakeMesh") != -1 ||
+ theCreationCmd->GetMethod() == "MakeBoundaryMesh" ||
+ theCreationCmd->GetMethod() == "MakeBoundaryElements" )
+ {
+ // this mesh depends on a source mesh
+ // (theCreationCmd is already Process()ed by _pyMeshEditor)
+ const _pyID& meshID = theCreationCmd->GetObject();
+ addFatherMesh( meshID );
+ }
+
+ // convert my creation command
+ Handle(_pyCommand) creationCmd = GetCreationCmd();
+ creationCmd->SetObject( SMESH_2smeshpy::SmeshpyName() );
+ theGen->SetAccessorMethod( meshId, _pyMesh::AccessorMethod() );
+}
+
+//================================================================================
+/*!
+ * \brief Convert an IDL API command of SMESH::SMESH_Mesh to a method call of python Mesh
+ * \param theCommand - Engine method called for this mesh
+ */
+//================================================================================
+
+void _pyMesh::Process( const Handle(_pyCommand)& theCommand )
+{
+ // some methods of SMESH_Mesh interface needs special conversion
+ // to methods of Mesh python class
+ //
+ // 1. GetSubMesh(geom, name) + AddHypothesis(geom, algo)
+ // --> in Mesh_Algorithm.Create(mesh, geom, hypo, so)
+ // 2. AddHypothesis(geom, hyp)
+ // --> in Mesh_Algorithm.Hypothesis(hyp, args, so)
+ // 3. CreateGroupFromGEOM(type, name, grp)
+ // --> in Mesh.Group(grp, name="")
+ // 4. ExportToMED(f, auto_groups, version)
+ // --> in Mesh.ExportMED( f, auto_groups, version )
+ // 5. etc
+
+ const TCollection_AsciiString& method = theCommand->GetMethod();
+ // ----------------------------------------------------------------------
+ if ( method == "Compute" ) // in snapshot mode, clear the previous Compute()
+ {
+ if ( !theGen->IsToKeepAllCommands() ) // !historical
+ {
+ list< Handle(_pyHypothesis) >::iterator hyp;
+ if ( !myLastComputeCmd.IsNull() )
+ {
+ for ( hyp = myHypos.begin(); hyp != myHypos.end(); ++hyp )
+ (*hyp)->ComputeDiscarded( myLastComputeCmd );
+
+ myLastComputeCmd->Clear();
+ }
+ myLastComputeCmd = theCommand;
+
+ for ( hyp = myHypos.begin(); hyp != myHypos.end(); ++hyp )
+ (*hyp)->MeshComputed( myLastComputeCmd );
+ }
+ Flush();
+ }
+ // ----------------------------------------------------------------------
+ else if ( method == "Clear" ) // in snapshot mode, clear all previous commands
+ {
+ if ( !theGen->IsToKeepAllCommands() ) // !historical
+ {
+ int untilCmdNb =
+ myChildMeshes.empty() ? 0 : myChildMeshes.back()->GetCreationCmd()->GetOrderNb();
+ // list< Handle(_pyCommand) >::reverse_iterator cmd = myProcessedCmds.rbegin();
+ // for ( ; cmd != myProcessedCmds.rend() && (*cmd)->GetOrderNb() > untilCmdNb; ++cmd )
+ // (*cmd)->Clear();
+ if ( !myLastComputeCmd.IsNull() )
+ {
+ list< Handle(_pyHypothesis) >::iterator hyp;
+ for ( hyp = myHypos.begin(); hyp != myHypos.end(); ++hyp )
+ (*hyp)->ComputeDiscarded( myLastComputeCmd );
+
+ myLastComputeCmd->Clear();
+ }
+
+ list< Handle(_pyMeshEditor)>::iterator e = myEditors.begin();
+ for ( ; e != myEditors.end(); ++e )
+ {
+ list< Handle(_pyCommand)>& cmds = (*e)->GetProcessedCmds();
+ list< Handle(_pyCommand) >::reverse_iterator cmd = cmds.rbegin();
+ for ( ; cmd != cmds.rend() && (*cmd)->GetOrderNb() > untilCmdNb; ++cmd )
+ if ( !(*cmd)->IsEmpty() )
+ {
+ if ( (*cmd)->GetStudyEntries( (*cmd)->GetResultValue() ).empty() ) // no object created
+ (*cmd)->Clear();
+ }
+ }
+ myLastComputeCmd = theCommand; // to clear Clear() the same way as Compute()
+ }
+ }
+ // ----------------------------------------------------------------------
+ else if ( method == "GetSubMesh" ) { // collect submeshes of the mesh
+ Handle(_pySubMesh) subMesh = theGen->FindSubMesh( theCommand->GetResultValue() );
+ if ( !subMesh.IsNull() ) {
+ subMesh->SetCreator( this );
+ mySubmeshes.push_back( subMesh );
+ }
+ }
+ // ----------------------------------------------------------------------
+ else if ( method == "AddHypothesis" ) { // mesh.AddHypothesis(geom, HYPO )
+ myAddHypCmds.push_back( theCommand );
+ // set mesh to hypo
+ const _pyID& hypID = theCommand->GetArg( 2 );
+ Handle(_pyHypothesis) hyp = theGen->FindHyp( hypID );
+ if ( !hyp.IsNull() ) {
+ myHypos.push_back( hyp );
+ if ( hyp->GetMesh().IsEmpty() )
+ hyp->SetMesh( this->GetID() );
+ }
+ }
+ // ----------------------------------------------------------------------
+ else if ( method == "CreateGroup" ||
+ method == "CreateGroupFromGEOM" ||
+ method == "CreateGroupFromFilter" )
+ {
+ Handle(_pyGroup) group = new _pyGroup( theCommand );
+ myGroups.push_back( group );
+ theGen->AddObject( group );
+ }
+ // update list of groups
+ else if ( method == "GetGroups" )
+ {
+ TCollection_AsciiString grIDs = theCommand->GetResultValue();
+ list< _pyID > idList = theCommand->GetStudyEntries( grIDs );
+ list< _pyID >::iterator grID = idList.begin();
+ for ( ; grID != idList.end(); ++grID )
+ {
+ Handle(_pyObject) obj = theGen->FindObject( *grID );
+ if ( obj.IsNull() )
+ {
+ Handle(_pyGroup) group = new _pyGroup( theCommand, *grID );
+ theGen->AddObject( group );
+ myGroups.push_back( group );
+ }
+ }
+ }
+ // notify a group about full removal
+ else if ( method == "RemoveGroupWithContents" )
+ {
+ if ( !theGen->IsToKeepAllCommands() ) { // snapshot mode
+ const _pyID groupID = theCommand->GetArg( 1 );
+ Handle(_pyGroup) grp = Handle(_pyGroup)::DownCast( theGen->FindObject( groupID ));
+ if ( !grp.IsNull() )
+ grp->RemovedWithContents();
+ }
+ }
+ // ----------------------------------------------------------------------
+ else if ( theCommand->MethodStartsFrom( "Export" ))
+ {
+ if ( method == "ExportToMED" || // ExportToMED() --> ExportMED()
+ method == "ExportToMEDX" ) { // ExportToMEDX() --> ExportMED()
+ theCommand->SetMethod( "ExportMED" );
+ }
+ else if ( method == "ExportCGNS" )
+ { // ExportCGNS(part, ...) -> ExportCGNS(..., part)
+ _pyID partID = theCommand->GetArg( 1 );
+ int nbArgs = theCommand->GetNbArgs();
+ for ( int i = 2; i <= nbArgs; ++i )
+ theCommand->SetArg( i-1, theCommand->GetArg( i ));
+ theCommand->SetArg( nbArgs, partID );
+ }
+ else if ( method == "ExportGMF" )
+ { // ExportGMF(part,file,bool) -> ExportCGNS(file, part)
+ _pyID partID = theCommand->GetArg( 1 );
+ _AString file = theCommand->GetArg( 2 );
+ theCommand->RemoveArgs();
+ theCommand->SetArg( 1, file );
+ theCommand->SetArg( 2, partID );
+ }
+ else if ( theCommand->MethodStartsFrom( "ExportPartTo" ))
+ { // ExportPartTo*(part, ...) -> Export*(..., part)
+ //
+ // remove "PartTo" from the method
+ TCollection_AsciiString newMethod = method;
+ newMethod.Remove( 7, 6 );
+ theCommand->SetMethod( newMethod );
+ // make the 1st arg be the last one
+ _pyID partID = theCommand->GetArg( 1 );
+ int nbArgs = theCommand->GetNbArgs();
+ for ( int i = 2; i <= nbArgs; ++i )
+ theCommand->SetArg( i-1, theCommand->GetArg( i ));
+ theCommand->SetArg( nbArgs, partID );
+ }
+ // remember file name
+ theGen->AddExportedMesh( theCommand->GetArg( 1 ),
+ ExportedMeshData( this, myLastComputeCmd ));
+ }
+ // ----------------------------------------------------------------------
+ else if ( method == "RemoveHypothesis" ) // (geom, hyp)
+ {
+ _pyID hypID = theCommand->GetArg( 2 );
+ _pyID geomID = theCommand->GetArg( 1 );
+ bool isLocal = ( geomID != GetGeom() );
+
+ // check if this mesh still has corresponding addition command
+ Handle(_pyCommand) addCmd;
+ list< Handle(_pyCommand) >::iterator cmd;
+ list< Handle(_pyCommand) >* addCmds[2] = { &myAddHypCmds, &myNotConvertedAddHypCmds };
+ for ( int i = 0; i < 2; ++i )
+ {
+ list< Handle(_pyCommand )> & addHypCmds = *(addCmds[i]);
+ for ( cmd = addHypCmds.begin(); cmd != addHypCmds.end(); )
+ {
+ bool sameHyp = true;
+ if ( hypID != (*cmd)->GetArg( 1 ) && hypID != (*cmd)->GetArg( 2 ))
+ sameHyp = false; // other hyp
+ if ( (*cmd)->GetNbArgs() == 2 &&
+ geomID != (*cmd)->GetArg( 1 ) && geomID != (*cmd)->GetArg( 2 ))
+ sameHyp = false; // other geom
+ if ( (*cmd)->GetNbArgs() == 1 && isLocal )
+ sameHyp = false; // other geom
+ if ( sameHyp )
+ {
+ addCmd = *cmd;
+ cmd = addHypCmds.erase( cmd );
+ if ( !theGen->IsToKeepAllCommands() ) {
+ addCmd->Clear();
+ theCommand->Clear();
+ }
+ }
+ else
+ {
+ ++cmd;
+ }
+ }
+ }
+ Handle(_pyHypothesis) hyp = theGen->FindHyp( hypID );
+ if ( !theCommand->IsEmpty() && !hypID.IsEmpty() ) {
+ // RemoveHypothesis(geom, hyp) --> RemoveHypothesis( hyp, geom=0 )
+ _pyID geom = theCommand->GetArg( 1 );
+ theCommand->RemoveArgs();
+ theCommand->SetArg( 1, hypID );
+ if ( geom != GetGeom() )
+ theCommand->SetArg( 2, geom );
+ }
+ // remove hyp from myHypos
+ myHypos.remove( hyp );
+ }
+ // check for SubMesh order commands
+ else if ( method == "GetMeshOrder" || method == "SetMeshOrder" )
+ {
+ // make commands GetSubMesh() returning sub-meshes be before using sub-meshes
+ // by GetMeshOrder() and SetMeshOrder(), since by defalut GetSubMesh()
+ // commands are moved at the end of the script
+ TCollection_AsciiString subIDs =
+ ( method == "SetMeshOrder" ) ? theCommand->GetArg(1) : theCommand->GetResultValue();
+ list< _pyID > idList = theCommand->GetStudyEntries( subIDs );
+ list< _pyID >::iterator subID = idList.begin();
+ for ( ; subID != idList.end(); ++subID )
+ {
+ Handle(_pySubMesh) subMesh = theGen->FindSubMesh( *subID );
+ if ( !subMesh.IsNull() )
+ subMesh->Process( theCommand ); // it moves GetSubMesh() before theCommand
+ }
+ }
+ // add accessor method if necessary
+ else
+ {
+ if ( NeedMeshAccess( theCommand ))
+ // apply theCommand to the mesh wrapped by smeshpy mesh
+ AddMeshAccess( theCommand );
+ }
+}
+
+//================================================================================
+/*!
+ * \brief Return True if addition of accesor method is needed
+ */
+//================================================================================
+
+bool _pyMesh::NeedMeshAccess( const Handle(_pyCommand)& theCommand )
+{
+ // names of SMESH_Mesh methods fully equal to methods of python class Mesh,
+ // so no conversion is needed for them at all:
+ static TStringSet sameMethods;
+ if ( sameMethods.empty() ) {
+ const char * names[] =
+ { "ExportDAT","ExportUNV","ExportSTL","ExportSAUV", "RemoveGroup","RemoveGroupWithContents",
+ "GetGroups","UnionGroups","IntersectGroups","CutGroups","GetLog","GetId","ClearLog",
+ "GetStudyId","HasDuplicatedGroupNamesMED","GetMEDMesh","NbNodes","NbElements",
+ "NbEdges","NbEdgesOfOrder","NbFaces","NbFacesOfOrder","NbTriangles",
+ "NbTrianglesOfOrder","NbQuadrangles","NbQuadranglesOfOrder","NbPolygons","NbVolumes",
+ "NbVolumesOfOrder","NbTetras","NbTetrasOfOrder","NbHexas","NbHexasOfOrder",
+ "NbPyramids","NbPyramidsOfOrder","NbPrisms","NbPrismsOfOrder","NbPolyhedrons",
+ "NbSubMesh","GetElementsId","GetElementsByType","GetNodesId","GetElementType",
+ "GetSubMeshElementsId","GetSubMeshNodesId","GetSubMeshElementType","Dump","GetNodeXYZ",
+ "GetNodeInverseElements","GetShapeID","GetShapeIDForElem","GetElemNbNodes",
+ "GetElemNode","IsMediumNode","IsMediumNodeOfAnyElem","ElemNbEdges","ElemNbFaces",
+ "IsPoly","IsQuadratic","BaryCenter","GetHypothesisList", "SetAutoColor", "GetAutoColor",
+ "Clear", "ConvertToStandalone", "GetMeshOrder", "SetMeshOrder"
+ ,"" }; // <- mark of end
+ sameMethods.Insert( names );
+ }
+
+ return !sameMethods.Contains( theCommand->GetMethod() );
+}
+
+//================================================================================
+/*!
+ * \brief Convert creation and addition of all algos and hypos
+ */
+//================================================================================
+
+void _pyMesh::Flush()
+{
+ {
+ // get the meshes this mesh depends on via hypotheses
+ list< Handle(_pyMesh) > fatherMeshes;
+ list< Handle(_pyHypothesis) >::iterator hyp = myHypos.begin();
+ for ( ; hyp != myHypos.end(); ++hyp )
+ if ( ! (*hyp)->GetReferredMeshesAndGeom( fatherMeshes ))
+ myGeomNotInStudy = true;
+
+ list< Handle(_pyMesh) >::iterator m = fatherMeshes.begin();
+ for ( ; m != fatherMeshes.end(); ++m )
+ addFatherMesh( *m );
+ // if ( removedGeom )
+ // SetRemovedFromStudy(); // as reffered geometry not in study
+ }
+ if ( myGeomNotInStudy )
+ return;
+
+ list < Handle(_pyCommand) >::iterator cmd;
+
+ // try to convert algo addition like this:
+ // mesh.AddHypothesis(geom, ALGO ) --> ALGO = mesh.Algo()
+ for ( cmd = myAddHypCmds.begin(); cmd != myAddHypCmds.end(); ++cmd )
+ {
+ Handle(_pyCommand) addCmd = *cmd;
+
+ _pyID algoID = addCmd->GetArg( 2 );
+ Handle(_pyHypothesis) algo = theGen->FindHyp( algoID );
+ if ( algo.IsNull() || !algo->IsAlgo() )
+ continue;
+
+ // check and create new algorithm instance if it is already wrapped
+ if ( algo->IsWrapped() ) {
+ _pyID localAlgoID = theGen->GenerateNewID( algoID );
+ TCollection_AsciiString aNewCmdStr = addCmd->GetIndentation() + localAlgoID +
+ TCollection_AsciiString( " = " ) + theGen->GetID() +
+ TCollection_AsciiString( ".CreateHypothesis( \"" ) + algo->GetAlgoType() +
+ TCollection_AsciiString( "\" )" );
+
+ Handle(_pyCommand) newCmd = theGen->AddCommand( aNewCmdStr );
+ Handle(_pyAlgorithm) newAlgo = Handle(_pyAlgorithm)::DownCast(theGen->FindHyp( localAlgoID ));
+ if ( !newAlgo.IsNull() ) {
+ newAlgo->Assign( algo, this->GetID() );
+ newAlgo->SetCreationCmd( newCmd );
+ algo = newAlgo;
+ // set algorithm creation
+ theGen->SetCommandBefore( newCmd, addCmd );
+ myHypos.push_back( newAlgo );
+ if ( !myLastComputeCmd.IsNull() &&
+ newCmd->GetOrderNb() == myLastComputeCmd->GetOrderNb() + 1)
+ newAlgo->MeshComputed( myLastComputeCmd );
+ }
+ else
+ newCmd->Clear();
+ }
+ _pyID geom = addCmd->GetArg( 1 );
+ bool isLocalAlgo = ( geom != GetGeom() );
+
+ // try to convert
+ if ( algo->Addition2Creation( addCmd, this->GetID() )) // OK
+ {
+ // wrapped algo is created after mesh creation
+ GetCreationCmd()->AddDependantCmd( addCmd );
+
+ if ( isLocalAlgo ) {
+ // mesh.AddHypothesis(geom, ALGO ) --> mesh.AlgoMethod(geom)
+ addCmd->SetArg( addCmd->GetNbArgs() + 1,
+ TCollection_AsciiString( "geom=" ) + geom );
+ // sm = mesh.GetSubMesh(geom, name) --> sm = ALGO.GetSubMesh()
+ list < Handle(_pySubMesh) >::iterator smIt;
+ for ( smIt = mySubmeshes.begin(); smIt != mySubmeshes.end(); ++smIt ) {
+ Handle(_pySubMesh) subMesh = *smIt;
+ Handle(_pyCommand) subCmd = subMesh->GetCreationCmd();
+ if ( geom == subCmd->GetArg( 1 )) {
+ subCmd->SetObject( algo->GetID() );
+ subCmd->RemoveArgs();
+ subMesh->SetCreator( algo );
+ }
+ }
+ }
+ }
+ else // KO - ALGO was already created
+ {
+ // mesh.AddHypothesis(geom, ALGO) --> mesh.AddHypothesis(ALGO, geom=0)
+ addCmd->RemoveArgs();
+ addCmd->SetArg( 1, algoID );
+ if ( isLocalAlgo )
+ addCmd->SetArg( 2, geom );
+ myNotConvertedAddHypCmds.push_back( addCmd );
+ }
+ }
+
+ // try to convert hypo addition like this:
+ // mesh.AddHypothesis(geom, HYPO ) --> HYPO = algo.Hypo()
+ for ( cmd = myAddHypCmds.begin(); cmd != myAddHypCmds.end(); ++cmd )
+ {
+ Handle(_pyCommand) addCmd = *cmd;
+ _pyID hypID = addCmd->GetArg( 2 );
+ Handle(_pyHypothesis) hyp = theGen->FindHyp( hypID );
+ if ( hyp.IsNull() || hyp->IsAlgo() )
+ continue;
+ bool converted = hyp->Addition2Creation( addCmd, this->GetID() );
+ if ( !converted ) {
+ // mesh.AddHypothesis(geom, HYP) --> mesh.AddHypothesis(HYP, geom=0)
+ _pyID geom = addCmd->GetArg( 1 );
+ addCmd->RemoveArgs();
+ addCmd->SetArg( 1, hypID );
+ if ( geom != GetGeom() )
+ addCmd->SetArg( 2, geom );
+ myNotConvertedAddHypCmds.push_back( addCmd );
+ }
+ }
+
+ myAddHypCmds.clear();
+ mySubmeshes.clear();
+
+ // flush hypotheses
+ list< Handle(_pyHypothesis) >::iterator hyp = myHypos.begin();
+ for ( hyp = myHypos.begin(); hyp != myHypos.end(); ++hyp )
+ (*hyp)->Flush();
+}
+
+//================================================================================
+/*!
+ * \brief Sets myIsPublished of me and of all objects depending on me.
+ */
+//================================================================================
+
+void _pyMesh::SetRemovedFromStudy(const bool isRemoved)
+{
+ _pyObject::SetRemovedFromStudy(isRemoved);
+
+ list< Handle(_pySubMesh) >::iterator sm = mySubmeshes.begin();
+ for ( ; sm != mySubmeshes.end(); ++sm )
+ (*sm)->SetRemovedFromStudy(isRemoved);
+
+ list< Handle(_pyGroup) >::iterator gr = myGroups.begin();
+ for ( ; gr != myGroups.end(); ++gr )
+ (*gr)->SetRemovedFromStudy(isRemoved);
+
+ list< Handle(_pyMesh) >::iterator m = myChildMeshes.begin();
+ for ( ; m != myChildMeshes.end(); ++m )
+ (*m)->SetRemovedFromStudy(isRemoved);
+
+ list< Handle(_pyMeshEditor)>::iterator e = myEditors.begin();
+ for ( ; e != myEditors.end(); ++e )
+ (*e)->SetRemovedFromStudy(isRemoved);
+}
+
+//================================================================================
+/*!
+ * \brief Return true if none of myChildMeshes is in study
+ */
+//================================================================================
+
+bool _pyMesh::CanClear()
+{
+ if ( IsInStudy() )
+ return false;
+
+ list< Handle(_pyMesh) >::iterator m = myChildMeshes.begin();
+ for ( ; m != myChildMeshes.end(); ++m )
+ if ( !(*m)->CanClear() )
+ return false;
+
+ return true;
+}
+
+//================================================================================
+/*!
+ * \brief Clear my commands and commands of mesh editor
+ */
+//================================================================================
+
+void _pyMesh::ClearCommands()
+{
+ if ( !CanClear() )
+ {
+ if ( !IsInStudy() )
+ {
+ // mark all sub-objects as not removed, except child meshes
+ list< Handle(_pyMesh) > children;
+ children.swap( myChildMeshes );
+ SetRemovedFromStudy( false );
+ children.swap( myChildMeshes );
+ }
+ return;
+ }
+ _pyObject::ClearCommands();
+
+ list< Handle(_pySubMesh) >::iterator sm = mySubmeshes.begin();
+ for ( ; sm != mySubmeshes.end(); ++sm )
+ (*sm)->ClearCommands();
+
+ list< Handle(_pyGroup) >::iterator gr = myGroups.begin();
+ for ( ; gr != myGroups.end(); ++gr )
+ (*gr)->ClearCommands();
+
+ list< Handle(_pyMeshEditor)>::iterator e = myEditors.begin();
+ for ( ; e != myEditors.end(); ++e )
+ (*e)->ClearCommands();
+}
+
+//================================================================================
+/*!
+ * \brief Add a father mesh by ID
+ */
+//================================================================================
+
+void _pyMesh::addFatherMesh( const _pyID& meshID )
+{
+ if ( !meshID.IsEmpty() && meshID != GetID() )
+ addFatherMesh( Handle(_pyMesh)::DownCast( theGen->FindObject( meshID )));
+}
+
+//================================================================================
+/*!
+ * \brief Add a father mesh
+ */
+//================================================================================
+
+void _pyMesh::addFatherMesh( const Handle(_pyMesh)& mesh )
+{
+ if ( !mesh.IsNull() && mesh->GetID() != GetID() )
+ {
+ //myFatherMeshes.push_back( mesh );
+ mesh->myChildMeshes.push_back( this );
+
+ // protect last Compute() from clearing by the next Compute()
+ mesh->myLastComputeCmd.Nullify();
+ }
+}
+
+//================================================================================
+/*!
+ * \brief MeshEditor convert its commands to ones of mesh
+ */
+//================================================================================
+
+_pyMeshEditor::_pyMeshEditor(const Handle(_pyCommand)& theCreationCmd):
+ _pyObject( theCreationCmd )
+{
+ myMesh = theCreationCmd->GetObject();
+ myCreationCmdStr = theCreationCmd->GetString();
+ theCreationCmd->Clear();
+
+ Handle(_pyMesh) mesh = ObjectToMesh( theGen->FindObject( myMesh ));
+ if ( !mesh.IsNull() )
+ mesh->AddEditor( this );
+}
+
+//================================================================================
+/*!
+ * \brief convert its commands to ones of mesh
+ */
+//================================================================================
+
+void _pyMeshEditor::Process( const Handle(_pyCommand)& theCommand)
+{
+ // names of SMESH_MeshEditor methods fully equal to methods of the python class Mesh, so
+ // commands calling this methods are converted to calls of Mesh methods
+ static TStringSet sameMethods;
+ if ( sameMethods.empty() ) {
+ const char * names[] = {
+ "RemoveElements","RemoveNodes","RemoveOrphanNodes","AddNode","Add0DElement","AddEdge","AddFace","AddPolygonalFace","AddBall",
+ "AddVolume","AddPolyhedralVolume","AddPolyhedralVolumeByFaces","MoveNode", "MoveClosestNodeToPoint",
+ "InverseDiag","DeleteDiag","Reorient","ReorientObject",
+ "TriToQuad","TriToQuadObject", "SplitQuad","SplitQuadObject",
+ "BestSplit","Smooth","SmoothObject","SmoothParametric","SmoothParametricObject",
+ "ConvertToQuadratic","ConvertFromQuadratic","RenumberNodes","RenumberElements",
+ "RotationSweep","RotationSweepObject","RotationSweepObject1D","RotationSweepObject2D",
+ "ExtrusionSweep","AdvancedExtrusion","ExtrusionSweepObject","ExtrusionSweepObject1D","ExtrusionSweepObject2D",
+ "ExtrusionAlongPath","ExtrusionAlongPathObject","ExtrusionAlongPathX",
+ "ExtrusionAlongPathObject1D","ExtrusionAlongPathObject2D",
+ "Mirror","MirrorObject","Translate","TranslateObject","Rotate","RotateObject",
+ "FindCoincidentNodes",/*"FindCoincidentNodesOnPart",*/"MergeNodes","FindEqualElements",
+ "MergeElements","MergeEqualElements","SewFreeBorders","SewConformFreeBorders",
+ "SewBorderToSide","SewSideElements","ChangeElemNodes","GetLastCreatedNodes",
+ "GetLastCreatedElems",
+ "MirrorMakeMesh","MirrorObjectMakeMesh","TranslateMakeMesh",
+ "TranslateObjectMakeMesh","RotateMakeMesh","RotateObjectMakeMesh","MakeBoundaryMesh",
+ "MakeBoundaryElements", "SplitVolumesIntoTetra"
+ ,"" }; // <- mark of the end
+ sameMethods.Insert( names );
+ }
+
+ // names of SMESH_MeshEditor commands in which only a method name must be replaced
+ TStringMap diffMethods;
+ if ( diffMethods.empty() ) {
+ const char * orig2newName[] = {
+ // original name --------------> new name
+ "ExtrusionAlongPathObjX" , "ExtrusionAlongPathX",
+ "FindCoincidentNodesOnPartBut", "FindCoincidentNodesOnPart",
+ "ConvertToQuadraticObject" , "ConvertToQuadratic",
+ "ConvertFromQuadraticObject" , "ConvertFromQuadratic",
+ "Create0DElementsOnAllNodes" , "Add0DElementsToAllNodes",
+ ""};// <- mark of the end
+ diffMethods.Insert( orig2newName );
+ }
+
+ // names of SMESH_MeshEditor methods which differ from methods of Mesh class
+ // only by last two arguments
+ static TStringSet diffLastTwoArgsMethods;
+ if (diffLastTwoArgsMethods.empty() ) {
+ const char * names[] = {
+ "MirrorMakeGroups","MirrorObjectMakeGroups",
+ "TranslateMakeGroups","TranslateObjectMakeGroups",
+ "RotateMakeGroups","RotateObjectMakeGroups",
+ ""};// <- mark of the end
+ diffLastTwoArgsMethods.Insert( names );
+ }
+
+ // only a method name is to change?
+ const TCollection_AsciiString & method = theCommand->GetMethod();
+ bool isPyMeshMethod = sameMethods.Contains( method );
+ if ( !isPyMeshMethod )
+ {
+ TCollection_AsciiString newMethod = diffMethods.Value( method );
+ if (( isPyMeshMethod = ( newMethod.Length() > 0 )))
+ theCommand->SetMethod( newMethod );
+ }
+ // ConvertToBiQuadratic(...) -> ConvertToQuadratic(...,True)
+ if ( !isPyMeshMethod && (method == "ConvertToBiQuadratic" || method == "ConvertToBiQuadraticObject") )
+ {
+ isPyMeshMethod = true;
+ theCommand->SetMethod( method.SubString( 1, 9) + method.SubString( 12, method.Length()));
+ theCommand->SetArg( theCommand->GetNbArgs() + 1, "True" );
+ }
+
+ if ( !isPyMeshMethod )
+ {
+ // Replace SMESH_MeshEditor "*MakeGroups" functions by the Mesh
+ // functions with the flag "theMakeGroups = True" like:
+ // SMESH_MeshEditor.CmdMakeGroups => Mesh.Cmd(...,True)
+ int pos = method.Search("MakeGroups");
+ if( pos != -1)
+ {
+ isPyMeshMethod = true;
+ bool is0DmethId = ( method == "ExtrusionSweepMakeGroups0D" );
+ bool is0DmethObj = ( method == "ExtrusionSweepObject0DMakeGroups");
+
+ // 1. Remove "MakeGroups" from the Command
+ TCollection_AsciiString aMethod = theCommand->GetMethod();
+ int nbArgsToAdd = diffLastTwoArgsMethods.Contains(aMethod) ? 2 : 1;
+
+ if(is0DmethObj)
+ pos = pos-2; //Remove "0D" from the Command too
+ aMethod.Trunc(pos-1);
+ theCommand->SetMethod(aMethod);
+
+ // 2. And add last "True" argument(s)
+ while(nbArgsToAdd--)
+ theCommand->SetArg(theCommand->GetNbArgs()+1,"True");
+ if( is0DmethId || is0DmethObj )
+ theCommand->SetArg(theCommand->GetNbArgs()+1,"True");
+ }
+ }
+
+ // ExtrusionSweep0D() -> ExtrusionSweep()
+ // ExtrusionSweepObject0D() -> ExtrusionSweepObject()
+ if ( !isPyMeshMethod && ( method == "ExtrusionSweep0D" ||
+ method == "ExtrusionSweepObject0D" ))
+ {
+ isPyMeshMethod = true;
+ theCommand->SetMethod( method.SubString( 1, method.Length()-2));
+ theCommand->SetArg(theCommand->GetNbArgs()+1,"False"); //sets flag "MakeGroups = False"
+ theCommand->SetArg(theCommand->GetNbArgs()+1,"True"); //sets flag "IsNode = True"
+ }
+
+ // DoubleNode...New(...) -> DoubleNode...(...,True)
+ if ( !isPyMeshMethod && ( method == "DoubleNodeElemGroupNew" ||
+ method == "DoubleNodeElemGroupsNew" ||
+ method == "DoubleNodeGroupNew" ||
+ method == "DoubleNodeGroupsNew" ||
+ method == "DoubleNodeElemGroup2New" ||
+ method == "DoubleNodeElemGroups2New"))
+ {
+ isPyMeshMethod = true;
+ const int excessLen = 3 + int( method.Value( method.Length()-3 ) == '2' );
+ theCommand->SetMethod( method.SubString( 1, method.Length()-excessLen));
+ if ( excessLen == 3 )
+ {
+ theCommand->SetArg(theCommand->GetNbArgs()+1,"True");
+ }
+ else if ( theCommand->GetArg(4) == "0" ||
+ theCommand->GetArg(5) == "0" )
+ {
+ // [ nothing, Group ] = DoubleNodeGroup2New(,,,False, True) ->
+ // Group = DoubleNodeGroup2New(,,,False, True)
+ _pyID groupID = theCommand->GetResultValue( 1 + int( theCommand->GetArg(4) == "0"));
+ theCommand->SetResultValue( groupID );
+ }
+ }
+ // FindAmongElementsByPoint(meshPart, x, y, z, elementType) ->
+ // FindElementsByPoint(x, y, z, elementType, meshPart)
+ if ( !isPyMeshMethod && method == "FindAmongElementsByPoint" )
+ {
+ isPyMeshMethod = true;
+ theCommand->SetMethod( "FindElementsByPoint" );
+ // make the 1st arg be the last one
+ _pyID partID = theCommand->GetArg( 1 );
+ int nbArgs = theCommand->GetNbArgs();
+ for ( int i = 2; i <= nbArgs; ++i )
+ theCommand->SetArg( i-1, theCommand->GetArg( i ));
+ theCommand->SetArg( nbArgs, partID );
+ }
+ // Reorient2D( mesh, dir, face, point ) -> Reorient2D( mesh, dir, faceORpoint )
+ if ( !isPyMeshMethod && method == "Reorient2D" )
+ {
+ isPyMeshMethod = true;
+ _AString mesh = theCommand->GetArg( 1 );
+ _AString dir = theCommand->GetArg( 2 );
+ _AString face = theCommand->GetArg( 3 );
+ _AString point = theCommand->GetArg( 4 );
+ theCommand->RemoveArgs();
+ theCommand->SetArg( 1, mesh );
+ theCommand->SetArg( 2, dir );
+ if ( face.Value(1) == '-' || face.Value(1) == '0' ) // invalid: face <= 0
+ theCommand->SetArg( 3, point );
+ else
+ theCommand->SetArg( 3, face );
+ }
+
+ if ( method == "QuadToTri" || method == "QuadToTriObject" )
+ {
+ isPyMeshMethod = true;
+ int crit_arg = theCommand->GetNbArgs();
+ const _AString& crit = theCommand->GetArg(crit_arg);
+ if (crit.Search("MaxElementLength2D") != -1)
+ theCommand->SetArg(crit_arg, "");
+ }
+
+ if ( isPyMeshMethod )
+ {
+ theCommand->SetObject( myMesh );
+ }
+ else
+ {
+ // editor creation command is needed only if any editor function is called
+ theGen->AddMeshAccessorMethod( theCommand ); // for *Object() methods
+ if ( !myCreationCmdStr.IsEmpty() ) {
+ GetCreationCmd()->GetString() = myCreationCmdStr;
+ myCreationCmdStr.Clear();
+ }
+ }
+}
+
+//================================================================================
+/*!
+ * \brief Return true if my mesh can be removed
+ */
+//================================================================================
+
+bool _pyMeshEditor::CanClear()
+{
+ Handle(_pyMesh) mesh = ObjectToMesh( theGen->FindObject( myMesh ));
+ return mesh.IsNull() ? true : mesh->CanClear();
+}
+
+//================================================================================
+/*!
+ * \brief _pyHypothesis constructor
+ * \param theCreationCmd -
+ */
+//================================================================================
+
+_pyHypothesis::_pyHypothesis(const Handle(_pyCommand)& theCreationCmd):
+ _pyObject( theCreationCmd ), myCurCrMethod(0)
+{
+ myIsAlgo = myIsWrapped = /*myIsConverted = myIsLocal = myDim = */false;
+}
+
+//================================================================================
+/*!
+ * \brief Creates algorithm or hypothesis
+ * \param theCreationCmd - The engine command creating a hypothesis
+ * \retval Handle(_pyHypothesis) - Result _pyHypothesis
+ */
+//================================================================================
+
+Handle(_pyHypothesis) _pyHypothesis::NewHypothesis( const Handle(_pyCommand)& theCreationCmd)
+{
+ // theCreationCmd: CreateHypothesis( "theHypType", "theLibName" )
+ ASSERT (( theCreationCmd->GetMethod() == "CreateHypothesis"));
+
+ Handle(_pyHypothesis) hyp, algo;
+
+ // "theHypType"
+ const TCollection_AsciiString & hypTypeQuoted = theCreationCmd->GetArg( 1 );
+ if ( hypTypeQuoted.IsEmpty() )
+ return hyp;
+ // theHypType
+ TCollection_AsciiString hypType =
+ hypTypeQuoted.SubString( 2, hypTypeQuoted.Length() - 1 );
+
+ algo = new _pyAlgorithm( theCreationCmd );
+ hyp = new _pyHypothesis( theCreationCmd );
+
+ if ( hypType == "NumberOfSegments" ) {
+ hyp = new _pyNumberOfSegmentsHyp( theCreationCmd );
+ hyp->SetConvMethodAndType( "NumberOfSegments", "Regular_1D");
+ // arg of SetNumberOfSegments() will become the 1-st arg of hyp creation command
+ hyp->AddArgMethod( "SetNumberOfSegments" );
+ // arg of SetScaleFactor() will become the 2-nd arg of hyp creation command
+ hyp->AddArgMethod( "SetScaleFactor" );
+ hyp->AddArgMethod( "SetReversedEdges" );
+ // same for ""CompositeSegment_1D:
+ hyp->SetConvMethodAndType( "NumberOfSegments", "CompositeSegment_1D");
+ hyp->AddArgMethod( "SetNumberOfSegments" );
+ hyp->AddArgMethod( "SetScaleFactor" );
+ hyp->AddArgMethod( "SetReversedEdges" );
+ }
+ else if ( hypType == "SegmentLengthAroundVertex" ) {
+ hyp = new _pySegmentLengthAroundVertexHyp( theCreationCmd );
+ hyp->SetConvMethodAndType( "LengthNearVertex", "Regular_1D" );
+ hyp->AddArgMethod( "SetLength" );
+ // same for ""CompositeSegment_1D:
+ hyp->SetConvMethodAndType( "LengthNearVertex", "CompositeSegment_1D");
+ hyp->AddArgMethod( "SetLength" );
+ }
+ else if ( hypType == "LayerDistribution2D" ) {
+ hyp = new _pyLayerDistributionHypo( theCreationCmd, "Get2DHypothesis" );
+ hyp->SetConvMethodAndType( "LayerDistribution", "RadialQuadrangle_1D2D");
+ }
+ else if ( hypType == "LayerDistribution" ) {
+ hyp = new _pyLayerDistributionHypo( theCreationCmd, "Get3DHypothesis" );
+ hyp->SetConvMethodAndType( "LayerDistribution", "RadialPrism_3D");
+ }
+ else if ( hypType == "CartesianParameters3D" ) {
+ hyp = new _pyComplexParamHypo( theCreationCmd );
+ hyp->SetConvMethodAndType( "SetGrid", "Cartesian_3D");
+ for ( int iArg = 0; iArg < 4; ++iArg )
+ hyp->setCreationArg( iArg+1, "[]");
+ }
+ else
+ {
+ hyp = theGen->GetHypothesisReader()->GetHypothesis( hypType, theCreationCmd );
+ }
+
+ return algo->IsValid() ? algo : hyp;
+}
+
+//================================================================================
+/*!
+ * \brief Returns true if addition of this hypothesis to a given mesh can be
+ * wrapped into hypothesis creation
+ */
+//================================================================================
+
+bool _pyHypothesis::IsWrappable(const _pyID& theMesh) const
+{
+ if ( !myIsWrapped && myMesh == theMesh && IsInStudy() )
+ {
+ Handle(_pyObject) pyMesh = theGen->FindObject( myMesh );
+ if ( !pyMesh.IsNull() && pyMesh->IsInStudy() )
+ return true;
+ }
+ return false;
+}
+
+//================================================================================
+/*!
+ * \brief Convert the command adding a hypothesis to mesh into a smesh command
+ * \param theCmd - The command like mesh.AddHypothesis( geom, hypo )
+ * \param theAlgo - The algo that can create this hypo
+ * \retval bool - false if the command cant be converted
+ */
+//================================================================================
+
+bool _pyHypothesis::Addition2Creation( const Handle(_pyCommand)& theCmd,
+ const _pyID& theMesh)
+{
+ ASSERT(( theCmd->GetMethod() == "AddHypothesis" ));
+
+ if ( !IsWrappable( theMesh ))
+ return false;
+
+ myGeom = theCmd->GetArg( 1 );
+
+ Handle(_pyHypothesis) algo;
+ if ( !IsAlgo() ) {
+ // find algo created on myGeom in theMesh
+ algo = theGen->FindAlgo( myGeom, theMesh, this );
+ if ( algo.IsNull() )
+ return false;
+ // attach hypothesis creation command to be after algo creation command
+ // because it can be new created instance of algorithm
+ algo->GetCreationCmd()->AddDependantCmd( theCmd );
+ }
+ myIsWrapped = true;
+
+ // mesh.AddHypothesis(geom,hyp) --> hyp = <theMesh or algo>.myCreationMethod(args)
+ theCmd->SetResultValue( GetID() );
+ theCmd->SetObject( IsAlgo() ? theMesh : algo->GetID());
+ theCmd->SetMethod( IsAlgo() ? GetAlgoCreationMethod() : GetCreationMethod( algo->GetAlgoType() ));
+ // set args (geom will be set by _pyMesh calling this method)
+ theCmd->RemoveArgs();
+ for ( size_t i = 0; i < myCurCrMethod->myArgs.size(); ++i ) {
+ if ( !myCurCrMethod->myArgs[ i ].IsEmpty() )
+ theCmd->SetArg( i+1, myCurCrMethod->myArgs[ i ]);
+ else
+ theCmd->SetArg( i+1, "[]");
+ }
+ // set a new creation command
+ GetCreationCmd()->Clear();
+ // replace creation command by wrapped instance
+ // please note, that hypothesis attaches to algo creation command (see upper)
+ SetCreationCmd( theCmd );
+
+
+ // clear commands setting arg values
+ list < Handle(_pyCommand) >::iterator argCmd = myArgCommands.begin();
+ for ( ; argCmd != myArgCommands.end(); ++argCmd )
+ (*argCmd)->Clear();
+
+ // set unknown arg commands after hypo creation
+ Handle(_pyCommand) afterCmd = myIsWrapped ? theCmd : GetCreationCmd();
+ list<Handle(_pyCommand)>::iterator cmd = myUnusedCommands.begin();
+ for ( ; cmd != myUnusedCommands.end(); ++cmd ) {
+ afterCmd->AddDependantCmd( *cmd );
+ }
+
+ return myIsWrapped;
+}
+
+//================================================================================
+/*!
+ * \brief Remember hypothesis parameter values
+ * \param theCommand - The called hypothesis method
+ */
+//================================================================================
+
+void _pyHypothesis::Process( const Handle(_pyCommand)& theCommand)
+{
+ ASSERT( !myIsAlgo );
+ if ( !theGen->IsToKeepAllCommands() )
+ rememberCmdOfParameter( theCommand );
+ // set args
+ bool usedCommand = false;
+ TType2CrMethod::iterator type2meth = myAlgoType2CreationMethod.begin();
+ for ( ; type2meth != myAlgoType2CreationMethod.end(); ++type2meth )
+ {
+ CreationMethod& crMethod = type2meth->second;
+ for ( size_t i = 0; i < crMethod.myArgMethods.size(); ++i ) {
+ if ( crMethod.myArgMethods[ i ] == theCommand->GetMethod() ) {
+ if ( !usedCommand )
+ myArgCommands.push_back( theCommand );
+ usedCommand = true;
+ while ( crMethod.myArgs.size() < i+1 )
+ crMethod.myArgs.push_back( "[]" );
+ crMethod.myArgs[ i ] = theCommand->GetArg( crMethod.myArgNb[i] );
+ }
+ }
+ }
+ if ( !usedCommand )
+ myUnusedCommands.push_back( theCommand );
+}
+
+//================================================================================
+/*!
+ * \brief Finish conversion
+ */
+//================================================================================
+
+void _pyHypothesis::Flush()
+{
+ if ( !IsAlgo() )
+ {
+ list < Handle(_pyCommand) >::iterator cmd = myArgCommands.begin();
+ for ( ; cmd != myArgCommands.end(); ++cmd ) {
+ // Add access to a wrapped mesh
+ theGen->AddMeshAccessorMethod( *cmd );
+ // Add access to a wrapped algorithm
+ theGen->AddAlgoAccessorMethod( *cmd );
+ }
+ cmd = myUnusedCommands.begin();
+ for ( ; cmd != myUnusedCommands.end(); ++cmd ) {
+ // Add access to a wrapped mesh
+ theGen->AddMeshAccessorMethod( *cmd );
+ // Add access to a wrapped algorithm
+ theGen->AddAlgoAccessorMethod( *cmd );
+ }
+ }
+ // forget previous hypothesis modifications
+ myArgCommands.clear();
+ myUnusedCommands.clear();
+}
+
+//================================================================================
+/*!
+ * \brief clear creation, arg and unkown commands
+ */
+//================================================================================
+
+void _pyHypothesis::ClearAllCommands()
+{
+ GetCreationCmd()->Clear();
+ list<Handle(_pyCommand)>::iterator cmd = myArgCommands.begin();
+ for ( ; cmd != myArgCommands.end(); ++cmd )
+ ( *cmd )->Clear();
+ cmd = myUnusedCommands.begin();
+ for ( ; cmd != myUnusedCommands.end(); ++cmd )
+ ( *cmd )->Clear();
+}
+
+
+//================================================================================
+/*!
+ * \brief Assign fields of theOther to me except myIsWrapped
+ */
+//================================================================================
+
+void _pyHypothesis::Assign( const Handle(_pyHypothesis)& theOther,
+ const _pyID& theMesh )
+{
+ // myCreationCmd = theOther->myCreationCmd;
+ myIsAlgo = theOther->myIsAlgo;
+ myIsWrapped = false;
+ myGeom = theOther->myGeom;
+ myMesh = theMesh;
+ myAlgoType2CreationMethod = theOther->myAlgoType2CreationMethod;
+ myAccumulativeMethods = theOther->myAccumulativeMethods;
+ //myUnusedCommands = theOther->myUnusedCommands;
+ // init myCurCrMethod
+ GetCreationMethod( theOther->GetAlgoType() );
+}
+
+//================================================================================
+/*!
+ * \brief Analyze my erasability depending on myReferredObjs
+ */
+//================================================================================
+
+bool _pyHypothesis::CanClear()
+{
+ if ( IsInStudy() )
+ {
+ list< Handle(_pyObject) >::iterator obj = myReferredObjs.begin();
+ for ( ; obj != myReferredObjs.end(); ++obj )
+ if ( (*obj)->CanClear() )
+ return true;
+ return false;
+ }
+ return true;
+}
+
+//================================================================================
+/*!
+ * \brief Clear my commands depending on usage by meshes
+ */
+//================================================================================
+
+void _pyHypothesis::ClearCommands()
+{
+ // if ( !theGen->IsToKeepAllCommands() )
+ // {
+ // bool isUsed = false;
+ // int lastComputeOrder = 0;
+ // list<Handle(_pyCommand) >::iterator cmd = myComputeCmds.begin();
+ // for ( ; cmd != myComputeCmds.end(); ++cmd )
+ // if ( ! (*cmd)->IsEmpty() )
+ // {
+ // isUsed = true;
+ // if ( (*cmd)->GetOrderNb() > lastComputeOrder )
+ // lastComputeOrder = (*cmd)->GetOrderNb();
+ // }
+ // if ( !isUsed )
+ // {
+ // SetRemovedFromStudy( true );
+ // }
+ // else
+ // {
+ // // clear my commands invoked after lastComputeOrder
+ // // map<TCollection_AsciiString, list< Handle(_pyCommand) > >::iterator m2c;
+ // // for ( m2c = myMeth2Commands.begin(); m2c != myMeth2Commands.end(); ++m2c )
+ // // {
+ // // list< Handle(_pyCommand)> & cmds = m2c->second;
+ // // if ( !cmds.empty() && cmds.back()->GetOrderNb() > lastComputeOrder )
+ // // cmds.back()->Clear();
+ // // }
+ // }
+ // }
+ _pyObject::ClearCommands();
+}
+
+//================================================================================
+/*!
+ * \brief Find arguments that are objects like mesh, group, geometry
+ * \param meshes - referred meshes (directly or indirrectly)
+ * \retval bool - false if a referred geometry is not in the study
+ */
+//================================================================================
+
+bool _pyHypothesis::GetReferredMeshesAndGeom( list< Handle(_pyMesh) >& meshes )
+{
+ if ( IsAlgo() ) return true;
+
+ bool geomPublished = true;
+ vector< _AString > args;
+ TType2CrMethod::iterator type2meth = myAlgoType2CreationMethod.begin();
+ for ( ; type2meth != myAlgoType2CreationMethod.end(); ++type2meth )
+ {
+ CreationMethod& crMethod = type2meth->second;
+ args.insert( args.end(), crMethod.myArgs.begin(), crMethod.myArgs.end());
+ }
+ list<Handle(_pyCommand)>::iterator cmd = myUnusedCommands.begin();
+ for ( ; cmd != myUnusedCommands.end(); ++cmd ) {
+ for ( int nb = (*cmd)->GetNbArgs(); nb; --nb )
+ args.push_back( (*cmd)->GetArg( nb ));
+ }
+
+ for ( size_t i = 0; i < args.size(); ++i )
+ {
+ list< _pyID > idList = _pyCommand::GetStudyEntries( args[ i ]);
+ if ( idList.empty() && !args[ i ].IsEmpty() )
+ idList.push_back( args[ i ]);
+ list< _pyID >::iterator id = idList.begin();
+ for ( ; id != idList.end(); ++id )
+ {
+ Handle(_pyObject) obj = theGen->FindObject( *id );
+ if ( obj.IsNull() ) obj = theGen->FindHyp( *id );
+ if ( obj.IsNull() )
+ {
+ if ( theGen->IsGeomObject( *id ) && theGen->IsNotPublished( *id ))
+ geomPublished = false;
+ }
+ else
+ {
+ myReferredObjs.push_back( obj );
+ Handle(_pyMesh) mesh = ObjectToMesh( obj );
+ if ( !mesh.IsNull() )
+ meshes.push_back( mesh );
+ // prevent clearing not published hyps referred e.g. by "LayerDistribution"
+ else if ( obj->IsKind( STANDARD_TYPE( _pyHypothesis )) && this->IsInStudy() )
+ obj->SetRemovedFromStudy( false );
+ }
+ }
+ }
+ return geomPublished;
+}
+
+//================================================================================
+/*!
+ * \brief Remember theCommand setting a parameter
+ */
+//================================================================================
+
+void _pyHypothesis::rememberCmdOfParameter( const Handle(_pyCommand) & theCommand )
+{
+ // parameters are discriminated by method name
+ _AString method = theCommand->GetMethod();
+ if ( myAccumulativeMethods.count( method ))
+ return; // this method adds values and not override the previus value
+
+ // discriminate commands setting different parameters via one method
+ // by passing parameter names like e.g. SetOption("size", "0.2")
+ if ( theCommand->GetString().FirstLocationInSet( "'\"", 1, theCommand->Length() ) &&
+ theCommand->GetNbArgs() > 1 )
+ {
+ // mangle method by appending a 1st textual arg
+ for ( int iArg = 1; iArg <= theCommand->GetNbArgs(); ++iArg )
+ {
+ const TCollection_AsciiString& arg = theCommand->GetArg( iArg );
+ if ( arg.Value(1) != '\"' && arg.Value(1) != '\'' ) continue;
+ if ( !isalpha( arg.Value(2))) continue;
+ method += arg;
+ break;
+ }
+ }
+ // parameters are discriminated by method name
+ list< Handle(_pyCommand)>& cmds = myMeth2Commands[ method /*theCommand->GetMethod()*/ ];
+ if ( !cmds.empty() && !isCmdUsedForCompute( cmds.back() ))
+ {
+ cmds.back()->Clear(); // previous parameter value has not been used
+ cmds.back() = theCommand;
+ }
+ else
+ {
+ cmds.push_back( theCommand );
+ }
+}
+
+//================================================================================
+/*!
+ * \brief Return true if a setting parameter command ha been used to compute mesh
+ */
+//================================================================================
+
+bool _pyHypothesis::isCmdUsedForCompute( const Handle(_pyCommand) & cmd,
+ _pyCommand::TAddr avoidComputeAddr ) const
+{
+ bool isUsed = false;
+ map< _pyCommand::TAddr, list<Handle(_pyCommand) > >::const_iterator addr2cmds =
+ myComputeAddr2Cmds.begin();
+ for ( ; addr2cmds != myComputeAddr2Cmds.end() && !isUsed; ++addr2cmds )
+ {
+ if ( addr2cmds->first == avoidComputeAddr ) continue;
+ const list<Handle(_pyCommand)> & cmds = addr2cmds->second;
+ isUsed = ( std::find( cmds.begin(), cmds.end(), cmd ) != cmds.end() );
+ }
+ return isUsed;
+}
+
+//================================================================================
+/*!
+ * \brief Save commands setting parameters as they are used for a mesh computation
+ */
+//================================================================================
+
+void _pyHypothesis::MeshComputed( const Handle(_pyCommand)& theComputeCmd )
+{
+ myComputeCmds.push_back( theComputeCmd );
+ list<Handle(_pyCommand)>& savedCmds = myComputeAddr2Cmds[ theComputeCmd->GetAddress() ];
+
+ map<TCollection_AsciiString, list< Handle(_pyCommand) > >::iterator m2c;
+ for ( m2c = myMeth2Commands.begin(); m2c != myMeth2Commands.end(); ++m2c )
+ savedCmds.push_back( m2c->second.back() );
+}
+
+//================================================================================
+/*!
+ * \brief Clear commands setting parameters as a mesh computed using them is cleared
+ */
+//================================================================================
+
+void _pyHypothesis::ComputeDiscarded( const Handle(_pyCommand)& theComputeCmd )
+{
+ list<Handle(_pyCommand)>& savedCmds = myComputeAddr2Cmds[ theComputeCmd->GetAddress() ];
+
+ list<Handle(_pyCommand)>::iterator cmd = savedCmds.begin();
+ for ( ; cmd != savedCmds.end(); ++cmd )
+ {
+ // check if a cmd has been used to compute another mesh
+ if ( isCmdUsedForCompute( *cmd, theComputeCmd->GetAddress() ))
+ continue;
+ // check if a cmd is a sole command setting its parameter;
+ // don't use method name for search as it can change
+ map<TCollection_AsciiString, list<Handle(_pyCommand)> >::iterator
+ m2cmds = myMeth2Commands.begin();
+ for ( ; m2cmds != myMeth2Commands.end(); ++m2cmds )
+ {
+ list< Handle(_pyCommand)>& cmds = m2cmds->second;
+ list< Handle(_pyCommand)>::iterator cmdIt = std::find( cmds.begin(), cmds.end(), *cmd );
+ if ( cmdIt != cmds.end() )
+ {
+ if ( cmds.back() != *cmd )
+ {
+ cmds.erase( cmdIt );
+ (*cmd)->Clear();
+ }
+ break;
+ }
+ }
+ }
+ myComputeAddr2Cmds.erase( theComputeCmd->GetAddress() );
+}
+
+//================================================================================
+/*!
+ * \brief Sets an argNb-th argument of current creation command
+ * \param argNb - argument index countered from 1
+ */
+//================================================================================
+
+void _pyHypothesis::setCreationArg( const int argNb, const _AString& arg )
+{
+ if ( myCurCrMethod )
+ {
+ while ( myCurCrMethod->myArgs.size() < argNb )
+ myCurCrMethod->myArgs.push_back( "None" );
+ if ( arg.IsEmpty() )
+ myCurCrMethod->myArgs[ argNb-1 ] = "None";
+ else
+ myCurCrMethod->myArgs[ argNb-1 ] = arg;
+ }
+}
+
+
+//================================================================================
+/*!
+ * \brief Remember hypothesis parameter values
+ * \param theCommand - The called hypothesis method
+ */
+//================================================================================
+
+void _pyComplexParamHypo::Process( const Handle(_pyCommand)& theCommand)
+{
+ if ( GetAlgoType() == "Cartesian_3D" )
+ {
+ // CartesianParameters3D hyp
+
+ if ( theCommand->GetMethod() == "SetSizeThreshold" )
+ {
+ setCreationArg( 4, theCommand->GetArg( 1 ));
+ myArgCommands.push_back( theCommand );
+ return;
+ }
+ if ( theCommand->GetMethod() == "SetGrid" ||
+ theCommand->GetMethod() == "SetGridSpacing" )
+ {
+ TCollection_AsciiString axis = theCommand->GetArg( theCommand->GetNbArgs() );
+ int iArg = axis.Value(1) - '0';
+ if ( theCommand->GetMethod() == "SetGrid" )
+ {
+ setCreationArg( 1+iArg, theCommand->GetArg( 1 ));
+ }
+ else
+ {
+ myCurCrMethod->myArgs[ iArg ] = "[ ";
+ myCurCrMethod->myArgs[ iArg ] += theCommand->GetArg( 1 );
+ myCurCrMethod->myArgs[ iArg ] += ", ";
+ myCurCrMethod->myArgs[ iArg ] += theCommand->GetArg( 2 );
+ myCurCrMethod->myArgs[ iArg ] += "]";
+ }
+ myArgCommands.push_back( theCommand );
+ rememberCmdOfParameter( theCommand );
+ return;
+ }
+ }
+
+ if( theCommand->GetMethod() == "SetLength" )
+ {
+ // NOW it is OBSOLETE
+ // ex: hyp.SetLength(start, 1)
+ // hyp.SetLength(end, 0)
+ ASSERT(( theCommand->GetArg( 2 ).IsIntegerValue() ));
+ int i = 1 - theCommand->GetArg( 2 ).IntegerValue();
+ TType2CrMethod::iterator type2meth = myAlgoType2CreationMethod.begin();
+ for ( ; type2meth != myAlgoType2CreationMethod.end(); ++type2meth )
+ {
+ CreationMethod& crMethod = type2meth->second;
+ while ( crMethod.myArgs.size() < i+1 )
+ crMethod.myArgs.push_back( "[]" );
+ crMethod.myArgs[ i ] = theCommand->GetArg( 1 ); // arg value
+ }
+ myArgCommands.push_back( theCommand );
+ }
+ else
+ {
+ _pyHypothesis::Process( theCommand );
+ }
+}
+//================================================================================
+/*!
+ * \brief Clear SetObjectEntry() as it is called by methods of Mesh_Segment
+ */
+//================================================================================
+
+void _pyComplexParamHypo::Flush()
+{
+ if ( IsWrapped() )
+ {
+ list < Handle(_pyCommand) >::iterator cmd = myUnusedCommands.begin();
+ for ( ; cmd != myUnusedCommands.end(); ++cmd )
+ if ((*cmd)->GetMethod() == "SetObjectEntry" )
+ (*cmd)->Clear();
+ }
+}
+
+//================================================================================
+/*!
+ * \brief Convert methods of 1D hypotheses to my own methods
+ * \param theCommand - The called hypothesis method
+ */
+//================================================================================
+
+void _pyLayerDistributionHypo::Process( const Handle(_pyCommand)& theCommand)
+{
+ if ( theCommand->GetMethod() != "SetLayerDistribution" )
+ return;
+
+ const _pyID& hyp1dID = theCommand->GetArg( 1 );
+ // Handle(_pyHypothesis) hyp1d = theGen->FindHyp( hyp1dID );
+ // if ( hyp1d.IsNull() && ! my1dHyp.IsNull()) // apparently hypId changed at study restoration
+ // {
+ // TCollection_AsciiString cmd =
+ // my1dHyp->GetCreationCmd()->GetIndentation() + hyp1dID + " = " + my1dHyp->GetID();
+ // Handle(_pyCommand) newCmd = theGen->AddCommand( cmd );
+ // theGen->SetCommandAfter( newCmd, my1dHyp->GetCreationCmd() );
+ // hyp1d = my1dHyp;
+ // }
+ // else if ( !my1dHyp.IsNull() && hyp1dID != my1dHyp->GetID() )
+ // {
+ // // 1D hypo is already set, so distribution changes and the old
+ // // 1D hypo is thrown away
+ // my1dHyp->ClearAllCommands();
+ // }
+ // my1dHyp = hyp1d;
+ // //my1dHyp->SetRemovedFromStudy( false );
+
+ // if ( !myArgCommands.empty() )
+ // myArgCommands.back()->Clear();
+ myCurCrMethod->myArgs.push_back( hyp1dID );
+ myArgCommands.push_back( theCommand );
+}
+
+//================================================================================
+/*!
+ * \brief
+ * \param theAdditionCmd - command to be converted
+ * \param theMesh - mesh instance
+ * \retval bool - status
+ */
+//================================================================================
+
+bool _pyLayerDistributionHypo::Addition2Creation( const Handle(_pyCommand)& theAdditionCmd,
+ const _pyID& theMesh)
+{
+ myIsWrapped = false;
+
+ if ( my1dHyp.IsNull() )
+ return false;
+
+ // set "SetLayerDistribution()" after addition cmd
+ theAdditionCmd->AddDependantCmd( myArgCommands.front() );
+
+ _pyID geom = theAdditionCmd->GetArg( 1 );
+
+ Handle(_pyHypothesis) algo = theGen->FindAlgo( geom, theMesh, this );
+ if ( !algo.IsNull() )
+ {
+ my1dHyp->SetMesh( theMesh );
+ my1dHyp->SetConvMethodAndType(my1dHyp->GetAlgoCreationMethod().ToCString(),
+ algo->GetAlgoType().ToCString());
+ if ( !my1dHyp->Addition2Creation( theAdditionCmd, theMesh ))
+ return false;
+
+ // clear "SetLayerDistribution()" cmd
+ myArgCommands.back()->Clear();
+
+ // Convert my creation => me = RadialPrismAlgo.Get3DHypothesis()
+
+ // find RadialPrism algo created on <geom> for theMesh
+ GetCreationCmd()->SetObject( algo->GetID() );
+ GetCreationCmd()->SetMethod( myAlgoMethod );
+ GetCreationCmd()->RemoveArgs();
+ theAdditionCmd->AddDependantCmd( GetCreationCmd() );
+ myIsWrapped = true;
+ }
+ return myIsWrapped;
+}
+
+//================================================================================
+/*!
+ * \brief
+ */
+//================================================================================
+
+void _pyLayerDistributionHypo::Flush()
+{
+ // as creation of 1D hyp was written later then it's edition,
+ // we need to find all it's edition calls and process them
+ list< Handle(_pyCommand) >::iterator cmd = myArgCommands.begin();
+ _pyID prevNewName;
+ for ( cmd = myArgCommands.begin(); cmd != myArgCommands.end(); ++cmd )
+ {
+ const _pyID& hyp1dID = (*cmd)->GetArg( 1 );
+ if ( hyp1dID.IsEmpty() ) continue;
+
+ Handle(_pyHypothesis) hyp1d = theGen->FindHyp( hyp1dID );
+
+ // make a new name for 1D hyp = "HypType" + "_Distribution"
+ _pyID newName;
+ if ( hyp1d.IsNull() ) // apparently hypId changed at study restoration
+ {
+ if ( prevNewName.IsEmpty() ) continue;
+ newName = prevNewName;
+ }
+ else
+ {
+ if ( hyp1d->IsWrapped() ) {
+ newName = hyp1d->GetCreationCmd()->GetMethod();
+ }
+ else {
+ TCollection_AsciiString hypTypeQuoted = hyp1d->GetCreationCmd()->GetArg(1);
+ newName = hypTypeQuoted.SubString( 2, hypTypeQuoted.Length() - 1 );
+ }
+ newName += "_Distribution";
+ prevNewName = newName;
+
+ hyp1d->GetCreationCmd()->SetResultValue( newName );
+ }
+ list< Handle(_pyCommand) >& cmds = theGen->GetCommands();
+ list< Handle(_pyCommand) >::iterator cmdIt = cmds.begin();
+ for ( ; cmdIt != cmds.end(); ++cmdIt ) {
+ const _pyID& objID = (*cmdIt)->GetObject();
+ if ( objID == hyp1dID ) {
+ if ( !hyp1d.IsNull() )
+ {
+ hyp1d->Process( *cmdIt );
+ hyp1d->GetCreationCmd()->AddDependantCmd( *cmdIt );
+ }
+ ( *cmdIt )->SetObject( newName );
+ }
+ }
+ // Set new hyp name to SetLayerDistribution(hyp1dID) cmd
+ (*cmd)->SetArg( 1, newName );
+ }
+}
+
+//================================================================================
+/*!
+ * \brief additionally to Addition2Creation, clears SetDistrType() command
+ * \param theCmd - AddHypothesis() command
+ * \param theMesh - mesh to which a hypothesis is added
+ * \retval bool - convertion result
+ */
+//================================================================================
+
+bool _pyNumberOfSegmentsHyp::Addition2Creation( const Handle(_pyCommand)& theCmd,
+ const _pyID& theMesh)
+{
+ if ( IsWrappable( theMesh ) && myCurCrMethod->myArgs.size() > 1 ) {
+ // scale factor (2-nd arg) is provided: clear SetDistrType(1) command
+ bool scaleDistrType = false;
+ list<Handle(_pyCommand)>::reverse_iterator cmd = myUnusedCommands.rbegin();
+ for ( ; cmd != myUnusedCommands.rend(); ++cmd ) {
+ if ( (*cmd)->GetMethod() == "SetDistrType" ) {
+ if ( (*cmd)->GetArg( 1 ) == "1" ) {
+ scaleDistrType = true;
+ (*cmd)->Clear();
+ }
+ else if ( !scaleDistrType ) {
+ // distribution type changed: remove scale factor from args
+ TType2CrMethod::iterator type2meth = myAlgoType2CreationMethod.begin();
+ for ( ; type2meth != myAlgoType2CreationMethod.end(); ++type2meth )
+ {
+ CreationMethod& crMethod = type2meth->second;
+ if ( crMethod.myArgs.size() == 2 )
+ crMethod.myArgs.pop_back();
+ }
+ break;
+ }
+ }
+ }
+ }
+ return _pyHypothesis::Addition2Creation( theCmd, theMesh );
+}
+
+//================================================================================
+/*!
+ * \brief remove repeated commands defining distribution
+ */
+//================================================================================
+
+void _pyNumberOfSegmentsHyp::Flush()
+{
+ // find number of the last SetDistrType() command
+ list<Handle(_pyCommand)>::reverse_iterator cmd = myUnusedCommands.rbegin();
+ int distrTypeNb = 0;
+ for ( ; !distrTypeNb && cmd != myUnusedCommands.rend(); ++cmd )
+ if ( (*cmd)->GetMethod() == "SetDistrType" ) {
+ if ( cmd != myUnusedCommands.rbegin() )
+ distrTypeNb = (*cmd)->GetOrderNb();
+ }
+ else if (IsWrapped() && (*cmd)->GetMethod() == "SetObjectEntry" ) {
+ (*cmd)->Clear();
+ }
+ // clear commands before the last SetDistrType()
+ list<Handle(_pyCommand)> * cmds[2] = { &myArgCommands, &myUnusedCommands };
+ set< int > treatedCmdNbs; // avoid treating same cmd twice
+ for ( int i = 0; i < 2; ++i ) {
+ set<TCollection_AsciiString> uniqueMethods;
+ list<Handle(_pyCommand)> & cmdList = *cmds[i];
+ for ( cmd = cmdList.rbegin(); cmd != cmdList.rend(); ++cmd )
+ {
+ if ( !treatedCmdNbs.insert( (*cmd)->GetOrderNb() ).second )
+ continue;// avoid treating same cmd twice
+ bool clear = ( (*cmd)->GetOrderNb() < distrTypeNb );
+ const TCollection_AsciiString& method = (*cmd)->GetMethod();
+ if ( !clear || method == "SetNumberOfSegments" ) {
+ bool isNewInSet = uniqueMethods.insert( method ).second;
+ clear = !isNewInSet;
+ }
+ if ( clear )
+ (*cmd)->Clear();
+ }
+ cmdList.clear();
+ }
+}
+
+//================================================================================
+/*!
+ * \brief Convert the command adding "SegmentLengthAroundVertex" to mesh
+ * into regular1D.LengthNearVertex( length, vertex )
+ * \param theCmd - The command like mesh.AddHypothesis( vertex, SegmentLengthAroundVertex )
+ * \param theMesh - The mesh needing this hypo
+ * \retval bool - false if the command cant be converted
+ */
+//================================================================================
+
+bool _pySegmentLengthAroundVertexHyp::Addition2Creation( const Handle(_pyCommand)& theCmd,
+ const _pyID& theMeshID)
+{
+ if ( IsWrappable( theMeshID )) {
+
+ _pyID vertex = theCmd->GetArg( 1 );
+
+ // the problem here is that segment algo will not be found
+ // by pyHypothesis::Addition2Creation() for <vertex>, so we try to find
+ // geometry where segment algorithm is assigned
+ Handle(_pyHypothesis) algo;
+ _pyID geom = vertex;
+ while ( algo.IsNull() && !geom.IsEmpty()) {
+ // try to find geom as a father of <vertex>
+ geom = FatherID( geom );
+ algo = theGen->FindAlgo( geom, theMeshID, this );
+ }
+ if ( algo.IsNull() )
+ return false; // also possible to find geom as brother of veretex...
+ // set geom instead of vertex
+ theCmd->SetArg( 1, geom );
+
+ // set vertex as a second arg
+ if ( myCurCrMethod->myArgs.size() < 1) setCreationArg( 1, "1" ); // :(
+ setCreationArg( 2, vertex );
+
+ // mesh.AddHypothesis(vertex, SegmentLengthAroundVertex) -->
+ // theMeshID.LengthNearVertex( length, vertex )
+ return _pyHypothesis::Addition2Creation( theCmd, theMeshID );
+ }
+ return false;
+}
+
+//================================================================================
+/*!
+ * \brief _pyAlgorithm constructor
+ * \param theCreationCmd - The command like "algo = smeshgen.CreateHypothesis(type,lib)"
+ */
+//================================================================================
+
+_pyAlgorithm::_pyAlgorithm(const Handle(_pyCommand)& theCreationCmd)
+ : _pyHypothesis( theCreationCmd )
+{
+ myIsAlgo = true;
+}
+
+//================================================================================
+/*!
+ * \brief Convert the command adding an algorithm to mesh
+ * \param theCmd - The command like mesh.AddHypothesis( geom, algo )
+ * \param theMesh - The mesh needing this algo
+ * \retval bool - false if the command cant be converted
+ */
+//================================================================================
+
+bool _pyAlgorithm::Addition2Creation( const Handle(_pyCommand)& theCmd,
+ const _pyID& theMeshID)
+{
+ // mesh.AddHypothesis(geom,algo) --> theMeshID.myCreationMethod()
+ if ( _pyHypothesis::Addition2Creation( theCmd, theMeshID )) {
+ theGen->SetAccessorMethod( GetID(), "GetAlgorithm()" );
+ return true;
+ }
+ return false;
+}
+
+//================================================================================
+/*!
+ * \brief Return starting position of a part of python command
+ * \param thePartIndex - The index of command part
+ * \retval int - Part position
+ */
+//================================================================================
+
+int _pyCommand::GetBegPos( int thePartIndex )
+{
+ if ( IsEmpty() )
+ return EMPTY;
+ if ( myBegPos.Length() < thePartIndex )
+ return UNKNOWN;
+ return myBegPos( thePartIndex );
+}
+
+//================================================================================
+/*!
+ * \brief Store starting position of a part of python command
+ * \param thePartIndex - The index of command part
+ * \param thePosition - Part position
+ */
+//================================================================================
+
+void _pyCommand::SetBegPos( int thePartIndex, int thePosition )
+{
+ while ( myBegPos.Length() < thePartIndex )
+ myBegPos.Append( UNKNOWN );
+ myBegPos( thePartIndex ) = thePosition;
+}
+
+//================================================================================
+/*!
+ * \brief Returns whitespace symbols at the line beginning
+ * \retval TCollection_AsciiString - result
+ */
+//================================================================================
+
+TCollection_AsciiString _pyCommand::GetIndentation()
+{
+ int end = 1;
+ if ( GetBegPos( RESULT_IND ) == UNKNOWN )
+ GetWord( myString, end, true );
+ else
+ end = GetBegPos( RESULT_IND );
+ return myString.SubString( 1, end - 1 );
+}
+
+//================================================================================
+/*!
+ * \brief Return substring of python command looking like ResultValue = Obj.Meth()
+ * \retval const TCollection_AsciiString & - ResultValue substring
+ */
+//================================================================================
+
+const TCollection_AsciiString & _pyCommand::GetResultValue()
+{
+ if ( GetBegPos( RESULT_IND ) == UNKNOWN )
+ {
+ SetBegPos( RESULT_IND, EMPTY );
+ int begPos, endPos = myString.Location( "=", 1, Length() );
+ if ( endPos )
+ {
+ begPos = 1;
+ while ( begPos < endPos && isspace( myString.Value( begPos ))) ++begPos;
+ if ( begPos < endPos )
+ {
+ SetBegPos( RESULT_IND, begPos );
+ --endPos;
+ while ( begPos < endPos && isspace( myString.Value( endPos ))) --endPos;
+ myRes = myString.SubString( begPos, endPos );
+ }
+ }
+ }
+ return myRes;
+}
+
+//================================================================================
+/*!
+ * \brief Return number of python command result value ResultValue = Obj.Meth()
+ */
+//================================================================================
+
+int _pyCommand::GetNbResultValues()
+{
+ int nb = 0;
+ int begPos = 1;
+ int endPos = myString.Location( "=", 1, Length() );
+ while ( begPos < endPos )
+ {
+ _AString str = GetWord( myString, begPos, true );
+ begPos = begPos+ str.Length();
+ nb++;
+ }
+ return (nb-1);
+}
+
+
+//================================================================================
+/*!
+ * \brief Return substring of python command looking like
+ * ResultValue1 , ResultValue2,... = Obj.Meth() with res index
+ * \retval const TCollection_AsciiString & - ResultValue with res index substring
+ */
+//================================================================================
+TCollection_AsciiString _pyCommand::GetResultValue(int res)
+{
+ int begPos = 1;
+ if ( SkipSpaces( myString, begPos ) && myString.Value( begPos ) == '[' )
+ ++begPos; // skip [, else the whole list is returned
+ int endPos = myString.Location( "=", 1, Length() );
+ int Nb=0;
+ while ( begPos < endPos) {
+ _AString result = GetWord( myString, begPos, true );
+ begPos = begPos + result.Length();
+ Nb++;
+ if(res == Nb) {
+ result.RemoveAll('[');
+ result.RemoveAll(']');
+ return result;
+ }
+ if(Nb>res)
+ break;
+ }
+ return theEmptyString;
+}
+
+//================================================================================
+/*!
+ * \brief Return substring of python command looking like ResVal = Object.Meth()
+ * \retval const TCollection_AsciiString & - Object substring
+ */
+//================================================================================
+
+const TCollection_AsciiString & _pyCommand::GetObject()
+{
+ if ( GetBegPos( OBJECT_IND ) == UNKNOWN )
+ {
+ // beginning
+ int begPos = GetBegPos( RESULT_IND ) + myRes.Length();
+ if ( begPos < 1 ) {
+ begPos = myString.Location( "=", 1, Length() ) + 1;
+ // is '=' in the string argument (for example, name) or not
+ int nb1 = 0; // number of ' character at the left of =
+ int nb2 = 0; // number of " character at the left of =
+ for ( int i = 1; i < begPos-1; i++ ) {
+ if ( myString.Value( i )=='\'' )
+ nb1 += 1;
+ else if ( myString.Value( i )=='"' )
+ nb2 += 1;
+ }
+ // if number of ' or " is not divisible by 2,
+ // then get an object at the start of the command
+ if ( nb1 % 2 != 0 || nb2 % 2 != 0 )
+ begPos = 1;
+ }
+ myObj = GetWord( myString, begPos, true );
+ // check if object is complex,
+ // so far consider case like "smesh.Method()"
+ if ( int bracketPos = myString.Location( "(", begPos, Length() )) {
+ //if ( bracketPos==0 ) bracketPos = Length();
+ int dotPos = begPos+myObj.Length();
+ while ( dotPos+1 < bracketPos ) {
+ if ( int pos = myString.Location( ".", dotPos+1, bracketPos ))
+ dotPos = pos;
+ else
+ break;
+ }
+ if ( dotPos > begPos+myObj.Length() )
+ myObj = myString.SubString( begPos, dotPos-1 );
+ }
+ // 1st word after '=' is an object
+ // else // no method -> no object
+ // {
+ // myObj.Clear();
+ // begPos = EMPTY;
+ // }
+ // store
+ SetBegPos( OBJECT_IND, begPos );
+ }
+ //SCRUTE(myObj);
+ return myObj;
+}
+
+//================================================================================
+/*!
+ * \brief Return substring of python command looking like ResVal = Obj.Method()
+ * \retval const TCollection_AsciiString & - Method substring
+ */
+//================================================================================
+
+const TCollection_AsciiString & _pyCommand::GetMethod()
+{
+ if ( GetBegPos( METHOD_IND ) == UNKNOWN )
+ {
+ // beginning
+ int begPos = GetBegPos( OBJECT_IND ) + myObj.Length();
+ bool forward = true;
+ if ( begPos < 1 ) {
+ begPos = myString.Location( "(", 1, Length() ) - 1;
+ forward = false;
+ }
+ // store
+ myMeth = GetWord( myString, begPos, forward );
+ SetBegPos( METHOD_IND, begPos );
+ }
+ //SCRUTE(myMeth);
+ return myMeth;
+}
+
+//================================================================================
+/*!
+ * \brief Return substring of python command looking like ResVal = Obj.Meth(Arg1,...)
+ * \retval const TCollection_AsciiString & - Arg<index> substring
+ */
+//================================================================================
+
+const TCollection_AsciiString & _pyCommand::GetArg( int index )
+{
+ if ( GetBegPos( ARG1_IND ) == UNKNOWN )
+ {
+ // Find all args
+
+ int pos = GetBegPos( METHOD_IND ) + myMeth.Length();
+ if ( pos < 1 )
+ pos = myString.Location( "(", 1, Length() );
+ else
+ --pos;
+
+ // we are at or before '(', skip it if present
+ if ( pos > 0 ) {
+ while ( pos <= Length() && myString.Value( pos ) != '(' ) ++pos;
+ if ( pos > Length() )
+ pos = 0;
+ }
+ if ( pos < 1 ) {
+ SetBegPos( ARG1_IND, 0 ); // even no '('
+ return theEmptyString;
+ }
+ ++pos;
+
+ list< TCollection_AsciiString > separatorStack( 1, ",)");
+ bool ignoreNesting = false;
+ int prevPos = pos;
+ while ( pos <= Length() )
+ {
+ const char chr = myString.Value( pos );
+
+ if ( separatorStack.back().Location( chr, 1, separatorStack.back().Length()))
+ {
+ if ( separatorStack.size() == 1 ) // a comma dividing args or a terminal ')' found
+ {
+ while ( pos-1 >= prevPos && isspace( myString.Value( prevPos )))
+ ++prevPos;
+ TCollection_AsciiString arg;
+ if ( pos-1 >= prevPos ) {
+ arg = myString.SubString( prevPos, pos-1 );
+ arg.RightAdjust(); // remove spaces
+ arg.LeftAdjust();
+ }
+ if ( !arg.IsEmpty() || chr == ',' )
+ {
+ SetBegPos( ARG1_IND + myArgs.Length(), prevPos );
+ myArgs.Append( arg );
+ }
+ if ( chr == ')' )
+ break;
+ prevPos = pos+1;
+ }
+ else // end of nesting args found
+ {
+ separatorStack.pop_back();
+ ignoreNesting = false;
+ }
+ }
+ else if ( !ignoreNesting )
+ {
+ switch ( chr ) {
+ case '(' : separatorStack.push_back(")"); break;
+ case '[' : separatorStack.push_back("]"); break;
+ case '\'': separatorStack.push_back("'"); ignoreNesting=true; break;
+ case '"' : separatorStack.push_back("\""); ignoreNesting=true; break;
+ default:;
+ }
+ }
+ ++pos;
+ }
+ }
+ if ( myArgs.Length() < index )
+ return theEmptyString;
+ return myArgs( index );
+}
+
+//================================================================================
+/*!
+ * \brief Check if char is a word part
+ * \param c - The character to check
+ * \retval bool - The check result
+ */
+//================================================================================
+
+static inline bool isWord(const char c, const bool dotIsWord)
+{
+ return
+ !isspace(c) && c != ',' && c != '=' && c != ')' && c != '(' && ( dotIsWord || c != '.');
+}
+
+//================================================================================
+/*!
+ * \brief Looks for a word in the string and returns word's beginning
+ * \param theString - The input string
+ * \param theStartPos - The position to start the search, returning word's beginning
+ * \param theForward - The search direction
+ * \retval TCollection_AsciiString - The found word
+ */
+//================================================================================
+
+TCollection_AsciiString _pyCommand::GetWord( const _AString & theString,
+ int & theStartPos,
+ const bool theForward,
+ const bool dotIsWord )
+{
+ int beg = theStartPos, end = theStartPos;
+ theStartPos = EMPTY;
+ if ( beg < 1 || beg > theString.Length() )
+ return theEmptyString;
+
+ if ( theForward ) { // search forward
+ // beg
+ while ( beg <= theString.Length() && !isWord( theString.Value( beg ), dotIsWord))
+ ++beg;
+ if ( beg > theString.Length() )
+ return theEmptyString; // no word found
+ // end
+ end = beg + 1;
+ char begChar = theString.Value( beg );
+ if ( begChar == '"' || begChar == '\'' || begChar == '[') {
+ char endChar = ( begChar == '[' ) ? ']' : begChar;
+ // end is at the corresponding quoting mark or bracket
+ while ( end < theString.Length() &&
+ ( theString.Value( end ) != endChar || theString.Value( end-1 ) == '\\'))
+ ++end;
+ }
+ else {
+ while ( end <= theString.Length() && isWord( theString.Value( end ), dotIsWord))
+ ++end;
+ --end;
+ }
+ }
+ else { // search backward
+ // end
+ while ( end > 0 && !isWord( theString.Value( end ), dotIsWord))
+ --end;
+ if ( end == 0 )
+ return theEmptyString; // no word found
+ beg = end - 1;
+ char endChar = theString.Value( end );
+ if ( endChar == '"' || endChar == '\'' || endChar == ']') {
+ char begChar = ( endChar == ']' ) ? '[' : endChar;
+ // beg is at the corresponding quoting mark
+ while ( beg > 1 &&
+ ( theString.Value( beg ) != begChar || theString.Value( beg-1 ) == '\\'))
+ --beg;
+ }
+ else {
+ while ( beg > 0 && isWord( theString.Value( beg ), dotIsWord))
+ --beg;
+ ++beg;
+ }
+ }
+ theStartPos = beg;
+ //cout << theString << " ---- " << beg << " - " << end << endl;
+ return theString.SubString( beg, end );
+}
+
+//================================================================================
+/*!
+ * \brief Returns true if the string looks like a study entry
+ */
+//================================================================================
+
+bool _pyCommand::IsStudyEntry( const TCollection_AsciiString& str )