Salome HOME
IPAL52867: Crash at attempt to create group by 'Color of Group' filter
[modules/smesh.git] / src / SMESH_I / SMESH_Filter_i.cxx
1 // Copyright (C) 2007-2015  CEA/DEN, EDF R&D, OPEN CASCADE
2 //
3 // Copyright (C) 2003-2007  OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
4 // CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
5 //
6 // This library is free software; you can redistribute it and/or
7 // modify it under the terms of the GNU Lesser General Public
8 // License as published by the Free Software Foundation; either
9 // version 2.1 of the License, or (at your option) any later version.
10 //
11 // This library is distributed in the hope that it will be useful,
12 // but WITHOUT ANY WARRANTY; without even the implied warranty of
13 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
14 // Lesser General Public License for more details.
15 //
16 // You should have received a copy of the GNU Lesser General Public
17 // License along with this library; if not, write to the Free Software
18 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
19 //
20 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
21 //
22
23 //  SMESH SMESH_I : idl implementation based on 'SMESH' unit's calsses
24 //  File   : SMESH_Filter_i.cxx
25 //  Author : Alexey Petrov, OCC
26 //  Module : SMESH
27
28 #include "SMESH_Filter_i.hxx"
29
30 #include "SMDS_ElemIterator.hxx"
31 #include "SMDS_Mesh.hxx"
32 #include "SMDS_MeshElement.hxx"
33 #include "SMDS_MeshNode.hxx"
34 #include "SMESHDS_GroupBase.hxx"
35 #include "SMESHDS_Mesh.hxx"
36 #include "SMESH_Gen_i.hxx"
37 #include "SMESH_Group_i.hxx"
38 #include "SMESH_PythonDump.hxx"
39
40 #include <SALOMEDS_wrap.hxx>
41 #include <GEOM_wrap.hxx>
42
43 #include <BRep_Tool.hxx>
44 #include <Geom_CylindricalSurface.hxx>
45 #include <Geom_Plane.hxx>
46 #include <LDOMParser.hxx>
47 #include <LDOMString.hxx>
48 #include <LDOM_Document.hxx>
49 #include <LDOM_Element.hxx>
50 #include <LDOM_Node.hxx>
51 #include <LDOM_XmlWriter.hxx>
52 #include <Precision.hxx>
53 #include <TColStd_ListIteratorOfListOfInteger.hxx>
54 #include <TColStd_ListIteratorOfListOfReal.hxx>
55 #include <TColStd_ListOfInteger.hxx>
56 #include <TColStd_ListOfReal.hxx>
57 #include <TColStd_SequenceOfHAsciiString.hxx>
58 #include <TCollection_HAsciiString.hxx>
59 #include <TopExp_Explorer.hxx>
60 #include <TopoDS.hxx>
61 #include <TopoDS_Shape.hxx>
62
63 using namespace SMESH;
64 using namespace SMESH::Controls;
65
66
67 namespace SMESH
68 {
69   Predicate_i*
70   GetPredicate( Predicate_ptr thePredicate )
71   {
72     return DownCast<Predicate_i*>(thePredicate);
73   }
74 }
75
76 /*
77                             AUXILIARY METHODS
78 */
79
80 inline
81 const SMDS_Mesh*
82 MeshPtr2SMDSMesh( SMESH_Mesh_ptr theMesh )
83 {
84   SMESH_Mesh_i* anImplPtr = DownCast<SMESH_Mesh_i*>(theMesh);
85   return anImplPtr ? anImplPtr->GetImpl().GetMeshDS() : 0;
86 }
87
88 inline
89 SMESH::long_array*
90 toArray( const TColStd_ListOfInteger& aList )
91 {
92   SMESH::long_array_var anArray = new SMESH::long_array;
93   anArray->length( aList.Extent() );
94   TColStd_ListIteratorOfListOfInteger anIter( aList );
95   int i = 0;
96   for( ; anIter.More(); anIter.Next() )
97     anArray[ i++ ] = anIter.Value();
98
99   return anArray._retn();
100 }
101
102 inline
103 SMESH::double_array*
104 toArray( const TColStd_ListOfReal& aList )
105 {
106   SMESH::double_array_var anArray = new SMESH::double_array;
107   anArray->length( aList.Extent() );
108   TColStd_ListIteratorOfListOfReal anIter( aList );
109   int i = 0;
110   for( ; anIter.More(); anIter.Next() )
111     anArray[ i++ ] = anIter.Value();
112
113   return anArray._retn();
114 }
115
116 static SMESH::Filter::Criterion createCriterion()
117 {
118   SMESH::Filter::Criterion aCriterion;
119
120   aCriterion.Type          = FT_Undefined;
121   aCriterion.Compare       = FT_Undefined;
122   aCriterion.Threshold     = 0;
123   aCriterion.UnaryOp       = FT_Undefined;
124   aCriterion.BinaryOp      = FT_Undefined;
125   aCriterion.ThresholdStr  = "";
126   aCriterion.ThresholdID   = "";
127   aCriterion.Tolerance     = Precision::Confusion();
128   aCriterion.TypeOfElement = SMESH::ALL;
129   aCriterion.Precision     = -1;
130
131   return aCriterion;
132 }
133
134 static TopoDS_Shape getShapeByName( const char* theName )
135 {
136   if ( theName != 0 )
137   {
138     SMESH_Gen_i* aSMESHGen     = SMESH_Gen_i::GetSMESHGen();
139     SALOMEDS::Study_var aStudy = aSMESHGen->GetCurrentStudy();
140     if ( !aStudy->_is_nil() )
141     {
142       SALOMEDS::Study::ListOfSObject_var aList = aStudy->FindObjectByName( theName, "GEOM" );
143       if ( aList->length() > 0 )
144       {
145         CORBA::Object_var        anObj = aList[ 0 ]->GetObject();
146         GEOM::GEOM_Object_var aGeomObj = GEOM::GEOM_Object::_narrow( anObj );
147         TopoDS_Shape             shape = aSMESHGen->GeomObjectToShape( aGeomObj );
148         SALOME::UnRegister( aList ); // UnRegister() objects in aList
149         return shape;
150       }
151     }
152   }
153   return TopoDS_Shape();
154 }
155
156 static TopoDS_Shape getShapeByID (const char* theID)
157 {
158   if ( theID && strlen( theID ) > 0 ) {
159     SMESH_Gen_i*     aSMESHGen = SMESH_Gen_i::GetSMESHGen();
160     SALOMEDS::Study_var aStudy = aSMESHGen->GetCurrentStudy();
161     if ( !aStudy->_is_nil() ) {
162       SALOMEDS::SObject_wrap aSObj = aStudy->FindObjectID(theID);
163       if ( !aSObj->_is_nil() ) {
164         CORBA::Object_var          obj = aSObj->GetObject();
165         GEOM::GEOM_Object_var aGeomObj = GEOM::GEOM_Object::_narrow(obj);
166         return aSMESHGen->GeomObjectToShape( aGeomObj );
167       }
168     }
169   }
170   return TopoDS_Shape();
171 }
172
173 static std::string getShapeNameByID (const char* theID)
174 {
175   if ( theID && strlen( theID ) > 0 ) {
176     SMESH_Gen_i*     aSMESHGen = SMESH_Gen_i::GetSMESHGen();
177     SALOMEDS::Study_var aStudy = aSMESHGen->GetCurrentStudy();
178     if ( !aStudy->_is_nil() ) {
179       SALOMEDS::SObject_wrap aSObj = aStudy->FindObjectID(theID);
180       if ( !aSObj->_is_nil() ) {
181         CORBA::String_var name = aSObj->GetName();
182         return name.in();
183       }
184     }
185   }
186   return "";
187 }
188
189 /*
190                                 FUNCTORS
191 */
192
193 /*
194   Class       : Functor_i
195   Description : An abstact class for all functors
196 */
197 Functor_i::Functor_i():
198   SALOME::GenericObj_i( SMESH_Gen_i::GetPOA() )
199 {
200   //Base class Salome_GenericObject do it inmplicitly by overriding PortableServer::POA_ptr _default_POA() method  
201   //PortableServer::ObjectId_var anObjectId =
202   //  SMESH_Gen_i::GetPOA()->activate_object( this );
203 }
204
205 Functor_i::~Functor_i()
206 {
207   //TPythonDump()<<this<<".UnRegister()";
208 }
209
210 void Functor_i::SetMesh( SMESH_Mesh_ptr theMesh )
211 {
212   myFunctorPtr->SetMesh( MeshPtr2SMDSMesh( theMesh ) );
213   TPythonDump()<<this<<".SetMesh("<<theMesh<<")";
214 }
215
216 ElementType Functor_i::GetElementType()
217 {
218   return ( ElementType )myFunctorPtr->GetType();
219 }
220
221
222 /*
223   Class       : NumericalFunctor_i
224   Description : Base class for numerical functors
225 */
226 CORBA::Double NumericalFunctor_i::GetValue( CORBA::Long theId )
227 {
228   return myNumericalFunctorPtr->GetValue( theId );
229 }
230
231 SMESH::Histogram* NumericalFunctor_i::GetHistogram(CORBA::Short nbIntervals, CORBA::Boolean isLogarithmic)
232 {
233   std::vector<int> nbEvents;
234   std::vector<double> funValues;
235   std::vector<int> elements;
236   myNumericalFunctorPtr->GetHistogram(nbIntervals,nbEvents,funValues,elements,0,isLogarithmic);
237
238   SMESH::Histogram_var histogram = new SMESH::Histogram;
239
240   nbIntervals = CORBA::Short( Min( int( nbEvents.size()),
241                                    int( funValues.size() - 1 )));
242   if ( nbIntervals > 0 )
243   {
244     histogram->length( nbIntervals );
245     for ( int i = 0; i < nbIntervals; ++i )
246     {
247       HistogramRectangle& rect = histogram[i];
248       rect.nbEvents = nbEvents[i];
249       rect.min = funValues[i];
250       rect.max = funValues[i+1];
251     }
252   }
253   return histogram._retn();
254 }
255
256 SMESH::Histogram* NumericalFunctor_i::GetLocalHistogram(CORBA::Short              nbIntervals,
257                                                         CORBA::Boolean            isLogarithmic,
258                                                         SMESH::SMESH_IDSource_ptr object)
259 {
260   SMESH::Histogram_var histogram = new SMESH::Histogram;
261
262   std::vector<int>    nbEvents;
263   std::vector<double> funValues;
264   std::vector<int>    elements;
265
266   SMDS_ElemIteratorPtr elemIt;
267   if ( SMESH::DownCast< SMESH_GroupOnFilter_i* >( object ) ||
268        SMESH::DownCast< SMESH::Filter_i* >( object ))
269   {
270     elemIt = SMESH_Mesh_i::GetElements( object, GetElementType() );
271     if ( !elemIt ) return histogram._retn();
272   }
273   else
274   {
275     SMESH::SMESH_Mesh_var        mesh = object->GetMesh();
276     SMESH::long_array_var  objNbElems = object->GetNbElementsByType();
277     SMESH::long_array_var meshNbElems = mesh->  GetNbElementsByType();
278     if ( meshNbElems[ GetElementType() ] !=
279          objNbElems [ GetElementType() ] )
280     {
281       elements.reserve( objNbElems[ GetElementType() ]);
282       elemIt = SMESH_Mesh_i::GetElements( object, GetElementType() );
283     }
284   }
285   if ( elemIt )
286   {
287     while ( elemIt->more() )
288       elements.push_back( elemIt->next()->GetID() );
289     if ( elements.empty() ) return histogram._retn();
290   }
291
292   myNumericalFunctorPtr->GetHistogram(nbIntervals,nbEvents,funValues,elements,0,isLogarithmic);
293
294   nbIntervals = CORBA::Short( Min( int( nbEvents.size()),
295                                    int( funValues.size() - 1 )));
296   if ( nbIntervals > 0 )
297   {
298     histogram->length( nbIntervals );
299     for ( int i = 0; i < nbIntervals; ++i )
300     {
301       HistogramRectangle& rect = histogram[i];
302       rect.nbEvents = nbEvents[i];
303       rect.min = funValues[i];
304       rect.max = funValues[i+1];
305     }
306   }
307   return histogram._retn();
308 }
309
310 void NumericalFunctor_i::SetPrecision( CORBA::Long thePrecision )
311 {
312   myNumericalFunctorPtr->SetPrecision( thePrecision );
313   TPythonDump()<<this<<".SetPrecision("<<thePrecision<<")";
314 }
315
316 CORBA::Long NumericalFunctor_i::GetPrecision()
317 {
318  return myNumericalFunctorPtr->GetPrecision();
319 }
320
321 Controls::NumericalFunctorPtr NumericalFunctor_i::GetNumericalFunctor()
322 {
323   return myNumericalFunctorPtr;
324 }
325
326
327 /*
328   Class       : SMESH_MinimumAngle
329   Description : Functor for calculation of minimum angle
330 */
331 MinimumAngle_i::MinimumAngle_i()
332 {
333   myNumericalFunctorPtr.reset( new Controls::MinimumAngle() );
334   myFunctorPtr = myNumericalFunctorPtr;
335 }
336
337 FunctorType MinimumAngle_i::GetFunctorType()
338 {
339   return SMESH::FT_MinimumAngle;
340 }
341
342
343 /*
344   Class       : AspectRatio
345   Description : Functor for calculating aspect ratio
346 */
347 AspectRatio_i::AspectRatio_i()
348 {
349   myNumericalFunctorPtr.reset( new Controls::AspectRatio() );
350   myFunctorPtr = myNumericalFunctorPtr;
351 }
352
353 FunctorType AspectRatio_i::GetFunctorType()
354 {
355   return SMESH::FT_AspectRatio;
356 }
357
358
359 /*
360   Class       : AspectRatio3D
361   Description : Functor for calculating aspect ratio 3D
362 */
363 AspectRatio3D_i::AspectRatio3D_i()
364 {
365   myNumericalFunctorPtr.reset( new Controls::AspectRatio3D() );
366   myFunctorPtr = myNumericalFunctorPtr;
367 }
368
369 FunctorType AspectRatio3D_i::GetFunctorType()
370 {
371   return SMESH::FT_AspectRatio3D;
372 }
373
374
375 /*
376   Class       : Warping_i
377   Description : Functor for calculating warping
378 */
379 Warping_i::Warping_i()
380 {
381   myNumericalFunctorPtr.reset( new Controls::Warping() );
382   myFunctorPtr = myNumericalFunctorPtr;
383 }
384
385 FunctorType Warping_i::GetFunctorType()
386 {
387   return SMESH::FT_Warping;
388 }
389
390
391 /*
392   Class       : Taper_i
393   Description : Functor for calculating taper
394 */
395 Taper_i::Taper_i()
396 {
397   myNumericalFunctorPtr.reset( new Controls::Taper() );
398   myFunctorPtr = myNumericalFunctorPtr;
399 }
400
401 FunctorType Taper_i::GetFunctorType()
402 {
403   return SMESH::FT_Taper;
404 }
405
406 /*
407   Class       : Skew_i
408   Description : Functor for calculating skew in degrees
409 */
410 Skew_i::Skew_i()
411 {
412   myNumericalFunctorPtr.reset( new Controls::Skew() );
413   myFunctorPtr = myNumericalFunctorPtr;
414 }
415
416 FunctorType Skew_i::GetFunctorType()
417 {
418   return SMESH::FT_Skew;
419 }
420
421 /*
422   Class       : Area_i
423   Description : Functor for calculating area
424 */
425 Area_i::Area_i()
426 {
427   myNumericalFunctorPtr.reset( new Controls::Area() );
428   myFunctorPtr = myNumericalFunctorPtr;
429 }
430
431 FunctorType Area_i::GetFunctorType()
432 {
433   return SMESH::FT_Area;
434 }
435
436 /*
437   Class       : Volume3D_i
438   Description : Functor for calculating volume of 3D element
439 */
440 Volume3D_i::Volume3D_i()
441 {
442   myNumericalFunctorPtr.reset( new Controls::Volume() );
443   myFunctorPtr = myNumericalFunctorPtr;
444 }
445
446 FunctorType Volume3D_i::GetFunctorType()
447 {
448   return SMESH::FT_Volume3D;
449 }
450
451 /*
452   Class       : MaxElementLength2D_i
453   Description : Functor for calculating maximum length of 2D element
454 */
455 MaxElementLength2D_i::MaxElementLength2D_i()
456 {
457   myNumericalFunctorPtr.reset( new Controls::MaxElementLength2D() );
458   myFunctorPtr = myNumericalFunctorPtr;
459 }
460
461 FunctorType MaxElementLength2D_i::GetFunctorType()
462 {
463   return SMESH::FT_MaxElementLength2D;
464 }
465
466 /*
467   Class       : MaxElementLength3D_i
468   Description : Functor for calculating maximum length of 3D element
469 */
470 MaxElementLength3D_i::MaxElementLength3D_i()
471 {
472   myNumericalFunctorPtr.reset( new Controls::MaxElementLength3D() );
473   myFunctorPtr = myNumericalFunctorPtr;
474 }
475
476 FunctorType MaxElementLength3D_i::GetFunctorType()
477 {
478   return SMESH::FT_MaxElementLength3D;
479 }
480
481 /*
482   Class       : Length_i
483   Description : Functor for calculating length off edge
484 */
485 Length_i::Length_i()
486 {
487   myNumericalFunctorPtr.reset( new Controls::Length() );
488   myFunctorPtr = myNumericalFunctorPtr;
489 }
490
491 FunctorType Length_i::GetFunctorType()
492 {
493   return SMESH::FT_Length;
494 }
495
496 /*
497   Class       : Length2D_i
498   Description : Functor for calculating length of edge
499 */
500 Length2D_i::Length2D_i()
501 {
502   myNumericalFunctorPtr.reset( new Controls::Length2D() );
503   myFunctorPtr = myNumericalFunctorPtr;
504 }
505
506 FunctorType Length2D_i::GetFunctorType()
507 {
508   return SMESH::FT_Length2D;
509 }
510
511 SMESH::Length2D::Values* Length2D_i::GetValues()
512 {
513   INFOS("Length2D_i::GetValues");
514   SMESH::Controls::Length2D::TValues aValues;
515   (dynamic_cast<SMESH::Controls::Length2D*>(myFunctorPtr.get()))->GetValues( aValues );
516
517   long i = 0, iEnd = aValues.size();
518
519   SMESH::Length2D::Values_var aResult = new SMESH::Length2D::Values(iEnd);
520   aResult->length(iEnd);
521
522   SMESH::Controls::Length2D::TValues::const_iterator anIter;
523   for ( anIter = aValues.begin() ; anIter != aValues.end(); anIter++, i++ )
524   {
525     const SMESH::Controls::Length2D::Value&  aVal = *anIter;
526     SMESH::Length2D::Value &aValue = aResult[ i ];
527
528     aValue.myLength = aVal.myLength;
529     aValue.myPnt1 = aVal.myPntId[ 0 ];
530     aValue.myPnt2 = aVal.myPntId[ 1 ];
531   }
532
533   INFOS("Length2D_i::GetValuess~");
534   return aResult._retn();
535 }
536
537 /*
538   Class       : MultiConnection_i
539   Description : Functor for calculating number of faces conneted to the edge
540 */
541 MultiConnection_i::MultiConnection_i()
542 {
543   myNumericalFunctorPtr.reset( new Controls::MultiConnection() );
544   myFunctorPtr = myNumericalFunctorPtr;
545 }
546
547 FunctorType MultiConnection_i::GetFunctorType()
548 {
549   return SMESH::FT_MultiConnection;
550 }
551
552 /*
553   Class       : BallDiameter_i
554   Description : Functor returning diameter of a ball element
555 */
556 BallDiameter_i::BallDiameter_i()
557 {
558   myNumericalFunctorPtr.reset( new Controls::BallDiameter() );
559   myFunctorPtr = myNumericalFunctorPtr;
560 }
561
562 FunctorType BallDiameter_i::GetFunctorType()
563 {
564   return SMESH::FT_BallDiameter;
565 }
566
567 /*
568   Class       : MultiConnection2D_i
569   Description : Functor for calculating number of faces conneted to the edge
570 */
571 MultiConnection2D_i::MultiConnection2D_i()
572 {
573   myNumericalFunctorPtr.reset( new Controls::MultiConnection2D() );
574   myFunctorPtr = myNumericalFunctorPtr;
575 }
576
577 FunctorType MultiConnection2D_i::GetFunctorType()
578 {
579   return SMESH::FT_MultiConnection2D;
580 }
581
582 SMESH::MultiConnection2D::Values* MultiConnection2D_i::GetValues()
583 {
584   INFOS("MultiConnection2D_i::GetValues");
585   SMESH::Controls::MultiConnection2D::MValues aValues;
586   (dynamic_cast<SMESH::Controls::MultiConnection2D*>(myFunctorPtr.get()))->GetValues( aValues );
587   
588   long i = 0, iEnd = aValues.size();
589
590   SMESH::MultiConnection2D::Values_var aResult = new SMESH::MultiConnection2D::Values(iEnd);
591   aResult->length(iEnd);
592
593   SMESH::Controls::MultiConnection2D::MValues::const_iterator anIter;
594   for ( anIter = aValues.begin() ; anIter != aValues.end(); anIter++, i++ )
595   {
596     const SMESH::Controls::MultiConnection2D::Value&  aVal = (*anIter).first;
597     SMESH::MultiConnection2D::Value &aValue = aResult[ i ];
598
599     aValue.myPnt1 = aVal.myPntId[ 0 ];
600     aValue.myPnt2 = aVal.myPntId[ 1 ];
601     aValue.myNbConnects = (*anIter).second;
602   }
603
604   INFOS("Multiconnection2D_i::GetValuess~");
605   return aResult._retn();
606 }
607
608 /*
609                             PREDICATES
610 */
611
612
613 /*
614   Class       : Predicate_i
615   Description : Base class for all predicates
616 */
617 CORBA::Boolean Predicate_i::IsSatisfy( CORBA::Long theId )
618 {
619   return myPredicatePtr->IsSatisfy( theId );
620 }
621
622 CORBA::Long Predicate_i::NbSatisfying( SMESH::SMESH_IDSource_ptr obj )
623 {
624   SMESH::SMESH_Mesh_var meshVar = obj->GetMesh();
625   const SMDS_Mesh*       meshDS = MeshPtr2SMDSMesh( meshVar );
626   if ( !meshDS )
627     return 0;
628   myPredicatePtr->SetMesh( meshDS );
629
630   SMDSAbs_ElementType elemType = SMDSAbs_ElementType( GetElementType() );
631
632   int nb = 0;
633   SMDS_ElemIteratorPtr elemIt =
634     SMESH::DownCast<SMESH_Mesh_i*>( meshVar )->GetElements( obj, GetElementType() );
635   if ( elemIt )
636     while ( elemIt->more() )
637     {
638       const SMDS_MeshElement* e = elemIt->next();
639       if ( e && e->GetType() == elemType )
640         nb += myPredicatePtr->IsSatisfy( e->GetID() );
641     }
642   return nb;
643 }
644
645 Controls::PredicatePtr Predicate_i::GetPredicate()
646 {
647   return myPredicatePtr;
648 }
649
650 /*
651   Class       : BadOrientedVolume_i
652   Description : Verify whether a mesh volume is incorrectly oriented from
653                 the point of view of MED convention
654 */
655 BadOrientedVolume_i::BadOrientedVolume_i()
656 {
657   Controls::PredicatePtr control( new Controls::BadOrientedVolume() );
658   myFunctorPtr = myPredicatePtr = control;
659 };
660
661 FunctorType BadOrientedVolume_i::GetFunctorType()
662 {
663   return SMESH::FT_BadOrientedVolume;
664 }
665
666 /*
667   Class       : BareBorderVolume_i
668   Description : Verify whether a mesh volume has a free facet without a face on it
669 */
670 BareBorderVolume_i::BareBorderVolume_i()
671 {
672   Controls::PredicatePtr control( new Controls::BareBorderVolume() );
673   myFunctorPtr = myPredicatePtr = control;
674 };
675
676 FunctorType BareBorderVolume_i::GetFunctorType()
677 {
678   return SMESH::FT_BareBorderVolume;
679 }
680
681 /*
682   Class       : BareBorderFace_i
683   Description : Verify whether a mesh face has a free border without an edge on it
684 */
685 BareBorderFace_i::BareBorderFace_i()
686 {
687   Controls::PredicatePtr control( new Controls::BareBorderFace() );
688   myFunctorPtr = myPredicatePtr = control;
689 };
690
691 FunctorType BareBorderFace_i::GetFunctorType()
692 {
693   return SMESH::FT_BareBorderFace;
694 }
695
696 /*
697   Class       : OverConstrainedVolume_i
698   Description : Verify whether a mesh volume has only one facet shared with other volumes
699 */
700 OverConstrainedVolume_i::OverConstrainedVolume_i()
701 {
702   Controls::PredicatePtr control( new Controls::OverConstrainedVolume() );
703   myFunctorPtr = myPredicatePtr = control;
704 };
705
706 FunctorType OverConstrainedVolume_i::GetFunctorType()
707 {
708   return SMESH::FT_OverConstrainedVolume;
709 }
710
711 /*
712   Class       : OverConstrainedFace_i
713   Description : Verify whether a mesh face has only one border shared with other faces
714 */
715 OverConstrainedFace_i::OverConstrainedFace_i()
716 {
717   Controls::PredicatePtr control( new Controls::OverConstrainedFace() );
718   myFunctorPtr = myPredicatePtr = control;
719 };
720
721 FunctorType OverConstrainedFace_i::GetFunctorType()
722 {
723   return SMESH::FT_OverConstrainedFace;
724 }
725
726 /*
727   Class       : BelongToMeshGroup_i
728   Description : Verify whether a mesh element is included into a mesh group
729 */
730 BelongToMeshGroup_i::BelongToMeshGroup_i()
731 {
732   myBelongToMeshGroup = Controls::BelongToMeshGroupPtr( new Controls::BelongToMeshGroup() );
733   myFunctorPtr = myPredicatePtr = myBelongToMeshGroup;
734 }
735
736 BelongToMeshGroup_i::~BelongToMeshGroup_i()
737 {
738   SetGroup( SMESH::SMESH_GroupBase::_nil() );
739 }
740
741 void BelongToMeshGroup_i::SetGroup( SMESH::SMESH_GroupBase_ptr theGroup )
742 {
743   if ( myGroup->_is_equivalent( theGroup ))
744     return;
745
746   if ( ! myGroup->_is_nil() )
747     myGroup->UnRegister();
748
749   myGroup = SMESH_GroupBase::_duplicate( theGroup );
750
751   myBelongToMeshGroup->SetGroup( 0 );
752   if ( SMESH_GroupBase_i* gr_i = SMESH::DownCast< SMESH_GroupBase_i* >( myGroup ))
753   {
754     myBelongToMeshGroup->SetGroup( gr_i->GetGroupDS() );
755     myGroup->Register();
756   }
757 }
758
759 void BelongToMeshGroup_i::SetGroupID( const char* theID ) // IOR or StoreName
760 {
761   myID = theID;
762   if ( strncmp( "IOR:", myID.c_str(), 4 ) == 0 ) // transient mode, no GUI
763   {
764     CORBA::Object_var obj = SMESH_Gen_i::GetORB()->string_to_object( myID.c_str() );
765     SetGroup( SMESH::SMESH_GroupBase::_narrow( obj ));
766   }
767   else if ( strncmp( "0:", myID.c_str(), 2 ) == 0 ) // transient mode + GUI
768   {
769     SMESH_Gen_i* aSMESHGen     = SMESH_Gen_i::GetSMESHGen();
770     SALOMEDS::Study_var aStudy = aSMESHGen->GetCurrentStudy();
771     if ( !aStudy->_is_nil() ) {
772       SALOMEDS::SObject_wrap aSObj = aStudy->FindObjectID( myID.c_str() );
773       if ( !aSObj->_is_nil() ) {
774         CORBA::Object_var obj = aSObj->GetObject();
775         SetGroup( SMESH::SMESH_GroupBase::_narrow( obj ));
776       }
777     }
778   }
779   else if ( !myID.empty() ) // persistent mode
780   {
781     myBelongToMeshGroup->SetStoreName( myID );
782   }
783 }
784
785 std::string BelongToMeshGroup_i::GetGroupID()
786 {
787   if ( myGroup->_is_nil() )
788     SMESH::SMESH_GroupBase_var( GetGroup() );
789
790   if ( !myGroup->_is_nil() )
791     myID = SMESH_Gen_i::GetORB()->object_to_string( myGroup );
792
793   return myID;
794 }
795
796 SMESH::SMESH_GroupBase_ptr BelongToMeshGroup_i::GetGroup()
797 {
798   if ( myGroup->_is_nil() && myBelongToMeshGroup->GetGroup() )
799   {
800     // search for a group in a current study
801     SMESH_Gen_i* aSMESHGen = SMESH_Gen_i::GetSMESHGen();
802     if ( StudyContext*  sc = aSMESHGen->GetCurrentStudyContext() )
803     {
804       int id = 1;
805       std::string ior;
806       while (true)
807       {
808         ior = sc->getIORbyId( id++ );
809         if ( ior.empty() ) break;
810         CORBA::Object_var obj = aSMESHGen->GetORB()->string_to_object( ior.c_str() );
811         if ( SMESH_GroupBase_i* g_i = SMESH::DownCast<SMESH_GroupBase_i*>( obj ))
812           if ( g_i->GetGroupDS() == myBelongToMeshGroup->GetGroup() )
813           {
814             SetGroup( g_i->_this() );
815             break;
816           }
817       }
818     }
819   }
820   return SMESH::SMESH_GroupBase::_duplicate( myGroup );
821 }
822
823 FunctorType BelongToMeshGroup_i::GetFunctorType()
824 {
825   return SMESH::FT_BelongToMeshGroup;
826 }
827
828 /*
829   Class       : BelongToGeom_i
830   Description : Predicate for selection on geometrical support
831 */
832 BelongToGeom_i::BelongToGeom_i()
833 {
834   myBelongToGeomPtr.reset( new Controls::BelongToGeom() );
835   myFunctorPtr = myPredicatePtr = myBelongToGeomPtr;
836   myShapeName = 0;
837   myShapeID   = 0;
838 }
839
840 BelongToGeom_i::~BelongToGeom_i()
841 {
842   CORBA::string_free( myShapeName );
843   CORBA::string_free( myShapeID );
844 }
845
846 void BelongToGeom_i::SetGeom( GEOM::GEOM_Object_ptr theGeom )
847 {
848   if ( theGeom->_is_nil() )
849     return;
850   SMESH_Gen_i* aSMESHGen = SMESH_Gen_i::GetSMESHGen();
851   GEOM::GEOM_Gen_ptr aGEOMGen = SMESH_Gen_i::GetGeomEngine();
852   TopoDS_Shape aLocShape = aSMESHGen->GetShapeReader()->GetShape( aGEOMGen, theGeom );
853   myBelongToGeomPtr->SetGeom( aLocShape );
854   TPythonDump()<<this<<".SetGeom("<<theGeom<<")";
855 }
856
857 void BelongToGeom_i::SetGeom( const TopoDS_Shape& theShape )
858 {
859   myBelongToGeomPtr->SetGeom( theShape );
860 }
861
862 void BelongToGeom_i::SetElementType(ElementType theType)
863 {
864   myBelongToGeomPtr->SetType(SMDSAbs_ElementType(theType));
865   TPythonDump()<<this<<".SetElementType("<<theType<<")";
866 }
867
868 FunctorType BelongToGeom_i::GetFunctorType()
869 {
870   return SMESH::FT_BelongToGeom;
871 }
872
873 void BelongToGeom_i::SetShapeName( const char* theName )
874 {
875   CORBA::string_free( myShapeName );
876   myShapeName = CORBA::string_dup( theName );
877   myBelongToGeomPtr->SetGeom( getShapeByName( myShapeName ) );
878   TPythonDump()<<this<<".SetShapeName('"<<theName<<"')";
879 }
880
881 void BelongToGeom_i::SetShape( const char* theID, const char* theName )
882 {
883   CORBA::string_free( myShapeName );
884   CORBA::string_free( myShapeID );
885   myShapeName = CORBA::string_dup( theName );
886   myShapeID   = CORBA::string_dup( theID );
887   bool hasName = ( theName && theName[0] );
888   bool hasID   = ( theID   && theID[0] );
889
890   TopoDS_Shape S;
891   if ( hasName && hasID )
892   {
893     S = getShapeByID( myShapeID );
894     if ( S.IsNull() )
895       S = getShapeByName( myShapeName );
896   }
897   else
898   {
899     S = hasID ? getShapeByID( myShapeID ) : getShapeByName( myShapeName );
900   }
901   myBelongToGeomPtr->SetGeom( S );
902 }
903
904 char* BelongToGeom_i::GetShapeName()
905 {
906   return CORBA::string_dup( myShapeName );
907 }
908
909 char* BelongToGeom_i::GetShapeID()
910 {
911   return CORBA::string_dup( myShapeID );
912 }
913
914 void BelongToGeom_i::SetTolerance( CORBA::Double theToler )
915 {
916   myBelongToGeomPtr->SetTolerance( theToler );
917   TPythonDump()<<this<<".SetTolerance("<<theToler<<")";
918 }
919
920 CORBA::Double BelongToGeom_i::GetTolerance()
921 {
922   return myBelongToGeomPtr->GetTolerance();
923 }
924
925 /*
926   Class       : BelongToSurface_i
927   Description : Predicate for selection on geometrical support
928 */
929 BelongToSurface_i::BelongToSurface_i( const Handle(Standard_Type)& theSurfaceType )
930 {
931   myElementsOnSurfacePtr.reset( new Controls::ElementsOnSurface() );
932   myFunctorPtr = myPredicatePtr = myElementsOnSurfacePtr;
933   myShapeName = 0;
934   myShapeID   = 0;
935   mySurfaceType = theSurfaceType;
936 }
937
938 BelongToSurface_i::~BelongToSurface_i()
939 {
940   CORBA::string_free( myShapeName );
941   CORBA::string_free( myShapeID );
942 }
943
944 void BelongToSurface_i::SetSurface( GEOM::GEOM_Object_ptr theGeom, ElementType theType )
945 {
946   if ( theGeom->_is_nil() )
947     return;
948   SMESH_Gen_i* aSMESHGen = SMESH_Gen_i::GetSMESHGen();
949   GEOM::GEOM_Gen_ptr aGEOMGen = SMESH_Gen_i::GetGeomEngine();
950   TopoDS_Shape aLocShape = aSMESHGen->GetShapeReader()->GetShape( aGEOMGen, theGeom );
951
952   if ( aLocShape.ShapeType() == TopAbs_FACE )
953   {
954     Handle(Geom_Surface) aSurf = BRep_Tool::Surface( TopoDS::Face( aLocShape ) );
955     if ( !aSurf.IsNull() && aSurf->DynamicType() == mySurfaceType )
956     {
957       myElementsOnSurfacePtr->SetSurface( aLocShape, (SMDSAbs_ElementType)theType );
958       return;
959     }
960   }
961
962   myElementsOnSurfacePtr->SetSurface( TopoDS_Shape(), (SMDSAbs_ElementType)theType );
963 }
964
965 void BelongToSurface_i::SetShapeName( const char* theName, ElementType theType )
966 {
967   CORBA::string_free( myShapeName );
968   myShapeName = CORBA::string_dup( theName );
969   myElementsOnSurfacePtr->SetSurface( getShapeByName( myShapeName ), (SMDSAbs_ElementType)theType );
970   TPythonDump()<<this<<".SetShapeName('"<<theName<<"',"<<theType<<")";
971 }
972
973 void BelongToSurface_i::SetShape( const char* theID,  const char* theName, ElementType theType )
974 {
975   CORBA::string_free( myShapeName );
976   CORBA::string_free( myShapeID );
977   myShapeName = CORBA::string_dup( theName );
978   myShapeID   = CORBA::string_dup( theID );
979   bool hasName = ( theName && theName[0] );
980   bool hasID   = ( theID   && theID[0] );
981
982   TopoDS_Shape S;
983   if ( hasName && hasID )
984   {
985     S = getShapeByID( myShapeID );
986     if ( S.IsNull() )
987       S = getShapeByName( myShapeName );
988   }
989   else
990   {
991     S = hasID ? getShapeByID( myShapeID ) : getShapeByName( myShapeName );
992   }
993   myElementsOnSurfacePtr->SetSurface( S, (SMDSAbs_ElementType)theType );
994 }
995
996 char* BelongToSurface_i::GetShapeName()
997 {
998   return CORBA::string_dup( myShapeName );
999 }
1000
1001 char* BelongToSurface_i::GetShapeID()
1002 {
1003   return CORBA::string_dup( myShapeID );
1004 }
1005
1006 void BelongToSurface_i::SetTolerance( CORBA::Double theToler )
1007 {
1008   myElementsOnSurfacePtr->SetTolerance( theToler );
1009   TPythonDump()<<this<<".SetTolerance("<<theToler<<")";
1010 }
1011
1012 CORBA::Double BelongToSurface_i::GetTolerance()
1013 {
1014   return myElementsOnSurfacePtr->GetTolerance();
1015 }
1016
1017 void BelongToSurface_i::SetUseBoundaries( CORBA::Boolean theUseBndRestrictions )
1018 {
1019   myElementsOnSurfacePtr->SetUseBoundaries( theUseBndRestrictions );
1020   TPythonDump()<<this<<".SetUseBoundaries( " << theUseBndRestrictions << " )";
1021 }
1022
1023 CORBA::Boolean BelongToSurface_i::GetUseBoundaries()
1024 {
1025   return myElementsOnSurfacePtr->GetUseBoundaries();
1026 }
1027
1028
1029 /*
1030   Class       : BelongToPlane_i
1031   Description : Verify whether mesh element lie in pointed Geom planar object
1032 */
1033
1034 BelongToPlane_i::BelongToPlane_i()
1035 : BelongToSurface_i( STANDARD_TYPE( Geom_Plane ) )
1036 {
1037 }
1038
1039 void BelongToPlane_i::SetPlane( GEOM::GEOM_Object_ptr theGeom, ElementType theType )
1040 {
1041   BelongToSurface_i::SetSurface( theGeom, theType );
1042   TPythonDump()<<this<<".SetPlane("<<theGeom<<","<<theType<<")";
1043 }
1044
1045 FunctorType BelongToPlane_i::GetFunctorType()
1046 {
1047   return FT_BelongToPlane;
1048 }
1049
1050 /*
1051   Class       : BelongToCylinder_i
1052   Description : Verify whether mesh element lie in pointed Geom planar object
1053 */
1054
1055 BelongToCylinder_i::BelongToCylinder_i()
1056 : BelongToSurface_i( STANDARD_TYPE( Geom_CylindricalSurface ) )
1057 {
1058 }
1059
1060 void BelongToCylinder_i::SetCylinder( GEOM::GEOM_Object_ptr theGeom, ElementType theType )
1061 {
1062   BelongToSurface_i::SetSurface( theGeom, theType );
1063   TPythonDump()<<this<<".SetCylinder("<<theGeom<<","<<theType<<")";
1064 }
1065
1066 FunctorType BelongToCylinder_i::GetFunctorType()
1067 {
1068   return FT_BelongToCylinder;
1069 }
1070
1071 /*
1072   Class       : BelongToGenSurface_i
1073   Description : Verify whether mesh element lie in pointed Geom planar object
1074 */
1075
1076 BelongToGenSurface_i::BelongToGenSurface_i()
1077 : BelongToSurface_i( STANDARD_TYPE( Geom_CylindricalSurface ) )
1078 {
1079 }
1080
1081 void BelongToGenSurface_i::SetSurface( GEOM::GEOM_Object_ptr theGeom, ElementType theType )
1082 {
1083   if ( theGeom->_is_nil() )
1084     return;
1085   TopoDS_Shape aLocShape = SMESH_Gen_i::GetSMESHGen()->GeomObjectToShape( theGeom );
1086   if ( !aLocShape.IsNull() && aLocShape.ShapeType() != TopAbs_FACE )
1087     aLocShape.Nullify();
1088   
1089   BelongToSurface_i::myElementsOnSurfacePtr->SetSurface( aLocShape, (SMDSAbs_ElementType)theType );
1090   TPythonDump()<<this<<".SetGenSurface("<<theGeom<<","<<theType<<")";
1091 }
1092
1093 FunctorType BelongToGenSurface_i::GetFunctorType()
1094 {
1095   return FT_BelongToGenSurface;
1096 }
1097
1098 /*
1099   Class       : LyingOnGeom_i
1100   Description : Predicate for selection on geometrical support
1101 */
1102 LyingOnGeom_i::LyingOnGeom_i()
1103 {
1104   myLyingOnGeomPtr.reset( new Controls::LyingOnGeom() );
1105   myFunctorPtr = myPredicatePtr = myLyingOnGeomPtr;
1106   myShapeName = 0;
1107   myShapeID = 0;
1108 }
1109
1110 LyingOnGeom_i::~LyingOnGeom_i()
1111 {
1112   CORBA::string_free( myShapeName );
1113   CORBA::string_free( myShapeID );
1114 }
1115
1116 void LyingOnGeom_i::SetGeom( GEOM::GEOM_Object_ptr theGeom )
1117 {
1118   if ( theGeom->_is_nil() )
1119     return;
1120   SMESH_Gen_i* aSMESHGen = SMESH_Gen_i::GetSMESHGen();
1121   GEOM::GEOM_Gen_ptr aGEOMGen = SMESH_Gen_i::GetGeomEngine();
1122   TopoDS_Shape aLocShape = aSMESHGen->GetShapeReader()->GetShape( aGEOMGen, theGeom );
1123   myLyingOnGeomPtr->SetGeom( aLocShape );
1124   TPythonDump()<<this<<".SetGeom("<<theGeom<<")";
1125 }
1126
1127 void LyingOnGeom_i::SetGeom( const TopoDS_Shape& theShape )
1128 {
1129   myLyingOnGeomPtr->SetGeom( theShape );
1130 }
1131
1132 void LyingOnGeom_i::SetElementType(ElementType theType){
1133   myLyingOnGeomPtr->SetType(SMDSAbs_ElementType(theType));
1134   TPythonDump()<<this<<".SetElementType("<<theType<<")";
1135 }
1136
1137 FunctorType LyingOnGeom_i::GetFunctorType()
1138 {
1139   return SMESH::FT_LyingOnGeom;
1140 }
1141
1142 void LyingOnGeom_i::SetShapeName( const char* theName )
1143 {
1144   CORBA::string_free( myShapeName );
1145   myShapeName = CORBA::string_dup( theName );
1146   myLyingOnGeomPtr->SetGeom( getShapeByName( myShapeName ) );
1147   TPythonDump()<<this<<".SetShapeName('"<<theName<<"')";
1148 }
1149
1150 void LyingOnGeom_i::SetShape( const char* theID, const char* theName )
1151 {
1152   CORBA::string_free( myShapeName );
1153   CORBA::string_free( myShapeID   );
1154   myShapeName = CORBA::string_dup( theName );
1155   myShapeID   = CORBA::string_dup( theID   );
1156   bool hasName = ( theName && theName[0] );
1157   bool hasID   = ( theID   && theID[0]   );
1158
1159   TopoDS_Shape S;
1160   if ( hasName && hasID )
1161   {
1162     S = getShapeByID( myShapeID );
1163     if ( S.IsNull() )
1164       S = getShapeByName( myShapeName );
1165   }
1166   else
1167   {
1168     S = hasID ? getShapeByID( myShapeID ) : getShapeByName( myShapeName );
1169   }
1170   myLyingOnGeomPtr->SetGeom( S );
1171 }
1172
1173 char* LyingOnGeom_i::GetShapeName()
1174 {
1175   return CORBA::string_dup( myShapeName );
1176 }
1177
1178 char* LyingOnGeom_i::GetShapeID()
1179 {
1180   return CORBA::string_dup( myShapeID );
1181 }
1182
1183 void LyingOnGeom_i::SetTolerance( CORBA::Double theToler )
1184 {
1185   myLyingOnGeomPtr->SetTolerance( theToler );
1186   TPythonDump()<<this<<".SetTolerance("<<theToler<<")";
1187 }
1188
1189 CORBA::Double LyingOnGeom_i::GetTolerance()
1190 {
1191   return myLyingOnGeomPtr->GetTolerance();
1192 }
1193
1194 /*
1195   Class       : FreeBorders_i
1196   Description : Predicate for free borders
1197 */
1198 FreeBorders_i::FreeBorders_i()
1199 {
1200   myPredicatePtr.reset(new Controls::FreeBorders());
1201   myFunctorPtr = myPredicatePtr;
1202 }
1203
1204 FunctorType FreeBorders_i::GetFunctorType()
1205 {
1206   return SMESH::FT_FreeBorders;
1207 }
1208
1209 /*
1210   Class       : FreeEdges_i
1211   Description : Predicate for free borders
1212 */
1213 FreeEdges_i::FreeEdges_i()
1214 : myFreeEdgesPtr( new Controls::FreeEdges() )
1215 {
1216   myFunctorPtr = myPredicatePtr = myFreeEdgesPtr;
1217 }
1218
1219 SMESH::FreeEdges::Borders* FreeEdges_i::GetBorders()
1220 {
1221   INFOS("FreeEdges_i::GetBorders");
1222   SMESH::Controls::FreeEdges::TBorders aBorders;
1223   myFreeEdgesPtr->GetBoreders( aBorders );
1224
1225   long i = 0, iEnd = aBorders.size();
1226
1227   SMESH::FreeEdges::Borders_var aResult = new SMESH::FreeEdges::Borders;
1228   aResult->length(iEnd);
1229
1230   SMESH::Controls::FreeEdges::TBorders::const_iterator anIter;
1231   for ( anIter = aBorders.begin() ; anIter != aBorders.end(); anIter++, i++ )
1232   {
1233     const SMESH::Controls::FreeEdges::Border&  aBord = *anIter;
1234     SMESH::FreeEdges::Border &aBorder = aResult[ i ];
1235
1236     aBorder.myElemId = aBord.myElemId;
1237     aBorder.myPnt1 = aBord.myPntId[ 0 ];
1238     aBorder.myPnt2 = aBord.myPntId[ 1 ];
1239   }
1240
1241   INFOS("FreeEdges_i::GetBorders~");
1242   return aResult._retn();
1243 }
1244
1245 FunctorType FreeEdges_i::GetFunctorType()
1246 {
1247   return SMESH::FT_FreeEdges;
1248 }
1249
1250 /*
1251   Class       : FreeFaces_i
1252   Description : Predicate for free faces
1253 */
1254 FreeFaces_i::FreeFaces_i()
1255 {
1256   myPredicatePtr.reset(new Controls::FreeFaces());
1257   myFunctorPtr = myPredicatePtr;
1258 }
1259
1260 FunctorType FreeFaces_i::GetFunctorType()
1261 {
1262   return SMESH::FT_FreeFaces;
1263 }
1264
1265 /*
1266   Class       : FreeNodes_i
1267   Description : Predicate for free nodes
1268 */
1269 FreeNodes_i::FreeNodes_i()
1270 {
1271   myPredicatePtr.reset(new Controls::FreeNodes());
1272   myFunctorPtr = myPredicatePtr;
1273 }
1274
1275 FunctorType FreeNodes_i::GetFunctorType()
1276 {
1277   return SMESH::FT_FreeNodes;
1278 }
1279
1280 /*
1281   Class       : EqualNodes_i
1282   Description : Predicate for Equal nodes
1283 */
1284 EqualNodes_i::EqualNodes_i()
1285 {
1286   myCoincidentNodesPtr.reset(new Controls::CoincidentNodes());
1287   myFunctorPtr = myPredicatePtr = myCoincidentNodesPtr;
1288 }
1289
1290 FunctorType EqualNodes_i::GetFunctorType()
1291 {
1292   return SMESH::FT_EqualNodes;
1293 }
1294
1295 void EqualNodes_i::SetTolerance( double tol )
1296 {
1297   myCoincidentNodesPtr->SetTolerance( tol );
1298 }
1299
1300 double EqualNodes_i::GetTolerance()
1301 {
1302   return myCoincidentNodesPtr->GetTolerance();
1303 }
1304
1305 /*
1306   Class       : EqualEdges_i
1307   Description : Predicate for Equal Edges
1308 */
1309 EqualEdges_i::EqualEdges_i()
1310 {
1311   myPredicatePtr.reset(new Controls::CoincidentElements1D());
1312   myFunctorPtr = myPredicatePtr;
1313 }
1314
1315 FunctorType EqualEdges_i::GetFunctorType()
1316 {
1317   return SMESH::FT_EqualEdges;
1318 }
1319
1320 /*
1321   Class       : EqualFaces_i
1322   Description : Predicate for Equal Faces
1323 */
1324 EqualFaces_i::EqualFaces_i()
1325 {
1326   myPredicatePtr.reset(new Controls::CoincidentElements2D());
1327   myFunctorPtr = myPredicatePtr;
1328 }
1329
1330 FunctorType EqualFaces_i::GetFunctorType()
1331 {
1332   return SMESH::FT_EqualFaces;
1333 }
1334
1335 /*
1336   Class       : EqualVolumes_i
1337   Description : Predicate for Equal Volumes
1338 */
1339 EqualVolumes_i::EqualVolumes_i()
1340 {
1341   myPredicatePtr.reset(new Controls::CoincidentElements3D());
1342   myFunctorPtr = myPredicatePtr;
1343 }
1344
1345 FunctorType EqualVolumes_i::GetFunctorType()
1346 {
1347   return SMESH::FT_EqualVolumes;
1348 }
1349
1350
1351 /*
1352   Class       : RangeOfIds_i
1353   Description : Predicate for Range of Ids.
1354                 Range may be specified with two ways.
1355                 1. Using AddToRange method
1356                 2. With SetRangeStr method. Parameter of this method is a string
1357                    like as "1,2,3,50-60,63,67,70-"
1358 */
1359
1360 RangeOfIds_i::RangeOfIds_i()
1361 {
1362   myRangeOfIdsPtr.reset( new Controls::RangeOfIds() );
1363   myFunctorPtr = myPredicatePtr = myRangeOfIdsPtr;
1364 }
1365
1366 void RangeOfIds_i::SetRange( const SMESH::long_array& theIds )
1367 {
1368   CORBA::Long iEnd = theIds.length();
1369   for ( CORBA::Long i = 0; i < iEnd; i++ )
1370     myRangeOfIdsPtr->AddToRange( theIds[ i ] );
1371   TPythonDump()<<this<<".SetRange("<<theIds<<")";
1372 }
1373
1374 CORBA::Boolean RangeOfIds_i::SetRangeStr( const char* theRange )
1375 {
1376   TPythonDump()<<this<<".SetRangeStr('"<<theRange<<"')";
1377   return myRangeOfIdsPtr->SetRangeStr(
1378     TCollection_AsciiString( (Standard_CString)theRange ) );
1379 }
1380
1381 char* RangeOfIds_i::GetRangeStr()
1382 {
1383   TCollection_AsciiString aStr;
1384   myRangeOfIdsPtr->GetRangeStr( aStr );
1385   return CORBA::string_dup( aStr.ToCString() );
1386 }
1387
1388 void RangeOfIds_i::SetElementType( ElementType theType )
1389 {
1390   myRangeOfIdsPtr->SetType( SMDSAbs_ElementType( theType ) );
1391   TPythonDump()<<this<<".SetElementType("<<theType<<")";
1392 }
1393
1394 FunctorType RangeOfIds_i::GetFunctorType()
1395 {
1396   return SMESH::FT_RangeOfIds;
1397 }
1398
1399 /*
1400   Class       : LinearOrQuadratic_i
1401   Description : Predicate to verify whether a mesh element is linear
1402 */
1403 LinearOrQuadratic_i::LinearOrQuadratic_i()
1404 {
1405   myLinearOrQuadraticPtr.reset(new Controls::LinearOrQuadratic());
1406   myFunctorPtr = myPredicatePtr = myLinearOrQuadraticPtr;
1407 }
1408
1409 void LinearOrQuadratic_i::SetElementType(ElementType theType)
1410 {
1411   myLinearOrQuadraticPtr->SetType(SMDSAbs_ElementType(theType));
1412   TPythonDump()<<this<<".SetElementType("<<theType<<")";
1413 }
1414
1415 FunctorType LinearOrQuadratic_i::GetFunctorType()
1416 {
1417   return SMESH::FT_LinearOrQuadratic;
1418 }
1419
1420 /*
1421   Class       : GroupColor_i
1422   Description : Functor for check color of group to whic mesh element belongs to
1423 */
1424 GroupColor_i::GroupColor_i()
1425 {
1426   myGroupColorPtr.reset(new Controls::GroupColor());
1427   myFunctorPtr = myPredicatePtr = myGroupColorPtr;
1428 }
1429
1430 FunctorType GroupColor_i::GetFunctorType()
1431 {
1432   return SMESH::FT_GroupColor;
1433 }
1434
1435 void GroupColor_i::SetColorStr( const char* theColor )
1436 {
1437   myGroupColorPtr->SetColorStr(
1438     TCollection_AsciiString( (Standard_CString)theColor ) );
1439   TPythonDump()<<this<<".SetColorStr('"<<theColor<<"')";
1440 }
1441
1442 char* GroupColor_i::GetColorStr()
1443 {
1444   TCollection_AsciiString aStr;
1445   myGroupColorPtr->GetColorStr( aStr );
1446   return CORBA::string_dup( aStr.ToCString() );
1447 }
1448
1449 void GroupColor_i::SetElementType(ElementType theType)
1450 {
1451   myGroupColorPtr->SetType(SMDSAbs_ElementType(theType));
1452   TPythonDump()<<this<<".SetElementType("<<theType<<")";
1453 }
1454
1455 /*
1456   Class       : ElemGeomType_i
1457   Description : Predicate check is element has indicated geometry type
1458 */
1459 ElemGeomType_i::ElemGeomType_i()
1460 {
1461   myElemGeomTypePtr.reset(new Controls::ElemGeomType());
1462   myFunctorPtr = myPredicatePtr = myElemGeomTypePtr;
1463 }
1464
1465 void ElemGeomType_i::SetElementType(ElementType theType)
1466 {
1467   myElemGeomTypePtr->SetType(SMDSAbs_ElementType(theType));
1468   TPythonDump()<<this<<".SetElementType("<<theType<<")";
1469 }
1470
1471 void ElemGeomType_i::SetGeometryType(GeometryType theType)
1472 {
1473   myElemGeomTypePtr->SetGeomType(SMDSAbs_GeometryType(theType));
1474   TPythonDump()<<this<<".SetGeometryType("<<theType<<")";
1475 }
1476
1477 GeometryType ElemGeomType_i::GetGeometryType() const
1478 {
1479   return (GeometryType)myElemGeomTypePtr->GetGeomType();
1480 }
1481
1482 FunctorType ElemGeomType_i::GetFunctorType()
1483 {
1484   return SMESH::FT_ElemGeomType;
1485 }
1486
1487 /*
1488   Class       : ElemEntityType_i
1489   Description : Predicate check is element has indicated entity type
1490 */
1491 ElemEntityType_i::ElemEntityType_i()
1492 {
1493   myElemEntityTypePtr.reset(new Controls::ElemEntityType());
1494   myFunctorPtr = myPredicatePtr = myElemEntityTypePtr;
1495 }
1496
1497 void ElemEntityType_i::SetElementType(ElementType theType)
1498 {
1499   myElemEntityTypePtr->SetType(SMDSAbs_ElementType(theType));
1500   TPythonDump()<<this<<".SetElementType("<<theType<<")";
1501 }
1502
1503 void ElemEntityType_i::SetEntityType(EntityType theEntityType)
1504 {
1505   myElemEntityTypePtr->SetElemEntityType(SMDSAbs_EntityType (theEntityType));
1506   TPythonDump()<<this<<".SetEntityType("<<theEntityType<<")";
1507 }
1508 EntityType ElemEntityType_i::GetEntityType() const
1509 {
1510  return (EntityType) myElemEntityTypePtr->GetElemEntityType();
1511 }
1512
1513 FunctorType ElemEntityType_i::GetFunctorType()
1514 {
1515   return SMESH::FT_EntityType;
1516 }
1517
1518 /*
1519   Class       : CoplanarFaces_i
1520   Description : Returns true if a mesh face is a coplanar neighbour to a given one
1521 */
1522 CoplanarFaces_i::CoplanarFaces_i()
1523 {
1524   myCoplanarFacesPtr.reset(new Controls::CoplanarFaces());
1525   myFunctorPtr = myPredicatePtr = myCoplanarFacesPtr;
1526 }
1527
1528 void CoplanarFaces_i::SetFace ( CORBA::Long theFaceID )
1529 {
1530   myCoplanarFacesPtr->SetFace(theFaceID);
1531   TPythonDump()<<this<<".SetFace("<<theFaceID<<")";
1532 }
1533
1534 void CoplanarFaces_i::SetTolerance( CORBA::Double theToler )
1535 {
1536   myCoplanarFacesPtr->SetTolerance(theToler);
1537   TPythonDump()<<this<<".SetTolerance("<<theToler<<")";
1538 }
1539
1540 CORBA::Long CoplanarFaces_i::GetFace () const
1541 {
1542   return myCoplanarFacesPtr->GetFace();
1543 }
1544
1545 char* CoplanarFaces_i::GetFaceAsString () const
1546 {
1547   TCollection_AsciiString str(Standard_Integer(myCoplanarFacesPtr->GetFace()));
1548   return CORBA::string_dup( str.ToCString() );
1549 }
1550
1551 CORBA::Double CoplanarFaces_i::GetTolerance() const
1552 {
1553   return myCoplanarFacesPtr->GetTolerance();
1554 }
1555
1556 FunctorType CoplanarFaces_i::GetFunctorType()
1557 {
1558   return SMESH::FT_CoplanarFaces;
1559 }
1560
1561 /*
1562  * Class       : ConnectedElements_i
1563  * Description : Returns true if an element is connected via other elements to the element
1564  *               located at a given point.
1565  */
1566 ConnectedElements_i::ConnectedElements_i()
1567 {
1568   myConnectedElementsPtr.reset(new Controls::ConnectedElements());
1569   myFunctorPtr = myPredicatePtr = myConnectedElementsPtr;
1570 }
1571
1572 FunctorType ConnectedElements_i::GetFunctorType()
1573 {
1574   return FT_ConnectedElements;
1575 }
1576
1577 void ConnectedElements_i::SetElementType( ElementType theType )
1578 {
1579   myConnectedElementsPtr->SetType( SMDSAbs_ElementType( theType ));
1580   TPythonDump() << this << ".SetElementType( " << theType << " )";
1581 }
1582
1583 void ConnectedElements_i::SetPoint( CORBA::Double x, CORBA::Double y, CORBA::Double z )
1584 {
1585   myConnectedElementsPtr->SetPoint( x,y,z );
1586   myVertexID.clear();
1587   TPythonDump() << this << ".SetPoint( " << x << ", " << y << ", " << z << " )";
1588 }
1589
1590 void ConnectedElements_i::SetVertex( GEOM::GEOM_Object_ptr vertex )
1591   throw (SALOME::SALOME_Exception)
1592 {
1593   TopoDS_Shape shape = SMESH_Gen_i::GetSMESHGen()->GeomObjectToShape( vertex );
1594   if ( shape.IsNull() )
1595     THROW_SALOME_CORBA_EXCEPTION( "ConnectedElements_i::SetVertex(): NULL Vertex",
1596                                   SALOME::BAD_PARAM );
1597
1598   TopExp_Explorer v( shape, TopAbs_VERTEX );
1599   if ( !v.More() )
1600     THROW_SALOME_CORBA_EXCEPTION( "ConnectedElements_i::SetVertex(): empty vertex",
1601                                   SALOME::BAD_PARAM );
1602
1603   gp_Pnt p = BRep_Tool::Pnt( TopoDS::Vertex( v.Current() ));
1604   myConnectedElementsPtr->SetPoint( p.X(), p.Y(), p.Z() );
1605   //
1606   CORBA::String_var id = vertex->GetStudyEntry();
1607   myVertexID = id.in();
1608
1609   TPythonDump() << this << ".SetVertex( " << vertex << " )";
1610 }
1611
1612 void ConnectedElements_i::SetNode ( CORBA::Long nodeID )
1613   throw (SALOME::SALOME_Exception)
1614 {
1615   if ( nodeID < 1 )
1616     THROW_SALOME_CORBA_EXCEPTION( "ConnectedElements_i::SetNode(): nodeID must be > 0",
1617                                   SALOME::BAD_PARAM );
1618
1619   myConnectedElementsPtr->SetNode( nodeID );
1620   myVertexID.clear();
1621   TPythonDump() << this << ".SetNode( " << nodeID << " )";
1622 }
1623
1624 /*!
1625  * \brief This is a comfort method for Filter dialog
1626  */
1627 void ConnectedElements_i::SetThreshold ( const char*                             threshold,
1628                                          SMESH::ConnectedElements::ThresholdType type )
1629   throw (SALOME::SALOME_Exception)
1630 {
1631   if ( !threshold )
1632     THROW_SALOME_CORBA_EXCEPTION( "ConnectedElements_i::SetThreshold(): NULL threshold",
1633                                   SALOME::BAD_PARAM );
1634   switch ( type )
1635   {
1636   case SMESH::ConnectedElements::POINT: // read 3 node coordinates ///////////////////
1637     {
1638       vector< double > xyz;
1639       char* endptr;
1640       do
1641       {
1642         // skip a separator
1643         while ( *threshold &&
1644                 *threshold != '+' &&
1645                 *threshold != '-' &&
1646                 !isdigit( *threshold ))
1647           ++threshold;
1648         if ( !*threshold )
1649           break;
1650         // read a coordinate
1651         xyz.push_back( strtod( threshold, &endptr ));
1652         if ( threshold == endptr )
1653         {
1654           xyz.resize( xyz.size() - 1 );
1655           break;
1656         }
1657         threshold = endptr;
1658       }
1659       while ( xyz.size() < 3 );
1660
1661       if ( xyz.size() < 3 )
1662         THROW_SALOME_CORBA_EXCEPTION
1663           ( "ConnectedElements_i::SetThreshold(): invalid point coordinates", SALOME::BAD_PARAM );
1664
1665       SetPoint( xyz[0], xyz[1], xyz[2] );
1666       break;
1667     }
1668   case SMESH::ConnectedElements::VERTEX: // get a VERTEX by its entry /////////////////
1669     {
1670       SALOMEDS::Study_var study = SMESH_Gen_i::GetSMESHGen()->GetCurrentStudy();
1671       if ( study->_is_nil() )
1672         THROW_SALOME_CORBA_EXCEPTION
1673           ( "ConnectedElements_i::SetThreshold(): NULL current study", SALOME::BAD_PARAM );
1674       SALOMEDS::SObject_wrap sobj = study->FindObjectID( threshold );
1675       if ( sobj->_is_nil() )
1676         THROW_SALOME_CORBA_EXCEPTION
1677           ( "ConnectedElements_i::SetThreshold(): invalid vertex study entry", SALOME::BAD_PARAM );
1678       CORBA::Object_var        obj = sobj->GetObject();
1679       GEOM::GEOM_Object_var vertex = GEOM::GEOM_Object::_narrow( obj );
1680       if ( vertex->_is_nil() )
1681         THROW_SALOME_CORBA_EXCEPTION
1682           ( "ConnectedElements_i::SetThreshold(): no GEOM_Object in SObject", SALOME::BAD_PARAM );
1683       SetVertex( vertex );
1684       break;
1685     }
1686   case SMESH::ConnectedElements::NODE: // read a node ID ////////////////////////////
1687     {
1688       char* endptr;
1689       int id = strtol( threshold, &endptr, 10 );
1690       if ( threshold == endptr )
1691         THROW_SALOME_CORBA_EXCEPTION
1692           ( "ConnectedElements_i::SetThreshold(): invalid node ID", SALOME::BAD_PARAM );
1693       SetNode( id );
1694       break;
1695     }
1696   default:
1697     THROW_SALOME_CORBA_EXCEPTION
1698       ( "ConnectedElements_i::SetThreshold(): invalid ThresholdType", SALOME::BAD_PARAM );
1699   }
1700 }
1701
1702 char* ConnectedElements_i::GetThreshold ( SMESH::ConnectedElements::ThresholdType& type )
1703 {
1704   std::string threshold;
1705   if ( !myVertexID.empty() )
1706   {
1707     threshold = myVertexID;
1708     type      = SMESH::ConnectedElements::VERTEX;
1709   }
1710   else
1711   {
1712     std::vector<double> xyz = myConnectedElementsPtr->GetPoint();
1713     if ( xyz.size() == 3 )
1714     {
1715       threshold = SMESH_Comment( xyz[0] ) << "; " << xyz[1] << "; " << xyz[2];
1716       type      = SMESH::ConnectedElements::POINT;
1717     }
1718     else
1719     {
1720       threshold = SMESH_Comment( myConnectedElementsPtr->GetNode() );
1721       type      = SMESH::ConnectedElements::NODE;
1722     }
1723   }
1724   return CORBA::string_dup( threshold.c_str() );
1725 }
1726
1727 /*
1728   Class       : Comparator_i
1729   Description : Base class for comparators
1730 */
1731 Comparator_i::Comparator_i():
1732   myNumericalFunctor( NULL )
1733 {}
1734
1735 Comparator_i::~Comparator_i()
1736 {
1737   if ( myNumericalFunctor )
1738     myNumericalFunctor->UnRegister();
1739 }
1740
1741 void Comparator_i::SetMargin( CORBA::Double theValue )
1742 {
1743   myComparatorPtr->SetMargin( theValue );
1744   TPythonDump()<<this<<".SetMargin("<<theValue<<")";
1745 }
1746
1747 CORBA::Double Comparator_i::GetMargin()
1748 {
1749   return myComparatorPtr->GetMargin();
1750 }
1751
1752 void Comparator_i::SetNumFunctor( NumericalFunctor_ptr theFunct )
1753 {
1754   if ( myNumericalFunctor )
1755     myNumericalFunctor->UnRegister();
1756
1757   myNumericalFunctor = DownCast<NumericalFunctor_i*>(theFunct);
1758
1759   if ( myNumericalFunctor )
1760   {
1761     myComparatorPtr->SetNumFunctor( myNumericalFunctor->GetNumericalFunctor() );
1762     myNumericalFunctor->Register();
1763     TPythonDump()<<this<<".SetNumFunctor("<<myNumericalFunctor<<")";
1764   }
1765 }
1766
1767 Controls::ComparatorPtr Comparator_i::GetComparator()
1768 {
1769   return myComparatorPtr;
1770 }
1771
1772 NumericalFunctor_i* Comparator_i::GetNumFunctor_i()
1773 {
1774   return myNumericalFunctor;
1775 }
1776
1777
1778 /*
1779   Class       : LessThan_i
1780   Description : Comparator "<"
1781 */
1782 LessThan_i::LessThan_i()
1783 {
1784   myComparatorPtr.reset( new Controls::LessThan() );
1785   myFunctorPtr = myPredicatePtr = myComparatorPtr;
1786 }
1787
1788 FunctorType LessThan_i::GetFunctorType()
1789 {
1790   return SMESH::FT_LessThan;
1791 }
1792
1793
1794 /*
1795   Class       : MoreThan_i
1796   Description : Comparator ">"
1797 */
1798 MoreThan_i::MoreThan_i()
1799 {
1800   myComparatorPtr.reset( new Controls::MoreThan() );
1801   myFunctorPtr = myPredicatePtr = myComparatorPtr;
1802 }
1803
1804 FunctorType MoreThan_i::GetFunctorType()
1805 {
1806   return SMESH::FT_MoreThan;
1807 }
1808
1809
1810 /*
1811   Class       : EqualTo_i
1812   Description : Comparator "="
1813 */
1814 EqualTo_i::EqualTo_i()
1815 : myEqualToPtr( new Controls::EqualTo() )
1816 {
1817   myFunctorPtr = myPredicatePtr = myComparatorPtr = myEqualToPtr;
1818 }
1819
1820 void EqualTo_i::SetTolerance( CORBA::Double theToler )
1821 {
1822   myEqualToPtr->SetTolerance( theToler );
1823   TPythonDump()<<this<<".SetTolerance("<<theToler<<")";
1824 }
1825
1826 CORBA::Double EqualTo_i::GetTolerance()
1827 {
1828   return myEqualToPtr->GetTolerance();
1829 }
1830
1831 FunctorType EqualTo_i::GetFunctorType()
1832 {
1833   return SMESH::FT_EqualTo;
1834 }
1835
1836 /*
1837   Class       : LogicalNOT_i
1838   Description : Logical NOT predicate
1839 */
1840 LogicalNOT_i::LogicalNOT_i()
1841 : myPredicate( NULL ),
1842   myLogicalNOTPtr( new Controls::LogicalNOT() )
1843 {
1844   myFunctorPtr = myPredicatePtr = myLogicalNOTPtr;
1845 }
1846
1847 LogicalNOT_i::~LogicalNOT_i()
1848 {
1849   if ( myPredicate )
1850     myPredicate->UnRegister();
1851 }
1852
1853 void LogicalNOT_i::SetPredicate( Predicate_ptr thePredicate )
1854 {
1855   if ( myPredicate )
1856     myPredicate->UnRegister();
1857
1858   myPredicate = SMESH::GetPredicate(thePredicate);
1859
1860   if ( myPredicate ){
1861     myLogicalNOTPtr->SetPredicate(myPredicate->GetPredicate());
1862     myPredicate->Register();
1863     TPythonDump()<<this<<".SetPredicate("<<myPredicate<<")";
1864   }
1865 }
1866
1867 FunctorType LogicalNOT_i::GetFunctorType()
1868 {
1869   return SMESH::FT_LogicalNOT;
1870 }
1871
1872 Predicate_i* LogicalNOT_i::GetPredicate_i()
1873 {
1874   return myPredicate;
1875 }
1876
1877
1878
1879 /*
1880   Class       : LogicalBinary_i
1881   Description : Base class for binary logical predicate
1882 */
1883 LogicalBinary_i::LogicalBinary_i()
1884 : myPredicate1( NULL ),
1885   myPredicate2( NULL )
1886 {}
1887
1888 LogicalBinary_i::~LogicalBinary_i()
1889 {
1890   if ( myPredicate1 )
1891     myPredicate1->UnRegister();
1892
1893   if ( myPredicate2 )
1894     myPredicate2->UnRegister();
1895 }
1896
1897 void LogicalBinary_i::SetMesh( SMESH_Mesh_ptr theMesh )
1898 {
1899   if ( myPredicate1 )
1900     myPredicate1->SetMesh( theMesh );
1901
1902   if ( myPredicate2 )
1903     myPredicate2->SetMesh( theMesh );
1904 }
1905
1906 void LogicalBinary_i::SetPredicate1( Predicate_ptr thePredicate )
1907 {
1908   if ( myPredicate1 )
1909     myPredicate1->UnRegister();
1910
1911   myPredicate1 = SMESH::GetPredicate(thePredicate);
1912
1913   if ( myPredicate1 ){
1914     myLogicalBinaryPtr->SetPredicate1(myPredicate1->GetPredicate());
1915     myPredicate1->Register();
1916     TPythonDump()<<this<<".SetPredicate1("<<myPredicate1<<")";
1917   }
1918 }
1919
1920 void LogicalBinary_i::SetPredicate2( Predicate_ptr thePredicate )
1921 {
1922   if ( myPredicate2 )
1923     myPredicate2->UnRegister();
1924
1925   myPredicate2 = SMESH::GetPredicate(thePredicate);
1926
1927   if ( myPredicate2 ){
1928     myLogicalBinaryPtr->SetPredicate2(myPredicate2->GetPredicate());
1929     myPredicate2->Register();
1930     TPythonDump()<<this<<".SetPredicate2("<<myPredicate2<<")";
1931   }
1932 }
1933
1934 Controls::LogicalBinaryPtr LogicalBinary_i::GetLogicalBinary()
1935 {
1936   return myLogicalBinaryPtr;
1937 }
1938
1939 Predicate_i* LogicalBinary_i::GetPredicate1_i()
1940 {
1941   return myPredicate1;
1942 }
1943 Predicate_i* LogicalBinary_i::GetPredicate2_i()
1944 {
1945   return myPredicate2;
1946 }
1947
1948
1949 /*
1950   Class       : LogicalAND_i
1951   Description : Logical AND
1952 */
1953 LogicalAND_i::LogicalAND_i()
1954 {
1955   myLogicalBinaryPtr.reset( new Controls::LogicalAND() );
1956   myFunctorPtr = myPredicatePtr = myLogicalBinaryPtr;
1957 }
1958
1959 FunctorType LogicalAND_i::GetFunctorType()
1960 {
1961   return SMESH::FT_LogicalAND;
1962 }
1963
1964
1965 /*
1966   Class       : LogicalOR_i
1967   Description : Logical OR
1968 */
1969 LogicalOR_i::LogicalOR_i()
1970 {
1971   myLogicalBinaryPtr.reset( new Controls::LogicalOR() );
1972   myFunctorPtr = myPredicatePtr = myLogicalBinaryPtr;
1973 }
1974
1975 FunctorType LogicalOR_i::GetFunctorType()
1976 {
1977   return SMESH::FT_LogicalOR;
1978 }
1979
1980
1981 /*
1982                             FILTER MANAGER
1983 */
1984
1985 FilterManager_i::FilterManager_i()
1986 : SALOME::GenericObj_i( SMESH_Gen_i::GetPOA() )
1987 {
1988   //Base class Salome_GenericObject do it inmplicitly by overriding PortableServer::POA_ptr _default_POA() method
1989   //PortableServer::ObjectId_var anObjectId =
1990   //  SMESH_Gen_i::GetPOA()->activate_object( this );
1991 }
1992
1993
1994 FilterManager_i::~FilterManager_i()
1995 {
1996   //TPythonDump()<<this<<".UnRegister()";
1997 }
1998
1999
2000 MinimumAngle_ptr FilterManager_i::CreateMinimumAngle()
2001 {
2002   SMESH::MinimumAngle_i* aServant = new SMESH::MinimumAngle_i();
2003   SMESH::MinimumAngle_var anObj = aServant->_this();
2004   TPythonDump()<<aServant<<" = "<<this<<".CreateMinimumAngle()";
2005   return anObj._retn();
2006 }
2007
2008
2009 AspectRatio_ptr FilterManager_i::CreateAspectRatio()
2010 {
2011   SMESH::AspectRatio_i* aServant = new SMESH::AspectRatio_i();
2012   SMESH::AspectRatio_var anObj = aServant->_this();
2013   TPythonDump()<<aServant<<" = "<<this<<".CreateAspectRatio()";
2014   return anObj._retn();
2015 }
2016
2017
2018 AspectRatio3D_ptr FilterManager_i::CreateAspectRatio3D()
2019 {
2020   SMESH::AspectRatio3D_i* aServant = new SMESH::AspectRatio3D_i();
2021   SMESH::AspectRatio3D_var anObj = aServant->_this();
2022   TPythonDump()<<aServant<<" = "<<this<<".CreateAspectRatio3D()";
2023   return anObj._retn();
2024 }
2025
2026
2027 Warping_ptr FilterManager_i::CreateWarping()
2028 {
2029   SMESH::Warping_i* aServant = new SMESH::Warping_i();
2030   SMESH::Warping_var anObj = aServant->_this();
2031   TPythonDump()<<aServant<<" = "<<this<<".CreateWarping()";
2032   return anObj._retn();
2033 }
2034
2035
2036 Taper_ptr FilterManager_i::CreateTaper()
2037 {
2038   SMESH::Taper_i* aServant = new SMESH::Taper_i();
2039   SMESH::Taper_var anObj = aServant->_this();
2040   TPythonDump()<<aServant<<" = "<<this<<".CreateTaper()";
2041   return anObj._retn();
2042 }
2043
2044
2045 Skew_ptr FilterManager_i::CreateSkew()
2046 {
2047   SMESH::Skew_i* aServant = new SMESH::Skew_i();
2048   SMESH::Skew_var anObj = aServant->_this();
2049   TPythonDump()<<aServant<<" = "<<this<<".CreateSkew()";
2050   return anObj._retn();
2051 }
2052
2053
2054 Area_ptr FilterManager_i::CreateArea()
2055 {
2056   SMESH::Area_i* aServant = new SMESH::Area_i();
2057   SMESH::Area_var anObj = aServant->_this();
2058   TPythonDump()<<aServant<<" = "<<this<<".CreateArea()";
2059   return anObj._retn();
2060 }
2061
2062
2063 Volume3D_ptr FilterManager_i::CreateVolume3D()
2064 {
2065   SMESH::Volume3D_i* aServant = new SMESH::Volume3D_i();
2066   SMESH::Volume3D_var anObj = aServant->_this();
2067   TPythonDump()<<aServant<<" = "<<this<<".CreateVolume3D()";
2068   return anObj._retn();
2069 }
2070
2071
2072 MaxElementLength2D_ptr FilterManager_i::CreateMaxElementLength2D()
2073 {
2074   SMESH::MaxElementLength2D_i* aServant = new SMESH::MaxElementLength2D_i();
2075   SMESH::MaxElementLength2D_var anObj = aServant->_this();
2076   TPythonDump()<<aServant<<" = "<<this<<".CreateMaxElementLength2D()";
2077   return anObj._retn();
2078 }
2079
2080
2081 MaxElementLength3D_ptr FilterManager_i::CreateMaxElementLength3D()
2082 {
2083   SMESH::MaxElementLength3D_i* aServant = new SMESH::MaxElementLength3D_i();
2084   SMESH::MaxElementLength3D_var anObj = aServant->_this();
2085   TPythonDump()<<aServant<<" = "<<this<<".CreateMaxElementLength3D()";
2086   return anObj._retn();
2087 }
2088
2089
2090 Length_ptr FilterManager_i::CreateLength()
2091 {
2092   SMESH::Length_i* aServant = new SMESH::Length_i();
2093   SMESH::Length_var anObj = aServant->_this();
2094   TPythonDump()<<aServant<<" = "<<this<<".CreateLength()";
2095   return anObj._retn();
2096 }
2097
2098 Length2D_ptr FilterManager_i::CreateLength2D()
2099 {
2100   SMESH::Length2D_i* aServant = new SMESH::Length2D_i();
2101   SMESH::Length2D_var anObj = aServant->_this();
2102   TPythonDump()<<aServant<<" = "<<this<<".CreateLength2D()";
2103   return anObj._retn();
2104 }
2105
2106 MultiConnection_ptr FilterManager_i::CreateMultiConnection()
2107 {
2108   SMESH::MultiConnection_i* aServant = new SMESH::MultiConnection_i();
2109   SMESH::MultiConnection_var anObj = aServant->_this();
2110   TPythonDump()<<aServant<<" = "<<this<<".CreateMultiConnection()";
2111   return anObj._retn();
2112 }
2113
2114 MultiConnection2D_ptr FilterManager_i::CreateMultiConnection2D()
2115 {
2116   SMESH::MultiConnection2D_i* aServant = new SMESH::MultiConnection2D_i();
2117   SMESH::MultiConnection2D_var anObj = aServant->_this();
2118   TPythonDump()<<aServant<<" = "<<this<<".CreateMultiConnection2D()";
2119   return anObj._retn();
2120 }
2121
2122 BallDiameter_ptr FilterManager_i::CreateBallDiameter()
2123 {
2124   SMESH::BallDiameter_i* aServant = new SMESH::BallDiameter_i();
2125   SMESH::BallDiameter_var anObj = aServant->_this();
2126   TPythonDump()<<aServant<<" = "<<this<<".CreateBallDiameter()";
2127   return anObj._retn();
2128 }
2129
2130 BelongToMeshGroup_ptr FilterManager_i::CreateBelongToMeshGroup()
2131 {
2132   SMESH::BelongToMeshGroup_i* aServant = new SMESH::BelongToMeshGroup_i();
2133   SMESH::BelongToMeshGroup_var anObj = aServant->_this();
2134   TPythonDump()<<aServant<<" = "<<this<<".CreateBelongToMeshGroup()";
2135   return anObj._retn();
2136 }
2137
2138 BelongToGeom_ptr FilterManager_i::CreateBelongToGeom()
2139 {
2140   SMESH::BelongToGeom_i* aServant = new SMESH::BelongToGeom_i();
2141   SMESH::BelongToGeom_var anObj = aServant->_this();
2142   TPythonDump()<<aServant<<" = "<<this<<".CreateBelongToGeom()";
2143   return anObj._retn();
2144 }
2145
2146 BelongToPlane_ptr FilterManager_i::CreateBelongToPlane()
2147 {
2148   SMESH::BelongToPlane_i* aServant = new SMESH::BelongToPlane_i();
2149   SMESH::BelongToPlane_var anObj = aServant->_this();
2150   TPythonDump()<<aServant<<" = "<<this<<".CreateBelongToPlane()";
2151   return anObj._retn();
2152 }
2153
2154 BelongToCylinder_ptr FilterManager_i::CreateBelongToCylinder()
2155 {
2156   SMESH::BelongToCylinder_i* aServant = new SMESH::BelongToCylinder_i();
2157   SMESH::BelongToCylinder_var anObj = aServant->_this();
2158   TPythonDump()<<aServant<<" = "<<this<<".CreateBelongToCylinder()";
2159   return anObj._retn();
2160 }
2161
2162 BelongToGenSurface_ptr FilterManager_i::CreateBelongToGenSurface()
2163 {
2164   SMESH::BelongToGenSurface_i* aServant = new SMESH::BelongToGenSurface_i();
2165   SMESH::BelongToGenSurface_var anObj = aServant->_this();
2166   TPythonDump()<<aServant<<" = "<<this<<".CreateBelongToGenSurface()";
2167   return anObj._retn();
2168 }
2169
2170 LyingOnGeom_ptr FilterManager_i::CreateLyingOnGeom()
2171 {
2172   SMESH::LyingOnGeom_i* aServant = new SMESH::LyingOnGeom_i();
2173   SMESH::LyingOnGeom_var anObj = aServant->_this();
2174   TPythonDump()<<aServant<<" = "<<this<<".CreateLyingOnGeom()";
2175   return anObj._retn();
2176 }
2177
2178 CoplanarFaces_ptr FilterManager_i::CreateCoplanarFaces()
2179 {
2180   SMESH::CoplanarFaces_i* aServant = new SMESH::CoplanarFaces_i();
2181   SMESH::CoplanarFaces_var anObj = aServant->_this();
2182   TPythonDump()<<aServant<<" = "<<this<<".CreateCoplanarFaces()";
2183   return anObj._retn();
2184 }
2185
2186 ConnectedElements_ptr FilterManager_i::CreateConnectedElements()
2187 {
2188   SMESH::ConnectedElements_i* aServant = new SMESH::ConnectedElements_i();
2189   SMESH::ConnectedElements_var anObj = aServant->_this();
2190   TPythonDump()<<aServant<<" = "<<this<<".CreateConnectedElements()";
2191   return anObj._retn();
2192 }
2193
2194 FreeBorders_ptr FilterManager_i::CreateFreeBorders()
2195 {
2196   SMESH::FreeBorders_i* aServant = new SMESH::FreeBorders_i();
2197   SMESH::FreeBorders_var anObj = aServant->_this();
2198   TPythonDump()<<aServant<<" = "<<this<<".CreateFreeBorders()";
2199   return anObj._retn();
2200 }
2201
2202 FreeEdges_ptr FilterManager_i::CreateFreeEdges()
2203 {
2204   SMESH::FreeEdges_i* aServant = new SMESH::FreeEdges_i();
2205   SMESH::FreeEdges_var anObj = aServant->_this();
2206   TPythonDump()<<aServant<<" = "<<this<<".CreateFreeEdges()";
2207   return anObj._retn();
2208 }
2209
2210 FreeFaces_ptr FilterManager_i::CreateFreeFaces()
2211 {
2212   SMESH::FreeFaces_i* aServant = new SMESH::FreeFaces_i();
2213   SMESH::FreeFaces_var anObj = aServant->_this();
2214   TPythonDump()<<aServant<<" = "<<this<<".CreateFreeFaces()";
2215   return anObj._retn();
2216 }
2217
2218 FreeNodes_ptr FilterManager_i::CreateFreeNodes()
2219 {
2220   SMESH::FreeNodes_i* aServant = new SMESH::FreeNodes_i();
2221   SMESH::FreeNodes_var anObj = aServant->_this();
2222   TPythonDump()<<aServant<<" = "<<this<<".CreateFreeNodes()";
2223   return anObj._retn();
2224 }
2225
2226 EqualNodes_ptr FilterManager_i::CreateEqualNodes()
2227 {
2228   SMESH::EqualNodes_i* aServant = new SMESH::EqualNodes_i();
2229   SMESH::EqualNodes_var anObj = aServant->_this();
2230   TPythonDump()<<aServant<<" = "<<this<<".CreateEqualNodes()";
2231   return anObj._retn();
2232 }
2233
2234 EqualEdges_ptr FilterManager_i::CreateEqualEdges()
2235 {
2236   SMESH::EqualEdges_i* aServant = new SMESH::EqualEdges_i();
2237   SMESH::EqualEdges_var anObj = aServant->_this();
2238   TPythonDump()<<aServant<<" = "<<this<<".CreateEqualEdges()";
2239   return anObj._retn();
2240 }
2241 EqualFaces_ptr FilterManager_i::CreateEqualFaces()
2242 {
2243   SMESH::EqualFaces_i* aServant = new SMESH::EqualFaces_i();
2244   SMESH::EqualFaces_var anObj = aServant->_this();
2245   TPythonDump()<<aServant<<" = "<<this<<".CreateEqualFaces()";
2246   return anObj._retn();
2247 }
2248 EqualVolumes_ptr FilterManager_i::CreateEqualVolumes()
2249 {
2250   SMESH::EqualVolumes_i* aServant = new SMESH::EqualVolumes_i();
2251   SMESH::EqualVolumes_var anObj = aServant->_this();
2252   TPythonDump()<<aServant<<" = "<<this<<".CreateEqualVolumes()";
2253   return anObj._retn();
2254 }
2255
2256 RangeOfIds_ptr FilterManager_i::CreateRangeOfIds()
2257 {
2258   SMESH::RangeOfIds_i* aServant = new SMESH::RangeOfIds_i();
2259   SMESH::RangeOfIds_var anObj = aServant->_this();
2260   TPythonDump()<<aServant<<" = "<<this<<".CreateRangeOfIds()";
2261   return anObj._retn();
2262 }
2263
2264 BadOrientedVolume_ptr FilterManager_i::CreateBadOrientedVolume()
2265 {
2266   SMESH::BadOrientedVolume_i* aServant = new SMESH::BadOrientedVolume_i();
2267   SMESH::BadOrientedVolume_var anObj = aServant->_this();
2268   TPythonDump()<<aServant<<" = "<<this<<".CreateBadOrientedVolume()";
2269   return anObj._retn();
2270 }
2271
2272 BareBorderVolume_ptr FilterManager_i::CreateBareBorderVolume()
2273 {
2274   SMESH::BareBorderVolume_i* aServant = new SMESH::BareBorderVolume_i();
2275   SMESH::BareBorderVolume_var anObj = aServant->_this();
2276   TPythonDump()<<aServant<<" = "<<this<<".CreateBareBorderVolume()";
2277   return anObj._retn();
2278 }
2279
2280 BareBorderFace_ptr FilterManager_i::CreateBareBorderFace()
2281 {
2282   SMESH::BareBorderFace_i* aServant = new SMESH::BareBorderFace_i();
2283   SMESH::BareBorderFace_var anObj = aServant->_this();
2284   TPythonDump()<<aServant<<" = "<<this<<".CreateBareBorderFace()";
2285   return anObj._retn();
2286 }
2287
2288 OverConstrainedVolume_ptr FilterManager_i::CreateOverConstrainedVolume()
2289 {
2290   SMESH::OverConstrainedVolume_i* aServant = new SMESH::OverConstrainedVolume_i();
2291   SMESH::OverConstrainedVolume_var anObj = aServant->_this();
2292   TPythonDump()<<aServant<<" = "<<this<<".CreateOverConstrainedVolume()";
2293   return anObj._retn();
2294 }
2295
2296 OverConstrainedFace_ptr FilterManager_i::CreateOverConstrainedFace()
2297 {
2298   SMESH::OverConstrainedFace_i* aServant = new SMESH::OverConstrainedFace_i();
2299   SMESH::OverConstrainedFace_var anObj = aServant->_this();
2300   TPythonDump()<<aServant<<" = "<<this<<".CreateOverConstrainedFace()";
2301   return anObj._retn();
2302 }
2303
2304 LessThan_ptr FilterManager_i::CreateLessThan()
2305 {
2306   SMESH::LessThan_i* aServant = new SMESH::LessThan_i();
2307   SMESH::LessThan_var anObj = aServant->_this();
2308   TPythonDump()<<aServant<<" = "<<this<<".CreateLessThan()";
2309   return anObj._retn();
2310 }
2311
2312 MoreThan_ptr FilterManager_i::CreateMoreThan()
2313 {
2314   SMESH::MoreThan_i* aServant = new SMESH::MoreThan_i();
2315   SMESH::MoreThan_var anObj = aServant->_this();
2316   TPythonDump()<<aServant<<" = "<<this<<".CreateMoreThan()";
2317   return anObj._retn();
2318 }
2319
2320 EqualTo_ptr FilterManager_i::CreateEqualTo()
2321 {
2322   SMESH::EqualTo_i* aServant = new SMESH::EqualTo_i();
2323   SMESH::EqualTo_var anObj = aServant->_this();
2324   TPythonDump()<<aServant<<" = "<<this<<".CreateEqualTo()";
2325   return anObj._retn();
2326 }
2327
2328 LogicalNOT_ptr FilterManager_i::CreateLogicalNOT()
2329 {
2330   SMESH::LogicalNOT_i* aServant = new SMESH::LogicalNOT_i();
2331   SMESH::LogicalNOT_var anObj = aServant->_this();
2332   TPythonDump()<<aServant<<" = "<<this<<".CreateLogicalNOT()";
2333   return anObj._retn();
2334 }
2335
2336 LogicalAND_ptr FilterManager_i::CreateLogicalAND()
2337 {
2338   SMESH::LogicalAND_i* aServant = new SMESH::LogicalAND_i();
2339   SMESH::LogicalAND_var anObj = aServant->_this();
2340   TPythonDump()<<aServant<<" = "<<this<<".CreateLogicalAND()";
2341   return anObj._retn();
2342 }
2343
2344 LogicalOR_ptr FilterManager_i::CreateLogicalOR()
2345 {
2346   SMESH::LogicalOR_i* aServant = new SMESH::LogicalOR_i();
2347   SMESH::LogicalOR_var anObj = aServant->_this();
2348   TPythonDump()<<aServant<<" = "<<this<<".CreateLogicalOR()";
2349   return anObj._retn();
2350 }
2351
2352 LinearOrQuadratic_ptr FilterManager_i::CreateLinearOrQuadratic()
2353 {
2354   SMESH::LinearOrQuadratic_i* aServant = new SMESH::LinearOrQuadratic_i();
2355   SMESH::LinearOrQuadratic_var anObj = aServant->_this();
2356   TPythonDump()<<aServant<<" = "<<this<<".CreateLinearOrQuadratic()";
2357   return anObj._retn();
2358 }
2359
2360 GroupColor_ptr FilterManager_i::CreateGroupColor()
2361 {
2362   SMESH::GroupColor_i* aServant = new SMESH::GroupColor_i();
2363   SMESH::GroupColor_var anObj = aServant->_this();
2364   TPythonDump()<<aServant<<" = "<<this<<".CreateGroupColor()";
2365   return anObj._retn();
2366 }
2367
2368 ElemGeomType_ptr FilterManager_i::CreateElemGeomType()
2369 {
2370   SMESH::ElemGeomType_i* aServant = new SMESH::ElemGeomType_i();
2371   SMESH::ElemGeomType_var anObj = aServant->_this();
2372   TPythonDump()<<aServant<<" = "<<this<<".CreateElemGeomType()";
2373   return anObj._retn();
2374 }
2375
2376 ElemEntityType_ptr FilterManager_i::CreateElemEntityType()
2377 {
2378   SMESH::ElemEntityType_i* aServant = new SMESH::ElemEntityType_i();
2379   SMESH::ElemEntityType_var anObj = aServant->_this();
2380   TPythonDump()<<aServant<<" = "<<this<<".CreateElemEntityType()";
2381   return anObj._retn();
2382 }
2383
2384 Filter_ptr FilterManager_i::CreateFilter()
2385 {
2386   SMESH::Filter_i* aServant = new SMESH::Filter_i();
2387   SMESH::Filter_var anObj = aServant->_this();
2388   TPythonDump()<<aServant<<" = "<<this<<".CreateFilter()";
2389   return anObj._retn();
2390 }
2391
2392 FilterLibrary_ptr FilterManager_i::LoadLibrary( const char* aFileName )
2393 {
2394   SMESH::FilterLibrary_i* aServant = new SMESH::FilterLibrary_i( aFileName );
2395   SMESH::FilterLibrary_var anObj = aServant->_this();
2396   TPythonDump()<<aServant<<" = "<<this<<".LoadLibrary('"<<aFileName<<"')";
2397   return anObj._retn();
2398 }
2399
2400 FilterLibrary_ptr FilterManager_i::CreateLibrary()
2401 {
2402   SMESH::FilterLibrary_i* aServant = new SMESH::FilterLibrary_i();
2403   SMESH::FilterLibrary_var anObj = aServant->_this();
2404   TPythonDump()<<aServant<<" = "<<this<<".CreateLibrary()";
2405   return anObj._retn();
2406 }
2407
2408 CORBA::Boolean FilterManager_i::DeleteLibrary( const char* aFileName )
2409 {
2410   TPythonDump()<<this<<".DeleteLibrary("<<aFileName<<")";
2411   return remove( aFileName ) ? false : true;
2412 }
2413
2414 //=============================================================================
2415 /*!
2416  *  SMESH_Gen_i::CreateFilterManager
2417  *
2418  *  Create filter manager
2419  */
2420 //=============================================================================
2421
2422 SMESH::FilterManager_ptr SMESH_Gen_i::CreateFilterManager()
2423 {
2424   SMESH::FilterManager_i* aFilter = new SMESH::FilterManager_i();
2425   SMESH::FilterManager_var anObj = aFilter->_this();
2426   return anObj._retn();
2427 }
2428
2429
2430 /*
2431                               FILTER
2432 */
2433
2434 //=======================================================================
2435 // name    : Filter_i::Filter_i
2436 // Purpose : Constructor
2437 //=======================================================================
2438 Filter_i::Filter_i()
2439 : myPredicate( NULL )
2440 {}
2441
2442 //=======================================================================
2443 // name    : Filter_i::~Filter_i
2444 // Purpose : Destructor
2445 //=======================================================================
2446 Filter_i::~Filter_i()
2447 {
2448   if ( myPredicate )
2449     myPredicate->UnRegister();
2450
2451   if(!CORBA::is_nil(myMesh))
2452     myMesh->UnRegister();
2453
2454   myPredicate = 0;
2455   FindBaseObjects();
2456 }
2457
2458 //=======================================================================
2459 // name    : Filter_i::SetPredicate
2460 // Purpose : Set predicate
2461 //=======================================================================
2462 void Filter_i::SetPredicate( Predicate_ptr thePredicate )
2463 {
2464   if ( myPredicate )
2465     myPredicate->UnRegister();
2466
2467   myPredicate = SMESH::GetPredicate(thePredicate);
2468
2469   if ( myPredicate )
2470   {
2471     myFilter.SetPredicate( myPredicate->GetPredicate() );
2472     myPredicate->Register();
2473     if ( const SMDS_Mesh* aMesh = MeshPtr2SMDSMesh(myMesh))
2474       myPredicate->GetPredicate()->SetMesh( aMesh );
2475     TPythonDump()<<this<<".SetPredicate("<<myPredicate<<")";
2476   }
2477   NotifyerAndWaiter::Modified();
2478 }
2479
2480 //=======================================================================
2481 // name    : Filter_i::GetElementType
2482 // Purpose : Get entity type
2483 //=======================================================================
2484 SMESH::ElementType Filter_i::GetElementType()
2485 {
2486   return myPredicate != 0 ? myPredicate->GetElementType() : SMESH::ALL;
2487 }
2488
2489 //=======================================================================
2490 // name    : Filter_i::SetMesh
2491 // Purpose : Set mesh
2492 //=======================================================================
2493 void
2494 Filter_i::
2495 SetMesh( SMESH_Mesh_ptr theMesh )
2496 {
2497   if(!CORBA::is_nil(theMesh))
2498     theMesh->Register();
2499
2500   if(!CORBA::is_nil(myMesh))
2501     myMesh->UnRegister();
2502
2503   myMesh = SMESH_Mesh::_duplicate( theMesh );
2504   TPythonDump()<<this<<".SetMesh("<<theMesh<<")";
2505
2506   if ( myPredicate )
2507     if ( const SMDS_Mesh* aMesh = MeshPtr2SMDSMesh(theMesh))
2508       myPredicate->GetPredicate()->SetMesh( aMesh );
2509 }
2510
2511 SMESH::long_array*
2512 Filter_i::
2513 GetIDs()
2514 {
2515   return GetElementsId(myMesh);
2516 }
2517
2518 //=======================================================================
2519 // name    : Filter_i::GetElementsId
2520 // Purpose : Get ids of entities
2521 //=======================================================================
2522 void
2523 Filter_i::
2524 GetElementsId( Predicate_i* thePredicate,
2525                const SMDS_Mesh* theMesh,
2526                Controls::Filter::TIdSequence& theSequence )
2527 {
2528   if (thePredicate)
2529     Controls::Filter::GetElementsId(theMesh,thePredicate->GetPredicate(),theSequence);
2530 }
2531
2532 void
2533 Filter_i::
2534 GetElementsId( Predicate_i* thePredicate,
2535                SMESH_Mesh_ptr theMesh,
2536                Controls::Filter::TIdSequence& theSequence )
2537 {
2538   if (thePredicate) 
2539     if(const SMDS_Mesh* aMesh = MeshPtr2SMDSMesh(theMesh))
2540       Controls::Filter::GetElementsId(aMesh,thePredicate->GetPredicate(),theSequence);
2541 }
2542
2543 SMESH::long_array*
2544 Filter_i::
2545 GetElementsId( SMESH_Mesh_ptr theMesh )
2546 {
2547   SMESH::long_array_var anArray = new SMESH::long_array;
2548   if(!CORBA::is_nil(theMesh) && myPredicate){
2549     theMesh->Load();
2550     Controls::Filter::TIdSequence aSequence;
2551     GetElementsId(myPredicate,theMesh,aSequence);
2552     long i = 0, iEnd = aSequence.size();
2553     anArray->length( iEnd );
2554     for ( ; i < iEnd; i++ )
2555       anArray[ i ] = aSequence[i];
2556   }
2557   return anArray._retn();
2558 }
2559
2560 //=============================================================================
2561 /*!
2562  * \brief Returns number of mesh elements per each \a EntityType
2563  */
2564 //=============================================================================
2565
2566 SMESH::long_array* ::Filter_i::GetMeshInfo()
2567 {
2568   SMESH::long_array_var aRes = new SMESH::long_array();
2569   aRes->length(SMESH::Entity_Last);
2570   for (int i = 0; i < SMESH::Entity_Last; i++)
2571     aRes[i] = 0;
2572
2573   if ( !CORBA::is_nil(myMesh) && myPredicate )
2574   {
2575     const SMDS_Mesh*  aMesh = MeshPtr2SMDSMesh(myMesh);
2576     SMDS_ElemIteratorPtr it = aMesh->elementsIterator( SMDSAbs_ElementType( GetElementType() ));
2577     while ( it->more() )
2578     {
2579       const SMDS_MeshElement* anElem = it->next();
2580       if ( myPredicate->IsSatisfy( anElem->GetID() ) )
2581         aRes[ anElem->GetEntityType() ]++;
2582     }
2583   }
2584
2585   return aRes._retn();  
2586 }
2587
2588 //=============================================================================
2589 /*!
2590  * \brief Returns number of mesh elements of each \a ElementType
2591  */
2592 //=============================================================================
2593
2594 SMESH::long_array* ::Filter_i::GetNbElementsByType()
2595 {
2596   SMESH::long_array_var aRes = new SMESH::long_array();
2597   aRes->length(SMESH::NB_ELEMENT_TYPES);
2598   for (int i = 0; i < SMESH::NB_ELEMENT_TYPES; i++)
2599     aRes[i] = 0;
2600
2601   if ( !CORBA::is_nil(myMesh) && myPredicate ) {
2602     const SMDS_Mesh*  aMesh = MeshPtr2SMDSMesh(myMesh);
2603     SMDS_ElemIteratorPtr it = aMesh->elementsIterator( SMDSAbs_ElementType( GetElementType() ));
2604     CORBA::Long& nbElems = aRes[ GetElementType() ];
2605     while ( it->more() )
2606     {
2607       const SMDS_MeshElement* anElem = it->next();
2608       if ( myPredicate->IsSatisfy( anElem->GetID() ) )
2609         nbElems++;
2610     }
2611   }
2612
2613   return aRes._retn();  
2614 }
2615
2616
2617 //================================================================================
2618 /*!
2619  * \brief Return GetElementType() within an array
2620  * Implement SMESH_IDSource interface
2621  */
2622 //================================================================================
2623
2624 SMESH::array_of_ElementType* Filter_i::GetTypes()
2625 {
2626   SMESH::array_of_ElementType_var types = new SMESH::array_of_ElementType;
2627
2628   // check if any element passes through the filter
2629   if ( !CORBA::is_nil(myMesh) && myPredicate )
2630   {
2631     const SMDS_Mesh* aMesh = MeshPtr2SMDSMesh(myMesh);
2632     SMDS_ElemIteratorPtr it = aMesh->elementsIterator( SMDSAbs_ElementType( GetElementType() ));
2633     bool satisfies = false;
2634     while ( !satisfies && it->more() )
2635       satisfies = myPredicate->IsSatisfy( it->next()->GetID() );
2636     if ( satisfies ) {
2637       types->length( 1 );
2638       types[0] = GetElementType();
2639     }
2640   }
2641   return types._retn();
2642 }
2643
2644 //=======================================================================
2645 //function : GetMesh
2646 //purpose  : Returns mesh
2647 //=======================================================================
2648
2649 SMESH::SMESH_Mesh_ptr Filter_i::GetMesh()
2650 {
2651   return SMESH_Mesh::_duplicate( myMesh );
2652 }
2653
2654 //=======================================================================
2655 //function : GetVtkUgStream
2656 //purpose  : Return data vtk unstructured grid (not implemented)
2657 //=======================================================================
2658
2659 SALOMEDS::TMPFile* Filter_i::GetVtkUgStream()
2660 {
2661   SALOMEDS::TMPFile_var SeqFile;
2662   return SeqFile._retn();
2663 }
2664 //=======================================================================
2665 // name    : getCriteria
2666 // Purpose : Retrieve criterions from predicate
2667 //=======================================================================
2668 static inline void getPrediacates( Predicate_i*                thePred,
2669                                    std::vector<Predicate_i*> & thePredVec )
2670 {
2671   const int aFType = thePred->GetFunctorType();
2672
2673   switch ( aFType )
2674   {
2675   case FT_LogicalNOT:
2676   {
2677     Predicate_i* aPred = ( dynamic_cast<LogicalNOT_i*>( thePred ) )->GetPredicate_i();
2678     getPrediacates( aPred, thePredVec );
2679     break;
2680   }
2681   case FT_LogicalAND:
2682   case FT_LogicalOR:
2683   {
2684     Predicate_i* aPred1 = ( dynamic_cast<LogicalBinary_i*>( thePred ) )->GetPredicate1_i();
2685     Predicate_i* aPred2 = ( dynamic_cast<LogicalBinary_i*>( thePred ) )->GetPredicate2_i();
2686     getPrediacates( aPred1, thePredVec );
2687     getPrediacates( aPred2, thePredVec );
2688     break;
2689    }
2690   default:;
2691   }
2692   thePredVec.push_back( thePred );
2693 }
2694
2695 //=======================================================================
2696 // name    : getCriteria
2697 // Purpose : Retrieve criterions from predicate
2698 //=======================================================================
2699 static inline bool getCriteria( Predicate_i*                thePred,
2700                                 SMESH::Filter::Criteria_out theCriteria )
2701 {
2702   const int aFType = thePred->GetFunctorType();
2703
2704   switch ( aFType )
2705   {
2706   case FT_LogicalNOT:
2707     {
2708       Predicate_i* aPred = ( dynamic_cast<LogicalNOT_i*>( thePred ) )->GetPredicate_i();
2709       getCriteria( aPred, theCriteria );
2710       theCriteria[ theCriteria->length() - 1 ].UnaryOp = FT_LogicalNOT;
2711     }
2712     return true;
2713
2714   case FT_LogicalAND:
2715   case FT_LogicalOR:
2716     {
2717       Predicate_i* aPred1 = ( dynamic_cast<LogicalBinary_i*>( thePred ) )->GetPredicate1_i();
2718       Predicate_i* aPred2 = ( dynamic_cast<LogicalBinary_i*>( thePred ) )->GetPredicate2_i();
2719       if ( !getCriteria( aPred1, theCriteria ) )
2720         return false;
2721       theCriteria[ theCriteria->length() - 1 ].BinaryOp = aFType;
2722       return getCriteria( aPred2, theCriteria );
2723     }
2724   case FT_Undefined:
2725     return false;
2726   }
2727
2728   // resize theCriteria
2729   CORBA::ULong i = theCriteria->length();
2730   theCriteria->length( i + 1 );
2731   theCriteria[ i ] = createCriterion();
2732
2733   // set members of the added Criterion
2734
2735   theCriteria[ i ].Type = aFType;
2736   theCriteria[ i ].TypeOfElement = thePred->GetElementType();
2737
2738   switch ( aFType )
2739   {
2740   case FT_FreeBorders:
2741   case FT_FreeEdges:
2742   case FT_FreeFaces:
2743   case FT_LinearOrQuadratic:
2744   case FT_FreeNodes:
2745   case FT_EqualEdges:
2746   case FT_EqualFaces:
2747   case FT_EqualVolumes:
2748   case FT_BadOrientedVolume:
2749   case FT_BareBorderVolume:
2750   case FT_BareBorderFace:
2751   case FT_OverConstrainedVolume:
2752   case FT_OverConstrainedFace:
2753     {
2754       return true;
2755     }
2756   case FT_BelongToMeshGroup:
2757     {
2758       BelongToMeshGroup_i* aPred = dynamic_cast<BelongToMeshGroup_i*>( thePred );
2759       SMESH::SMESH_GroupBase_var grp = aPred->GetGroup();
2760       if ( !grp->_is_nil() )
2761       {
2762         theCriteria[ i ].ThresholdStr = grp->GetName();
2763         theCriteria[ i ].ThresholdID  = aPred->GetGroupID().c_str();
2764       }
2765       return true;
2766     }
2767   case FT_BelongToGeom:
2768     {
2769       BelongToGeom_i* aPred = dynamic_cast<BelongToGeom_i*>( thePred );
2770       theCriteria[ i ].ThresholdStr  = aPred->GetShapeName();
2771       theCriteria[ i ].ThresholdID   = aPred->GetShapeID();
2772       theCriteria[ i ].Tolerance     = aPred->GetTolerance();
2773       return true;
2774     }
2775   case FT_BelongToPlane:
2776   case FT_BelongToCylinder:
2777   case FT_BelongToGenSurface:
2778     {
2779       BelongToSurface_i* aPred = dynamic_cast<BelongToSurface_i*>( thePred );
2780       theCriteria[ i ].ThresholdStr  = aPred->GetShapeName();
2781       theCriteria[ i ].ThresholdID   = aPred->GetShapeID();
2782       theCriteria[ i ].Tolerance     = aPred->GetTolerance();
2783       return true;
2784     }
2785   case FT_LyingOnGeom:
2786     {
2787       LyingOnGeom_i* aPred = dynamic_cast<LyingOnGeom_i*>( thePred );
2788       theCriteria[ i ].ThresholdStr  = aPred->GetShapeName();
2789       theCriteria[ i ].ThresholdID   = aPred->GetShapeID();
2790       theCriteria[ i ].Tolerance     = aPred->GetTolerance();
2791       return true;
2792     }
2793   case FT_CoplanarFaces:
2794     {
2795       CoplanarFaces_i* aPred = dynamic_cast<CoplanarFaces_i*>( thePred );
2796       theCriteria[ i ].ThresholdID   = aPred->GetFaceAsString();
2797       theCriteria[ i ].Tolerance     = aPred->GetTolerance();
2798       return true;
2799     }
2800   case FT_ConnectedElements:
2801     {
2802       ConnectedElements_i* aPred = dynamic_cast<ConnectedElements_i*>( thePred );
2803       SMESH::ConnectedElements::ThresholdType type;
2804       CORBA::String_var threshold = aPred->GetThreshold( type );
2805       switch ( type ) {
2806       case SMESH::ConnectedElements::POINT:
2807         theCriteria[ i ].ThresholdStr = threshold; break;
2808       case SMESH::ConnectedElements::VERTEX:
2809         theCriteria[ i ].ThresholdID = threshold; break;
2810       case SMESH::ConnectedElements::NODE:
2811         theCriteria[ i ].Threshold = atof( threshold.in() ); break;
2812       default:;
2813       }
2814       return true;
2815     }
2816   case FT_EqualNodes:
2817     {
2818       EqualNodes_i* aPred = dynamic_cast<EqualNodes_i*>( thePred );
2819       theCriteria[ i ].Tolerance  = aPred->GetTolerance();
2820       return true;
2821     }
2822   case FT_RangeOfIds:
2823     {
2824       RangeOfIds_i* aPred = dynamic_cast<RangeOfIds_i*>( thePred );
2825       theCriteria[ i ].ThresholdStr  = aPred->GetRangeStr();
2826       return true;
2827     }
2828   case FT_LessThan:
2829   case FT_MoreThan:
2830   case FT_EqualTo:
2831     {
2832       Comparator_i* aCompar = dynamic_cast<Comparator_i*>( thePred );
2833       theCriteria[ i ].Type      = aCompar->GetNumFunctor_i()->GetFunctorType();
2834       theCriteria[ i ].Compare   = aFType;
2835       theCriteria[ i ].Threshold = aCompar->GetMargin();
2836       if ( aFType == FT_EqualTo )
2837       {
2838         EqualTo_i* aCompar = dynamic_cast<EqualTo_i*>( thePred );
2839         theCriteria[ i ].Tolerance = aCompar->GetTolerance();
2840       }
2841       return true;
2842     }
2843   case FT_GroupColor:
2844     {
2845       GroupColor_i* aPred = dynamic_cast<GroupColor_i*>( thePred );
2846       theCriteria[ i ].ThresholdStr = aPred->GetColorStr();
2847       return true;
2848     }
2849   case FT_ElemGeomType:
2850     {
2851       ElemGeomType_i* aPred = dynamic_cast<ElemGeomType_i*>( thePred );
2852       theCriteria[ i ].Threshold = (double)aPred->GetGeometryType();
2853       return true;
2854     }
2855   case FT_EntityType:
2856     {
2857       ElemEntityType_i* aPred = dynamic_cast<ElemEntityType_i*>( thePred );
2858       theCriteria[ i ].Threshold = (double)aPred->GetEntityType();
2859       return true;
2860     }
2861   default:
2862     return false;
2863   }
2864 }
2865
2866 //=======================================================================
2867 // name    : Filter_i::GetCriteria
2868 // Purpose : Retrieve criterions from predicate
2869 //=======================================================================
2870 CORBA::Boolean Filter_i::GetCriteria( SMESH::Filter::Criteria_out theCriteria )
2871 {
2872   theCriteria = new SMESH::Filter::Criteria;
2873   return myPredicate != 0 ? getCriteria( myPredicate, theCriteria ) : true;
2874 }
2875
2876 //=======================================================================
2877 // name    : Filter_i::SetCriteria
2878 // Purpose : Create new predicate and set criterions in it
2879 //=======================================================================
2880 CORBA::Boolean Filter_i::SetCriteria( const SMESH::Filter::Criteria& theCriteria )
2881 {
2882   SetPredicate( SMESH::Predicate::_nil() );
2883
2884   SMESH::FilterManager_i* aFilter = new SMESH::FilterManager_i();
2885   FilterManager_ptr aFilterMgr = aFilter->_this();
2886
2887   // CREATE two lists ( PREDICATES  and LOG OP )
2888
2889   // Criterion
2890   TPythonDump()<<"aCriteria = []";
2891   std::list<SMESH::Predicate_ptr> aPredicates;
2892   std::list<int>                  aBinaries;
2893   for ( int i = 0, n = theCriteria.length(); i < n; i++ )
2894   {
2895     int         aCriterion    = theCriteria[ i ].Type;
2896     int         aCompare      = theCriteria[ i ].Compare;
2897     double      aThreshold    = theCriteria[ i ].Threshold;
2898     const char* aThresholdStr = theCriteria[ i ].ThresholdStr;
2899     const char* aThresholdID  = theCriteria[ i ].ThresholdID;
2900     int         aUnary        = theCriteria[ i ].UnaryOp;
2901     int         aBinary       = theCriteria[ i ].BinaryOp;
2902     double      aTolerance    = theCriteria[ i ].Tolerance;
2903     ElementType aTypeOfElem   = theCriteria[ i ].TypeOfElement;
2904     long        aPrecision    = theCriteria[ i ].Precision;
2905
2906     {
2907       TPythonDump pd;
2908       pd << "aCriterion = SMESH.Filter.Criterion("
2909          << aCriterion    << ", "
2910          << aCompare      << ", "
2911          << aThreshold    << ", '"
2912          << aThresholdStr << "', '"
2913          << aThresholdID  << "', "
2914          << aUnary        << ", "
2915          << aBinary       << ", "
2916          << aTolerance    << ", "
2917          << aTypeOfElem   << ", "
2918          << aPrecision    << ")";
2919     }
2920     TPythonDump pd;
2921
2922     SMESH::Predicate_ptr aPredicate = SMESH::Predicate::_nil();
2923     SMESH::NumericalFunctor_ptr aFunctor = SMESH::NumericalFunctor::_nil();
2924
2925     switch ( aCriterion )
2926     {
2927       // Functors
2928
2929       case SMESH::FT_MultiConnection:
2930         aFunctor = aFilterMgr->CreateMultiConnection();
2931         break;
2932       case SMESH::FT_MultiConnection2D:
2933         aFunctor = aFilterMgr->CreateMultiConnection2D();
2934         break;
2935       case SMESH::FT_Length:
2936         aFunctor = aFilterMgr->CreateLength();
2937         break;
2938       case SMESH::FT_Length2D:
2939         aFunctor = aFilterMgr->CreateLength2D();
2940         break;
2941       case SMESH::FT_AspectRatio:
2942         aFunctor = aFilterMgr->CreateAspectRatio();
2943         break;
2944       case SMESH::FT_AspectRatio3D:
2945         aFunctor = aFilterMgr->CreateAspectRatio3D();
2946         break;
2947       case SMESH::FT_Warping:
2948         aFunctor = aFilterMgr->CreateWarping();
2949         break;
2950       case SMESH::FT_MinimumAngle:
2951         aFunctor = aFilterMgr->CreateMinimumAngle();
2952         break;
2953       case SMESH::FT_Taper:
2954         aFunctor = aFilterMgr->CreateTaper();
2955         break;
2956       case SMESH::FT_Skew:
2957         aFunctor = aFilterMgr->CreateSkew();
2958         break;
2959       case SMESH::FT_Area:
2960         aFunctor = aFilterMgr->CreateArea();
2961         break;
2962       case SMESH::FT_Volume3D:
2963         aFunctor = aFilterMgr->CreateVolume3D();
2964         break;
2965       case SMESH::FT_MaxElementLength2D:
2966         aFunctor = aFilterMgr->CreateMaxElementLength2D();
2967         break;
2968       case SMESH::FT_MaxElementLength3D:
2969         aFunctor = aFilterMgr->CreateMaxElementLength3D();
2970         break;
2971       case SMESH::FT_BallDiameter:
2972         aFunctor = aFilterMgr->CreateBallDiameter();
2973         break;
2974
2975       // Predicates
2976
2977       case SMESH::FT_FreeBorders:
2978         aPredicate = aFilterMgr->CreateFreeBorders();
2979         break;
2980       case SMESH::FT_FreeEdges:
2981         aPredicate = aFilterMgr->CreateFreeEdges();
2982         break;
2983       case SMESH::FT_FreeFaces:
2984         aPredicate = aFilterMgr->CreateFreeFaces();
2985         break;
2986       case SMESH::FT_FreeNodes:
2987         aPredicate = aFilterMgr->CreateFreeNodes();
2988         break;
2989       case SMESH::FT_EqualNodes:
2990         {
2991           SMESH::EqualNodes_ptr pred = aFilterMgr->CreateEqualNodes();
2992           pred->SetTolerance( aTolerance );
2993           aPredicate = pred;
2994           break;
2995         }
2996       case SMESH::FT_EqualEdges:
2997         aPredicate = aFilterMgr->CreateEqualEdges();
2998         break;
2999       case SMESH::FT_EqualFaces:
3000         aPredicate = aFilterMgr->CreateEqualFaces();
3001         break;
3002       case SMESH::FT_EqualVolumes:
3003         aPredicate = aFilterMgr->CreateEqualVolumes();
3004         break;
3005       case SMESH::FT_BelongToMeshGroup:
3006         {
3007           SMESH::BelongToMeshGroup_ptr tmpPred = aFilterMgr->CreateBelongToMeshGroup();
3008           tmpPred->SetGroupID( aThresholdID );
3009           aPredicate = tmpPred;
3010         }
3011         break;
3012       case SMESH::FT_BelongToGeom:
3013         {
3014           SMESH::BelongToGeom_ptr tmpPred = aFilterMgr->CreateBelongToGeom();
3015           tmpPred->SetElementType( aTypeOfElem );
3016           tmpPred->SetShape( aThresholdID, aThresholdStr );
3017           tmpPred->SetTolerance( aTolerance );
3018           aPredicate = tmpPred;
3019         }
3020         break;
3021       case SMESH::FT_BelongToPlane:
3022       case SMESH::FT_BelongToCylinder:
3023       case SMESH::FT_BelongToGenSurface:
3024         {
3025           SMESH::BelongToSurface_ptr tmpPred;
3026           switch ( aCriterion ) {
3027           case SMESH::FT_BelongToPlane:
3028             tmpPred = aFilterMgr->CreateBelongToPlane(); break;
3029           case SMESH::FT_BelongToCylinder:
3030             tmpPred = aFilterMgr->CreateBelongToCylinder(); break;
3031           default:
3032             tmpPred = aFilterMgr->CreateBelongToGenSurface();
3033           }
3034           tmpPred->SetShape( aThresholdID, aThresholdStr, aTypeOfElem );
3035           tmpPred->SetTolerance( aTolerance );
3036           aPredicate = tmpPred;
3037         }
3038         break;
3039       case SMESH::FT_LyingOnGeom:
3040         {
3041           SMESH::LyingOnGeom_ptr tmpPred = aFilterMgr->CreateLyingOnGeom();
3042           tmpPred->SetElementType( aTypeOfElem );
3043           tmpPred->SetShape( aThresholdID, aThresholdStr );
3044           tmpPred->SetTolerance( aTolerance );
3045           aPredicate = tmpPred;
3046         }
3047         break;
3048       case SMESH::FT_RangeOfIds:
3049         {
3050           SMESH::RangeOfIds_ptr tmpPred = aFilterMgr->CreateRangeOfIds();
3051           tmpPred->SetRangeStr( aThresholdStr );
3052           tmpPred->SetElementType( aTypeOfElem );
3053           aPredicate = tmpPred;
3054         }
3055         break;
3056       case SMESH::FT_BadOrientedVolume:
3057         {
3058           aPredicate = aFilterMgr->CreateBadOrientedVolume();
3059         }
3060         break;
3061       case SMESH::FT_BareBorderVolume:
3062         {
3063           aPredicate = aFilterMgr->CreateBareBorderVolume();
3064         }
3065         break;
3066       case SMESH::FT_BareBorderFace:
3067         {
3068           aPredicate = aFilterMgr->CreateBareBorderFace();
3069         }
3070         break;
3071       case SMESH::FT_OverConstrainedVolume:
3072         {
3073           aPredicate = aFilterMgr->CreateOverConstrainedVolume();
3074         }
3075         break;
3076       case SMESH::FT_OverConstrainedFace:
3077         {
3078           aPredicate = aFilterMgr->CreateOverConstrainedFace();
3079         }
3080         break;
3081       case SMESH::FT_LinearOrQuadratic:
3082         {
3083           SMESH::LinearOrQuadratic_ptr tmpPred = aFilterMgr->CreateLinearOrQuadratic();
3084           tmpPred->SetElementType( aTypeOfElem );
3085           aPredicate = tmpPred;
3086           break;
3087         }
3088       case SMESH::FT_GroupColor:
3089         {
3090           SMESH::GroupColor_ptr tmpPred = aFilterMgr->CreateGroupColor();
3091           tmpPred->SetElementType( aTypeOfElem );
3092           tmpPred->SetColorStr( aThresholdStr );
3093           aPredicate = tmpPred;
3094           break;
3095         }
3096       case SMESH::FT_ElemGeomType:
3097         {
3098           SMESH::ElemGeomType_ptr tmpPred = aFilterMgr->CreateElemGeomType();
3099           tmpPred->SetElementType( aTypeOfElem );
3100           tmpPred->SetGeometryType( (GeometryType)(int)(aThreshold + 0.5) );
3101           aPredicate = tmpPred;
3102           break;
3103         }
3104       case SMESH::FT_EntityType:
3105         {
3106           SMESH::ElemEntityType_ptr tmpPred = aFilterMgr->CreateElemEntityType();
3107           tmpPred->SetElementType( aTypeOfElem );
3108           tmpPred->SetEntityType( EntityType( (int (aThreshold + 0.5))));
3109           aPredicate = tmpPred;
3110           break;
3111         }
3112       case SMESH::FT_CoplanarFaces:
3113         {
3114           SMESH::CoplanarFaces_ptr tmpPred = aFilterMgr->CreateCoplanarFaces();
3115           tmpPred->SetFace( atol (aThresholdID ));
3116           tmpPred->SetTolerance( aTolerance );
3117           aPredicate = tmpPred;
3118           break;
3119         }
3120       case SMESH::FT_ConnectedElements:
3121         {
3122           SMESH::ConnectedElements_ptr tmpPred = aFilterMgr->CreateConnectedElements();
3123           if ( strlen( aThresholdID ) > 0 ) // shape ID
3124             tmpPred->SetThreshold( aThresholdID, SMESH::ConnectedElements::VERTEX );
3125           else if ( strlen( aThresholdStr ) > 0 ) // point coords
3126             tmpPred->SetThreshold( aThresholdStr, SMESH::ConnectedElements::POINT );
3127           else if ( aThreshold >= 1 )
3128             tmpPred->SetNode( (CORBA::Long) aThreshold ); // node ID
3129           tmpPred->SetElementType( aTypeOfElem );
3130           aPredicate = tmpPred;
3131           break;
3132         }
3133
3134       default:
3135         continue;
3136     }
3137
3138     // Comparator
3139     if ( !aFunctor->_is_nil() && aPredicate->_is_nil() )
3140     {
3141       SMESH::Comparator_ptr aComparator = SMESH::Comparator::_nil();
3142
3143       if ( aCompare == SMESH::FT_LessThan )
3144         aComparator = aFilterMgr->CreateLessThan();
3145       else if ( aCompare == SMESH::FT_MoreThan )
3146         aComparator = aFilterMgr->CreateMoreThan();
3147       else if ( aCompare == SMESH::FT_EqualTo )
3148         aComparator = aFilterMgr->CreateEqualTo();
3149       else
3150         continue;
3151
3152       aComparator->SetNumFunctor( aFunctor );
3153       aComparator->SetMargin( aThreshold );
3154
3155       if ( aCompare == FT_EqualTo )
3156       {
3157         SMESH::EqualTo_var anEqualTo = SMESH::EqualTo::_narrow( aComparator );
3158         anEqualTo->SetTolerance( aTolerance );
3159       }
3160
3161       aPredicate = aComparator;
3162
3163       aFunctor->SetPrecision( aPrecision );
3164     }
3165
3166     // Logical not
3167     if ( aUnary == FT_LogicalNOT )
3168     {
3169       SMESH::LogicalNOT_ptr aNotPred = aFilterMgr->CreateLogicalNOT();
3170       aNotPred->SetPredicate( aPredicate );
3171       aPredicate = aNotPred;
3172     }
3173
3174     // logical op
3175     aPredicates.push_back( aPredicate );
3176     aBinaries.push_back( aBinary );
3177     pd <<"aCriteria.append(aCriterion)";
3178
3179   } // end of for
3180   TPythonDump pd; pd<<this<<".SetCriteria(aCriteria)";
3181
3182   // CREATE ONE PREDICATE FROM PREVIOUSLY CREATED MAP
3183
3184   // combine all "AND" operations
3185
3186   std::list<SMESH::Predicate_ptr> aResList;
3187
3188   std::list<SMESH::Predicate_ptr>::iterator aPredIter;
3189   std::list<int>::iterator                  aBinaryIter;
3190
3191   SMESH::Predicate_ptr aPrevPredicate = SMESH::Predicate::_nil();
3192   int aPrevBinary = SMESH::FT_Undefined;
3193   if ( !aBinaries.empty() )
3194     aBinaries.back() = SMESH::FT_Undefined;
3195
3196   for ( aPredIter = aPredicates.begin(), aBinaryIter = aBinaries.begin();
3197         aPredIter != aPredicates.end() && aBinaryIter != aBinaries.end();
3198         ++aPredIter, ++aBinaryIter )
3199   {
3200     int aCurrBinary = *aBinaryIter;
3201
3202     SMESH::Predicate_ptr aCurrPred = SMESH::Predicate::_nil();
3203
3204     if ( aPrevBinary == SMESH::FT_LogicalAND )
3205     {
3206
3207       SMESH::LogicalBinary_ptr aBinaryPred = aFilterMgr->CreateLogicalAND();
3208       aBinaryPred->SetPredicate1( aPrevPredicate );
3209       aBinaryPred->SetPredicate2( *aPredIter );
3210       aCurrPred = aBinaryPred;
3211     }
3212     else
3213       aCurrPred = *aPredIter;
3214
3215     if ( aCurrBinary != SMESH::FT_LogicalAND )
3216       aResList.push_back( aCurrPred );
3217
3218     aPrevPredicate = aCurrPred;
3219     aPrevBinary = aCurrBinary;
3220   }
3221
3222   // combine all "OR" operations
3223
3224   SMESH::Predicate_ptr aResPredicate = SMESH::Predicate::_nil();
3225
3226   if ( aResList.size() == 1 )
3227     aResPredicate = *aResList.begin();
3228   else if ( aResList.size() > 1 )
3229   {
3230     std::list<SMESH::Predicate_ptr>::iterator anIter = aResList.begin();
3231     aResPredicate = *anIter;
3232     anIter++;
3233     for ( ; anIter != aResList.end(); ++anIter )
3234     {
3235       SMESH::LogicalBinary_ptr aBinaryPred = aFilterMgr->CreateLogicalOR();
3236       aBinaryPred->SetPredicate1( aResPredicate );
3237       aBinaryPred->SetPredicate2( *anIter );
3238       aResPredicate = aBinaryPred;
3239     }
3240   }
3241
3242   SetPredicate( aResPredicate );
3243   if ( !aResPredicate->_is_nil() )
3244     aResPredicate->UnRegister();
3245
3246   return !aResPredicate->_is_nil();
3247 }
3248
3249 //=======================================================================
3250 // name    : Filter_i::GetPredicate_i
3251 // Purpose : Get implementation of predicate
3252 //=======================================================================
3253 Predicate_i* Filter_i::GetPredicate_i()
3254 {
3255   return myPredicate;
3256 }
3257
3258 //=======================================================================
3259 // name    : Filter_i::GetPredicate
3260 // Purpose : Get predicate
3261 //=======================================================================
3262 Predicate_ptr Filter_i::GetPredicate()
3263 {
3264   if ( myPredicate == 0 )
3265     return SMESH::Predicate::_nil();
3266   else
3267   {
3268     SMESH::Predicate_var anObj = myPredicate->_this();
3269     // if ( SMESH::Functor_i* fun = SMESH::DownCast<SMESH::Functor_i*>( anObj ))
3270     //   TPythonDump() << fun << " = " << this << ".GetPredicate()";
3271     return anObj._retn();
3272   }
3273 }
3274
3275 //================================================================================
3276 /*!
3277  * \brief Find groups it depends on
3278  */
3279 //================================================================================
3280
3281 void Filter_i::FindBaseObjects()
3282 {
3283   // release current groups
3284   for ( size_t i = 0; i < myBaseGroups.size(); ++i )
3285     if ( myBaseGroups[i] )
3286     {
3287       myBaseGroups[i]->RemoveModifWaiter( this );
3288       myBaseGroups[i]->UnRegister();
3289     }
3290
3291   // remember new groups
3292   myBaseGroups.clear();
3293   if ( myPredicate )
3294   {
3295     std::vector<Predicate_i*> predicates;
3296     getPrediacates( myPredicate, predicates );
3297     for ( size_t i = 0; i < predicates.size(); ++i )
3298       if ( BelongToMeshGroup_i* bmg = dynamic_cast< BelongToMeshGroup_i* >( predicates[i] ))
3299       {
3300         SMESH::SMESH_GroupBase_var g = bmg->GetGroup();
3301         SMESH_GroupBase_i* g_i = SMESH::DownCast< SMESH_GroupBase_i*>( g );
3302         if ( g_i )
3303         {
3304           g_i->AddModifWaiter( this );
3305           g_i->Register();
3306           myBaseGroups.push_back( g_i );
3307         }
3308       }
3309   }
3310 }
3311
3312 //================================================================================
3313 /*!
3314  * \brief When notified on removal of myBaseGroups[i], remove a reference to a
3315  *        group from a predicate
3316  */
3317 //================================================================================
3318
3319 void Filter_i::OnBaseObjModified(NotifyerAndWaiter* group, bool removed)
3320 {
3321   if ( !removed )
3322     return; // a GroupOnFilter holding this filter is notified automatically
3323
3324   if ( myPredicate )
3325   {
3326     std::vector<Predicate_i*> predicates;
3327     getPrediacates( myPredicate, predicates );
3328     for ( size_t i = 0; i < predicates.size(); ++i )
3329       if ( BelongToMeshGroup_i* bmg = dynamic_cast< BelongToMeshGroup_i* >( predicates[i] ))
3330       {
3331         SMESH::SMESH_GroupBase_var g = bmg->GetGroup();
3332         SMESH_GroupBase_i* g_i = SMESH::DownCast< SMESH_GroupBase_i*>( g );
3333         if ( g_i == group )
3334         {
3335           bmg->SetGroup( SMESH::SMESH_GroupBase::_nil() );
3336           bmg->SetGroupID( "" );
3337         }
3338       }
3339   }
3340
3341   FindBaseObjects(); // release and update myBaseGroups;
3342 }
3343
3344 /*
3345                             FILTER LIBRARY
3346 */
3347
3348 #define ATTR_TYPE          "type"
3349 #define ATTR_COMPARE       "compare"
3350 #define ATTR_THRESHOLD     "threshold"
3351 #define ATTR_UNARY         "unary"
3352 #define ATTR_BINARY        "binary"