Salome HOME
c9c4fdc7f60569175f305433e685c0e75a2189aa
[modules/smesh.git] / src / SMESHClient / SMESH_Client.cxx
1 // Copyright (C) 2007-2016  CEA/DEN, EDF R&D, OPEN CASCADE
2 //
3 // Copyright (C) 2003-2007  OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
4 // CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
5 //
6 // This library is free software; you can redistribute it and/or
7 // modify it under the terms of the GNU Lesser General Public
8 // License as published by the Free Software Foundation; either
9 // version 2.1 of the License, or (at your option) any later version.
10 //
11 // This library is distributed in the hope that it will be useful,
12 // but WITHOUT ANY WARRANTY; without even the implied warranty of
13 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
14 // Lesser General Public License for more details.
15 //
16 // You should have received a copy of the GNU Lesser General Public
17 // License along with this library; if not, write to the Free Software
18 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
19 //
20 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
21 //
22 //  File   : SMESH_Client.cxx
23 //  Author : Pavel TELKOV
24 //  Module : SMESH
25
26 #include "SMESH_Client.hxx"
27 #include "SMESH_Mesh.hxx"
28 #include "SMESHDS_Script.hxx"
29
30 #include "SALOME_NamingService.hxx"
31 #include "SALOME_LifeCycleCORBA.hxx"
32
33 #include <SALOMEconfig.h>
34 #include CORBA_SERVER_HEADER(SALOME_Component)
35 #include CORBA_SERVER_HEADER(SALOME_Exception)
36
37 #include "Basics_Utils.hxx"
38 #include "utilities.h"
39
40 #ifdef WIN32
41 #include <process.h>
42 #else
43 #include <unistd.h>
44 #endif
45
46 #include <stdexcept>
47
48 #ifndef EXCEPTION
49 #define EXCEPTION(TYPE, MSG) {\
50   std::ostringstream aStream;\
51   aStream<<__FILE__<<"["<<__LINE__<<"]::"<<MSG;\
52   throw TYPE(aStream.str());\
53 }
54 #endif
55
56 #ifdef _DEBUG_
57 static int MYDEBUG = 0;
58 #else
59 static int MYDEBUG = 0;
60 #endif
61
62 namespace
63 {
64   using std::runtime_error;
65
66   //=======================================================================
67   //function : FindNode
68   //=======================================================================
69   inline const SMDS_MeshNode* FindNode(const SMDS_Mesh* theMesh, int theId){
70     if(const SMDS_MeshNode* anElem = theMesh->FindNode(theId)) return anElem;
71     EXCEPTION(runtime_error,"SMDS_Mesh::FindNode - cannot find a SMDS_MeshNode for ID = "<<theId);
72   }
73
74
75   //=======================================================================
76   //function : FindElement
77   //=======================================================================
78   inline const SMDS_MeshElement* FindElement(const SMDS_Mesh* theMesh, int theId){
79     if(const SMDS_MeshElement* anElem = theMesh->FindElement(theId)) return anElem;
80     EXCEPTION(runtime_error,"SMDS_Mesh::FindElement - cannot find a SMDS_MeshElement for ID = "<<theId);
81   }
82
83
84   //=======================================================================
85   //function : AddNodesWithID
86   //=======================================================================
87   inline void AddNodesWithID(SMDS_Mesh* theMesh,
88                              SMESH::log_array_var& theSeq,
89                              CORBA::Long theId)
90   {
91     const SMESH::double_array& aCoords = theSeq[theId].coords;
92     const SMESH::long_array& anIndexes = theSeq[theId].indexes;
93     CORBA::Long anElemId = 0, aNbElems = theSeq[theId].number;
94     if(3*aNbElems != (CORBA::Long) aCoords.length())
95       EXCEPTION(runtime_error,"AddNodesWithID - 3*aNbElems != aCoords.length()");
96     for(CORBA::Long aCoordId = 0; anElemId < aNbElems; anElemId++, aCoordId+=3){
97       SMDS_MeshElement* anElem = theMesh->AddNodeWithID(aCoords[aCoordId],
98                                                         aCoords[aCoordId+1],
99                                                         aCoords[aCoordId+2],
100                                                         anIndexes[anElemId]);
101       if(!anElem)
102         EXCEPTION(runtime_error,"SMDS_Mesh::FindElement - cannot AddNodeWithID for ID = "<<anElemId);
103     }
104   }
105
106
107   //=======================================================================
108   //function : Add0DElementsWithID
109   //=======================================================================
110   inline void Add0DElementsWithID(SMDS_Mesh* theMesh,
111                                   SMESH::log_array_var& theSeq,
112                                   CORBA::Long theId)
113   {
114     const SMESH::long_array& anIndexes = theSeq[theId].indexes;
115     CORBA::Long anElemId = 0, aNbElems = theSeq[theId].number;
116     if (2*aNbElems != (CORBA::Long) anIndexes.length())
117       EXCEPTION(runtime_error,"AddEdgeWithID - 2*aNbElems != aCoords.length()");
118     CORBA::Long anIndexId = 0;
119     for (; anElemId < aNbElems; anElemId++, anIndexId+=2)
120     {
121       SMDS_MeshElement* anElem = theMesh->Add0DElementWithID(anIndexes[anIndexId+1],
122                                                              anIndexes[anIndexId]);
123       if (!anElem)
124         EXCEPTION(runtime_error,"SMDS_Mesh::FindElement - cannot Add0DElementWithID for ID = "<<anElemId);
125     }
126   }
127
128
129   //=======================================================================
130   //function : AddBallsWithID
131   //=======================================================================
132   inline void AddBallsWithID(SMDS_Mesh*            theMesh,
133                              SMESH::log_array_var& theSeq,
134                              CORBA::Long           theId)
135   {
136     const SMESH::double_array& aDiameter = theSeq[theId].coords;
137     const SMESH::long_array& anIndexes   = theSeq[theId].indexes;
138     CORBA::Long anElemId = 0, aNbElems   = theSeq[theId].number;
139     if (2*aNbElems != (CORBA::Long) anIndexes.length() )
140       EXCEPTION(runtime_error,"AddEdgeWithID - 2*aNbElems != anIndexes.length()");
141     if (aNbElems != (CORBA::Long) aDiameter.length())
142       EXCEPTION(runtime_error,"AddEdgeWithID - aNbElems != aDiameter.length()");
143     CORBA::Long anIndexId = 0;
144     for (; anElemId < aNbElems; anElemId++, anIndexId+=2)
145     {
146       SMDS_MeshElement* anElem = theMesh->AddBallWithID(anIndexes[anIndexId+1],
147                                                         aDiameter[anElemId],
148                                                         anIndexes[anIndexId]);
149       if (!anElem)
150         EXCEPTION(runtime_error,"cannot SMDS_Mesh::AddBallsWithID for ID = "<<anElemId);
151     }
152   }
153
154
155   //=======================================================================
156   //function : AddEdgesWithID
157   //=======================================================================
158   inline void AddEdgesWithID(SMDS_Mesh* theMesh,
159                              SMESH::log_array_var& theSeq,
160                              CORBA::Long theId)
161   {
162     const SMESH::long_array& anIndexes = theSeq[theId].indexes;
163     CORBA::Long anElemId = 0, aNbElems = theSeq[theId].number;
164     if(3*aNbElems != (CORBA::Long) anIndexes.length())
165       EXCEPTION(runtime_error,"AddEdgeWithID - 3*aNbElems != aCoords.length()");
166     for(CORBA::Long anIndexId = 0; anElemId < aNbElems; anElemId++, anIndexId+=3){
167       SMDS_MeshElement* anElem = theMesh->AddEdgeWithID(anIndexes[anIndexId+1],
168                                                         anIndexes[anIndexId+2],
169                                                         anIndexes[anIndexId]);
170       if(!anElem)
171         EXCEPTION(runtime_error,"SMDS_Mesh::FindElement - cannot AddEdgeWithID for ID = "<<anElemId);
172     }
173   }
174
175
176   //=======================================================================
177   //function : AddTriasWithID
178   //=======================================================================
179   inline void AddTriasWithID(SMDS_Mesh* theMesh,
180                              SMESH::log_array_var& theSeq,
181                              CORBA::Long theId)
182   {
183     const SMESH::long_array& anIndexes = theSeq[theId].indexes;
184     CORBA::Long anElemId = 0, aNbElems = theSeq[theId].number;
185     if(4*aNbElems != (CORBA::Long) anIndexes.length())
186       EXCEPTION(runtime_error,"AddTriasWithID - 4*aNbElems != anIndexes.length()");
187     for(CORBA::Long anIndexId = 0; anElemId < aNbElems; anElemId++, anIndexId+=4){
188       SMDS_MeshElement* anElem = theMesh->AddFaceWithID(anIndexes[anIndexId+1],
189                                                         anIndexes[anIndexId+2],
190                                                         anIndexes[anIndexId+3],
191                                                         anIndexes[anIndexId]);
192       if(!anElem)
193         EXCEPTION(runtime_error,"SMDS_Mesh::FindElement - cannot AddFaceWithID for ID = "<<anElemId);
194     }
195   }
196
197
198   //=======================================================================
199   //function : AddQuadsWithID
200   //=======================================================================
201   inline void AddQuadsWithID(SMDS_Mesh* theMesh,
202                              SMESH::log_array_var theSeq,
203                              CORBA::Long theId)
204   {
205     const SMESH::long_array& anIndexes = theSeq[theId].indexes;
206     CORBA::Long anElemId = 0, aNbElems = theSeq[theId].number;
207     if(5*aNbElems != (CORBA::Long) anIndexes.length())
208       EXCEPTION(runtime_error,"AddQuadsWithID - 4*aNbElems != anIndexes.length()");
209     for(CORBA::Long anIndexId = 0; anElemId < aNbElems; anElemId++, anIndexId+=5){
210       SMDS_MeshElement* anElem = theMesh->AddFaceWithID(anIndexes[anIndexId+1],
211                                                         anIndexes[anIndexId+2],
212                                                         anIndexes[anIndexId+3],
213                                                         anIndexes[anIndexId+4],
214                                                         anIndexes[anIndexId]);
215       if(!anElem)
216         EXCEPTION(runtime_error,"SMDS_Mesh::FindElement - cannot AddFaceWithID for ID = "<<anElemId);
217     }
218   }
219
220
221   //=======================================================================
222   //function : AddPolygonsWithID
223   //=======================================================================
224   inline void AddPolygonsWithID(SMDS_Mesh* theMesh,
225                                 SMESH::log_array_var& theSeq,
226                                 CORBA::Long theId)
227   {
228     const SMESH::long_array& anIndexes = theSeq[theId].indexes;
229     CORBA::Long anIndexId = 0, aNbElems = theSeq[theId].number;
230
231     for (CORBA::Long anElemId = 0; anElemId < aNbElems; anElemId++) {
232       int aFaceId = anIndexes[anIndexId++];
233
234       int aNbNodes = anIndexes[anIndexId++];
235       std::vector<int> nodes_ids (aNbNodes);
236       for (int i = 0; i < aNbNodes; i++) {
237         nodes_ids[i] = anIndexes[anIndexId++];
238       }
239
240       SMDS_MeshElement* anElem = theMesh->AddPolygonalFaceWithID(nodes_ids, aFaceId);
241       if (!anElem)
242         EXCEPTION(runtime_error, "SMDS_Mesh::FindElement - cannot AddPolygonalFaceWithID for ID = "
243                   << anElemId);
244     }
245   }
246
247
248   //=======================================================================
249   //function : AddQaudPolygonsWithID
250   //=======================================================================
251   inline void AddQuadPolygonsWithID(SMDS_Mesh* theMesh,
252                                     SMESH::log_array_var& theSeq,
253                                     CORBA::Long theId)
254   {
255     const SMESH::long_array& anIndexes = theSeq[theId].indexes;
256     CORBA::Long anIndexId = 0, aNbElems = theSeq[theId].number;
257
258     for (CORBA::Long anElemId = 0; anElemId < aNbElems; anElemId++) {
259       int aFaceId = anIndexes[anIndexId++];
260
261       int aNbNodes = anIndexes[anIndexId++];
262       std::vector<int> nodes_ids (aNbNodes);
263       for (int i = 0; i < aNbNodes; i++) {
264         nodes_ids[i] = anIndexes[anIndexId++];
265       }
266
267       SMDS_MeshElement* anElem = theMesh->AddQuadPolygonalFaceWithID(nodes_ids, aFaceId);
268       if (!anElem)
269         EXCEPTION(runtime_error, "SMDS_Mesh::FindElement - cannot AddQaudPolygonalFaceWithID for ID = "
270                   << anElemId);
271     }
272   }
273
274
275   //=======================================================================
276   //function : AddTetrasWithID
277   //=======================================================================
278   inline void AddTetrasWithID(SMDS_Mesh* theMesh,
279                               SMESH::log_array_var& theSeq,
280                               CORBA::Long theId)
281   {
282     const SMESH::long_array& anIndexes = theSeq[theId].indexes;
283     CORBA::Long anElemId = 0, aNbElems = theSeq[theId].number;
284     if(5*aNbElems != (CORBA::Long) anIndexes.length())
285       EXCEPTION(runtime_error,"AddTetrasWithID - 5*aNbElems != anIndexes.length()");
286     for(CORBA::Long anIndexId = 0; anElemId < aNbElems; anElemId++, anIndexId+=5){
287       SMDS_MeshElement* anElem = theMesh->AddVolumeWithID(anIndexes[anIndexId+1],
288                                                           anIndexes[anIndexId+2],
289                                                           anIndexes[anIndexId+3],
290                                                           anIndexes[anIndexId+4],
291                                                           anIndexes[anIndexId]);
292       if(!anElem)
293         EXCEPTION(runtime_error,"SMDS_Mesh::FindElement - cannot AddVolumeWithID for ID = "<<anElemId);
294     }
295   }
296
297
298   //=======================================================================
299   //function : AddPiramidsWithID
300   //=======================================================================
301   inline void AddPiramidsWithID(SMDS_Mesh* theMesh,
302                                 SMESH::log_array_var& theSeq,
303                                 CORBA::Long theId)
304   {
305     const SMESH::long_array& anIndexes = theSeq[theId].indexes;
306     CORBA::Long anElemId = 0, aNbElems = theSeq[theId].number;
307     if(6*aNbElems != (CORBA::Long) anIndexes.length())
308       EXCEPTION(runtime_error,"AddPiramidsWithID - 6*aNbElems != anIndexes.length()");
309     for(CORBA::Long anIndexId = 0; anElemId < aNbElems; anElemId++, anIndexId+=6){
310       SMDS_MeshElement* anElem = theMesh->AddVolumeWithID(anIndexes[anIndexId+1],
311                                                           anIndexes[anIndexId+2],
312                                                           anIndexes[anIndexId+3],
313                                                           anIndexes[anIndexId+4],
314                                                           anIndexes[anIndexId+5],
315                                                           anIndexes[anIndexId]);
316       if(!anElem)
317         EXCEPTION(runtime_error,"SMDS_Mesh::FindElement - cannot AddVolumeWithID for ID = "<<anElemId);
318     }
319   }
320
321
322   //=======================================================================
323   //function : AddPrismsWithID
324   //=======================================================================
325   inline void AddPrismsWithID(SMDS_Mesh* theMesh,
326                               SMESH::log_array_var& theSeq,
327                               CORBA::Long theId)
328   {
329     const SMESH::long_array& anIndexes = theSeq[theId].indexes;
330     CORBA::Long anElemId = 0, aNbElems = theSeq[theId].number;
331     if(7*aNbElems != (CORBA::Long) anIndexes.length())
332       EXCEPTION(runtime_error,"AddPrismsWithID - 7*aNbElems != anIndexes.length()");
333     for(CORBA::Long anIndexId = 0; anElemId < aNbElems; anElemId++, anIndexId+=7){
334       SMDS_MeshElement* anElem = theMesh->AddVolumeWithID(anIndexes[anIndexId+1],
335                                                           anIndexes[anIndexId+2],
336                                                           anIndexes[anIndexId+3],
337                                                           anIndexes[anIndexId+4],
338                                                           anIndexes[anIndexId+5],
339                                                           anIndexes[anIndexId+6],
340                                                           anIndexes[anIndexId]);
341       if(!anElem)
342         EXCEPTION(runtime_error,"SMDS_Mesh::FindElement - cannot AddVolumeWithID for ID = "<<anElemId);
343     }
344   }
345
346
347   //=======================================================================
348   //function : AddHexasWithID
349   //=======================================================================
350   inline void AddHexasWithID(SMDS_Mesh* theMesh,
351                              SMESH::log_array_var& theSeq,
352                              CORBA::Long theId)
353   {
354     const SMESH::long_array& anIndexes = theSeq[theId].indexes;
355     CORBA::Long anElemId = 0, aNbElems = theSeq[theId].number;
356     if(9*aNbElems != (CORBA::Long) anIndexes.length())
357       EXCEPTION(runtime_error,"AddHexasWithID - 9*aNbElems != anIndexes.length()");
358     for(CORBA::Long anIndexId = 0; anElemId < aNbElems; anElemId++, anIndexId+=9){
359       SMDS_MeshElement* anElem = theMesh->AddVolumeWithID(anIndexes[anIndexId+1],
360                                                           anIndexes[anIndexId+2],
361                                                           anIndexes[anIndexId+3],
362                                                           anIndexes[anIndexId+4],
363                                                           anIndexes[anIndexId+5],
364                                                           anIndexes[anIndexId+6],
365                                                           anIndexes[anIndexId+7],
366                                                           anIndexes[anIndexId+8],
367                                                           anIndexes[anIndexId]);
368       if(!anElem)
369         EXCEPTION(runtime_error,"SMDS_Mesh::FindElement - cannot AddVolumeWithID for ID = "<<anElemId);
370     }
371   }
372
373   //=======================================================================
374   //function : AddHexPrismWithID
375   //=======================================================================
376   inline void AddHexPrismWithID(SMDS_Mesh* theMesh,
377                                 SMESH::log_array_var& theSeq,
378                                 CORBA::Long theId)
379   {
380     const SMESH::long_array& anIndexes = theSeq[theId].indexes;
381     CORBA::Long anElemId = 0, aNbElems = theSeq[theId].number;
382     if(13*aNbElems != (CORBA::Long) anIndexes.length())
383       EXCEPTION(runtime_error,"AddHexPrismWithID - 13*aNbElems != anIndexes.length()");
384     for(CORBA::Long anIndexId = 0; anElemId < aNbElems; anElemId++, anIndexId+=13){
385       SMDS_MeshElement* anElem = theMesh->AddVolumeWithID(anIndexes[anIndexId+1],
386                                                           anIndexes[anIndexId+2],
387                                                           anIndexes[anIndexId+3],
388                                                           anIndexes[anIndexId+4],
389                                                           anIndexes[anIndexId+5],
390                                                           anIndexes[anIndexId+6],
391                                                           anIndexes[anIndexId+7],
392                                                           anIndexes[anIndexId+8],
393                                                           anIndexes[anIndexId+9],
394                                                           anIndexes[anIndexId+10],
395                                                           anIndexes[anIndexId+11],
396                                                           anIndexes[anIndexId+12],
397                                                           anIndexes[anIndexId]);
398       if(!anElem)
399         EXCEPTION(runtime_error,"AddHexPrismWithID - cannot AddVolumeWithID for ID = "<<anElemId);
400     }
401   }
402
403
404   //=======================================================================
405   //function : AddPolyhedronsWithID
406   //=======================================================================
407   inline void AddPolyhedronsWithID (SMDS_Mesh* theMesh,
408                                     SMESH::log_array_var& theSeq,
409                                     CORBA::Long theId)
410   {
411     const SMESH::long_array& anIndexes = theSeq[theId].indexes;
412     CORBA::Long anIndexId = 0, aNbElems = theSeq[theId].number;
413
414     for (CORBA::Long anElemId = 0; anElemId < aNbElems; anElemId++) {
415       int aFaceId = anIndexes[anIndexId++];
416
417       int aNbNodes = anIndexes[anIndexId++];
418       std::vector<int> nodes_ids (aNbNodes);
419       for (int i = 0; i < aNbNodes; i++) {
420         nodes_ids[i] = anIndexes[anIndexId++];
421       }
422
423       int aNbFaces = anIndexes[anIndexId++];
424       std::vector<int> quantities (aNbFaces);
425       for (int i = 0; i < aNbFaces; i++) {
426         quantities[i] = anIndexes[anIndexId++];
427       }
428
429       SMDS_MeshElement* anElem =
430         theMesh->AddPolyhedralVolumeWithID(nodes_ids, quantities, aFaceId);
431       if (!anElem)
432         EXCEPTION(runtime_error, "SMDS_Mesh::FindElement - cannot AddPolyhedralVolumeWithID for ID = "
433                   << anElemId);
434     }
435   }
436
437
438   //=======================================================================
439   //function : AddQuadEdgesWithID
440   //=======================================================================
441   inline void AddQuadEdgesWithID(SMDS_Mesh* theMesh,
442                                  SMESH::log_array_var& theSeq,
443                                  CORBA::Long theId)
444   {
445     const SMESH::long_array& anIndexes = theSeq[theId].indexes;
446     CORBA::Long anElemId = 0, aNbElems = theSeq[theId].number;
447     if(4*aNbElems != (CORBA::Long) anIndexes.length())
448       EXCEPTION(runtime_error,"AddQuadEdgeWithID - 4*aNbElems != aCoords.length()");
449     for(CORBA::Long anIndexId = 0; anElemId < aNbElems; anElemId++, anIndexId+=4){
450       SMDS_MeshElement* anElem = theMesh->AddEdgeWithID(anIndexes[anIndexId+1],
451                                                         anIndexes[anIndexId+2],
452                                                         anIndexes[anIndexId+3],
453                                                         anIndexes[anIndexId]);
454       if(!anElem)
455         EXCEPTION(runtime_error,"SMDS_Mesh::FindElement - cannot AddEdgeWithID for ID = "<<anElemId);
456     }
457   }
458
459
460   //=======================================================================
461   //function : AddQuadTriasWithID
462   //=======================================================================
463   inline void AddQuadTriasWithID(SMDS_Mesh* theMesh,
464                                  SMESH::log_array_var& theSeq,
465                                  CORBA::Long theId)
466   {
467     const SMESH::long_array& anIndexes = theSeq[theId].indexes;
468     CORBA::Long anElemId = 0, aNbElems = theSeq[theId].number;
469     if(7*aNbElems != (CORBA::Long) anIndexes.length())
470       EXCEPTION(runtime_error,"AddQuadTriasWithID - 7*aNbElems != anIndexes.length()");
471     for(CORBA::Long anIndexId = 0; anElemId < aNbElems; anElemId++, anIndexId+=7){
472       SMDS_MeshElement* anElem = theMesh->AddFaceWithID(anIndexes[anIndexId+1],
473                                                         anIndexes[anIndexId+2],
474                                                         anIndexes[anIndexId+3],
475                                                         anIndexes[anIndexId+4],
476                                                         anIndexes[anIndexId+5],
477                                                         anIndexes[anIndexId+6],
478                                                         anIndexes[anIndexId]);
479       if(!anElem)
480         EXCEPTION(runtime_error,"SMDS_Mesh::FindElement - cannot AddFaceWithID for ID = "<<anElemId);
481     }
482   }
483
484
485   //=======================================================================
486   //function : AddQuadQuadsWithID
487   //=======================================================================
488   inline void AddQuadQuadsWithID(SMDS_Mesh* theMesh,
489                                  SMESH::log_array_var theSeq,
490                                  CORBA::Long theId)
491   {
492     const SMESH::long_array& anIndexes = theSeq[theId].indexes;
493     CORBA::Long anElemId = 0, aNbElems = theSeq[theId].number;
494     if(9*aNbElems != (CORBA::Long) anIndexes.length())
495       EXCEPTION(runtime_error,"AddQuadQuadsWithID - 9*aNbElems != anIndexes.length()");
496     for(CORBA::Long anIndexId = 0; anElemId < aNbElems; anElemId++, anIndexId+=9){
497       SMDS_MeshElement* anElem = theMesh->AddFaceWithID(anIndexes[anIndexId+1],
498                                                         anIndexes[anIndexId+2],
499                                                         anIndexes[anIndexId+3],
500                                                         anIndexes[anIndexId+4],
501                                                         anIndexes[anIndexId+5],
502                                                         anIndexes[anIndexId+6],
503                                                         anIndexes[anIndexId+7],
504                                                         anIndexes[anIndexId+8],
505                                                         anIndexes[anIndexId]);
506       if(!anElem)
507         EXCEPTION(runtime_error,"SMDS_Mesh::FindElement - cannot AddFaceWithID for ID = "<<anElemId);
508     }
509   }
510
511   //=======================================================================
512   //function : AddBiQuadQuadsWithID
513   //=======================================================================
514   inline void AddBiQuadQuadsWithID(SMDS_Mesh* theMesh,
515                                    SMESH::log_array_var theSeq,
516                                    CORBA::Long theId)
517   {
518     const SMESH::long_array& anIndexes = theSeq[theId].indexes;
519     CORBA::Long anElemId = 0, aNbElems = theSeq[theId].number;
520     if(10*aNbElems != (CORBA::Long) anIndexes.length())
521       EXCEPTION(runtime_error,"AddBiQuadQuadsWithID - 10*aNbElems != anIndexes.length()");
522     for(CORBA::Long anIndexId = 0; anElemId < aNbElems; anElemId++, anIndexId+=10){
523       SMDS_MeshElement* anElem = theMesh->AddFaceWithID(anIndexes[anIndexId+1],
524                                                         anIndexes[anIndexId+2],
525                                                         anIndexes[anIndexId+3],
526                                                         anIndexes[anIndexId+4],
527                                                         anIndexes[anIndexId+5],
528                                                         anIndexes[anIndexId+6],
529                                                         anIndexes[anIndexId+7],
530                                                         anIndexes[anIndexId+8],
531                                                         anIndexes[anIndexId+9],
532                                                         anIndexes[anIndexId]);
533       if(!anElem)
534         EXCEPTION(runtime_error,"AddBiQuadQuadsWithID() - cannot AddFaceWithID for ID = "<<anElemId);
535     }
536   }
537
538   //=======================================================================
539   //function : AddBiQuadTriasWithID
540   //=======================================================================
541   inline void AddBiQuadTriasWithID(SMDS_Mesh* theMesh,
542                                    SMESH::log_array_var theSeq,
543                                    CORBA::Long theId)
544   {
545     const SMESH::long_array& anIndexes = theSeq[theId].indexes;
546     CORBA::Long anElemId = 0, aNbElems = theSeq[theId].number;
547     if(8*aNbElems != (CORBA::Long) anIndexes.length())
548       EXCEPTION(runtime_error,"AddBiQuadTriasWithID - 8*aNbElems != anIndexes.length()");
549     for(CORBA::Long anIndexId = 0; anElemId < aNbElems; anElemId++, anIndexId+=8){
550       SMDS_MeshElement* anElem = theMesh->AddFaceWithID(anIndexes[anIndexId+1],
551                                                         anIndexes[anIndexId+2],
552                                                         anIndexes[anIndexId+3],
553                                                         anIndexes[anIndexId+4],
554                                                         anIndexes[anIndexId+5],
555                                                         anIndexes[anIndexId+6],
556                                                         anIndexes[anIndexId+7],
557                                                         anIndexes[anIndexId]);
558       if(!anElem)
559         EXCEPTION(runtime_error,"AddBiQuadTriasWithID() - cannot AddFaceWithID for ID = "<<anElemId);
560     }
561   }
562
563
564   //=======================================================================
565   //function : AddQuadTetrasWithID
566   //=======================================================================
567   inline void AddQuadTetrasWithID(SMDS_Mesh* theMesh,
568                                   SMESH::log_array_var& theSeq,
569                                   CORBA::Long theId)
570   {
571     const SMESH::long_array& anIndexes = theSeq[theId].indexes;
572     CORBA::Long anElemId = 0, aNbElems = theSeq[theId].number;
573     if(11*aNbElems != (CORBA::Long) anIndexes.length())
574       EXCEPTION(runtime_error,"AddQuadTetrasWithID - 11*aNbElems != anIndexes.length()");
575     for(CORBA::Long anIndexId = 0; anElemId < aNbElems; anElemId++, anIndexId+=11){
576       SMDS_MeshElement* anElem = theMesh->AddVolumeWithID(anIndexes[anIndexId+1],
577                                                           anIndexes[anIndexId+2],
578                                                           anIndexes[anIndexId+3],
579                                                           anIndexes[anIndexId+4],
580                                                           anIndexes[anIndexId+5],
581                                                           anIndexes[anIndexId+6],
582                                                           anIndexes[anIndexId+7],
583                                                           anIndexes[anIndexId+8],
584                                                           anIndexes[anIndexId+9],
585                                                           anIndexes[anIndexId+10],
586                                                           anIndexes[anIndexId]);
587       if(!anElem)
588         EXCEPTION(runtime_error,"SMDS_Mesh::FindElement - cannot AddVolumeWithID for ID = "<<anElemId);
589     }
590   }
591
592
593   //=======================================================================
594   //function : AddQuadPiramidsWithID
595   //=======================================================================
596   inline void AddQuadPiramidsWithID(SMDS_Mesh* theMesh,
597                                     SMESH::log_array_var& theSeq,
598                                     CORBA::Long theId)
599   {
600     const SMESH::long_array& anIndexes = theSeq[theId].indexes;
601     CORBA::Long anElemId = 0, aNbElems = theSeq[theId].number;
602     if(14*aNbElems != (CORBA::Long) anIndexes.length())
603       EXCEPTION(runtime_error,"AddQuadPiramidsWithID - 14*aNbElems != anIndexes.length()");
604     for(CORBA::Long anIndexId = 0; anElemId < aNbElems; anElemId++, anIndexId+=14){
605       SMDS_MeshElement* anElem = theMesh->AddVolumeWithID(anIndexes[anIndexId+1],
606                                                           anIndexes[anIndexId+2],
607                                                           anIndexes[anIndexId+3],
608                                                           anIndexes[anIndexId+4],
609                                                           anIndexes[anIndexId+5],
610                                                           anIndexes[anIndexId+6],
611                                                           anIndexes[anIndexId+7],
612                                                           anIndexes[anIndexId+8],
613                                                           anIndexes[anIndexId+9],
614                                                           anIndexes[anIndexId+10],
615                                                           anIndexes[anIndexId+11],
616                                                           anIndexes[anIndexId+12],
617                                                           anIndexes[anIndexId+13],
618                                                           anIndexes[anIndexId]);
619       if(!anElem)
620         EXCEPTION(runtime_error,"SMDS_Mesh::FindElement - cannot AddVolumeWithID for ID = "<<anElemId);
621     }
622   }
623
624
625   //=======================================================================
626   //function : AddQuadPentasWithID
627   //=======================================================================
628   inline void AddQuadPentasWithID(SMDS_Mesh* theMesh,
629                                   SMESH::log_array_var& theSeq,
630                                   CORBA::Long theId)
631   {
632     const SMESH::long_array& anIndexes = theSeq[theId].indexes;
633     CORBA::Long anElemId = 0, aNbElems = theSeq[theId].number;
634     if(16*aNbElems != (CORBA::Long) anIndexes.length())
635       EXCEPTION(runtime_error,"AddQuadPentasWithID - 16*aNbElems != anIndexes.length()");
636     for(CORBA::Long anIndexId = 0; anElemId < aNbElems; anElemId++, anIndexId+=16){
637       SMDS_MeshElement* anElem = theMesh->AddVolumeWithID(anIndexes[anIndexId+1],
638                                                           anIndexes[anIndexId+2],
639                                                           anIndexes[anIndexId+3],
640                                                           anIndexes[anIndexId+4],
641                                                           anIndexes[anIndexId+5],
642                                                           anIndexes[anIndexId+6],
643                                                           anIndexes[anIndexId+7],
644                                                           anIndexes[anIndexId+8],
645                                                           anIndexes[anIndexId+9],
646                                                           anIndexes[anIndexId+10],
647                                                           anIndexes[anIndexId+11],
648                                                           anIndexes[anIndexId+12],
649                                                           anIndexes[anIndexId+13],
650                                                           anIndexes[anIndexId+14],
651                                                           anIndexes[anIndexId+15],
652                                                           anIndexes[anIndexId]);
653       if(!anElem)
654         EXCEPTION(runtime_error,"SMDS_Mesh::FindElement - cannot AddVolumeWithID for ID = "<<anElemId);
655     }
656   }
657
658
659   //=======================================================================
660   //function : AddQuadHexasWithID
661   //=======================================================================
662   inline void AddQuadHexasWithID(SMDS_Mesh* theMesh,
663                                  SMESH::log_array_var& theSeq,
664                                  CORBA::Long theId)
665   {
666     const SMESH::long_array& anIndexes = theSeq[theId].indexes;
667     CORBA::Long anElemId = 0, aNbElems = theSeq[theId].number;
668     if(21*aNbElems != (CORBA::Long) anIndexes.length())
669       EXCEPTION(runtime_error,"AddQuadHexasWithID - 21*aNbElems != anIndexes.length()");
670     for(CORBA::Long anIndexId = 0; anElemId < aNbElems; anElemId++, anIndexId+=21){
671       SMDS_MeshElement* anElem = theMesh->AddVolumeWithID(anIndexes[anIndexId+1],
672                                                           anIndexes[anIndexId+2],
673                                                           anIndexes[anIndexId+3],
674                                                           anIndexes[anIndexId+4],
675                                                           anIndexes[anIndexId+5],
676                                                           anIndexes[anIndexId+6],
677                                                           anIndexes[anIndexId+7],
678                                                           anIndexes[anIndexId+8],
679                                                           anIndexes[anIndexId+9],
680                                                           anIndexes[anIndexId+10],
681                                                           anIndexes[anIndexId+11],
682                                                           anIndexes[anIndexId+12],
683                                                           anIndexes[anIndexId+13],
684                                                           anIndexes[anIndexId+14],
685                                                           anIndexes[anIndexId+15],
686                                                           anIndexes[anIndexId+16],
687                                                           anIndexes[anIndexId+17],
688                                                           anIndexes[anIndexId+18],
689                                                           anIndexes[anIndexId+19],
690                                                           anIndexes[anIndexId+20],
691                                                           anIndexes[anIndexId]);
692       if(!anElem)
693         EXCEPTION(runtime_error,"SMDS_Mesh::FindElement - cannot AddVolumeWithID for ID = "<<anElemId);
694     }
695   }
696
697   //=======================================================================
698   //function : AddTriQuadHexasWithID
699   //=======================================================================
700   inline void AddTriQuadHexasWithID(SMDS_Mesh* theMesh,
701                                     SMESH::log_array_var& theSeq,
702                                     CORBA::Long theId)
703   {
704     const SMESH::long_array& anIndexes = theSeq[theId].indexes;
705     CORBA::Long anElemId = 0, aNbElems = theSeq[theId].number;
706     if(28*aNbElems != (CORBA::Long) anIndexes.length())
707       EXCEPTION(runtime_error,"AddTriQuadHexasWithID - 28*aNbElems != anIndexes.length()");
708     for(CORBA::Long anIndexId = 0; anElemId < aNbElems; anElemId++, anIndexId+=28){
709       SMDS_MeshElement* anElem = theMesh->AddVolumeWithID(anIndexes[anIndexId+1],
710                                                           anIndexes[anIndexId+2],
711                                                           anIndexes[anIndexId+3],
712                                                           anIndexes[anIndexId+4],
713                                                           anIndexes[anIndexId+5],
714                                                           anIndexes[anIndexId+6],
715                                                           anIndexes[anIndexId+7],
716                                                           anIndexes[anIndexId+8],
717                                                           anIndexes[anIndexId+9],
718                                                           anIndexes[anIndexId+10],
719                                                           anIndexes[anIndexId+11],
720                                                           anIndexes[anIndexId+12],
721                                                           anIndexes[anIndexId+13],
722                                                           anIndexes[anIndexId+14],
723                                                           anIndexes[anIndexId+15],
724                                                           anIndexes[anIndexId+16],
725                                                           anIndexes[anIndexId+17],
726                                                           anIndexes[anIndexId+18],
727                                                           anIndexes[anIndexId+19],
728                                                           anIndexes[anIndexId+20],
729                                                           anIndexes[anIndexId+21],
730                                                           anIndexes[anIndexId+22],
731                                                           anIndexes[anIndexId+23],
732                                                           anIndexes[anIndexId+24],
733                                                           anIndexes[anIndexId+25],
734                                                           anIndexes[anIndexId+26],
735                                                           anIndexes[anIndexId+27],
736                                                           anIndexes[anIndexId]);
737       if(!anElem)
738         EXCEPTION(runtime_error,"AddTriQuadHexasWithID() - cannot AddVolumeWithID for ID = "<<anElemId);
739     }
740   }
741
742
743   //=======================================================================
744   //function : ChangePolyhedronNodes
745   //=======================================================================
746   inline void ChangePolyhedronNodes (SMDS_Mesh* theMesh,
747                                      SMESH::log_array_var& theSeq,
748                                      CORBA::Long theId)
749   {
750     const SMESH::long_array& anIndexes = theSeq[theId].indexes;
751     CORBA::Long iind = 0, aNbElems = theSeq[theId].number;
752
753     for (CORBA::Long anElemId = 0; anElemId < aNbElems; anElemId++)
754     {
755       // find element
756       const SMDS_MeshElement* elem = FindElement(theMesh, anIndexes[iind++]);
757       // nb nodes
758       int nbNodes = anIndexes[iind++];
759       // nodes
760       std::vector<const SMDS_MeshNode*> aNodes (nbNodes);
761       for (int iNode = 0; iNode < nbNodes; iNode++) {
762         aNodes[iNode] = FindNode(theMesh, anIndexes[iind++]);
763       }
764       // nb faces
765       int nbFaces = anIndexes[iind++];
766       // quantities
767       std::vector<int> quantities (nbFaces);
768       for (int iFace = 0; iFace < nbFaces; iFace++) {
769         quantities[iFace] = anIndexes[iind++];
770       }
771       // change
772       theMesh->ChangePolyhedronNodes(elem, aNodes, quantities);
773     }
774   }
775 }
776
777 //=======================================================================
778 SMESH::SMESH_Gen_var
779 SMESH_Client::GetSMESHGen(CORBA::ORB_ptr theORB,
780                           CORBA::Boolean& theIsEmbeddedMode)
781 {
782   static SMESH::SMESH_Gen_var aMeshGen;
783
784   if(CORBA::is_nil(aMeshGen.in())){
785 #ifdef WIN32
786     long aClientPID = (long)_getpid();
787 #else
788     long aClientPID =  (long)getpid();
789 #endif
790
791     SALOME_NamingService aNamingService(theORB);
792     SALOME_LifeCycleCORBA aLifeCycleCORBA(&aNamingService);
793     Engines::EngineComponent_var aComponent = aLifeCycleCORBA.FindOrLoad_Component("FactoryServer","SMESH");
794     aMeshGen = SMESH::SMESH_Gen::_narrow(aComponent);
795
796     std::string aClientHostName = Kernel_Utils::GetHostname();
797     Engines::Container_var aServerContainer = aMeshGen->GetContainerRef();
798     CORBA::String_var aServerHostName = aServerContainer->getHostName();
799     CORBA::Long aServerPID = aServerContainer->getPID();
800     aMeshGen->SetEmbeddedMode((aClientPID == aServerPID) && (aClientHostName == aServerHostName.in()));
801   }
802   theIsEmbeddedMode = aMeshGen->IsEmbeddedMode();
803
804   return aMeshGen;
805 }
806
807
808 //=======================================================================
809 // function : Create()
810 // purpose  :
811 //=======================================================================
812 SMESH_Client::SMESH_Client(CORBA::ORB_ptr theORB,
813                            SMESH::SMESH_Mesh_ptr theMesh):
814   myMeshServer(SMESH::SMESH_Mesh::_duplicate(theMesh)),
815   mySMESHDSMesh(NULL),
816   mySMDSMesh(NULL)
817 {
818   MESSAGE("SMESH_Client::SMESH_Client");
819   myMeshServer->Register();
820
821   CORBA::Boolean anIsEmbeddedMode;
822   GetSMESHGen(theORB,anIsEmbeddedMode);
823   if(anIsEmbeddedMode){
824     if ( MYDEBUG )
825       MESSAGE("Info: The same process, update mesh by pointer ");
826     // just set client mesh pointer to server mesh pointer
827     //SMESH_Mesh* aMesh = reinterpret_cast<SMESH_Mesh*>(theMesh->GetMeshPtr());
828     CORBA::LongLong pointeur = theMesh->GetMeshPtr();
829     if( MYDEBUG )
830       MESSAGE("SMESH_Client::SMESH_Client pointeur "<<pointeur);
831     SMESH_Mesh* aMesh = reinterpret_cast<SMESH_Mesh*> (pointeur);
832     if ( MYDEBUG )
833       MESSAGE("SMESH_Client::SMESH_Client aMesh "<<aMesh);
834     //if(aMesh->GetMeshDS()->IsEmbeddedMode()){
835     if(anIsEmbeddedMode){
836       mySMESHDSMesh = aMesh->GetMeshDS();
837       mySMDSMesh = mySMESHDSMesh;
838     }
839   }
840   if(!mySMDSMesh)
841     mySMDSMesh = new SMDS_Mesh();
842 }
843
844
845 //=================================================================================
846 // function : ~SMESH_Client
847 // purpose  : Destructor
848 //=================================================================================
849 SMESH_Client::~SMESH_Client()
850 {
851   myMeshServer->UnRegister();
852   if(!mySMESHDSMesh)
853     delete mySMDSMesh;
854 }
855
856
857 //=================================================================================
858 SMDS_Mesh*
859 SMESH_Client::GetMesh() const
860 {
861   return mySMDSMesh;
862 }
863
864
865 //=================================================================================
866 SMDS_Mesh*
867 SMESH_Client::operator->() const
868 {
869   return GetMesh();
870 }
871
872
873 //=================================================================================
874 SMESH::SMESH_Mesh_ptr
875 SMESH_Client::GetMeshServer()
876 {
877   return myMeshServer.in();
878 }
879
880
881 //=================================================================================
882 // function : SMESH_Client
883 // purpose  : Update mesh
884 //=================================================================================
885 bool
886 SMESH_Client::Update(bool theIsClear)
887 {
888   bool anIsModified = true;
889   if(mySMESHDSMesh)
890   {
891     if ( MYDEBUG ) MESSAGE("Update mySMESHDSMesh");
892     SMESHDS_Script* aScript = mySMESHDSMesh->GetScript();
893     anIsModified = aScript->IsModified();
894     aScript->SetModified(false);
895   }
896   else
897   {
898     if ( MYDEBUG ) MESSAGE("Update CORBA");
899     SMESH::log_array_var aSeq = myMeshServer->GetLog( theIsClear );
900     CORBA::Long aLength = aSeq->length();
901     anIsModified = aLength > 0;
902     if ( MYDEBUG ) MESSAGE( "Update: length of the script is "<<aLength );
903
904     if ( !anIsModified )
905       return false;
906
907     // update client mesh structure by logged changes commands
908     try
909     {
910       for ( CORBA::Long anId = 0; anId < aLength; anId++)
911       {
912         const SMESH::double_array& aCoords = aSeq[anId].coords;
913         const SMESH::long_array& anIndexes = aSeq[anId].indexes;
914         CORBA::Long anElemId = 0, aNbElems = aSeq[anId].number;
915         CORBA::Long aCommand = aSeq[anId].commandType;
916
917         switch(aCommand)
918         {
919         case SMESH::ADD_NODE             : AddNodesWithID      ( mySMDSMesh, aSeq, anId ); break;
920         case SMESH::ADD_EDGE             : AddEdgesWithID      ( mySMDSMesh, aSeq, anId ); break;
921         case SMESH::ADD_TRIANGLE         : AddTriasWithID      ( mySMDSMesh, aSeq, anId ); break;
922         case SMESH::ADD_QUADRANGLE       : AddQuadsWithID      ( mySMDSMesh, aSeq, anId ); break;
923         case SMESH::ADD_POLYGON          : AddPolygonsWithID   ( mySMDSMesh, aSeq, anId ); break;
924         case SMESH::ADD_TETRAHEDRON      : AddTetrasWithID     ( mySMDSMesh, aSeq, anId ); break;
925         case SMESH::ADD_PYRAMID          : AddPiramidsWithID   ( mySMDSMesh, aSeq, anId ); break;
926         case SMESH::ADD_PRISM            : AddPrismsWithID     ( mySMDSMesh, aSeq, anId ); break;
927         case SMESH::ADD_HEXAHEDRON       : AddHexasWithID      ( mySMDSMesh, aSeq, anId ); break;
928         case SMESH::ADD_HEXAGONAL_PRISM  : AddHexPrismWithID   ( mySMDSMesh, aSeq, anId ); break;
929         case SMESH::ADD_POLYHEDRON       : AddPolyhedronsWithID( mySMDSMesh, aSeq, anId ); break;
930         case SMESH::ADD_ELEM0D           : Add0DElementsWithID ( mySMDSMesh, aSeq, anId ); break;
931         case SMESH::ADD_BALL             : AddBallsWithID      ( mySMDSMesh, aSeq, anId ); break;
932
933         case SMESH::ADD_QUADEDGE         : AddQuadEdgesWithID   ( mySMDSMesh, aSeq, anId ); break;
934         case SMESH::ADD_QUADTRIANGLE     : AddQuadTriasWithID   ( mySMDSMesh, aSeq, anId ); break;
935         case SMESH::ADD_QUADQUADRANGLE   : AddQuadQuadsWithID   ( mySMDSMesh, aSeq, anId ); break;
936         case SMESH::ADD_QUADPOLYGON      : AddQuadPolygonsWithID( mySMDSMesh, aSeq, anId ); break;
937         case SMESH::ADD_QUADTETRAHEDRON  : AddQuadTetrasWithID  ( mySMDSMesh, aSeq, anId ); break;
938         case SMESH::ADD_QUADPYRAMID      : AddQuadPiramidsWithID( mySMDSMesh, aSeq, anId ); break;
939         case SMESH::ADD_QUADPENTAHEDRON  : AddQuadPentasWithID  ( mySMDSMesh, aSeq, anId ); break;
940         case SMESH::ADD_QUADHEXAHEDRON   : AddQuadHexasWithID   ( mySMDSMesh, aSeq, anId ); break;
941
942         case SMESH::ADD_BIQUAD_QUADRANGLE: AddBiQuadQuadsWithID ( mySMDSMesh, aSeq, anId ); break;
943         case SMESH::ADD_BIQUAD_TRIANGLE  : AddBiQuadTriasWithID ( mySMDSMesh, aSeq, anId ); break;
944         case SMESH::ADD_TRIQUAD_HEXA     : AddTriQuadHexasWithID( mySMDSMesh, aSeq, anId ); break;
945
946         case SMESH::CLEAR_MESH:
947           mySMDSMesh->Clear();
948           break;
949
950         case SMESH::REMOVE_NODE:
951           for( ; anElemId < aNbElems; anElemId++ )
952             mySMDSMesh->RemoveNode( FindNode( mySMDSMesh, anIndexes[anElemId] ) );
953         break;
954
955         case SMESH::REMOVE_ELEMENT:
956           for( ; anElemId < aNbElems; anElemId++ )
957             mySMDSMesh->RemoveElement( FindElement( mySMDSMesh, anIndexes[anElemId] ) );
958         break;
959
960         case SMESH::MOVE_NODE:
961           for(CORBA::Long aCoordId=0; anElemId < aNbElems; anElemId++, aCoordId+=3)
962           {
963             SMDS_MeshNode* node =
964               const_cast<SMDS_MeshNode*>( FindNode( mySMDSMesh, anIndexes[anElemId] ));
965             node->setXYZ( aCoords[aCoordId], aCoords[aCoordId+1], aCoords[aCoordId+2] );
966           }
967         break;
968
969         case SMESH::CHANGE_ELEMENT_NODES:
970           for ( CORBA::Long i = 0; anElemId < aNbElems; anElemId++ )
971           {
972             // find element
973             const SMDS_MeshElement* elem = FindElement( mySMDSMesh, anIndexes[i++] );
974             // nb nodes
975             int nbNodes = anIndexes[i++];
976             // nodes
977             //ASSERT( nbNodes < 9 );
978             std::vector<const SMDS_MeshNode*> aNodes( nbNodes );
979             for ( int iNode = 0; iNode < nbNodes; iNode++ )
980               aNodes[ iNode ] = FindNode( mySMDSMesh, anIndexes[i++] );
981             // change
982             mySMDSMesh->ChangeElementNodes( elem, &aNodes[0], nbNodes );
983           }
984           break;
985
986         case SMESH::CHANGE_POLYHEDRON_NODES:
987           ChangePolyhedronNodes(mySMDSMesh, aSeq, anId);
988           break;
989         case SMESH::RENUMBER:
990           for(CORBA::Long i=0; anElemId < aNbElems; anElemId++, i+=3)
991           {
992             mySMDSMesh->Renumber( anIndexes[i], anIndexes[i+1], anIndexes[i+2] );
993           }
994           break;
995
996         default:;
997         }
998       }
999     }
1000     catch ( SALOME::SALOME_Exception& exc )
1001     {
1002       INFOS("Following exception was cought:\n\t"<<exc.details.text);
1003     }
1004     catch( const std::exception& exc)
1005     {
1006       INFOS("Following exception was cought:\n\t"<<exc.what());
1007     }
1008     catch(...)
1009     {
1010       INFOS("Unknown exception was cought !!!");
1011     }
1012
1013     if ( MYDEBUG && mySMDSMesh )
1014     {
1015       MESSAGE("Update - mySMDSMesh->NbNodes() = "<<mySMDSMesh->NbNodes());
1016       MESSAGE("Update - mySMDSMesh->Nb0DElements() = "<<mySMDSMesh->Nb0DElements());
1017       MESSAGE("Update - mySMDSMesh->NbEdges() = "<<mySMDSMesh->NbEdges());
1018       MESSAGE("Update - mySMDSMesh->NbFaces() = "<<mySMDSMesh->NbFaces());
1019       MESSAGE("Update - mySMDSMesh->NbVolumes() = "<<mySMDSMesh->NbVolumes());
1020     }
1021   } // end of update mesh by log script
1022
1023   return anIsModified;
1024 }