Salome HOME
877874c338b3c8bb2b42ee2fb553570a7991d9af
[modules/smesh.git] / src / SMESH_I / SMESH_NoteBook.cxx
1 // Copyright (C) 2008  OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
2 // CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
3 //
4 // This library is free software; you can redistribute it and/or
5 // modify it under the terms of the GNU Lesser General Public
6 // License as published by the Free Software Foundation; either
7 // version 2.1 of the License.
8 //
9 // This library is distributed in the hope that it will be useful
10 // but WITHOUT ANY WARRANTY; without even the implied warranty of
11 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
12 // Lesser General Public License for more details.
13 //
14 // You should have received a copy of the GNU Lesser General Public
15 // License along with this library; if not, write to the Free Software
16 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
17 //
18 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
19 //
20 // File      : SMESH_NoteBook.cxx
21 // Author    : Roman NIKOLAEV
22
23 #include "SMESH_2smeshpy.hxx"
24 #include "SMESH_NoteBook.hxx"
25 #include "SMESH_Gen_i.hxx"
26 #include "SMESH_PythonDump.hxx"
27
28 #include <Resource_DataMapOfAsciiStringAsciiString.hxx>
29 #include <TColStd_SequenceOfAsciiString.hxx>
30 #include <TColStd_HSequenceOfInteger.hxx>
31
32 #include <vector>
33 #include <string>
34
35 #ifdef _DEBUG_
36 static int MYDEBUG = 0;
37 #else
38 static int MYDEBUG = 0;
39 #endif
40
41 using namespace std;
42
43
44
45 //================================================================================
46 /*!
47  * \brief Constructor
48  */
49 //================================================================================
50 ObjectStates::ObjectStates(TCollection_AsciiString theType)
51 {
52   _type = theType;
53   _dumpstate = 0;
54 }
55
56 //================================================================================
57 /*!
58  * \brief Destructor
59  */
60 //================================================================================
61 ObjectStates::~ObjectStates()
62 {
63 }
64
65 //================================================================================
66 /*!
67  * \brief Add new object state 
68  * \param theState - Object state (vector of notebook variable)
69  */
70 //================================================================================
71 void ObjectStates::AddState(const TState &theState)
72 {
73   _states.push_back(theState);
74 }
75
76 //================================================================================
77 /*!
78  * \brief Return current object state
79  * \\retval state - Object state (vector of notebook variable)
80  */
81 //================================================================================
82 TState ObjectStates::GetCurrectState() const
83 {
84   if(_states.size() > _dumpstate)
85     return _states[_dumpstate];
86   TState empty;
87   return empty;
88 }
89
90
91 //================================================================================
92 /*!
93  *
94  */
95 //================================================================================
96 TAllStates ObjectStates::GetAllStates() const
97 {
98   return _states;
99 }
100
101 //================================================================================
102 /*!
103  *
104  */
105 //================================================================================
106 void ObjectStates::IncrementState()
107 {
108   _dumpstate++;
109 }
110
111 //================================================================================
112 /*!
113  *
114  */
115 //================================================================================
116 TCollection_AsciiString ObjectStates::GetObjectType() const{
117   return _type;
118 }
119
120
121 //================================================================================
122 /*!
123  * \brief Constructor
124  */
125 //================================================================================
126 LayerDistributionStates::LayerDistributionStates():
127   ObjectStates("LayerDistribution")
128 {
129 }
130 //================================================================================
131 /*!
132  * \brief Destructor
133  */
134 //================================================================================
135 LayerDistributionStates::~LayerDistributionStates()
136 {
137 }
138
139
140 //================================================================================
141 /*!
142  * \brief AddDistribution
143  */
144 //================================================================================
145 void LayerDistributionStates::AddDistribution(const TCollection_AsciiString& theDistribution)
146 {
147   _distributions.insert(pair<TCollection_AsciiString,TCollection_AsciiString>(theDistribution,""));
148 }
149
150 //================================================================================
151 /*!
152  * \brief HasDistribution
153  */
154 //================================================================================
155 bool LayerDistributionStates::HasDistribution(const TCollection_AsciiString& theDistribution) const
156 {
157   return _distributions.find(theDistribution) != _distributions.end();
158 }
159
160 //================================================================================
161 /*!
162  * \brief SetDistributionType
163  */
164 //================================================================================
165 bool LayerDistributionStates::SetDistributionType(const TCollection_AsciiString& theDistribution,
166                                                   const TCollection_AsciiString& theType)
167 {
168   TDistributionMap::iterator it = _distributions.find(theDistribution);
169   if(it == _distributions.end())
170     return false;
171   (*it).second = theType;
172   return true;
173 }
174
175 //================================================================================
176 /*!
177  * \brief GetDistributionType
178  */
179 //================================================================================
180 TCollection_AsciiString LayerDistributionStates::
181 GetDistributionType(const TCollection_AsciiString& theDistribution) const
182 {
183   TDistributionMap::const_iterator it = _distributions.find(theDistribution);
184   return (it == _distributions.end()) ? TCollection_AsciiString() : (*it).second;
185 }
186
187 //================================================================================
188 /*!
189  * \brief Constructor
190  */
191 //================================================================================
192 SMESH_NoteBook::SMESH_NoteBook()
193 {
194   InitObjectMap();
195 }
196
197 //================================================================================
198 /*!
199  * \brief Destructor
200  */
201 //================================================================================
202 SMESH_NoteBook::~SMESH_NoteBook()
203 {
204   TVariablesMap::const_iterator it = _objectMap.begin();
205   for(;it!=_objectMap.end();it++) {
206     if((*it).second)
207       delete (*it).second;
208   }
209 }
210
211 //================================================================================
212 /*!
213  * \brief Replace parameters of the functions on the Salome NoteBook Variables
214  * \param theString - Input string
215  * \retval TCollection_AsciiString - Convertion result
216  */
217 //================================================================================
218 void SMESH_NoteBook::ReplaceVariables()
219 {
220   for(int i=0;i<_commands.size();i++) {
221     Handle(_pyCommand) aCmd = _commands[i];
222     TCollection_AsciiString aMethod = aCmd->GetMethod();
223     TCollection_AsciiString aObject = aCmd->GetObject();
224     TCollection_AsciiString aResultValue = aCmd->GetResultValue();
225     if(MYDEBUG) {
226       cout<<"Command before : "<< aCmd->GetString()<<endl;
227       cout<<"Method : "<< aMethod<<endl;
228       cout<<"Object : "<< aObject<<endl;
229       cout<<"Result : "<< aResultValue<<endl;
230     }
231     
232     // check if method modifies the object itself
233     TVariablesMap::const_iterator it = _objectMap.find(aObject);
234     if(it == _objectMap.end()) // check if method returns a new object
235       it = _objectMap.find(aResultValue);
236     
237     if(it == _objectMap.end()) { // check if method modifies a mesh using mesh editor
238       TMeshEditorMap::const_iterator meIt = myMeshEditors.find(aObject);
239       if(meIt != myMeshEditors.end()) {
240         TCollection_AsciiString aMesh = (*meIt).second;
241         it = _objectMap.find(aMesh);
242       }
243     }
244     
245     if(it != _objectMap.end()) {
246       ObjectStates *aStates = (*it).second;
247       // Case for LocalLength hypothesis
248       if(aStates->GetObjectType().IsEqual("LocalLength") && aStates->GetCurrectState().size() >= 2) {
249         if(aMethod.IsEqual("SetLength")) {
250           if(!aStates->GetCurrectState().at(0).IsEmpty() )
251             aCmd->SetArg(1,aStates->GetCurrectState().at(0));
252           aStates->IncrementState();
253         }
254         else if(aMethod.IsEqual("SetPrecision")) {
255           if(!aStates->GetCurrectState().at(1).IsEmpty() )
256             aCmd->SetArg(1,aStates->GetCurrectState().at(1));
257           aStates->IncrementState();
258         }
259       }
260       
261       // Case for SegmentLengthAroundVertex hypothesis
262       else if(aStates->GetObjectType().IsEqual("SegmentLengthAroundVertex")
263               && aStates->GetCurrectState().size() >= 1) {
264         if(aMethod == "SetLength") {
265           if(!aStates->GetCurrectState().at(0).IsEmpty() )
266             aCmd->SetArg(1,aStates->GetCurrectState().at(0));
267           aStates->IncrementState();
268         }
269       }
270
271       // Case for Arithmetic1D and StartEndLength hypothesis
272       else if(aStates->GetObjectType().IsEqual("Arithmetic1D") || 
273               aStates->GetObjectType().IsEqual("StartEndLength")) {
274         if(aMethod == "SetLength" &&
275            aStates->GetCurrectState().size() >= 2) {
276           if(aCmd->GetArg(2) == "1" && !aStates->GetCurrectState().at(0).IsEmpty())
277             aCmd->SetArg(1,aStates->GetCurrectState().at(0));
278           else if(!aStates->GetCurrectState().at(1).IsEmpty())
279             aCmd->SetArg(1,aStates->GetCurrectState().at(1));
280           aStates->IncrementState();
281         }
282       }
283       
284       //Case for Deflection1D hypothesis
285       else if(aStates->GetObjectType().IsEqual("Deflection1D")){
286         if(aMethod == "SetDeflection" && aStates->GetCurrectState().size() >= 1) {
287           if(!aStates->GetCurrectState().at(0).IsEmpty() )
288             aCmd->SetArg(1,aStates->GetCurrectState().at(0));
289           aStates->IncrementState();
290         }
291       }
292       
293       // Case for LayerDistribution hypothesis (not finished yet)
294       else if(aStates->GetObjectType() == "LayerDistribution") {
295         if(aMethod == "SetLayerDistribution"){
296           LayerDistributionStates* aLDStates = (LayerDistributionStates*)(aStates);
297           aLDStates->AddDistribution(aCmd->GetArg(1));
298         }
299       }
300       
301       // Case for MaxElementArea hypothesis
302       else if(aStates->GetObjectType().IsEqual("MaxElementArea")){
303         if(aMethod == "SetMaxElementArea" && aStates->GetCurrectState().size() >= 1) {
304           if(!aStates->GetCurrectState().at(0).IsEmpty() )
305             aCmd->SetArg(1,aStates->GetCurrectState().at(0));
306           aStates->IncrementState();
307         }
308       }
309
310       // Case for NumberOfLayers hypothesis
311       else if(aStates->GetObjectType().IsEqual("NumberOfLayers")){
312         if(aMethod == "SetNumberOfLayers" && aStates->GetCurrectState().size() >= 1) {
313           if(!aStates->GetCurrectState().at(0).IsEmpty() )
314             aCmd->SetArg(1,aStates->GetCurrectState().at(0));
315           aStates->IncrementState();
316         }
317       }
318
319       // Case for NumberOfSegments hypothesis
320       else if(aStates->GetObjectType().IsEqual("NumberOfSegments")){
321         if(aMethod == "SetNumberOfSegments" && aStates->GetCurrectState().size() >= 1) {
322           if(!aStates->GetCurrectState().at(0).IsEmpty() )
323             aCmd->SetArg(1,aStates->GetCurrectState().at(0));
324           if(aStates->GetCurrectState().size()==1)
325             aStates->IncrementState();
326         }
327         else if (aMethod == "SetScaleFactor" && aStates->GetCurrectState().size() >= 2) {
328           if(!aStates->GetCurrectState().at(1).IsEmpty() )
329             aCmd->SetArg(1,aStates->GetCurrectState().at(1));
330           aStates->IncrementState();
331         }
332       }
333       
334       else if(aStates->GetObjectType().IsEqual("Mesh")) {
335         TState aCurrentState = aStates->GetCurrectState();
336         int aCurrentStateSize = aCurrentState.size();
337         if(aMethod.IsEqual("Translate") ||
338            aMethod.IsEqual("TranslateMakeGroups") ||
339            aMethod.IsEqual("TranslateMakeMesh")) {
340           bool isVariableFound = false;
341           int anArgIndex = 0;
342           for(int i = 1, n = aCmd->GetNbArgs(); i <= n; i++) {
343             if(aCmd->GetArg(i).IsEqual("SMESH.PointStruct")) {
344               anArgIndex = i+1;
345               break;
346             }
347           }
348           if(anArgIndex > 0) {
349             if(aCurrentStateSize == 3) { // translation by dx, dy, dz
350               for(int j = 0; j < aCurrentStateSize; j++) {
351                 if(!aCurrentState.at(j).IsEmpty()) {
352                   isVariableFound = true;
353                   aCmd->SetArg(anArgIndex+j, aCurrentState.at(j));
354                 }
355               }
356             }
357             else if(aCurrentStateSize == 6) { // translation by x1, x2, y1, y2, z1, z2
358               isVariableFound = true;
359               for(int j = 0; j < 3; j++) {
360                 TCollection_AsciiString anArg = aCmd->GetArg(anArgIndex+j);
361                 TCollection_AsciiString aValue1 = aCurrentState.at(2*j), aValue2 = aCurrentState.at(2*j+1);
362                 bool aV1 = !aValue1.IsEmpty();
363                 bool aV2 = !aValue2.IsEmpty();
364                 double aValue, aCurrentValue = anArg.IsRealValue() ? anArg.RealValue() : 0;
365                 if(aV1 && !aV2) {
366                   if(!GetReal(aValue1, aValue))
367                     aValue = 0;
368                   aValue2 = TCollection_AsciiString( aValue + aCurrentValue );
369                 }
370                 else if(!aV1 && aV2) {
371                   if(!GetReal(aValue2, aValue))
372                     aValue = 0;
373                   aValue1 = TCollection_AsciiString( aValue - aCurrentValue );
374                 }
375                 else if(!aV1 && !aV2) {
376                   aValue1 = TCollection_AsciiString( 0 );
377                   aValue2 = TCollection_AsciiString( aCurrentValue );
378                 }
379                 aCmd->SetArg(anArgIndex+j, aValue1 + ", " + aValue2 );
380               }
381             }
382           }
383           if(isVariableFound) {
384             TCollection_AsciiString aDim;
385             if(aCurrentStateSize == 6)
386               aDim = "6";
387             aCmd->SetArg(anArgIndex - 1, TCollection_AsciiString(SMESH_2smeshpy::SmeshpyName())+".PointStructStr"+aDim);
388             aCmd->SetArg(anArgIndex - 2, TCollection_AsciiString(SMESH_2smeshpy::SmeshpyName())+".DirStructStr");
389           }
390           aStates->IncrementState();
391         }
392         else if(aMethod.IsEqual("Rotate") ||
393                 aMethod.IsEqual("RotateMakeGroups") ||
394                 aMethod.IsEqual("RotateMakeMesh") ||
395                 aMethod.IsEqual("RotationSweep") ||
396                 aMethod.IsEqual("RotationSweepMakeGroups") ||
397                 aMethod.IsEqual("Mirror") ||
398                 aMethod.IsEqual("MirrorMakeGroups")) {
399           bool isSubstitute = false;
400           int anArgIndex = 0;
401           for(int i = 1, n = aCmd->GetNbArgs(); i <= n; i++) {
402             if(aCmd->GetArg(i).IsEqual("SMESH.AxisStruct")) {
403               anArgIndex = i+1;
404               break;
405             }
406           }
407           if(anArgIndex > 0) {
408             for(int j = 0; j < aCurrentStateSize; j++) {
409               if(!aCurrentState.at(j).IsEmpty()) {
410                 if(j < 6) // 0-5 - axis struct, 6 - angle (rotation & sweep), 7-8 - nbSteps and tolerance (sweep)
411                   isSubstitute = true;
412                 aCmd->SetArg(anArgIndex+j, aCurrentState.at(j));
413               }
414             }
415           }
416           if(isSubstitute)
417             aCmd->SetArg(anArgIndex - 1, TCollection_AsciiString(SMESH_2smeshpy::SmeshpyName())+".AxisStructStr");
418           aStates->IncrementState();
419         }
420         else if(aMethod.IsEqual("AddNode") ||
421                 aMethod.IsEqual("MoveClosestNodeToPoint")) {
422           for(int j = 0; j < aCurrentStateSize; j++) {
423             if(!aCurrentState.at(j).IsEmpty())
424               aCmd->SetArg(j+1, aCurrentState.at(j));
425           }
426           aStates->IncrementState();
427         }
428         else if(aMethod.IsEqual("MoveNode")) {
429           for(int j = 0; j < aCurrentStateSize; j++) {
430             if(!aCurrentState.at(j).IsEmpty())
431               aCmd->SetArg(j+2, aCurrentState.at(j));
432           }
433           aStates->IncrementState();
434         }
435         else if(aMethod.IsEqual("ExtrusionSweep") ||
436                 aMethod.IsEqual("ExtrusionSweepMakeGroups")) {
437           bool isSubstitute = false;
438           int anArgIndex = 0;
439           for(int i = 1, n = aCmd->GetNbArgs(); i <= n; i++) {
440             if(aCmd->GetArg(i).IsEqual("SMESH.PointStruct")) {
441               anArgIndex = i+1;
442               break;
443             }
444           }
445           if(anArgIndex > 0) {
446             for(int j = 0; j < aCurrentStateSize; j++) {
447               if(!aCurrentState.at(j).IsEmpty()) {
448                 if(j < 3) // 0-2 - dir struct, 3 - number of steps
449                   isSubstitute = true;
450                 aCmd->SetArg(anArgIndex+j, aCurrentState.at(j));
451               }
452             }
453           }
454           if(isSubstitute) {
455             aCmd->SetArg(anArgIndex - 1, TCollection_AsciiString(SMESH_2smeshpy::SmeshpyName())+".PointStructStr");
456             aCmd->SetArg(anArgIndex - 2, TCollection_AsciiString(SMESH_2smeshpy::SmeshpyName())+".DirStructStr");
457           }
458           aStates->IncrementState();
459         }
460         else if(aMethod.IsEqual("ExtrusionAlongPath") ||
461                 aMethod.IsEqual("ExtrusionAlongPathMakeGroups") ||
462                 /* workaround for a bug in the command parsing algorithm */
463                 aCmd->GetString().Search("ExtrusionAlongPathMakeGroups") != -1) {
464           int aNbAngles = aCurrentStateSize-3; // State looks like "Angle1:...:AngleN:X:Y:Z"
465           bool isSubstitute = false;
466           int anArgIndex = 0;
467           for(int i = 1, n = aCmd->GetNbArgs(); i <= n; i++) {
468             if(aCmd->GetArg(i).IsEqual("SMESH.PointStruct")) {
469               anArgIndex = i-1-aNbAngles;
470               break;
471             }
472           }
473           if(anArgIndex > 0) {
474             int j = 0;
475             for(; j < aNbAngles; j++) {
476               if(!aCurrentState.at(j).IsEmpty()) {
477                 aCmd->SetArg(anArgIndex+j-1, aCurrentState.at(j));
478               }
479             }
480             for(; j < aNbAngles+3; j++) {
481               if(!aCurrentState.at(j).IsEmpty()) {
482                 isSubstitute = true;
483                 aCmd->SetArg(anArgIndex+j+2, aCurrentState.at(j));
484               }
485             }
486           }
487           if(isSubstitute)
488             aCmd->SetArg(anArgIndex + aNbAngles + 1,
489                          TCollection_AsciiString(SMESH_2smeshpy::SmeshpyName())+".PointStructStr");
490           aStates->IncrementState();
491         }
492         else if(aMethod.IsEqual("TriToQuad") ||
493                 aMethod.IsEqual("Concatenate") ||
494                 aMethod.IsEqual("ConcatenateWithGroups")) {
495           if(aCurrentStateSize && !aCurrentState.at(0).IsEmpty())
496             aCmd->SetArg(aCmd->GetNbArgs(), aCurrentState.at(0));
497           aStates->IncrementState();
498         }
499         else if(aMethod.IsEqual("Smooth") ||
500                 aMethod.IsEqual("SmoothParametric")) {
501           int anArgIndex = aCmd->GetNbArgs() - 2;
502           for(int j = 0; j < aCurrentStateSize; j++) {
503             if(!aCurrentState.at(j).IsEmpty())
504               aCmd->SetArg(anArgIndex+j, aCurrentState.at(j));
505           }
506           aStates->IncrementState();
507         }
508       }
509     }
510     if(MYDEBUG) {
511       cout<<"Command after: "<< aCmd->GetString()<<endl;
512     }
513   }
514   //  ProcessLayerDistribution();
515 }
516 //================================================================================
517 /*!
518  * \brief Private method
519  */
520 //================================================================================
521 void SMESH_NoteBook::InitObjectMap()
522 {
523   SMESH_Gen_i *aGen = SMESH_Gen_i::GetSMESHGen();
524   if(!aGen)
525     return;
526   
527   SALOMEDS::Study_ptr aStudy = aGen->GetCurrentStudy();
528   if(aStudy->_is_nil())
529     return;
530   
531   SALOMEDS::SObject_var aSO = aStudy->FindComponent(aGen->ComponentDataType());
532   if(CORBA::is_nil(aSO))
533     return;
534   
535   SALOMEDS::ChildIterator_var Itr = aStudy->NewChildIterator(aSO);
536   char* aParameters;
537   for(Itr->InitEx(true); Itr->More(); Itr->Next()) {
538     SALOMEDS::SObject_var aSObject = Itr->Value();
539     SALOMEDS::GenericAttribute_var anAttr;
540     if ( aSObject->FindAttribute(anAttr, "AttributeString")) {
541       aParameters = SALOMEDS::AttributeString::_narrow(anAttr)->Value();
542       SALOMEDS::ListOfListOfStrings_var aSections = aStudy->ParseVariables(aParameters);
543       if(MYDEBUG) {
544         cout<<"Entry : "<< aSObject->GetID()<<endl;
545         cout<<"aParameters : "<<aParameters<<endl;
546       }      
547       TCollection_AsciiString anObjType;
548       CORBA::Object_var anObject = SMESH_Gen_i::SObjectToObject(aSObject);
549       SMESH::SMESH_Hypothesis_var aHyp = SMESH::SMESH_Hypothesis::_narrow(anObject);
550       if(!aHyp->_is_nil()) {
551         anObjType = TCollection_AsciiString(aHyp->GetName());
552       }
553       else if(SMESH::SMESH_Mesh_var aMesh = SMESH::SMESH_Mesh::_narrow(anObject)) {
554         anObjType = TCollection_AsciiString("Mesh");
555       }
556       if(MYDEBUG)
557         cout<<"The object Type : "<<anObjType<<endl;
558       ObjectStates *aState = NULL;
559       if(anObjType == "LayerDistribution")
560         aState = new LayerDistributionStates();
561       else
562         aState = new  ObjectStates(anObjType);
563       
564       for(int i = 0; i < aSections->length(); i++) {
565         TState aVars;
566         SALOMEDS::ListOfStrings aListOfVars = aSections[i];
567         for(int j = 0;j<aListOfVars.length();j++) {
568           TCollection_AsciiString aVar(aListOfVars[j].in());
569           if(!aVar.IsEmpty() && aStudy->IsVariable(aVar.ToCString())) {
570             aVar.InsertBefore(1,"\"");
571             aVar.InsertAfter(aVar.Length(),"\"");
572           }
573           aVars.push_back(aVar);
574           if(MYDEBUG) {
575             cout<<"Variable: '"<<aVar<<"'"<<endl;
576           }
577         }
578         aState->AddState(aVars);
579       }
580       _objectMap.insert(pair<TCollection_AsciiString,ObjectStates*>(TCollection_AsciiString(aSObject->GetID()),aState));
581     }
582   }
583 }
584
585 //================================================================================
586 /*!
587  * 
588  */
589 //================================================================================
590 void SMESH_NoteBook::AddCommand(const TCollection_AsciiString& theString)
591 {
592   if(MYDEBUG)
593     cout<<theString<<endl;
594   Handle(_pyCommand) aCommand = new _pyCommand( theString, -1);
595   _commands.push_back(aCommand);
596
597   if ( aCommand->GetMethod() == "GetMeshEditor" ) { // MeshEditor creation
598     myMeshEditors.insert( make_pair( aCommand->GetResultValue(),
599                                      aCommand->GetObject() ) );
600   }
601 }
602
603 //================================================================================
604 /*!
605  * 
606  */
607 //================================================================================
608 void SMESH_NoteBook::ProcessLayerDistribution()
609 {
610   // 1) Find all LayerDistribution states
611   vector<LayerDistributionStates*> aLDS;
612   TVariablesMap::const_iterator it = _objectMap.begin();
613   for(;it != _objectMap.end();it++)
614     if(LayerDistributionStates* aLDStates = (LayerDistributionStates*)((*it).second)) {
615       aLDS.push_back(aLDStates);
616     }
617   
618   // 2) Initialize all type of 1D Distribution hypothesis
619   for(int i=0;i<_commands.size();i++){
620     for(int j =0;j < aLDS.size();j++){
621       TCollection_AsciiString aResultValue = _commands[i]->GetResultValue();
622       if(_commands[i]->GetMethod() == "CreateHypothesis" &&
623          aLDS[j]->HasDistribution(aResultValue)){
624         TCollection_AsciiString aType = _commands[i]->GetArg(1);
625         aType.RemoveAll('\'');
626         aLDS[j]->SetDistributionType(aResultValue,aType);
627       }
628     }
629   }
630   // 3) ... and replase variables ...
631
632   for(int i=0;i<_commands.size();i++){
633     for(int j =0;j < aLDS.size();j++){
634       TCollection_AsciiString anObject = _commands[i]->GetObject();
635
636       if(aLDS[j]->HasDistribution(anObject)) {
637         TCollection_AsciiString aType = aLDS[j]->GetDistributionType(anObject);
638         TCollection_AsciiString aMethod = _commands[i]->GetMethod();
639         if(aType == "LocalLength") {
640           if(aMethod == "SetLength") {
641             if(!aLDS[j]->GetCurrectState().at(0).IsEmpty() )
642               _commands[i]->SetArg(1,aLDS[j]->GetCurrectState().at(0));
643             aLDS[j]->IncrementState();
644           }
645           else if(aMethod == "SetPrecision") {
646             if(!aLDS[j]->GetCurrectState().at(1).IsEmpty() )
647               _commands[i]->SetArg(1,aLDS[j]->GetCurrectState().at(1));
648             aLDS[j]->IncrementState();
649           }
650         }
651       }
652     }
653   }
654 }
655 //================================================================================
656 /*!
657  *  \brief Return result script
658  */
659 //================================================================================
660 TCollection_AsciiString SMESH_NoteBook::GetResultScript() const
661 {
662   TCollection_AsciiString aResult;
663   for(int i=0;i<_commands.size();i++)
664     aResult+=_commands[i]->GetString()+"\n";
665   return aResult;
666 }
667
668 //================================================================================
669 /*!
670  *  \brief Return value of the variable
671  */
672 //================================================================================
673 bool SMESH_NoteBook::GetReal(const TCollection_AsciiString& theVarName, double& theValue)
674 {
675   bool ok = false;
676
677   SMESH_Gen_i *aGen = SMESH_Gen_i::GetSMESHGen();
678   if(!aGen)
679     return ok;
680
681   SALOMEDS::Study_ptr aStudy = aGen->GetCurrentStudy();
682   if(aStudy->_is_nil())
683     return ok;
684
685   TCollection_AsciiString aVarName = theVarName;
686   aVarName.RemoveAll('\"');
687
688   if(aVarName.IsEmpty())
689     return ok;
690
691   const char* aName = aVarName.ToCString();
692   if(aStudy->IsVariable(aName) && (aStudy->IsReal(aName) || aStudy->IsInteger(aName))) {
693     theValue = aStudy->GetReal(aVarName.ToCString());
694     ok = true;
695   }
696
697   return ok;
698 }