src/SMESH_I/SMESH_Filter_i.cxx

   1 // Copyright (C) 2007-2016  CEA/DEN, EDF R&D, OPEN CASCADE
   2 //
   3 // Copyright (C) 2003-2007  OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
   4 // CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
   5 //
   6 // This library is free software; you can redistribute it and/or
   7 // modify it under the terms of the GNU Lesser General Public
   8 // License as published by the Free Software Foundation; either
   9 // version 2.1 of the License, or (at your option) any later version.
  10 //
  11 // This library is distributed in the hope that it will be useful,
  12 // but WITHOUT ANY WARRANTY; without even the implied warranty of
  13 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  14 // Lesser General Public License for more details.
  15 //
  16 // You should have received a copy of the GNU Lesser General Public
  17 // License along with this library; if not, write to the Free Software
  18 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
  19 //
  20 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
  21 //
  22
  23 //  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 <Basics_OCCTVersion.hxx>
  44
  45 #include <BRep_Tool.hxx>
  46 #include <Geom_CylindricalSurface.hxx>
  47 #include <Geom_Plane.hxx>
  48 #include <LDOMParser.hxx>
  49 #include <LDOMString.hxx>
  50 #include <LDOM_Document.hxx>
  51 #include <LDOM_Element.hxx>
  52 #include <LDOM_Node.hxx>
  53 #include <LDOM_XmlWriter.hxx>
  54 #include <Precision.hxx>
  55 #include <TColStd_ListIteratorOfListOfInteger.hxx>
  56 #include <TColStd_ListIteratorOfListOfReal.hxx>
  57 #include <TColStd_ListOfInteger.hxx>
  58 #include <TColStd_ListOfReal.hxx>
  59 #include <TColStd_SequenceOfHAsciiString.hxx>
  60 #include <TCollection_HAsciiString.hxx>
  61 #include <TopExp_Explorer.hxx>
  62 #include <TopoDS.hxx>
  63 #include <TopoDS_Shape.hxx>
  64
  65 using namespace SMESH;
  66 using namespace SMESH::Controls;
  67
  68
  69 namespace SMESH
  70 {
  71   Predicate_i*
  72   GetPredicate( Predicate_ptr thePredicate )
  73   {
  74     return DownCast<Predicate_i*>(thePredicate);
  75   }
  76 }
  77
  78 /*
  79                             AUXILIARY METHODS
  80 */
  81
  82 inline
  83 const SMDS_Mesh*
  84 MeshPtr2SMDSMesh( SMESH_Mesh_ptr theMesh )
  85 {
  86   SMESH_Mesh_i* anImplPtr = DownCast<SMESH_Mesh_i*>(theMesh);
  87   return anImplPtr ? anImplPtr->GetImpl().GetMeshDS() : 0;
  88 }
  89
  90 inline
  91 SMESH::long_array*
  92 toArray( const TColStd_ListOfInteger& aList )
  93 {
  94   SMESH::long_array_var anArray = new SMESH::long_array;
  95   anArray->length( aList.Extent() );
  96   TColStd_ListIteratorOfListOfInteger anIter( aList );
  97   int i = 0;
  98   for( ; anIter.More(); anIter.Next() )
  99     anArray[ i++ ] = anIter.Value();
 100
 101   return anArray._retn();
 102 }
 103
 104 inline
 105 SMESH::double_array*
 106 toArray( const TColStd_ListOfReal& aList )
 107 {
 108   SMESH::double_array_var anArray = new SMESH::double_array;
 109   anArray->length( aList.Extent() );
 110   TColStd_ListIteratorOfListOfReal anIter( aList );
 111   int i = 0;
 112   for( ; anIter.More(); anIter.Next() )
 113     anArray[ i++ ] = anIter.Value();
 114
 115   return anArray._retn();
 116 }
 117
 118 static SMESH::Filter::Criterion createCriterion()
 119 {
 120   SMESH::Filter::Criterion aCriterion;
 121
 122   aCriterion.Type          = FT_Undefined;
 123   aCriterion.Compare       = FT_Undefined;
 124   aCriterion.Threshold     = 0;
 125   aCriterion.UnaryOp       = FT_Undefined;
 126   aCriterion.BinaryOp      = FT_Undefined;
 127   aCriterion.ThresholdStr  = "";
 128   aCriterion.ThresholdID   = "";
 129   aCriterion.Tolerance     = Precision::Confusion();
 130   aCriterion.TypeOfElement = SMESH::ALL;
 131   aCriterion.Precision     = -1;
 132
 133   return aCriterion;
 134 }
 135
 136 static TopoDS_Shape getShapeByName( const char* theName )
 137 {
 138   if ( theName != 0 )
 139   {
 140     SMESH_Gen_i* aSMESHGen     = SMESH_Gen_i::GetSMESHGen();
 141     SALOMEDS::Study_var aStudy = aSMESHGen->GetCurrentStudy();
 142     if ( !aStudy->_is_nil() )
 143     {
 144       SALOMEDS::Study::ListOfSObject_var aList = aStudy->FindObjectByName( theName, "GEOM" );
 145       if ( aList->length() > 0 )
 146       {
 147         CORBA::Object_var        anObj = aList[ 0 ]->GetObject();
 148         GEOM::GEOM_Object_var aGeomObj = GEOM::GEOM_Object::_narrow( anObj );
 149         TopoDS_Shape             shape = aSMESHGen->GeomObjectToShape( aGeomObj );
 150         SALOME::UnRegister( aList ); // UnRegister() objects in aList
 151         return shape;
 152       }
 153     }
 154   }
 155   return TopoDS_Shape();
 156 }
 157
 158 static TopoDS_Shape getShapeByID (const char* theID)
 159 {
 160   if ( theID && strlen( theID ) > 0 ) {
 161     SMESH_Gen_i*     aSMESHGen = SMESH_Gen_i::GetSMESHGen();
 162     SALOMEDS::Study_var aStudy = aSMESHGen->GetCurrentStudy();
 163     if ( !aStudy->_is_nil() ) {
 164       SALOMEDS::SObject_wrap aSObj = aStudy->FindObjectID(theID);
 165       if ( !aSObj->_is_nil() ) {
 166         CORBA::Object_var          obj = aSObj->GetObject();
 167         GEOM::GEOM_Object_var aGeomObj = GEOM::GEOM_Object::_narrow(obj);
 168         return aSMESHGen->GeomObjectToShape( aGeomObj );
 169       }
 170     }
 171   }
 172   return TopoDS_Shape();
 173 }
 174
 175 // static std::string getShapeNameByID (const char* theID)
 176 // {
 177 //   if ( theID && strlen( theID ) > 0 ) {
 178 //     SMESH_Gen_i*     aSMESHGen = SMESH_Gen_i::GetSMESHGen();
 179 //     SALOMEDS::Study_var aStudy = aSMESHGen->GetCurrentStudy();
 180 //     if ( !aStudy->_is_nil() ) {
 181 //       SALOMEDS::SObject_wrap aSObj = aStudy->FindObjectID(theID);
 182 //       if ( !aSObj->_is_nil() ) {
 183 //         CORBA::String_var name = aSObj->GetName();
 184 //         return name.in();
 185 //       }
 186 //     }
 187 //   }
 188 //   return "";
 189 // }
 190
 191 /*
 192                                 FUNCTORS
 193 */
 194
 195 /*
 196   Class       : Functor_i
 197   Description : An abstact class for all functors
 198 */
 199 Functor_i::Functor_i():
 200   SALOME::GenericObj_i( SMESH_Gen_i::GetPOA() )
 201 {
 202   //Base class Salome_GenericObject do it inmplicitly by overriding PortableServer::POA_ptr _default_POA() method
 203   //PortableServer::ObjectId_var anObjectId =
 204   //  SMESH_Gen_i::GetPOA()->activate_object( this );
 205 }
 206
 207 Functor_i::~Functor_i()
 208 {
 209   //TPythonDump()<<this<<".UnRegister()";
 210 }
 211
 212 void Functor_i::SetMesh( SMESH_Mesh_ptr theMesh )
 213 {
 214   myFunctorPtr->SetMesh( MeshPtr2SMDSMesh( theMesh ) );
 215   TPythonDump()<<this<<".SetMesh("<<theMesh<<")";
 216 }
 217
 218 ElementType Functor_i::GetElementType()
 219 {
 220   return ( ElementType )myFunctorPtr->GetType();
 221 }
 222
 223
 224 /*
 225   Class       : NumericalFunctor_i
 226   Description : Base class for numerical functors
 227 */
 228 CORBA::Double NumericalFunctor_i::GetValue( CORBA::Long theId )
 229 {
 230   return myNumericalFunctorPtr->GetValue( theId );
 231 }
 232
 233 SMESH::Histogram* NumericalFunctor_i::GetHistogram(CORBA::Short nbIntervals, CORBA::Boolean isLogarithmic)
 234 {
 235   std::vector<int> nbEvents;
 236   std::vector<double> funValues;
 237   std::vector<int> elements;
 238   myNumericalFunctorPtr->GetHistogram(nbIntervals,nbEvents,funValues,elements,0,isLogarithmic);
 239
 240   SMESH::Histogram_var histogram = new SMESH::Histogram;
 241
 242   nbIntervals = CORBA::Short( Min( int( nbEvents.size()),
 243                                    int( funValues.size() - 1 )));
 244   if ( nbIntervals > 0 )
 245   {
 246     histogram->length( nbIntervals );
 247     for ( int i = 0; i < nbIntervals; ++i )
 248     {
 249       HistogramRectangle& rect = histogram[i];
 250       rect.nbEvents = nbEvents[i];
 251       rect.min = funValues[i];
 252       rect.max = funValues[i+1];
 253     }
 254   }
 255   return histogram._retn();
 256 }
 257
 258 SMESH::Histogram* NumericalFunctor_i::GetLocalHistogram(CORBA::Short              nbIntervals,
 259                                                         CORBA::Boolean            isLogarithmic,
 260                                                         SMESH::SMESH_IDSource_ptr object)
 261 {
 262   SMESH::Histogram_var histogram = new SMESH::Histogram;
 263
 264   std::vector<int>    nbEvents;
 265   std::vector<double> funValues;
 266   std::vector<int>    elements;
 267
 268   SMDS_ElemIteratorPtr elemIt;
 269   if ( SMESH::DownCast< SMESH_GroupOnFilter_i* >( object ) ||
 270        SMESH::DownCast< SMESH::Filter_i* >( object ))
 271   {
 272     elemIt = SMESH_Mesh_i::GetElements( object, GetElementType() );
 273     if ( !elemIt ) return histogram._retn();
 274   }
 275   else
 276   {
 277     SMESH::SMESH_Mesh_var        mesh = object->GetMesh();
 278     SMESH::long_array_var  objNbElems = object->GetNbElementsByType();
 279     SMESH::long_array_var meshNbElems = mesh->  GetNbElementsByType();
 280     if ( meshNbElems[ GetElementType() ] !=
 281          objNbElems [ GetElementType() ] )
 282     {
 283       elements.reserve( objNbElems[ GetElementType() ]);
 284       elemIt = SMESH_Mesh_i::GetElements( object, GetElementType() );
 285     }
 286   }
 287   if ( elemIt )
 288   {
 289     while ( elemIt->more() )
 290       elements.push_back( elemIt->next()->GetID() );
 291     if ( elements.empty() ) return histogram._retn();
 292   }
 293
 294   myNumericalFunctorPtr->GetHistogram(nbIntervals,nbEvents,funValues,elements,0,isLogarithmic);
 295
 296   nbIntervals = CORBA::Short( Min( int( nbEvents.size()),
 297                                    int( funValues.size() - 1 )));
 298   if ( nbIntervals > 0 )
 299   {
 300     histogram->length( nbIntervals );
 301     for ( int i = 0; i < nbIntervals; ++i )
 302     {
 303       HistogramRectangle& rect = histogram[i];
 304       rect.nbEvents = nbEvents[i];
 305       rect.min = funValues[i];
 306       rect.max = funValues[i+1];
 307     }
 308   }
 309   return histogram._retn();
 310 }
 311
 312 void NumericalFunctor_i::SetPrecision( CORBA::Long thePrecision )
 313 {
 314   myNumericalFunctorPtr->SetPrecision( thePrecision );
 315   TPythonDump()<<this<<".SetPrecision("<<thePrecision<<")";
 316 }
 317
 318 CORBA::Long NumericalFunctor_i::GetPrecision()
 319 {
 320  return myNumericalFunctorPtr->GetPrecision();
 321 }
 322
 323 Controls::NumericalFunctorPtr NumericalFunctor_i::GetNumericalFunctor()
 324 {
 325   return myNumericalFunctorPtr;
 326 }
 327
 328
 329 /*
 330   Class       : SMESH_MinimumAngle
 331   Description : Functor for calculation of minimum angle
 332 */
 333 MinimumAngle_i::MinimumAngle_i()
 334 {
 335   myNumericalFunctorPtr.reset( new Controls::MinimumAngle() );
 336   myFunctorPtr = myNumericalFunctorPtr;
 337 }
 338
 339 FunctorType MinimumAngle_i::GetFunctorType()
 340 {
 341   return SMESH::FT_MinimumAngle;
 342 }
 343
 344
 345 /*
 346   Class       : AspectRatio
 347   Description : Functor for calculating aspect ratio
 348 */
 349 AspectRatio_i::AspectRatio_i()
 350 {
 351   myNumericalFunctorPtr.reset( new Controls::AspectRatio() );
 352   myFunctorPtr = myNumericalFunctorPtr;
 353 }
 354
 355 FunctorType AspectRatio_i::GetFunctorType()
 356 {
 357   return SMESH::FT_AspectRatio;
 358 }
 359
 360
 361 /*
 362   Class       : AspectRatio3D
 363   Description : Functor for calculating aspect ratio 3D
 364 */
 365 AspectRatio3D_i::AspectRatio3D_i()
 366 {
 367   myNumericalFunctorPtr.reset( new Controls::AspectRatio3D() );
 368   myFunctorPtr = myNumericalFunctorPtr;
 369 }
 370
 371 FunctorType AspectRatio3D_i::GetFunctorType()
 372 {
 373   return SMESH::FT_AspectRatio3D;
 374 }
 375
 376
 377 /*
 378   Class       : Warping_i
 379   Description : Functor for calculating warping
 380 */
 381 Warping_i::Warping_i()
 382 {
 383   myNumericalFunctorPtr.reset( new Controls::Warping() );
 384   myFunctorPtr = myNumericalFunctorPtr;
 385 }
 386
 387 FunctorType Warping_i::GetFunctorType()
 388 {
 389   return SMESH::FT_Warping;
 390 }
 391
 392
 393 /*
 394   Class       : Taper_i
 395   Description : Functor for calculating taper
 396 */
 397 Taper_i::Taper_i()
 398 {
 399   myNumericalFunctorPtr.reset( new Controls::Taper() );
 400   myFunctorPtr = myNumericalFunctorPtr;
 401 }
 402
 403 FunctorType Taper_i::GetFunctorType()
 404 {
 405   return SMESH::FT_Taper;
 406 }
 407
 408 /*
 409   Class       : Skew_i
 410   Description : Functor for calculating skew in degrees
 411 */
 412 Skew_i::Skew_i()
 413 {
 414   myNumericalFunctorPtr.reset( new Controls::Skew() );
 415   myFunctorPtr = myNumericalFunctorPtr;
 416 }
 417
 418 FunctorType Skew_i::GetFunctorType()
 419 {
 420   return SMESH::FT_Skew;
 421 }
 422
 423 /*
 424   Class       : Area_i
 425   Description : Functor for calculating area
 426 */
 427 Area_i::Area_i()
 428 {
 429   myNumericalFunctorPtr.reset( new Controls::Area() );
 430   myFunctorPtr = myNumericalFunctorPtr;
 431 }
 432
 433 FunctorType Area_i::GetFunctorType()
 434 {
 435   return SMESH::FT_Area;
 436 }
 437
 438 /*
 439   Class       : Volume3D_i
 440   Description : Functor for calculating volume of 3D element
 441 */
 442 Volume3D_i::Volume3D_i()
 443 {
 444   myNumericalFunctorPtr.reset( new Controls::Volume() );
 445   myFunctorPtr = myNumericalFunctorPtr;
 446 }
 447
 448 FunctorType Volume3D_i::GetFunctorType()
 449 {
 450   return SMESH::FT_Volume3D;
 451 }
 452
 453 /*
 454   Class       : MaxElementLength2D_i
 455   Description : Functor for calculating maximum length of 2D element
 456 */
 457 MaxElementLength2D_i::MaxElementLength2D_i()
 458 {
 459   myNumericalFunctorPtr.reset( new Controls::MaxElementLength2D() );
 460   myFunctorPtr = myNumericalFunctorPtr;
 461 }
 462
 463 FunctorType MaxElementLength2D_i::GetFunctorType()
 464 {
 465   return SMESH::FT_MaxElementLength2D;
 466 }
 467
 468 /*
 469   Class       : MaxElementLength3D_i
 470   Description : Functor for calculating maximum length of 3D element
 471 */
 472 MaxElementLength3D_i::MaxElementLength3D_i()
 473 {
 474   myNumericalFunctorPtr.reset( new Controls::MaxElementLength3D() );
 475   myFunctorPtr = myNumericalFunctorPtr;
 476 }
 477
 478 FunctorType MaxElementLength3D_i::GetFunctorType()
 479 {
 480   return SMESH::FT_MaxElementLength3D;
 481 }
 482
 483 /*
 484   Class       : Length_i
 485   Description : Functor for calculating length off edge
 486 */
 487 Length_i::Length_i()
 488 {
 489   myNumericalFunctorPtr.reset( new Controls::Length() );
 490   myFunctorPtr = myNumericalFunctorPtr;
 491 }
 492
 493 FunctorType Length_i::GetFunctorType()
 494 {
 495   return SMESH::FT_Length;
 496 }
 497
 498 /*
 499   Class       : Length2D_i
 500   Description : Functor for calculating length of edge
 501 */
 502 Length2D_i::Length2D_i()
 503 {
 504   myNumericalFunctorPtr.reset( new Controls::Length2D() );
 505   myFunctorPtr = myNumericalFunctorPtr;
 506 }
 507
 508 FunctorType Length2D_i::GetFunctorType()
 509 {
 510   return SMESH::FT_Length2D;
 511 }
 512
 513 SMESH::Length2D::Values* Length2D_i::GetValues()
 514 {
 515   SMESH::Controls::Length2D::TValues aValues;
 516   (dynamic_cast<SMESH::Controls::Length2D*>(myFunctorPtr.get()))->GetValues( aValues );
 517
 518   long i = 0, iEnd = aValues.size();
 519
 520   SMESH::Length2D::Values_var aResult = new SMESH::Length2D::Values(iEnd);
 521   aResult->length(iEnd);
 522
 523   SMESH::Controls::Length2D::TValues::const_iterator anIter;
 524   for ( anIter = aValues.begin() ; anIter != aValues.end(); anIter++, i++ )
 525   {
 526     const SMESH::Controls::Length2D::Value&  aVal = *anIter;
 527     SMESH::Length2D::Value &aValue = aResult[ i ];
 528
 529     aValue.myLength = aVal.myLength;
 530     aValue.myPnt1 = aVal.myPntId[ 0 ];
 531     aValue.myPnt2 = aVal.myPntId[ 1 ];
 532   }
 533
 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       : NodeConnectivityNumber_i
 569   Description : Functor returning diameter of a ball element
 570 */
 571 NodeConnectivityNumber_i::NodeConnectivityNumber_i()
 572 {
 573   myNumericalFunctorPtr.reset( new Controls::NodeConnectivityNumber() );
 574   myFunctorPtr = myNumericalFunctorPtr;
 575 }
 576
 577 FunctorType NodeConnectivityNumber_i::GetFunctorType()
 578 {
 579   return SMESH::FT_NodeConnectivityNumber;
 580 }
 581
 582 /*
 583   Class       : MultiConnection2D_i
 584   Description : Functor for calculating number of faces conneted to the edge
 585 */
 586 MultiConnection2D_i::MultiConnection2D_i()
 587 {
 588   myNumericalFunctorPtr.reset( new Controls::MultiConnection2D() );
 589   myFunctorPtr = myNumericalFunctorPtr;
 590 }
 591
 592 FunctorType MultiConnection2D_i::GetFunctorType()
 593 {
 594   return SMESH::FT_MultiConnection2D;
 595 }
 596
 597 SMESH::MultiConnection2D::Values* MultiConnection2D_i::GetValues()
 598 {
 599   SMESH::Controls::MultiConnection2D::MValues aValues;
 600   (dynamic_cast<SMESH::Controls::MultiConnection2D*>(myFunctorPtr.get()))->GetValues( aValues );
 601
 602   long i = 0, iEnd = aValues.size();
 603
 604   SMESH::MultiConnection2D::Values_var aResult = new SMESH::MultiConnection2D::Values(iEnd);
 605   aResult->length(iEnd);
 606
 607   SMESH::Controls::MultiConnection2D::MValues::const_iterator anIter;
 608   for ( anIter = aValues.begin() ; anIter != aValues.end(); anIter++, i++ )
 609   {
 610     const SMESH::Controls::MultiConnection2D::Value&  aVal = (*anIter).first;
 611     SMESH::MultiConnection2D::Value &aValue = aResult[ i ];
 612
 613     aValue.myPnt1 = aVal.myPntId[ 0 ];
 614     aValue.myPnt2 = aVal.myPntId[ 1 ];
 615     aValue.myNbConnects = (*anIter).second;
 616   }
 617
 618   return aResult._retn();
 619 }
 620
 621 /*
 622                             PREDICATES
 623 */
 624
 625
 626 /*
 627   Class       : Predicate_i
 628   Description : Base class for all predicates
 629 */
 630 CORBA::Boolean Predicate_i::IsSatisfy( CORBA::Long theId )
 631 {
 632   return myPredicatePtr->IsSatisfy( theId );
 633 }
 634
 635 CORBA::Long Predicate_i::NbSatisfying( SMESH::SMESH_IDSource_ptr obj )
 636 {
 637   SMESH::SMESH_Mesh_var meshVar = obj->GetMesh();
 638   const SMDS_Mesh*       meshDS = MeshPtr2SMDSMesh( meshVar );
 639   if ( !meshDS )
 640     return 0;
 641   myPredicatePtr->SetMesh( meshDS );
 642
 643   SMDSAbs_ElementType elemType = SMDSAbs_ElementType( GetElementType() );
 644
 645   int nb = 0;
 646   SMDS_ElemIteratorPtr elemIt =
 647     SMESH::DownCast<SMESH_Mesh_i*>( meshVar )->GetElements( obj, GetElementType() );
 648   if ( elemIt )
 649     while ( elemIt->more() )
 650     {
 651       const SMDS_MeshElement* e = elemIt->next();
 652       if ( e && e->GetType() == elemType )
 653         nb += myPredicatePtr->IsSatisfy( e->GetID() );
 654     }
 655   return nb;
 656 }
 657
 658 Controls::PredicatePtr Predicate_i::GetPredicate()
 659 {
 660   return myPredicatePtr;
 661 }
 662
 663 /*
 664   Class       : BadOrientedVolume_i
 665   Description : Verify whether a mesh volume is incorrectly oriented from
 666                 the point of view of MED convention
 667 */
 668 BadOrientedVolume_i::BadOrientedVolume_i()
 669 {
 670   Controls::PredicatePtr control( new Controls::BadOrientedVolume() );
 671   myFunctorPtr = myPredicatePtr = control;
 672 };
 673
 674 FunctorType BadOrientedVolume_i::GetFunctorType()
 675 {
 676   return SMESH::FT_BadOrientedVolume;
 677 }
 678
 679 /*
 680   Class       : BareBorderVolume_i
 681   Description : Verify whether a mesh volume has a free facet without a face on it
 682 */
 683 BareBorderVolume_i::BareBorderVolume_i()
 684 {
 685   Controls::PredicatePtr control( new Controls::BareBorderVolume() );
 686   myFunctorPtr = myPredicatePtr = control;
 687 };
 688
 689 FunctorType BareBorderVolume_i::GetFunctorType()
 690 {
 691   return SMESH::FT_BareBorderVolume;
 692 }
 693
 694 /*
 695   Class       : BareBorderFace_i
 696   Description : Verify whether a mesh face has a free border without an edge on it
 697 */
 698 BareBorderFace_i::BareBorderFace_i()
 699 {
 700   Controls::PredicatePtr control( new Controls::BareBorderFace() );
 701   myFunctorPtr = myPredicatePtr = control;
 702 };
 703
 704 FunctorType BareBorderFace_i::GetFunctorType()
 705 {
 706   return SMESH::FT_BareBorderFace;
 707 }
 708
 709 /*
 710   Class       : OverConstrainedVolume_i
 711   Description : Verify whether a mesh volume has only one facet shared with other volumes
 712 */
 713 OverConstrainedVolume_i::OverConstrainedVolume_i()
 714 {
 715   Controls::PredicatePtr control( new Controls::OverConstrainedVolume() );
 716   myFunctorPtr = myPredicatePtr = control;
 717 };
 718
 719 FunctorType OverConstrainedVolume_i::GetFunctorType()
 720 {
 721   return SMESH::FT_OverConstrainedVolume;
 722 }
 723
 724 /*
 725   Class       : OverConstrainedFace_i
 726   Description : Verify whether a mesh face has only one border shared with other faces
 727 */
 728 OverConstrainedFace_i::OverConstrainedFace_i()
 729 {
 730   Controls::PredicatePtr control( new Controls::OverConstrainedFace() );
 731   myFunctorPtr = myPredicatePtr = control;
 732 };
 733
 734 FunctorType OverConstrainedFace_i::GetFunctorType()
 735 {
 736   return SMESH::FT_OverConstrainedFace;
 737 }
 738
 739 /*
 740   Class       : BelongToMeshGroup_i
 741   Description : Verify whether a mesh element is included into a mesh group
 742 */
 743 BelongToMeshGroup_i::BelongToMeshGroup_i()
 744 {
 745   myBelongToMeshGroup = Controls::BelongToMeshGroupPtr( new Controls::BelongToMeshGroup() );
 746   myFunctorPtr = myPredicatePtr = myBelongToMeshGroup;
 747 }
 748
 749 BelongToMeshGroup_i::~BelongToMeshGroup_i()
 750 {
 751   SetGroup( SMESH::SMESH_GroupBase::_nil() );
 752 }
 753
 754 void BelongToMeshGroup_i::SetGroup( SMESH::SMESH_GroupBase_ptr theGroup )
 755 {
 756   if ( myGroup->_is_equivalent( theGroup ))
 757     return;
 758
 759   if ( ! myGroup->_is_nil() )
 760     myGroup->UnRegister();
 761
 762   myGroup = SMESH_GroupBase::_duplicate( theGroup );
 763
 764   myBelongToMeshGroup->SetGroup( 0 );
 765   if ( SMESH_GroupBase_i* gr_i = SMESH::DownCast< SMESH_GroupBase_i* >( myGroup ))
 766   {
 767     myBelongToMeshGroup->SetGroup( gr_i->GetGroupDS() );
 768     myGroup->Register();
 769   }
 770 }
 771
 772 void BelongToMeshGroup_i::SetGroupID( const char* theID ) // IOR or StoreName
 773 {
 774   myID = theID;
 775   if ( strncmp( "IOR:", myID.c_str(), 4 ) == 0 ) // transient mode, no GUI
 776   {
 777     CORBA::Object_var obj = SMESH_Gen_i::GetORB()->string_to_object( myID.c_str() );
 778     SetGroup( SMESH::SMESH_GroupBase::_narrow( obj ));
 779   }
 780   else if ( strncmp( "0:", myID.c_str(), 2 ) == 0 ) // transient mode + GUI
 781   {
 782     SMESH_Gen_i* aSMESHGen     = SMESH_Gen_i::GetSMESHGen();
 783     SALOMEDS::Study_var aStudy = aSMESHGen->GetCurrentStudy();
 784     if ( !aStudy->_is_nil() ) {
 785       SALOMEDS::SObject_wrap aSObj = aStudy->FindObjectID( myID.c_str() );
 786       if ( !aSObj->_is_nil() ) {
 787         CORBA::Object_var obj = aSObj->GetObject();
 788         SetGroup( SMESH::SMESH_GroupBase::_narrow( obj ));
 789       }
 790     }
 791   }
 792   else if ( !myID.empty() ) // persistent mode
 793   {
 794     myBelongToMeshGroup->SetStoreName( myID );
 795   }
 796 }
 797
 798 std::string BelongToMeshGroup_i::GetGroupID()
 799 {
 800   if ( myGroup->_is_nil() )
 801     SMESH::SMESH_GroupBase_var( GetGroup() );
 802
 803   if ( !myGroup->_is_nil() )
 804     myID = SMESH_Gen_i::GetORB()->object_to_string( myGroup );
 805
 806   return myID;
 807 }
 808
 809 SMESH::SMESH_GroupBase_ptr BelongToMeshGroup_i::GetGroup()
 810 {
 811   if ( myGroup->_is_nil() && myBelongToMeshGroup->GetGroup() )
 812   {
 813     // search for a group in a current study
 814     SMESH_Gen_i* aSMESHGen = SMESH_Gen_i::GetSMESHGen();
 815     if ( StudyContext*  sc = aSMESHGen->GetCurrentStudyContext() )
 816     {
 817       int id = 1;
 818       std::string ior;
 819       while (true)
 820       {
 821         ior = sc->getIORbyId( id++ );
 822         if ( ior.empty() ) break;
 823         CORBA::Object_var obj = aSMESHGen->GetORB()->string_to_object( ior.c_str() );
 824         if ( SMESH_GroupBase_i* g_i = SMESH::DownCast<SMESH_GroupBase_i*>( obj ))
 825           if ( g_i->GetGroupDS() == myBelongToMeshGroup->GetGroup() )
 826           {
 827             SetGroup( g_i->_this() );
 828             break;
 829           }
 830       }
 831     }
 832   }
 833   return SMESH::SMESH_GroupBase::_duplicate( myGroup );
 834 }
 835
 836 FunctorType BelongToMeshGroup_i::GetFunctorType()
 837 {
 838   return SMESH::FT_BelongToMeshGroup;
 839 }
 840
 841 /*
 842   Class       : BelongToGeom_i
 843   Description : Predicate for selection on geometrical support
 844 */
 845 BelongToGeom_i::BelongToGeom_i()
 846 {
 847   myBelongToGeomPtr.reset( new Controls::BelongToGeom() );
 848   myFunctorPtr = myPredicatePtr = myBelongToGeomPtr;
 849   myShapeName = 0;
 850   myShapeID   = 0;
 851 }
 852
 853 BelongToGeom_i::~BelongToGeom_i()
 854 {
 855   CORBA::string_free( myShapeName );
 856   CORBA::string_free( myShapeID );
 857 }
 858
 859 void BelongToGeom_i::SetGeom( GEOM::GEOM_Object_ptr theGeom )
 860 {
 861   if ( theGeom->_is_nil() )
 862     return;
 863   SMESH_Gen_i* aSMESHGen = SMESH_Gen_i::GetSMESHGen();
 864   GEOM::GEOM_Gen_ptr aGEOMGen = SMESH_Gen_i::GetGeomEngine();
 865   TopoDS_Shape aLocShape = aSMESHGen->GetShapeReader()->GetShape( aGEOMGen, theGeom );
 866   myBelongToGeomPtr->SetGeom( aLocShape );
 867   TPythonDump()<<this<<".SetGeom("<<theGeom<<")";
 868 }
 869
 870 void BelongToGeom_i::SetGeom( const TopoDS_Shape& theShape )
 871 {
 872   myBelongToGeomPtr->SetGeom( theShape );
 873 }
 874
 875 void BelongToGeom_i::SetElementType(ElementType theType)
 876 {
 877   myBelongToGeomPtr->SetType(SMDSAbs_ElementType(theType));
 878   TPythonDump()<<this<<".SetElementType("<<theType<<")";
 879 }
 880
 881 FunctorType BelongToGeom_i::GetFunctorType()
 882 {
 883   return SMESH::FT_BelongToGeom;
 884 }
 885
 886 void BelongToGeom_i::SetShapeName( const char* theName )
 887 {
 888   CORBA::string_free( myShapeName );
 889   myShapeName = CORBA::string_dup( theName );
 890   myBelongToGeomPtr->SetGeom( getShapeByName( myShapeName ) );
 891   TPythonDump()<<this<<".SetShapeName('"<<theName<<"')";
 892 }
 893
 894 void BelongToGeom_i::SetShape( const char* theID, const char* theName )
 895 {
 896   CORBA::string_free( myShapeName );
 897   CORBA::string_free( myShapeID );
 898   myShapeName = CORBA::string_dup( theName );
 899   myShapeID   = CORBA::string_dup( theID );
 900   bool hasName = ( theName && theName[0] );
 901   bool hasID   = ( theID   && theID[0] );
 902
 903   TopoDS_Shape S;
 904   if ( hasName && hasID )
 905   {
 906     S = getShapeByID( myShapeID );
 907     if ( S.IsNull() )
 908       S = getShapeByName( myShapeName );
 909   }
 910   else
 911   {
 912     S = hasID ? getShapeByID( myShapeID ) : getShapeByName( myShapeName );
 913   }
 914   myBelongToGeomPtr->SetGeom( S );
 915 }
 916
 917 char* BelongToGeom_i::GetShapeName()
 918 {
 919   return CORBA::string_dup( myShapeName );
 920 }
 921
 922 char* BelongToGeom_i::GetShapeID()
 923 {
 924   return CORBA::string_dup( myShapeID );
 925 }
 926
 927 void BelongToGeom_i::SetTolerance( CORBA::Double theToler )
 928 {
 929   myBelongToGeomPtr->SetTolerance( theToler );
 930   TPythonDump()<<this<<".SetTolerance("<<theToler<<")";
 931 }
 932
 933 CORBA::Double BelongToGeom_i::GetTolerance()
 934 {
 935   return myBelongToGeomPtr->GetTolerance();
 936 }
 937
 938 /*
 939   Class       : BelongToSurface_i
 940   Description : Predicate for selection on geometrical support
 941 */
 942 BelongToSurface_i::BelongToSurface_i( const Handle(Standard_Type)& theSurfaceType )
 943 {
 944   myElementsOnSurfacePtr.reset( new Controls::ElementsOnSurface() );
 945   myFunctorPtr = myPredicatePtr = myElementsOnSurfacePtr;
 946   myShapeName = 0;
 947   myShapeID   = 0;
 948   mySurfaceType = theSurfaceType;
 949 }
 950
 951 BelongToSurface_i::~BelongToSurface_i()
 952 {
 953   CORBA::string_free( myShapeName );
 954   CORBA::string_free( myShapeID );
 955 }
 956
 957 void BelongToSurface_i::SetSurface( GEOM::GEOM_Object_ptr theGeom, ElementType theType )
 958 {
 959   if ( theGeom->_is_nil() )
 960     return;
 961   SMESH_Gen_i* aSMESHGen = SMESH_Gen_i::GetSMESHGen();
 962   GEOM::GEOM_Gen_ptr aGEOMGen = SMESH_Gen_i::GetGeomEngine();
 963   TopoDS_Shape aLocShape = aSMESHGen->GetShapeReader()->GetShape( aGEOMGen, theGeom );
 964
 965   if ( aLocShape.ShapeType() == TopAbs_FACE )
 966   {
 967     Handle(Geom_Surface) aSurf = BRep_Tool::Surface( TopoDS::Face( aLocShape ) );
 968     if ( !aSurf.IsNull() && aSurf->DynamicType() == mySurfaceType )
 969     {
 970       myElementsOnSurfacePtr->SetSurface( aLocShape, (SMDSAbs_ElementType)theType );
 971       return;
 972     }
 973   }
 974
 975   myElementsOnSurfacePtr->SetSurface( TopoDS_Shape(), (SMDSAbs_ElementType)theType );
 976 }
 977
 978 void BelongToSurface_i::SetShapeName( const char* theName, ElementType theType )
 979 {
 980   CORBA::string_free( myShapeName );
 981   myShapeName = CORBA::string_dup( theName );
 982   myElementsOnSurfacePtr->SetSurface( getShapeByName( myShapeName ), (SMDSAbs_ElementType)theType );
 983   TPythonDump()<<this<<".SetShapeName('"<<theName<<"',"<<theType<<")";
 984 }
 985
 986 void BelongToSurface_i::SetShape( const char* theID,  const char* theName, ElementType theType )
 987 {
 988   CORBA::string_free( myShapeName );
 989   CORBA::string_free( myShapeID );
 990   myShapeName = CORBA::string_dup( theName );
 991   myShapeID   = CORBA::string_dup( theID );
 992   bool hasName = ( theName && theName[0] );
 993   bool hasID   = ( theID   && theID[0] );
 994
 995   TopoDS_Shape S;
 996   if ( hasName && hasID )
 997   {
 998     S = getShapeByID( myShapeID );
 999     if ( S.IsNull() )
1000       S = getShapeByName( myShapeName );
1001   }
1002   else
1003   {
1004     S = hasID ? getShapeByID( myShapeID ) : getShapeByName( myShapeName );
1005   }
1006   myElementsOnSurfacePtr->SetSurface( S, (SMDSAbs_ElementType)theType );
1007 }
1008
1009 char* BelongToSurface_i::GetShapeName()
1010 {
1011   return CORBA::string_dup( myShapeName );
1012 }
1013
1014 char* BelongToSurface_i::GetShapeID()
1015 {
1016   return CORBA::string_dup( myShapeID );
1017 }
1018
1019 void BelongToSurface_i::SetTolerance( CORBA::Double theToler )
1020 {
1021   myElementsOnSurfacePtr->SetTolerance( theToler );
1022   TPythonDump()<<this<<".SetTolerance("<<theToler<<")";
1023 }
1024
1025 CORBA::Double BelongToSurface_i::GetTolerance()
1026 {
1027   return myElementsOnSurfacePtr->GetTolerance();
1028 }
1029
1030 void BelongToSurface_i::SetUseBoundaries( CORBA::Boolean theUseBndRestrictions )
1031 {
1032   myElementsOnSurfacePtr->SetUseBoundaries( theUseBndRestrictions );
1033   TPythonDump()<<this<<".SetUseBoundaries( " << theUseBndRestrictions << " )";
1034 }
1035
1036 CORBA::Boolean BelongToSurface_i::GetUseBoundaries()
1037 {
1038   return myElementsOnSurfacePtr->GetUseBoundaries();
1039 }
1040
1041
1042 /*
1043   Class       : BelongToPlane_i
1044   Description : Verify whether mesh element lie in pointed Geom planar object
1045 */
1046
1047 BelongToPlane_i::BelongToPlane_i()
1048 : BelongToSurface_i( STANDARD_TYPE( Geom_Plane ) )
1049 {
1050 }
1051
1052 void BelongToPlane_i::SetPlane( GEOM::GEOM_Object_ptr theGeom, ElementType theType )
1053 {
1054   BelongToSurface_i::SetSurface( theGeom, theType );
1055   TPythonDump()<<this<<".SetPlane("<<theGeom<<","<<theType<<")";
1056 }
1057
1058 FunctorType BelongToPlane_i::GetFunctorType()
1059 {
1060   return FT_BelongToPlane;
1061 }
1062
1063 /*
1064   Class       : BelongToCylinder_i
1065   Description : Verify whether mesh element lie in pointed Geom planar object
1066 */
1067
1068 BelongToCylinder_i::BelongToCylinder_i()
1069 : BelongToSurface_i( STANDARD_TYPE( Geom_CylindricalSurface ) )
1070 {
1071 }
1072
1073 void BelongToCylinder_i::SetCylinder( GEOM::GEOM_Object_ptr theGeom, ElementType theType )
1074 {
1075   BelongToSurface_i::SetSurface( theGeom, theType );
1076   TPythonDump()<<this<<".SetCylinder("<<theGeom<<","<<theType<<")";
1077 }
1078
1079 FunctorType BelongToCylinder_i::GetFunctorType()
1080 {
1081   return FT_BelongToCylinder;
1082 }
1083
1084 /*
1085   Class       : BelongToGenSurface_i
1086   Description : Verify whether mesh element lie in pointed Geom planar object
1087 */
1088
1089 BelongToGenSurface_i::BelongToGenSurface_i()
1090 : BelongToSurface_i( STANDARD_TYPE( Geom_CylindricalSurface ) )
1091 {
1092 }
1093
1094 void BelongToGenSurface_i::SetSurface( GEOM::GEOM_Object_ptr theGeom, ElementType theType )
1095 {
1096   if ( theGeom->_is_nil() )
1097     return;
1098   TopoDS_Shape aLocShape = SMESH_Gen_i::GetSMESHGen()->GeomObjectToShape( theGeom );
1099   if ( !aLocShape.IsNull() && aLocShape.ShapeType() != TopAbs_FACE )
1100     aLocShape.Nullify();
1101
1102   BelongToSurface_i::myElementsOnSurfacePtr->SetSurface( aLocShape, (SMDSAbs_ElementType)theType );
1103   TPythonDump()<<this<<".SetGenSurface("<<theGeom<<","<<theType<<")";
1104 }
1105
1106 FunctorType BelongToGenSurface_i::GetFunctorType()
1107 {
1108   return FT_BelongToGenSurface;
1109 }
1110
1111 /*
1112   Class       : LyingOnGeom_i
1113   Description : Predicate for selection on geometrical support
1114 */
1115 LyingOnGeom_i::LyingOnGeom_i()
1116 {
1117   myLyingOnGeomPtr.reset( new Controls::LyingOnGeom() );
1118   myFunctorPtr = myPredicatePtr = myLyingOnGeomPtr;
1119   myShapeName = 0;
1120   myShapeID = 0;
1121 }
1122
1123 LyingOnGeom_i::~LyingOnGeom_i()
1124 {
1125   CORBA::string_free( myShapeName );
1126   CORBA::string_free( myShapeID );
1127 }
1128
1129 void LyingOnGeom_i::SetGeom( GEOM::GEOM_Object_ptr theGeom )
1130 {
1131   if ( theGeom->_is_nil() )
1132     return;
1133   SMESH_Gen_i* aSMESHGen = SMESH_Gen_i::GetSMESHGen();
1134   GEOM::GEOM_Gen_ptr aGEOMGen = SMESH_Gen_i::GetGeomEngine();
1135   TopoDS_Shape aLocShape = aSMESHGen->GetShapeReader()->GetShape( aGEOMGen, theGeom );
1136   myLyingOnGeomPtr->SetGeom( aLocShape );
1137   TPythonDump()<<this<<".SetGeom("<<theGeom<<")";
1138 }
1139
1140 void LyingOnGeom_i::SetGeom( const TopoDS_Shape& theShape )
1141 {
1142   myLyingOnGeomPtr->SetGeom( theShape );
1143 }
1144
1145 void LyingOnGeom_i::SetElementType(ElementType theType){
1146   myLyingOnGeomPtr->SetType(SMDSAbs_ElementType(theType));
1147   TPythonDump()<<this<<".SetElementType("<<theType<<")";
1148 }
1149
1150 FunctorType LyingOnGeom_i::GetFunctorType()
1151 {
1152   return SMESH::FT_LyingOnGeom;
1153 }
1154
1155 void LyingOnGeom_i::SetShapeName( const char* theName )
1156 {
1157   CORBA::string_free( myShapeName );
1158   myShapeName = CORBA::string_dup( theName );
1159   myLyingOnGeomPtr->SetGeom( getShapeByName( myShapeName ) );
1160   TPythonDump()<<this<<".SetShapeName('"<<theName<<"')";
1161 }
1162
1163 void LyingOnGeom_i::SetShape( const char* theID, const char* theName )
1164 {
1165   CORBA::string_free( myShapeName );
1166   CORBA::string_free( myShapeID   );
1167   myShapeName = CORBA::string_dup( theName );
1168   myShapeID   = CORBA::string_dup( theID   );
1169   bool hasName = ( theName && theName[0] );
1170   bool hasID   = ( theID   && theID[0]   );
1171
1172   TopoDS_Shape S;
1173   if ( hasName && hasID )
1174   {
1175     S = getShapeByID( myShapeID );
1176     if ( S.IsNull() )
1177       S = getShapeByName( myShapeName );
1178   }
1179   else
1180   {
1181     S = hasID ? getShapeByID( myShapeID ) : getShapeByName( myShapeName );
1182   }
1183   myLyingOnGeomPtr->SetGeom( S );
1184 }
1185
1186 char* LyingOnGeom_i::GetShapeName()
1187 {
1188   return CORBA::string_dup( myShapeName );
1189 }
1190
1191 char* LyingOnGeom_i::GetShapeID()
1192 {
1193   return CORBA::string_dup( myShapeID );
1194 }
1195
1196 void LyingOnGeom_i::SetTolerance( CORBA::Double theToler )
1197 {
1198   myLyingOnGeomPtr->SetTolerance( theToler );
1199   TPythonDump()<<this<<".SetTolerance("<<theToler<<")";
1200 }
1201
1202 CORBA::Double LyingOnGeom_i::GetTolerance()
1203 {
1204   return myLyingOnGeomPtr->GetTolerance();
1205 }
1206
1207 /*
1208   Class       : FreeBorders_i
1209   Description : Predicate for free borders
1210 */
1211 FreeBorders_i::FreeBorders_i()
1212 {
1213   myPredicatePtr.reset(new Controls::FreeBorders());
1214   myFunctorPtr = myPredicatePtr;
1215 }
1216
1217 FunctorType FreeBorders_i::GetFunctorType()
1218 {
1219   return SMESH::FT_FreeBorders;
1220 }
1221
1222 /*
1223   Class       : FreeEdges_i
1224   Description : Predicate for free borders
1225 */
1226 FreeEdges_i::FreeEdges_i()
1227 : myFreeEdgesPtr( new Controls::FreeEdges() )
1228 {
1229   myFunctorPtr = myPredicatePtr = myFreeEdgesPtr;
1230 }
1231
1232 SMESH::FreeEdges::Borders* FreeEdges_i::GetBorders()
1233 {
1234   SMESH::Controls::FreeEdges::TBorders aBorders;
1235   myFreeEdgesPtr->GetBoreders( aBorders );
1236
1237   long i = 0, iEnd = aBorders.size();
1238
1239   SMESH::FreeEdges::Borders_var aResult = new SMESH::FreeEdges::Borders;
1240   aResult->length(iEnd);
1241
1242   SMESH::Controls::FreeEdges::TBorders::const_iterator anIter;
1243   for ( anIter = aBorders.begin() ; anIter != aBorders.end(); anIter++, i++ )
1244   {
1245     const SMESH::Controls::FreeEdges::Border&  aBord = *anIter;
1246     SMESH::FreeEdges::Border &aBorder = aResult[ i ];
1247
1248     aBorder.myElemId = aBord.myElemId;
1249     aBorder.myPnt1 = aBord.myPntId[ 0 ];
1250     aBorder.myPnt2 = aBord.myPntId[ 1 ];
1251   }
1252   return aResult._retn();
1253 }
1254
1255 FunctorType FreeEdges_i::GetFunctorType()
1256 {
1257   return SMESH::FT_FreeEdges;
1258 }
1259
1260 /*
1261   Class       : FreeFaces_i
1262   Description : Predicate for free faces
1263 */
1264 FreeFaces_i::FreeFaces_i()
1265 {
1266   myPredicatePtr.reset(new Controls::FreeFaces());
1267   myFunctorPtr = myPredicatePtr;
1268 }
1269
1270 FunctorType FreeFaces_i::GetFunctorType()
1271 {
1272   return SMESH::FT_FreeFaces;
1273 }
1274
1275 /*
1276   Class       : FreeNodes_i
1277   Description : Predicate for free nodes
1278 */
1279 FreeNodes_i::FreeNodes_i()
1280 {
1281   myPredicatePtr.reset(new Controls::FreeNodes());
1282   myFunctorPtr = myPredicatePtr;
1283 }
1284
1285 FunctorType FreeNodes_i::GetFunctorType()
1286 {
1287   return SMESH::FT_FreeNodes;
1288 }
1289
1290 /*
1291   Class       : EqualNodes_i
1292   Description : Predicate for Equal nodes
1293 */
1294 EqualNodes_i::EqualNodes_i()
1295 {
1296   myCoincidentNodesPtr.reset(new Controls::CoincidentNodes());
1297   myFunctorPtr = myPredicatePtr = myCoincidentNodesPtr;
1298 }
1299
1300 FunctorType EqualNodes_i::GetFunctorType()
1301 {
1302   return SMESH::FT_EqualNodes;
1303 }
1304
1305 void EqualNodes_i::SetTolerance( double tol )
1306 {
1307   myCoincidentNodesPtr->SetTolerance( tol );
1308 }
1309
1310 double EqualNodes_i::GetTolerance()
1311 {
1312   return myCoincidentNodesPtr->GetTolerance();
1313 }
1314
1315 /*
1316   Class       : EqualEdges_i
1317   Description : Predicate for Equal Edges
1318 */
1319 EqualEdges_i::EqualEdges_i()
1320 {
1321   myPredicatePtr.reset(new Controls::CoincidentElements1D());
1322   myFunctorPtr = myPredicatePtr;
1323 }
1324
1325 FunctorType EqualEdges_i::GetFunctorType()
1326 {
1327   return SMESH::FT_EqualEdges;
1328 }
1329
1330 /*
1331   Class       : EqualFaces_i
1332   Description : Predicate for Equal Faces
1333 */
1334 EqualFaces_i::EqualFaces_i()
1335 {
1336   myPredicatePtr.reset(new Controls::CoincidentElements2D());
1337   myFunctorPtr = myPredicatePtr;
1338 }
1339
1340 FunctorType EqualFaces_i::GetFunctorType()
1341 {
1342   return SMESH::FT_EqualFaces;
1343 }
1344
1345 /*
1346   Class       : EqualVolumes_i
1347   Description : Predicate for Equal Volumes
1348 */
1349 EqualVolumes_i::EqualVolumes_i()
1350 {
1351   myPredicatePtr.reset(new Controls::CoincidentElements3D());
1352   myFunctorPtr = myPredicatePtr;
1353 }
1354
1355 FunctorType EqualVolumes_i::GetFunctorType()
1356 {
1357   return SMESH::FT_EqualVolumes;
1358 }
1359
1360
1361 /*
1362   Class       : RangeOfIds_i
1363   Description : Predicate for Range of Ids.
1364                 Range may be specified with two ways.
1365                 1. Using AddToRange method
1366                 2. With SetRangeStr method. Parameter of this method is a string
1367                    like as "1,2,3,50-60,63,67,70-"
1368 */
1369
1370 RangeOfIds_i::RangeOfIds_i()
1371 {
1372   myRangeOfIdsPtr.reset( new Controls::RangeOfIds() );
1373   myFunctorPtr = myPredicatePtr = myRangeOfIdsPtr;
1374 }
1375
1376 void RangeOfIds_i::SetRange( const SMESH::long_array& theIds )
1377 {
1378   CORBA::Long iEnd = theIds.length();
1379   for ( CORBA::Long i = 0; i < iEnd; i++ )
1380     myRangeOfIdsPtr->AddToRange( theIds[ i ] );
1381   TPythonDump()<<this<<".SetRange("<<theIds<<")";
1382 }
1383
1384 CORBA::Boolean RangeOfIds_i::SetRangeStr( const char* theRange )
1385 {
1386   TPythonDump()<<this<<".SetRangeStr('"<<theRange<<"')";
1387   return myRangeOfIdsPtr->SetRangeStr(
1388     TCollection_AsciiString( (Standard_CString)theRange ) );
1389 }
1390
1391 char* RangeOfIds_i::GetRangeStr()
1392 {
1393   TCollection_AsciiString aStr;
1394   myRangeOfIdsPtr->GetRangeStr( aStr );
1395   return CORBA::string_dup( aStr.ToCString() );
1396 }
1397
1398 void RangeOfIds_i::SetElementType( ElementType theType )
1399 {
1400   myRangeOfIdsPtr->SetType( SMDSAbs_ElementType( theType ) );
1401   TPythonDump()<<this<<".SetElementType("<<theType<<")";
1402 }
1403
1404 FunctorType RangeOfIds_i::GetFunctorType()
1405 {
1406   return SMESH::FT_RangeOfIds;
1407 }
1408
1409 /*
1410   Class       : LinearOrQuadratic_i
1411   Description : Predicate to verify whether a mesh element is linear
1412 */
1413 LinearOrQuadratic_i::LinearOrQuadratic_i()
1414 {
1415   myLinearOrQuadraticPtr.reset(new Controls::LinearOrQuadratic());
1416   myFunctorPtr = myPredicatePtr = myLinearOrQuadraticPtr;
1417 }
1418
1419 void LinearOrQuadratic_i::SetElementType(ElementType theType)
1420 {
1421   myLinearOrQuadraticPtr->SetType(SMDSAbs_ElementType(theType));
1422   TPythonDump()<<this<<".SetElementType("<<theType<<")";
1423 }
1424
1425 FunctorType LinearOrQuadratic_i::GetFunctorType()
1426 {
1427   return SMESH::FT_LinearOrQuadratic;
1428 }
1429
1430 /*
1431   Class       : GroupColor_i
1432   Description : Functor for check color of group to which mesh element belongs to
1433 */
1434 GroupColor_i::GroupColor_i()
1435 {
1436   myGroupColorPtr.reset(new Controls::GroupColor());
1437   myFunctorPtr = myPredicatePtr = myGroupColorPtr;
1438 }
1439
1440 FunctorType GroupColor_i::GetFunctorType()
1441 {
1442   return SMESH::FT_GroupColor;
1443 }
1444
1445 void GroupColor_i::SetColorStr( const char* theColor )
1446 {
1447   myGroupColorPtr->SetColorStr(
1448     TCollection_AsciiString( (Standard_CString)theColor ) );
1449   TPythonDump()<<this<<".SetColorStr('"<<theColor<<"')";
1450 }
1451
1452 char* GroupColor_i::GetColorStr()
1453 {
1454   TCollection_AsciiString aStr;
1455   myGroupColorPtr->GetColorStr( aStr );
1456   return CORBA::string_dup( aStr.ToCString() );
1457 }
1458
1459 void GroupColor_i::SetElementType(ElementType theType)
1460 {
1461   myGroupColorPtr->SetType(SMDSAbs_ElementType(theType));
1462   TPythonDump()<<this<<".SetElementType("<<theType<<")";
1463 }
1464
1465 /*
1466   Class       : ElemGeomType_i
1467   Description : Predicate check is element has indicated geometry type
1468 */
1469 ElemGeomType_i::ElemGeomType_i()
1470 {
1471   myElemGeomTypePtr.reset(new Controls::ElemGeomType());
1472   myFunctorPtr = myPredicatePtr = myElemGeomTypePtr;
1473 }
1474
1475 void ElemGeomType_i::SetElementType(ElementType theType)
1476 {
1477   myElemGeomTypePtr->SetType(SMDSAbs_ElementType(theType));
1478   TPythonDump()<<this<<".SetElementType("<<theType<<")";
1479 }
1480
1481 void ElemGeomType_i::SetGeometryType(GeometryType theType)
1482 {
1483   myElemGeomTypePtr->SetGeomType(SMDSAbs_GeometryType(theType));
1484   TPythonDump()<<this<<".SetGeometryType("<<theType<<")";
1485 }
1486
1487 GeometryType ElemGeomType_i::GetGeometryType() const
1488 {
1489   return (GeometryType)myElemGeomTypePtr->GetGeomType();
1490 }
1491
1492 FunctorType ElemGeomType_i::GetFunctorType()
1493 {
1494   return SMESH::FT_ElemGeomType;
1495 }
1496
1497 /*
1498   Class       : ElemEntityType_i
1499   Description : Predicate check is element has indicated entity type
1500 */
1501 ElemEntityType_i::ElemEntityType_i()
1502 {
1503   myElemEntityTypePtr.reset(new Controls::ElemEntityType());
1504   myFunctorPtr = myPredicatePtr = myElemEntityTypePtr;
1505 }
1506
1507 void ElemEntityType_i::SetElementType(ElementType theType)
1508 {
1509   myElemEntityTypePtr->SetType(SMDSAbs_ElementType(theType));
1510   TPythonDump()<<this<<".SetElementType("<<theType<<")";
1511 }
1512
1513 void ElemEntityType_i::SetEntityType(EntityType theEntityType)
1514 {
1515   myElemEntityTypePtr->SetElemEntityType(SMDSAbs_EntityType (theEntityType));
1516   TPythonDump()<<this<<".SetEntityType("<<theEntityType<<")";
1517 }
1518 EntityType ElemEntityType_i::GetEntityType() const
1519 {
1520  return (EntityType) myElemEntityTypePtr->GetElemEntityType();
1521 }
1522
1523 FunctorType ElemEntityType_i::GetFunctorType()
1524 {
1525   return SMESH::FT_EntityType;
1526 }
1527
1528 /*
1529   Class       : CoplanarFaces_i
1530   Description : Returns true if a mesh face is a coplanar neighbour to a given one
1531 */
1532 CoplanarFaces_i::CoplanarFaces_i()
1533 {
1534   myCoplanarFacesPtr.reset(new Controls::CoplanarFaces());
1535   myFunctorPtr = myPredicatePtr = myCoplanarFacesPtr;
1536 }
1537
1538 void CoplanarFaces_i::SetFace ( CORBA::Long theFaceID )
1539 {
1540   myCoplanarFacesPtr->SetFace(theFaceID);
1541   TPythonDump()<<this<<".SetFace("<<theFaceID<<")";
1542 }
1543
1544 void CoplanarFaces_i::SetTolerance( CORBA::Double theToler )
1545 {
1546   myCoplanarFacesPtr->SetTolerance(theToler);
1547   TPythonDump()<<this<<".SetTolerance("<<theToler<<")";
1548 }
1549
1550 CORBA::Long CoplanarFaces_i::GetFace () const
1551 {
1552   return myCoplanarFacesPtr->GetFace();
1553 }
1554
1555 char* CoplanarFaces_i::GetFaceAsString () const
1556 {
1557   TCollection_AsciiString str(Standard_Integer(myCoplanarFacesPtr->GetFace()));
1558   return CORBA::string_dup( str.ToCString() );
1559 }
1560
1561 CORBA::Double CoplanarFaces_i::GetTolerance() const
1562 {
1563   return myCoplanarFacesPtr->GetTolerance();
1564 }
1565
1566 FunctorType CoplanarFaces_i::GetFunctorType()
1567 {
1568   return SMESH::FT_CoplanarFaces;
1569 }
1570
1571 /*
1572  * Class       : ConnectedElements_i
1573  * Description : Returns true if an element is connected via other elements to the element
1574  *               located at a given point.
1575  */
1576 ConnectedElements_i::ConnectedElements_i()
1577 {
1578   myConnectedElementsPtr.reset(new Controls::ConnectedElements());
1579   myFunctorPtr = myPredicatePtr = myConnectedElementsPtr;
1580 }
1581
1582 FunctorType ConnectedElements_i::GetFunctorType()
1583 {
1584   return FT_ConnectedElements;
1585 }
1586
1587 void ConnectedElements_i::SetElementType( ElementType theType )
1588 {
1589   myConnectedElementsPtr->SetType( SMDSAbs_ElementType( theType ));
1590   TPythonDump() << this << ".SetElementType( " << theType << " )";
1591 }
1592
1593 void ConnectedElements_i::SetPoint( CORBA::Double x, CORBA::Double y, CORBA::Double z )
1594 {
1595   myConnectedElementsPtr->SetPoint( x,y,z );
1596   myVertexID.clear();
1597   TPythonDump() << this << ".SetPoint( " << x << ", " << y << ", " << z << " )";
1598 }
1599
1600 void ConnectedElements_i::SetVertex( GEOM::GEOM_Object_ptr vertex )
1601   throw (SALOME::SALOME_Exception)
1602 {
1603   TopoDS_Shape shape = SMESH_Gen_i::GetSMESHGen()->GeomObjectToShape( vertex );
1604   if ( shape.IsNull() )
1605     THROW_SALOME_CORBA_EXCEPTION( "ConnectedElements_i::SetVertex(): NULL Vertex",
1606                                   SALOME::BAD_PARAM );
1607
1608   TopExp_Explorer v( shape, TopAbs_VERTEX );
1609   if ( !v.More() )
1610     THROW_SALOME_CORBA_EXCEPTION( "ConnectedElements_i::SetVertex(): empty vertex",
1611                                   SALOME::BAD_PARAM );
1612
1613   gp_Pnt p = BRep_Tool::Pnt( TopoDS::Vertex( v.Current() ));
1614   myConnectedElementsPtr->SetPoint( p.X(), p.Y(), p.Z() );
1615   //
1616   CORBA::String_var id = vertex->GetStudyEntry();
1617   myVertexID = id.in();
1618
1619   TPythonDump() << this << ".SetVertex( " << vertex << " )";
1620 }
1621
1622 void ConnectedElements_i::SetNode ( CORBA::Long nodeID )
1623   throw (SALOME::SALOME_Exception)
1624 {
1625   if ( nodeID < 1 )
1626     THROW_SALOME_CORBA_EXCEPTION( "ConnectedElements_i::SetNode(): nodeID must be > 0",
1627                                   SALOME::BAD_PARAM );
1628
1629   myConnectedElementsPtr->SetNode( nodeID );
1630   myVertexID.clear();
1631   TPythonDump() << this << ".SetNode( " << nodeID << " )";
1632 }
1633
1634 /*!
1635  * \brief This is a comfort method for Filter dialog
1636  */
1637 void ConnectedElements_i::SetThreshold ( const char*                             threshold,
1638                                          SMESH::ConnectedElements::ThresholdType type )
1639   throw (SALOME::SALOME_Exception)
1640 {
1641   if ( !threshold )
1642     THROW_SALOME_CORBA_EXCEPTION( "ConnectedElements_i::SetThreshold(): NULL threshold",
1643                                   SALOME::BAD_PARAM );
1644   switch ( type )
1645   {
1646   case SMESH::ConnectedElements::POINT: // read 3 node coordinates ///////////////////
1647     {
1648       std::vector< double > xyz;
1649       char* endptr;
1650       do
1651       {
1652         // skip a separator
1653         while ( *threshold &&
1654                 *threshold != '+' &&
1655                 *threshold != '-' &&
1656                 !isdigit( *threshold ))
1657           ++threshold;
1658         if ( !*threshold )
1659           break;
1660         // read a coordinate
1661         xyz.push_back( strtod( threshold, &endptr ));
1662         if ( threshold == endptr )
1663         {
1664           xyz.resize( xyz.size() - 1 );
1665           break;
1666         }
1667         threshold = endptr;
1668       }
1669       while ( xyz.size() < 3 );
1670
1671       if ( xyz.size() < 3 )
1672         THROW_SALOME_CORBA_EXCEPTION
1673           ( "ConnectedElements_i::SetThreshold(): invalid point coordinates", SALOME::BAD_PARAM );
1674
1675       SetPoint( xyz[0], xyz[1], xyz[2] );
1676       break;
1677     }
1678   case SMESH::ConnectedElements::VERTEX: // get a VERTEX by its entry /////////////////
1679     {
1680       SALOMEDS::Study_var study = SMESH_Gen_i::GetSMESHGen()->GetCurrentStudy();
1681       if ( study->_is_nil() )
1682         THROW_SALOME_CORBA_EXCEPTION
1683           ( "ConnectedElements_i::SetThreshold(): NULL current study", SALOME::BAD_PARAM );
1684       SALOMEDS::SObject_wrap sobj = study->FindObjectID( threshold );
1685       if ( sobj->_is_nil() )
1686         THROW_SALOME_CORBA_EXCEPTION
1687           ( "ConnectedElements_i::SetThreshold(): invalid vertex study entry", SALOME::BAD_PARAM );
1688       CORBA::Object_var        obj = sobj->GetObject();
1689       GEOM::GEOM_Object_var vertex = GEOM::GEOM_Object::_narrow( obj );
1690       if ( vertex->_is_nil() )
1691         THROW_SALOME_CORBA_EXCEPTION
1692           ( "ConnectedElements_i::SetThreshold(): no GEOM_Object in SObject", SALOME::BAD_PARAM );
1693       SetVertex( vertex );
1694       break;
1695     }
1696   case SMESH::ConnectedElements::NODE: // read a node ID ////////////////////////////
1697     {
1698       char* endptr;
1699       int id = strtol( threshold, &endptr, 10 );
1700       if ( threshold == endptr )
1701         THROW_SALOME_CORBA_EXCEPTION
1702           ( "ConnectedElements_i::SetThreshold(): invalid node ID", SALOME::BAD_PARAM );
1703       SetNode( id );
1704       break;
1705     }
1706   default:
1707     THROW_SALOME_CORBA_EXCEPTION
1708       ( "ConnectedElements_i::SetThreshold(): invalid ThresholdType", SALOME::BAD_PARAM );
1709   }
1710 }
1711
1712 char* ConnectedElements_i::GetThreshold ( SMESH::ConnectedElements::ThresholdType& type )
1713 {
1714   std::string threshold;
1715   if ( !myVertexID.empty() )
1716   {
1717     threshold = myVertexID;
1718     type      = SMESH::ConnectedElements::VERTEX;
1719   }
1720   else
1721   {
1722     std::vector<double> xyz = myConnectedElementsPtr->GetPoint();
1723     if ( xyz.size() == 3 )
1724     {
1725       threshold = SMESH_Comment( xyz[0] ) << "; " << xyz[1] << "; " << xyz[2];
1726       type      = SMESH::ConnectedElements::POINT;
1727     }
1728     else
1729     {
1730       threshold = SMESH_Comment( myConnectedElementsPtr->GetNode() );
1731       type      = SMESH::ConnectedElements::NODE;
1732     }
1733   }
1734   return CORBA::string_dup( threshold.c_str() );
1735 }
1736
1737 /*
1738   Class       : Comparator_i
1739   Description : Base class for comparators
1740 */
1741 Comparator_i::Comparator_i():
1742   myNumericalFunctor( NULL )
1743 {}
1744
1745 Comparator_i::~Comparator_i()
1746 {
1747   if ( myNumericalFunctor )
1748     myNumericalFunctor->UnRegister();
1749 }
1750
1751 void Comparator_i::SetMargin( CORBA::Double theValue )
1752 {
1753   myComparatorPtr->SetMargin( theValue );
1754   TPythonDump()<<this<<".SetMargin("<<theValue<<")";
1755 }
1756
1757 CORBA::Double Comparator_i::GetMargin()
1758 {
1759   return myComparatorPtr->GetMargin();
1760 }
1761
1762 void Comparator_i::SetNumFunctor( NumericalFunctor_ptr theFunct )
1763 {
1764   if ( myNumericalFunctor )
1765     myNumericalFunctor->UnRegister();
1766
1767   myNumericalFunctor = DownCast<NumericalFunctor_i*>(theFunct);
1768
1769   if ( myNumericalFunctor )
1770   {
1771     myComparatorPtr->SetNumFunctor( myNumericalFunctor->GetNumericalFunctor() );
1772     myNumericalFunctor->Register();
1773     TPythonDump()<<this<<".SetNumFunctor("<<myNumericalFunctor<<")";
1774   }
1775 }
1776
1777 Controls::ComparatorPtr Comparator_i::GetComparator()
1778 {
1779   return myComparatorPtr;
1780 }
1781
1782 NumericalFunctor_i* Comparator_i::GetNumFunctor_i()
1783 {
1784   return myNumericalFunctor;
1785 }
1786
1787
1788 /*
1789   Class       : LessThan_i
1790   Description : Comparator "<"
1791 */
1792 LessThan_i::LessThan_i()
1793 {
1794   myComparatorPtr.reset( new Controls::LessThan() );
1795   myFunctorPtr = myPredicatePtr = myComparatorPtr;
1796 }
1797
1798 FunctorType LessThan_i::GetFunctorType()
1799 {
1800   return SMESH::FT_LessThan;
1801 }
1802
1803
1804 /*
1805   Class       : MoreThan_i
1806   Description : Comparator ">"
1807 */
1808 MoreThan_i::MoreThan_i()
1809 {
1810   myComparatorPtr.reset( new Controls::MoreThan() );
1811   myFunctorPtr = myPredicatePtr = myComparatorPtr;
1812 }
1813
1814 FunctorType MoreThan_i::GetFunctorType()
1815 {
1816   return SMESH::FT_MoreThan;
1817 }
1818
1819
1820 /*
1821   Class       : EqualTo_i
1822   Description : Comparator "="
1823 */
1824 EqualTo_i::EqualTo_i()
1825 : myEqualToPtr( new Controls::EqualTo() )
1826 {
1827   myFunctorPtr = myPredicatePtr = myComparatorPtr = myEqualToPtr;
1828 }
1829
1830 void EqualTo_i::SetTolerance( CORBA::Double theToler )
1831 {
1832   myEqualToPtr->SetTolerance( theToler );
1833   TPythonDump()<<this<<".SetTolerance("<<theToler<<")";
1834 }
1835
1836 CORBA::Double EqualTo_i::GetTolerance()
1837 {
1838   return myEqualToPtr->GetTolerance();
1839 }
1840
1841 FunctorType EqualTo_i::GetFunctorType()
1842 {
1843   return SMESH::FT_EqualTo;
1844 }
1845
1846 /*
1847   Class       : LogicalNOT_i
1848   Description : Logical NOT predicate
1849 */
1850 LogicalNOT_i::LogicalNOT_i():
1851   myLogicalNOTPtr( new Controls::LogicalNOT() ),
1852   myPredicate( NULL )
1853 {
1854   myFunctorPtr = myPredicatePtr = myLogicalNOTPtr;
1855 }
1856
1857 LogicalNOT_i::~LogicalNOT_i()
1858 {
1859   if ( myPredicate )
1860     myPredicate->UnRegister();
1861 }
1862
1863 void LogicalNOT_i::SetPredicate( Predicate_ptr thePredicate )
1864 {
1865   if ( myPredicate )
1866     myPredicate->UnRegister();
1867
1868   myPredicate = SMESH::GetPredicate(thePredicate);
1869
1870   if ( myPredicate ){
1871     myLogicalNOTPtr->SetPredicate(myPredicate->GetPredicate());
1872     myPredicate->Register();
1873     TPythonDump()<<this<<".SetPredicate("<<myPredicate<<")";
1874   }
1875 }
1876
1877 FunctorType LogicalNOT_i::GetFunctorType()
1878 {
1879   return SMESH::FT_LogicalNOT;
1880 }
1881
1882 Predicate_i* LogicalNOT_i::GetPredicate_i()
1883 {
1884   return myPredicate;
1885 }
1886
1887
1888
1889 /*
1890   Class       : LogicalBinary_i
1891   Description : Base class for binary logical predicate
1892 */
1893 LogicalBinary_i::LogicalBinary_i()
1894 : myPredicate1( NULL ),
1895   myPredicate2( NULL )
1896 {}
1897
1898 LogicalBinary_i::~LogicalBinary_i()
1899 {
1900   if ( myPredicate1 )
1901     myPredicate1->UnRegister();
1902
1903   if ( myPredicate2 )
1904     myPredicate2->UnRegister();
1905 }
1906
1907 void LogicalBinary_i::SetMesh( SMESH_Mesh_ptr theMesh )
1908 {
1909   if ( myPredicate1 )
1910     myPredicate1->SetMesh( theMesh );
1911
1912   if ( myPredicate2 )
1913     myPredicate2->SetMesh( theMesh );
1914 }
1915
1916 void LogicalBinary_i::SetPredicate1( Predicate_ptr thePredicate )
1917 {
1918   if ( myPredicate1 )
1919     myPredicate1->UnRegister();
1920
1921   myPredicate1 = SMESH::GetPredicate(thePredicate);
1922
1923   if ( myPredicate1 ){
1924     myLogicalBinaryPtr->SetPredicate1(myPredicate1->GetPredicate());
1925     myPredicate1->Register();
1926     TPythonDump()<<this<<".SetPredicate1("<<myPredicate1<<")";
1927   }
1928 }
1929
1930 void LogicalBinary_i::SetPredicate2( Predicate_ptr thePredicate )
1931 {
1932   if ( myPredicate2 )
1933     myPredicate2->UnRegister();
1934
1935   myPredicate2 = SMESH::GetPredicate(thePredicate);
1936
1937   if ( myPredicate2 ){
1938     myLogicalBinaryPtr->SetPredicate2(myPredicate2->GetPredicate());
1939     myPredicate2->Register();
1940     TPythonDump()<<this<<".SetPredicate2("<<myPredicate2<<")";
1941   }
1942 }
1943
1944 Controls::LogicalBinaryPtr LogicalBinary_i::GetLogicalBinary()
1945 {
1946   return myLogicalBinaryPtr;
1947 }
1948
1949 Predicate_i* LogicalBinary_i::GetPredicate1_i()
1950 {
1951   return myPredicate1;
1952 }
1953 Predicate_i* LogicalBinary_i::GetPredicate2_i()
1954 {
1955   return myPredicate2;
1956 }
1957
1958
1959 /*
1960   Class       : LogicalAND_i
1961   Description : Logical AND
1962 */
1963 LogicalAND_i::LogicalAND_i()
1964 {
1965   myLogicalBinaryPtr.reset( new Controls::LogicalAND() );
1966   myFunctorPtr = myPredicatePtr = myLogicalBinaryPtr;
1967 }
1968
1969 FunctorType LogicalAND_i::GetFunctorType()
1970 {
1971   return SMESH::FT_LogicalAND;
1972 }
1973
1974
1975 /*
1976   Class       : LogicalOR_i
1977   Description : Logical OR
1978 */
1979 LogicalOR_i::LogicalOR_i()
1980 {
1981   myLogicalBinaryPtr.reset( new Controls::LogicalOR() );
1982   myFunctorPtr = myPredicatePtr = myLogicalBinaryPtr;
1983 }
1984
1985 FunctorType LogicalOR_i::GetFunctorType()
1986 {
1987   return SMESH::FT_LogicalOR;
1988 }
1989
1990
1991 /*
1992                             FILTER MANAGER
1993 */
1994
1995 FilterManager_i::FilterManager_i()
1996 : SALOME::GenericObj_i( SMESH_Gen_i::GetPOA() )
1997 {
1998   //Base class Salome_GenericObject do it inmplicitly by overriding PortableServer::POA_ptr _default_POA() method
1999   //PortableServer::ObjectId_var anObjectId =
2000   //  SMESH_Gen_i::GetPOA()->activate_object( this );
2001 }
2002
2003
2004 FilterManager_i::~FilterManager_i()
2005 {
2006   //TPythonDump()<<this<<".UnRegister()";
2007 }
2008
2009
2010 MinimumAngle_ptr FilterManager_i::CreateMinimumAngle()
2011 {
2012   SMESH::MinimumAngle_i* aServant = new SMESH::MinimumAngle_i();
2013   SMESH::MinimumAngle_var anObj = aServant->_this();
2014   TPythonDump()<<aServant<<" = "<<this<<".CreateMinimumAngle()";
2015   return anObj._retn();
2016 }
2017
2018
2019 AspectRatio_ptr FilterManager_i::CreateAspectRatio()
2020 {
2021   SMESH::AspectRatio_i* aServant = new SMESH::AspectRatio_i();
2022   SMESH::AspectRatio_var anObj = aServant->_this();
2023   TPythonDump()<<aServant<<" = "<<this<<".CreateAspectRatio()";
2024   return anObj._retn();
2025 }
2026
2027
2028 AspectRatio3D_ptr FilterManager_i::CreateAspectRatio3D()
2029 {
2030   SMESH::AspectRatio3D_i* aServant = new SMESH::AspectRatio3D_i();
2031   SMESH::AspectRatio3D_var anObj = aServant->_this();
2032   TPythonDump()<<aServant<<" = "<<this<<".CreateAspectRatio3D()";
2033   return anObj._retn();
2034 }
2035
2036
2037 Warping_ptr FilterManager_i::CreateWarping()
2038 {
2039   SMESH::Warping_i* aServant = new SMESH::Warping_i();
2040   SMESH::Warping_var anObj = aServant->_this();
2041   TPythonDump()<<aServant<<" = "<<this<<".CreateWarping()";
2042   return anObj._retn();
2043 }
2044
2045
2046 Taper_ptr FilterManager_i::CreateTaper()
2047 {
2048   SMESH::Taper_i* aServant = new SMESH::Taper_i();
2049   SMESH::Taper_var anObj = aServant->_this();
2050   TPythonDump()<<aServant<<" = "<<this<<".CreateTaper()";
2051   return anObj._retn();
2052 }
2053
2054
2055 Skew_ptr FilterManager_i::CreateSkew()
2056 {
2057   SMESH::Skew_i* aServant = new SMESH::Skew_i();
2058   SMESH::Skew_var anObj = aServant->_this();
2059   TPythonDump()<<aServant<<" = "<<this<<".CreateSkew()";
2060   return anObj._retn();
2061 }
2062
2063
2064 Area_ptr FilterManager_i::CreateArea()
2065 {
2066   SMESH::Area_i* aServant = new SMESH::Area_i();
2067   SMESH::Area_var anObj = aServant->_this();
2068   TPythonDump()<<aServant<<" = "<<this<<".CreateArea()";
2069   return anObj._retn();
2070 }
2071
2072
2073 Volume3D_ptr FilterManager_i::CreateVolume3D()
2074 {
2075   SMESH::Volume3D_i* aServant = new SMESH::Volume3D_i();
2076   SMESH::Volume3D_var anObj = aServant->_this();
2077   TPythonDump()<<aServant<<" = "<<this<<".CreateVolume3D()";
2078   return anObj._retn();
2079 }
2080
2081
2082 MaxElementLength2D_ptr FilterManager_i::CreateMaxElementLength2D()
2083 {
2084   SMESH::MaxElementLength2D_i* aServant = new SMESH::MaxElementLength2D_i();
2085   SMESH::MaxElementLength2D_var anObj = aServant->_this();
2086   TPythonDump()<<aServant<<" = "<<this<<".CreateMaxElementLength2D()";
2087   return anObj._retn();
2088 }
2089
2090
2091 MaxElementLength3D_ptr FilterManager_i::CreateMaxElementLength3D()
2092 {
2093   SMESH::MaxElementLength3D_i* aServant = new SMESH::MaxElementLength3D_i();
2094   SMESH::MaxElementLength3D_var anObj = aServant->_this();
2095   TPythonDump()<<aServant<<" = "<<this<<".CreateMaxElementLength3D()";
2096   return anObj._retn();
2097 }
2098
2099
2100 Length_ptr FilterManager_i::CreateLength()
2101 {
2102   SMESH::Length_i* aServant = new SMESH::Length_i();
2103   SMESH::Length_var anObj = aServant->_this();
2104   TPythonDump()<<aServant<<" = "<<this<<".CreateLength()";
2105   return anObj._retn();
2106 }
2107
2108 Length2D_ptr FilterManager_i::CreateLength2D()
2109 {
2110   SMESH::Length2D_i* aServant = new SMESH::Length2D_i();
2111   SMESH::Length2D_var anObj = aServant->_this();
2112   TPythonDump()<<aServant<<" = "<<this<<".CreateLength2D()";
2113   return anObj._retn();
2114 }
2115
2116 MultiConnection_ptr FilterManager_i::CreateMultiConnection()
2117 {
2118   SMESH::MultiConnection_i* aServant = new SMESH::MultiConnection_i();
2119   SMESH::MultiConnection_var anObj = aServant->_this();
2120   TPythonDump()<<aServant<<" = "<<this<<".CreateMultiConnection()";
2121   return anObj._retn();
2122 }
2123
2124 MultiConnection2D_ptr FilterManager_i::CreateMultiConnection2D()
2125 {
2126   SMESH::MultiConnection2D_i* aServant = new SMESH::MultiConnection2D_i();
2127   SMESH::MultiConnection2D_var anObj = aServant->_this();
2128   TPythonDump()<<aServant<<" = "<<this<<".CreateMultiConnection2D()";
2129   return anObj._retn();
2130 }
2131
2132 BallDiameter_ptr FilterManager_i::CreateBallDiameter()
2133 {
2134   SMESH::BallDiameter_i* aServant = new SMESH::BallDiameter_i();
2135   SMESH::BallDiameter_var anObj = aServant->_this();
2136   TPythonDump()<<aServant<<" = "<<this<<".CreateBallDiameter()";
2137   return anObj._retn();
2138 }
2139
2140 NodeConnectivityNumber_ptr FilterManager_i::CreateNodeConnectivityNumber()
2141 {
2142   SMESH::NodeConnectivityNumber_i* aServant = new SMESH::NodeConnectivityNumber_i();
2143   SMESH::NodeConnectivityNumber_var anObj = aServant->_this();
2144   TPythonDump()<<aServant<<" = "<<this<<".CreateNodeConnectivityNumber()";
2145   return anObj._retn();
2146 }
2147
2148 BelongToMeshGroup_ptr FilterManager_i::CreateBelongToMeshGroup()
2149 {
2150   SMESH::BelongToMeshGroup_i* aServant = new SMESH::BelongToMeshGroup_i();
2151   SMESH::BelongToMeshGroup_var anObj = aServant->_this();
2152   TPythonDump()<<aServant<<" = "<<this<<".CreateBelongToMeshGroup()";
2153   return anObj._retn();
2154 }
2155
2156 BelongToGeom_ptr FilterManager_i::CreateBelongToGeom()
2157 {
2158   SMESH::BelongToGeom_i* aServant = new SMESH::BelongToGeom_i();
2159   SMESH::BelongToGeom_var anObj = aServant->_this();
2160   TPythonDump()<<aServant<<" = "<<this<<".CreateBelongToGeom()";
2161   return anObj._retn();
2162 }
2163
2164 BelongToPlane_ptr FilterManager_i::CreateBelongToPlane()
2165 {
2166   SMESH::BelongToPlane_i* aServant = new SMESH::BelongToPlane_i();
2167   SMESH::BelongToPlane_var anObj = aServant->_this();
2168   TPythonDump()<<aServant<<" = "<<this<<".CreateBelongToPlane()";
2169   return anObj._retn();
2170 }
2171
2172 BelongToCylinder_ptr FilterManager_i::CreateBelongToCylinder()
2173 {
2174   SMESH::BelongToCylinder_i* aServant = new SMESH::BelongToCylinder_i();
2175   SMESH::BelongToCylinder_var anObj = aServant->_this();
2176   TPythonDump()<<aServant<<" = "<<this<<".CreateBelongToCylinder()";
2177   return anObj._retn();
2178 }
2179
2180 BelongToGenSurface_ptr FilterManager_i::CreateBelongToGenSurface()
2181 {
2182   SMESH::BelongToGenSurface_i* aServant = new SMESH::BelongToGenSurface_i();
2183   SMESH::BelongToGenSurface_var anObj = aServant->_this();
2184   TPythonDump()<<aServant<<" = "<<this<<".CreateBelongToGenSurface()";
2185   return anObj._retn();
2186 }
2187
2188 LyingOnGeom_ptr FilterManager_i::CreateLyingOnGeom()
2189 {
2190   SMESH::LyingOnGeom_i* aServant = new SMESH::LyingOnGeom_i();
2191   SMESH::LyingOnGeom_var anObj = aServant->_this();
2192   TPythonDump()<<aServant<<" = "<<this<<".CreateLyingOnGeom()";
2193   return anObj._retn();
2194 }
2195
2196 CoplanarFaces_ptr FilterManager_i::CreateCoplanarFaces()
2197 {
2198   SMESH::CoplanarFaces_i* aServant = new SMESH::CoplanarFaces_i();
2199   SMESH::CoplanarFaces_var anObj = aServant->_this();
2200   TPythonDump()<<aServant<<" = "<<this<<".CreateCoplanarFaces()";
2201   return anObj._retn();
2202 }
2203
2204 ConnectedElements_ptr FilterManager_i::CreateConnectedElements()
2205 {
2206   SMESH::ConnectedElements_i* aServant = new SMESH::ConnectedElements_i();
2207   SMESH::ConnectedElements_var anObj = aServant->_this();
2208   TPythonDump()<<aServant<<" = "<<this<<".CreateConnectedElements()";
2209   return anObj._retn();
2210 }
2211
2212 FreeBorders_ptr FilterManager_i::CreateFreeBorders()
2213 {
2214   SMESH::FreeBorders_i* aServant = new SMESH::FreeBorders_i();
2215   SMESH::FreeBorders_var anObj = aServant->_this();
2216   TPythonDump()<<aServant<<" = "<<this<<".CreateFreeBorders()";
2217   return anObj._retn();
2218 }
2219
2220 FreeEdges_ptr FilterManager_i::CreateFreeEdges()
2221 {
2222   SMESH::FreeEdges_i* aServant = new SMESH::FreeEdges_i();
2223   SMESH::FreeEdges_var anObj = aServant->_this();
2224   TPythonDump()<<aServant<<" = "<<this<<".CreateFreeEdges()";
2225   return anObj._retn();
2226 }
2227
2228 FreeFaces_ptr FilterManager_i::CreateFreeFaces()
2229 {
2230   SMESH::FreeFaces_i* aServant = new SMESH::FreeFaces_i();
2231   SMESH::FreeFaces_var anObj = aServant->_this();
2232   TPythonDump()<<aServant<<" = "<<this<<".CreateFreeFaces()";
2233   return anObj._retn();
2234 }
2235
2236 FreeNodes_ptr FilterManager_i::CreateFreeNodes()
2237 {
2238   SMESH::FreeNodes_i* aServant = new SMESH::FreeNodes_i();
2239   SMESH::FreeNodes_var anObj = aServant->_this();
2240   TPythonDump()<<aServant<<" = "<<this<<".CreateFreeNodes()";
2241   return anObj._retn();
2242 }
2243
2244 EqualNodes_ptr FilterManager_i::CreateEqualNodes()
2245 {
2246   SMESH::EqualNodes_i* aServant = new SMESH::EqualNodes_i();
2247   SMESH::EqualNodes_var anObj = aServant->_this();
2248   TPythonDump()<<aServant<<" = "<<this<<".CreateEqualNodes()";
2249   return anObj._retn();
2250 }
2251
2252 EqualEdges_ptr FilterManager_i::CreateEqualEdges()
2253 {
2254   SMESH::EqualEdges_i* aServant = new SMESH::EqualEdges_i();
2255   SMESH::EqualEdges_var anObj = aServant->_this();
2256   TPythonDump()<<aServant<<" = "<<this<<".CreateEqualEdges()";
2257   return anObj._retn();
2258 }
2259 EqualFaces_ptr FilterManager_i::CreateEqualFaces()
2260 {
2261   SMESH::EqualFaces_i* aServant = new SMESH::EqualFaces_i();
2262   SMESH::EqualFaces_var anObj = aServant->_this();
2263   TPythonDump()<<aServant<<" = "<<this<<".CreateEqualFaces()";
2264   return anObj._retn();
2265 }
2266 EqualVolumes_ptr FilterManager_i::CreateEqualVolumes()
2267 {
2268   SMESH::EqualVolumes_i* aServant = new SMESH::EqualVolumes_i();
2269   SMESH::EqualVolumes_var anObj = aServant->_this();
2270   TPythonDump()<<aServant<<" = "<<this<<".CreateEqualVolumes()";
2271   return anObj._retn();
2272 }
2273
2274 RangeOfIds_ptr FilterManager_i::CreateRangeOfIds()
2275 {
2276   SMESH::RangeOfIds_i* aServant = new SMESH::RangeOfIds_i();
2277   SMESH::RangeOfIds_var anObj = aServant->_this();
2278   TPythonDump()<<aServant<<" = "<<this<<".CreateRangeOfIds()";
2279   return anObj._retn();
2280 }
2281
2282 BadOrientedVolume_ptr FilterManager_i::CreateBadOrientedVolume()
2283 {
2284   SMESH::BadOrientedVolume_i* aServant = new SMESH::BadOrientedVolume_i();
2285   SMESH::BadOrientedVolume_var anObj = aServant->_this();
2286   TPythonDump()<<aServant<<" = "<<this<<".CreateBadOrientedVolume()";
2287   return anObj._retn();
2288 }
2289
2290 BareBorderVolume_ptr FilterManager_i::CreateBareBorderVolume()
2291 {
2292   SMESH::BareBorderVolume_i* aServant = new SMESH::BareBorderVolume_i();
2293   SMESH::BareBorderVolume_var anObj = aServant->_this();
2294   TPythonDump()<<aServant<<" = "<<this<<".CreateBareBorderVolume()";
2295   return anObj._retn();
2296 }
2297
2298 BareBorderFace_ptr FilterManager_i::CreateBareBorderFace()
2299 {
2300   SMESH::BareBorderFace_i* aServant = new SMESH::BareBorderFace_i();
2301   SMESH::BareBorderFace_var anObj = aServant->_this();
2302   TPythonDump()<<aServant<<" = "<<this<<".CreateBareBorderFace()";
2303   return anObj._retn();
2304 }
2305
2306 OverConstrainedVolume_ptr FilterManager_i::CreateOverConstrainedVolume()
2307 {
2308   SMESH::OverConstrainedVolume_i* aServant = new SMESH::OverConstrainedVolume_i();
2309   SMESH::OverConstrainedVolume_var anObj = aServant->_this();
2310   TPythonDump()<<aServant<<" = "<<this<<".CreateOverConstrainedVolume()";
2311   return anObj._retn();
2312 }
2313
2314 OverConstrainedFace_ptr FilterManager_i::CreateOverConstrainedFace()
2315 {
2316   SMESH::OverConstrainedFace_i* aServant = new SMESH::OverConstrainedFace_i();
2317   SMESH::OverConstrainedFace_var anObj = aServant->_this();
2318   TPythonDump()<<aServant<<" = "<<this<<".CreateOverConstrainedFace()";
2319   return anObj._retn();
2320 }
2321
2322 LessThan_ptr FilterManager_i::CreateLessThan()
2323 {
2324   SMESH::LessThan_i* aServant = new SMESH::LessThan_i();
2325   SMESH::LessThan_var anObj = aServant->_this();
2326   TPythonDump()<<aServant<<" = "<<this<<".CreateLessThan()";
2327   return anObj._retn();
2328 }
2329
2330 MoreThan_ptr FilterManager_i::CreateMoreThan()
2331 {
2332   SMESH::MoreThan_i* aServant = new SMESH::MoreThan_i();
2333   SMESH::MoreThan_var anObj = aServant->_this();
2334   TPythonDump()<<aServant<<" = "<<this<<".CreateMoreThan()";
2335   return anObj._retn();
2336 }
2337
2338 EqualTo_ptr FilterManager_i::CreateEqualTo()
2339 {
2340   SMESH::EqualTo_i* aServant = new SMESH::EqualTo_i();
2341   SMESH::EqualTo_var anObj = aServant->_this();
2342   TPythonDump()<<aServant<<" = "<<this<<".CreateEqualTo()";
2343   return anObj._retn();
2344 }
2345
2346 LogicalNOT_ptr FilterManager_i::CreateLogicalNOT()
2347 {
2348   SMESH::LogicalNOT_i* aServant = new SMESH::LogicalNOT_i();
2349   SMESH::LogicalNOT_var anObj = aServant->_this();
2350   TPythonDump()<<aServant<<" = "<<this<<".CreateLogicalNOT()";
2351   return anObj._retn();
2352 }
2353
2354 LogicalAND_ptr FilterManager_i::CreateLogicalAND()
2355 {
2356   SMESH::LogicalAND_i* aServant = new SMESH::LogicalAND_i();
2357   SMESH::LogicalAND_var anObj = aServant->_this();
2358   TPythonDump()<<aServant<<" = "<<this<<".CreateLogicalAND()";
2359   return anObj._retn();
2360 }
2361
2362 LogicalOR_ptr FilterManager_i::CreateLogicalOR()
2363 {
2364   SMESH::LogicalOR_i* aServant = new SMESH::LogicalOR_i();
2365   SMESH::LogicalOR_var anObj = aServant->_this();
2366   TPythonDump()<<aServant<<" = "<<this<<".CreateLogicalOR()";
2367   return anObj._retn();
2368 }
2369
2370 LinearOrQuadratic_ptr FilterManager_i::CreateLinearOrQuadratic()
2371 {
2372   SMESH::LinearOrQuadratic_i* aServant = new SMESH::LinearOrQuadratic_i();
2373   SMESH::LinearOrQuadratic_var anObj = aServant->_this();
2374   TPythonDump()<<aServant<<" = "<<this<<".CreateLinearOrQuadratic()";
2375   return anObj._retn();
2376 }
2377
2378 GroupColor_ptr FilterManager_i::CreateGroupColor()
2379 {
2380   SMESH::GroupColor_i* aServant = new SMESH::GroupColor_i();
2381   SMESH::GroupColor_var anObj = aServant->_this();
2382   TPythonDump()<<aServant<<" = "<<this<<".CreateGroupColor()";
2383   return anObj._retn();
2384 }
2385
2386 ElemGeomType_ptr FilterManager_i::CreateElemGeomType()
2387 {
2388   SMESH::ElemGeomType_i* aServant = new SMESH::ElemGeomType_i();
2389   SMESH::ElemGeomType_var anObj = aServant->_this();
2390   TPythonDump()<<aServant<<" = "<<this<<".CreateElemGeomType()";
2391   return anObj._retn();
2392 }
2393
2394 ElemEntityType_ptr FilterManager_i::CreateElemEntityType()
2395 {
2396   SMESH::ElemEntityType_i* aServant = new SMESH::ElemEntityType_i();
2397   SMESH::ElemEntityType_var anObj = aServant->_this();
2398   TPythonDump()<<aServant<<" = "<<this<<".CreateElemEntityType()";
2399   return anObj._retn();
2400 }
2401
2402 Filter_ptr FilterManager_i::CreateFilter()
2403 {
2404   SMESH::Filter_i* aServant = new SMESH::Filter_i();
2405   SMESH::Filter_var anObj = aServant->_this();
2406   TPythonDump()<<aServant<<" = "<<this<<".CreateFilter()";
2407   return anObj._retn();
2408 }
2409
2410 FilterLibrary_ptr FilterManager_i::LoadLibrary( const char* aFileName )
2411 {
2412   SMESH::FilterLibrary_i* aServant = new SMESH::FilterLibrary_i( aFileName );
2413   SMESH::FilterLibrary_var anObj = aServant->_this();
2414   TPythonDump()<<aServant<<" = "<<this<<".LoadLibrary('"<<aFileName<<"')";
2415   return anObj._retn();
2416 }
2417
2418 FilterLibrary_ptr FilterManager_i::CreateLibrary()
2419 {
2420   SMESH::FilterLibrary_i* aServant = new SMESH::FilterLibrary_i();
2421   SMESH::FilterLibrary_var anObj = aServant->_this();
2422   TPythonDump()<<aServant<<" = "<<this<<".CreateLibrary()";
2423   return anObj._retn();
2424 }
2425
2426 CORBA::Boolean FilterManager_i::DeleteLibrary( const char* aFileName )
2427 {
2428   TPythonDump()<<this<<".DeleteLibrary("<<aFileName<<")";
2429   return remove( aFileName ) ? false : true;
2430 }
2431
2432 //=============================================================================
2433 /*!
2434  *  SMESH_Gen_i::CreateFilterManager
2435  *
2436  *  Create filter manager
2437  */
2438 //=============================================================================
2439
2440 SMESH::FilterManager_ptr SMESH_Gen_i::CreateFilterManager()
2441 {
2442   SMESH::FilterManager_i* aFilter = new SMESH::FilterManager_i();
2443   SMESH::FilterManager_var anObj = aFilter->_this();
2444   return anObj._retn();
2445 }
2446
2447
2448 /*
2449                               FILTER
2450 */
2451
2452 //=======================================================================
2453 // name    : Filter_i::Filter_i
2454 // Purpose : Constructor
2455 //=======================================================================
2456 Filter_i::Filter_i()
2457 : myPredicate( NULL )
2458 {}
2459
2460 //=======================================================================
2461 // name    : Filter_i::~Filter_i
2462 // Purpose : Destructor
2463 //=======================================================================
2464 Filter_i::~Filter_i()
2465 {
2466   if ( myPredicate )
2467     myPredicate->UnRegister();
2468
2469   if(!CORBA::is_nil(myMesh))
2470     myMesh->UnRegister();
2471
2472   myPredicate = 0;
2473   FindBaseObjects();
2474 }
2475
2476 //=======================================================================
2477 // name    : Filter_i::SetPredicate
2478 // Purpose : Set predicate
2479 //=======================================================================
2480 void Filter_i::SetPredicate( Predicate_ptr thePredicate )
2481 {
2482   if ( myPredicate )
2483     myPredicate->UnRegister();
2484
2485   myPredicate = SMESH::GetPredicate(thePredicate);
2486
2487   if ( myPredicate )
2488   {
2489     myFilter.SetPredicate( myPredicate->GetPredicate() );
2490     myPredicate->Register();
2491     if ( const SMDS_Mesh* aMesh = MeshPtr2SMDSMesh(myMesh))
2492       myPredicate->GetPredicate()->SetMesh( aMesh );
2493     TPythonDump()<<this<<".SetPredicate("<<myPredicate<<")";
2494   }
2495   NotifyerAndWaiter::Modified();
2496 }
2497
2498 //=======================================================================
2499 // name    : Filter_i::GetElementType
2500 // Purpose : Get entity type
2501 //=======================================================================
2502 SMESH::ElementType Filter_i::GetElementType()
2503 {
2504   return myPredicate != 0 ? myPredicate->GetElementType() : SMESH::ALL;
2505 }
2506
2507 //=======================================================================
2508 // name    : Filter_i::SetMesh
2509 // Purpose : Set mesh
2510 //=======================================================================
2511 void
2512 Filter_i::
2513 SetMesh( SMESH_Mesh_ptr theMesh )
2514 {
2515   if(!CORBA::is_nil(theMesh))
2516     theMesh->Register();
2517
2518   if(!CORBA::is_nil(myMesh))
2519     myMesh->UnRegister();
2520
2521   myMesh = SMESH_Mesh::_duplicate( theMesh );
2522   TPythonDump()<<this<<".SetMesh("<<theMesh<<")";
2523
2524   if ( myPredicate )
2525     if ( const SMDS_Mesh* aMesh = MeshPtr2SMDSMesh(theMesh))
2526       myPredicate->GetPredicate()->SetMesh( aMesh );
2527 }
2528
2529 SMESH::long_array*
2530 Filter_i::
2531 GetIDs()
2532 {
2533   return GetElementsId(myMesh);
2534 }
2535
2536 //=======================================================================
2537 // name    : Filter_i::GetElementsId
2538 // Purpose : Get ids of entities
2539 //=======================================================================
2540 void
2541 Filter_i::
2542 GetElementsId( Predicate_i* thePredicate,
2543                const SMDS_Mesh* theMesh,
2544                Controls::Filter::TIdSequence& theSequence )
2545 {
2546   if (thePredicate)
2547     Controls::Filter::GetElementsId(theMesh,thePredicate->GetPredicate(),theSequence);
2548 }
2549
2550 void
2551 Filter_i::
2552 GetElementsId( Predicate_i* thePredicate,
2553                SMESH_Mesh_ptr theMesh,
2554                Controls::Filter::TIdSequence& theSequence )
2555 {
2556   if (thePredicate)
2557     if(const SMDS_Mesh* aMesh = MeshPtr2SMDSMesh(theMesh))
2558       Controls::Filter::GetElementsId(aMesh,thePredicate->GetPredicate(),theSequence);
2559 }
2560
2561 SMESH::long_array*
2562 Filter_i::
2563 GetElementsId( SMESH_Mesh_ptr theMesh )
2564 {
2565   SMESH::long_array_var anArray = new SMESH::long_array;
2566   if(!CORBA::is_nil(theMesh) && myPredicate){
2567     theMesh->Load();
2568     Controls::Filter::TIdSequence aSequence;
2569     GetElementsId(myPredicate,theMesh,aSequence);
2570     long i = 0, iEnd = aSequence.size();
2571     anArray->length( iEnd );
2572     for ( ; i < iEnd; i++ )
2573       anArray[ i ] = aSequence[i];
2574   }
2575   return anArray._retn();
2576 }
2577
2578 //=============================================================================
2579 /*!
2580  * \brief Returns number of mesh elements per each \a EntityType
2581  */
2582 //=============================================================================
2583
2584 SMESH::long_array* ::Filter_i::GetMeshInfo()
2585 {
2586   SMESH::long_array_var aRes = new SMESH::long_array();
2587   aRes->length(SMESH::Entity_Last);
2588   for (int i = 0; i < SMESH::Entity_Last; i++)
2589     aRes[i] = 0;
2590
2591   if ( !CORBA::is_nil(myMesh) && myPredicate )
2592   {
2593     const SMDS_Mesh*  aMesh = MeshPtr2SMDSMesh(myMesh);
2594     SMDS_ElemIteratorPtr it = aMesh->elementsIterator( SMDSAbs_ElementType( GetElementType() ));
2595     while ( it->more() )
2596     {
2597       const SMDS_MeshElement* anElem = it->next();
2598       if ( myPredicate->IsSatisfy( anElem->GetID() ) )
2599         aRes[ anElem->GetEntityType() ]++;
2600     }
2601   }
2602
2603   return aRes._retn();
2604 }
2605
2606 //=============================================================================
2607 /*!
2608  * \brief Returns number of mesh elements of each \a ElementType
2609  */
2610 //=============================================================================
2611
2612 SMESH::long_array* ::Filter_i::GetNbElementsByType()
2613 {
2614   SMESH::long_array_var aRes = new SMESH::long_array();
2615   aRes->length(SMESH::NB_ELEMENT_TYPES);
2616   for (int i = 0; i < SMESH::NB_ELEMENT_TYPES; i++)
2617     aRes[i] = 0;
2618
2619   if ( !CORBA::is_nil(myMesh) && myPredicate ) {
2620     const SMDS_Mesh*  aMesh = MeshPtr2SMDSMesh(myMesh);
2621     SMDS_ElemIteratorPtr it = aMesh->elementsIterator( SMDSAbs_ElementType( GetElementType() ));
2622     CORBA::Long& nbElems = aRes[ GetElementType() ];
2623     while ( it->more() )
2624     {
2625       const SMDS_MeshElement* anElem = it->next();
2626       if ( myPredicate->IsSatisfy( anElem->GetID() ) )
2627         nbElems++;
2628     }
2629   }
2630
2631   return aRes._retn();
2632 }
2633
2634
2635 //================================================================================
2636 /*!
2637  * \brief Return GetElementType() within an array
2638  * Implement SMESH_IDSource interface
2639  */
2640 //================================================================================
2641
2642 SMESH::array_of_ElementType* Filter_i::GetTypes()
2643 {
2644   SMESH::array_of_ElementType_var types = new SMESH::array_of_ElementType;
2645
2646   // check if any element passes through the filter
2647   if ( !CORBA::is_nil(myMesh) && myPredicate )
2648   {
2649     const SMDS_Mesh* aMesh = MeshPtr2SMDSMesh(myMesh);
2650     SMDS_ElemIteratorPtr it = aMesh->elementsIterator( SMDSAbs_ElementType( GetElementType() ));
2651     bool satisfies = false;
2652     while ( !satisfies && it->more() )
2653       satisfies = myPredicate->IsSatisfy( it->next()->GetID() );
2654     if ( satisfies ) {
2655       types->length( 1 );
2656       types[0] = GetElementType();
2657     }
2658   }
2659   return types._retn();
2660 }
2661
2662 //=======================================================================
2663 //function : GetMesh
2664 //purpose  : Returns mesh
2665 //=======================================================================
2666
2667 SMESH::SMESH_Mesh_ptr Filter_i::GetMesh()
2668 {
2669   return SMESH_Mesh::_duplicate( myMesh );
2670 }
2671
2672 //=======================================================================
2673 //function : GetVtkUgStream
2674 //purpose  : Return data vtk unstructured grid (not implemented)
2675 //=======================================================================
2676
2677 SALOMEDS::TMPFile* Filter_i::GetVtkUgStream()
2678 {
2679   SALOMEDS::TMPFile_var SeqFile;
2680   return SeqFile._retn();
2681 }
2682 //=======================================================================
2683 // name    : getCriteria
2684 // Purpose : Retrieve criterions from predicate
2685 //=======================================================================
2686 static inline void getPrediacates( Predicate_i*                thePred,
2687                                    std::vector<Predicate_i*> & thePredVec )
2688 {
2689   const int aFType = thePred->GetFunctorType();
2690
2691   switch ( aFType )
2692   {
2693   case FT_LogicalNOT:
2694   {
2695     Predicate_i* aPred = ( dynamic_cast<LogicalNOT_i*>( thePred ) )->GetPredicate_i();
2696     getPrediacates( aPred, thePredVec );
2697     break;
2698   }
2699   case FT_LogicalAND:
2700   case FT_LogicalOR:
2701   {
2702     Predicate_i* aPred1 = ( dynamic_cast<LogicalBinary_i*>( thePred ) )->GetPredicate1_i();
2703     Predicate_i* aPred2 = ( dynamic_cast<LogicalBinary_i*>( thePred ) )->GetPredicate2_i();
2704     getPrediacates( aPred1, thePredVec );
2705     getPrediacates( aPred2, thePredVec );
2706     break;
2707    }
2708   default:;
2709   }
2710   thePredVec.push_back( thePred );
2711 }
2712
2713 //=======================================================================
2714 // name    : getCriteria
2715 // Purpose : Retrieve criterions from predicate
2716 //=======================================================================
2717 static inline bool getCriteria( Predicate_i*                thePred,
2718                                 SMESH::Filter::Criteria_out theCriteria )
2719 {
2720   const int aFType = thePred->GetFunctorType();
2721
2722   switch ( aFType )
2723   {
2724   case FT_LogicalNOT:
2725     {
2726       Predicate_i* aPred = ( dynamic_cast<LogicalNOT_i*>( thePred ) )->GetPredicate_i();
2727       getCriteria( aPred, theCriteria );
2728       theCriteria[ theCriteria->length() - 1 ].UnaryOp = FT_LogicalNOT;
2729     }
2730     return true;
2731
2732   case FT_LogicalAND:
2733   case FT_LogicalOR:
2734     {
2735       Predicate_i* aPred1 = ( dynamic_cast<LogicalBinary_i*>( thePred ) )->GetPredicate1_i();
2736       Predicate_i* aPred2 = ( dynamic_cast<LogicalBinary_i*>( thePred ) )->GetPredicate2_i();
2737       if ( !getCriteria( aPred1, theCriteria ) )
2738         return false;
2739       theCriteria[ theCriteria->length() - 1 ].BinaryOp = aFType;
2740       return getCriteria( aPred2, theCriteria );
2741     }
2742   case FT_Undefined:
2743     return false;
2744   }
2745
2746   // resize theCriteria
2747   CORBA::ULong i = theCriteria->length();
2748   theCriteria->length( i + 1 );
2749   theCriteria[ i ] = createCriterion();
2750
2751   // set members of the added Criterion
2752
2753   theCriteria[ i ].Type = aFType;
2754   theCriteria[ i ].TypeOfElement = thePred->GetElementType();
2755
2756   switch ( aFType )
2757   {
2758   case FT_FreeBorders:
2759   case FT_FreeEdges:
2760   case FT_FreeFaces:
2761   case FT_LinearOrQuadratic:
2762   case FT_FreeNodes:
2763   case FT_EqualEdges:
2764   case FT_EqualFaces:
2765   case FT_EqualVolumes:
2766   case FT_BadOrientedVolume:
2767   case FT_BareBorderVolume:
2768   case FT_BareBorderFace:
2769   case FT_OverConstrainedVolume:
2770   case FT_OverConstrainedFace:
2771     {
2772       return true;
2773     }
2774   case FT_BelongToMeshGroup:
2775     {
2776       BelongToMeshGroup_i* aPred = dynamic_cast<BelongToMeshGroup_i*>( thePred );
2777       SMESH::SMESH_GroupBase_var grp = aPred->GetGroup();
2778       if ( !grp->_is_nil() )
2779       {
2780         theCriteria[ i ].ThresholdStr = grp->GetName();
2781         theCriteria[ i ].ThresholdID  = aPred->GetGroupID().c_str();
2782       }
2783       return true;
2784     }
2785   case FT_BelongToGeom:
2786     {
2787       BelongToGeom_i* aPred = dynamic_cast<BelongToGeom_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_BelongToPlane:
2794   case FT_BelongToCylinder:
2795   case FT_BelongToGenSurface:
2796     {
2797       BelongToSurface_i* aPred = dynamic_cast<BelongToSurface_i*>( thePred );
2798       theCriteria[ i ].ThresholdStr  = aPred->GetShapeName();
2799       theCriteria[ i ].ThresholdID   = aPred->GetShapeID();
2800       theCriteria[ i ].Tolerance     = aPred->GetTolerance();
2801       return true;
2802     }
2803   case FT_LyingOnGeom:
2804     {
2805       LyingOnGeom_i* aPred = dynamic_cast<LyingOnGeom_i*>( thePred );
2806       theCriteria[ i ].ThresholdStr  = aPred->GetShapeName();
2807       theCriteria[ i ].ThresholdID   = aPred->GetShapeID();
2808       theCriteria[ i ].Tolerance     = aPred->GetTolerance();
2809       return true;
2810     }
2811   case FT_CoplanarFaces:
2812     {
2813       CoplanarFaces_i* aPred = dynamic_cast<CoplanarFaces_i*>( thePred );
2814       theCriteria[ i ].ThresholdID   = aPred->GetFaceAsString();
2815       theCriteria[ i ].Tolerance     = aPred->GetTolerance();
2816       return true;
2817     }
2818   case FT_ConnectedElements:
2819     {
2820       ConnectedElements_i* aPred = dynamic_cast<ConnectedElements_i*>( thePred );
2821       SMESH::ConnectedElements::ThresholdType type;
2822       CORBA::String_var threshold = aPred->GetThreshold( type );
2823       switch ( type ) {
2824       case SMESH::ConnectedElements::POINT:
2825         theCriteria[ i ].ThresholdStr = threshold; break;
2826       case SMESH::ConnectedElements::VERTEX:
2827         theCriteria[ i ].ThresholdID = threshold; break;
2828       case SMESH::ConnectedElements::NODE:
2829         theCriteria[ i ].Threshold = atof( threshold.in() ); break;
2830       default:;
2831       }
2832       return true;
2833     }
2834   case FT_EqualNodes:
2835     {
2836       EqualNodes_i* aPred = dynamic_cast<EqualNodes_i*>( thePred );
2837       theCriteria[ i ].Tolerance  = aPred->GetTolerance();
2838       return true;
2839     }
2840   case FT_RangeOfIds:
2841     {
2842       RangeOfIds_i* aPred = dynamic_cast<RangeOfIds_i*>( thePred );
2843       theCriteria[ i ].ThresholdStr  = aPred->GetRangeStr();
2844       return true;
2845     }
2846   case FT_LessThan:
2847   case FT_MoreThan:
2848   case FT_EqualTo:
2849     {
2850       Comparator_i* aCompar = dynamic_cast<Comparator_i*>( thePred );
2851       theCriteria[ i ].Type      = aCompar->GetNumFunctor_i()->GetFunctorType();
2852       theCriteria[ i ].Compare   = aFType;
2853       theCriteria[ i ].Threshold = aCompar->GetMargin();
2854       if ( aFType == FT_EqualTo )
2855       {
2856         EqualTo_i* aCompar = dynamic_cast<EqualTo_i*>( thePred );
2857         theCriteria[ i ].Tolerance = aCompar->GetTolerance();
2858       }
2859       return true;
2860     }
2861   case FT_GroupColor:
2862     {
2863       GroupColor_i* aPred = dynamic_cast<GroupColor_i*>( thePred );
2864       theCriteria[ i ].ThresholdStr = aPred->GetColorStr();
2865       return true;
2866     }
2867   case FT_ElemGeomType:
2868     {
2869       ElemGeomType_i* aPred = dynamic_cast<ElemGeomType_i*>( thePred );
2870       theCriteria[ i ].Threshold = (double)aPred->GetGeometryType();
2871       return true;
2872     }
2873   case FT_EntityType:
2874     {
2875       ElemEntityType_i* aPred = dynamic_cast<ElemEntityType_i*>( thePred );
2876       theCriteria[ i ].Threshold = (double)aPred->GetEntityType();
2877       return true;
2878     }
2879   default:
2880     return false;
2881   }
2882 }
2883
2884 //=======================================================================
2885 // name    : Filter_i::GetCriteria
2886 // Purpose : Retrieve criterions from predicate
2887 //=======================================================================
2888 CORBA::Boolean Filter_i::GetCriteria( SMESH::Filter::Criteria_out theCriteria )
2889 {
2890   theCriteria = new SMESH::Filter::Criteria;
2891   return myPredicate != 0 ? getCriteria( myPredicate, theCriteria ) : true;
2892 }
2893
2894 //=======================================================================
2895 // name    : Filter_i::SetCriteria
2896 // Purpose : Create new predicate and set criterions in it
2897 //=======================================================================
2898 CORBA::Boolean Filter_i::SetCriteria( const SMESH::Filter::Criteria& theCriteria )
2899 {
2900   SetPredicate( SMESH::Predicate::_nil() );
2901
2902   SMESH::FilterManager_i* aFilter = new SMESH::FilterManager_i();
2903   FilterManager_ptr aFilterMgr = aFilter->_this();
2904
2905   // CREATE two lists ( PREDICATES  and LOG OP )
2906
2907   // Criterion
2908   TPythonDump()<<"aCriteria = []";
2909   std::list<SMESH::Predicate_ptr> aPredicates;
2910   std::list<int>                  aBinaries;
2911   for ( int i = 0, n = theCriteria.length(); i < n; i++ )
2912   {
2913     int         aCriterion    = theCriteria[ i ].Type;
2914     int         aCompare      = theCriteria[ i ].Compare;
2915     double      aThreshold    = theCriteria[ i ].Threshold;
2916     const char* aThresholdStr = theCriteria[ i ].ThresholdStr;
2917     const char* aThresholdID  = theCriteria[ i ].ThresholdID;
2918     int         aUnary        = theCriteria[ i ].UnaryOp;
2919     int         aBinary       = theCriteria[ i ].BinaryOp;
2920     double      aTolerance    = theCriteria[ i ].Tolerance;
2921     ElementType aTypeOfElem   = theCriteria[ i ].TypeOfElement;
2922     long        aPrecision    = theCriteria[ i ].Precision;
2923
2924     {
2925       TPythonDump pd;
2926       pd << "aCriterion = SMESH.Filter.Criterion("
2927          << aCriterion    << ", "
2928          << aCompare      << ", "
2929          << aThreshold    << ", '"
2930          << aThresholdStr << "', '"
2931          << aThresholdID  << "', "
2932          << aUnary        << ", "
2933          << aBinary       << ", "
2934          << aTolerance    << ", "
2935          << aTypeOfElem   << ", "
2936          << aPrecision    << ")";
2937     }
2938     TPythonDump pd;
2939
2940     SMESH::Predicate_ptr aPredicate = SMESH::Predicate::_nil();
2941     SMESH::NumericalFunctor_ptr aFunctor = SMESH::NumericalFunctor::_nil();
2942
2943     switch ( aCriterion )
2944     {
2945       // Functors
2946
2947       case SMESH::FT_MultiConnection:
2948         aFunctor = aFilterMgr->CreateMultiConnection();
2949         break;
2950       case SMESH::FT_MultiConnection2D:
2951         aFunctor = aFilterMgr->CreateMultiConnection2D();
2952         break;
2953       case SMESH::FT_Length:
2954         aFunctor = aFilterMgr->CreateLength();
2955         break;
2956       case SMESH::FT_Length2D:
2957         aFunctor = aFilterMgr->CreateLength2D();
2958         break;
2959       case SMESH::FT_AspectRatio:
2960         aFunctor = aFilterMgr->CreateAspectRatio();
2961         break;
2962       case SMESH::FT_AspectRatio3D:
2963         aFunctor = aFilterMgr->CreateAspectRatio3D();
2964         break;
2965       case SMESH::FT_Warping:
2966         aFunctor = aFilterMgr->CreateWarping();
2967         break;
2968       case SMESH::FT_MinimumAngle:
2969         aFunctor = aFilterMgr->CreateMinimumAngle();
2970         break;
2971       case SMESH::FT_Taper:
2972         aFunctor = aFilterMgr->CreateTaper();
2973         break;
2974       case SMESH::FT_Skew:
2975         aFunctor = aFilterMgr->CreateSkew();
2976         break;
2977       case SMESH::FT_Area:
2978         aFunctor = aFilterMgr->CreateArea();
2979         break;
2980       case SMESH::FT_Volume3D:
2981         aFunctor = aFilterMgr->CreateVolume3D();
2982         break;
2983       case SMESH::FT_MaxElementLength2D:
2984         aFunctor = aFilterMgr->CreateMaxElementLength2D();
2985         break;
2986       case SMESH::FT_MaxElementLength3D:
2987         aFunctor = aFilterMgr->CreateMaxElementLength3D();
2988         break;
2989       case SMESH::FT_BallDiameter:
2990         aFunctor = aFilterMgr->CreateBallDiameter();
2991         break;
2992       case SMESH::FT_NodeConnectivityNumber:
2993         aFunctor = aFilterMgr->CreateNodeConnectivityNumber();
2994         break;
2995
2996       // Predicates
2997
2998       case SMESH::FT_FreeBorders:
2999         aPredicate = aFilterMgr->CreateFreeBorders();
3000         break;
3001       case SMESH::FT_FreeEdges:
3002         aPredicate = aFilterMgr->CreateFreeEdges();
3003         break;
3004       case SMESH::FT_FreeFaces:
3005         aPredicate = aFilterMgr->CreateFreeFaces();
3006         break;
3007       case SMESH::FT_FreeNodes:
3008         aPredicate = aFilterMgr->CreateFreeNodes();
3009         break;
3010       case SMESH::FT_EqualNodes:
3011         {
3012           SMESH::EqualNodes_ptr pred = aFilterMgr->CreateEqualNodes();
3013           pred->SetTolerance( aTolerance );
3014           aPredicate = pred;
3015           break;
3016         }
3017       case SMESH::FT_EqualEdges:
3018         aPredicate = aFilterMgr->CreateEqualEdges();
3019         break;
3020       case SMESH::FT_EqualFaces:
3021         aPredicate = aFilterMgr->CreateEqualFaces();
3022         break;
3023       case SMESH::FT_EqualVolumes:
3024         aPredicate = aFilterMgr->CreateEqualVolumes();
3025         break;
3026       case SMESH::FT_BelongToMeshGroup:
3027         {
3028           SMESH::BelongToMeshGroup_ptr tmpPred = aFilterMgr->CreateBelongToMeshGroup();
3029           tmpPred->SetGroupID( aThresholdID );
3030           aPredicate = tmpPred;
3031         }
3032         break;
3033       case SMESH::FT_BelongToGeom:
3034         {
3035           SMESH::BelongToGeom_ptr tmpPred = aFilterMgr->CreateBelongToGeom();
3036           tmpPred->SetElementType( aTypeOfElem );
3037           tmpPred->SetShape( aThresholdID, aThresholdStr );
3038           tmpPred->SetTolerance( aTolerance );
3039           aPredicate = tmpPred;
3040         }
3041         break;
3042       case SMESH::FT_BelongToPlane:
3043       case SMESH::FT_BelongToCylinder:
3044       case SMESH::FT_BelongToGenSurface:
3045         {
3046           SMESH::BelongToSurface_ptr tmpPred;
3047           switch ( aCriterion ) {
3048           case SMESH::FT_BelongToPlane:
3049             tmpPred = aFilterMgr->CreateBelongToPlane(); break;
3050           case SMESH::FT_BelongToCylinder:
3051             tmpPred = aFilterMgr->CreateBelongToCylinder(); break;
3052           default:
3053             tmpPred = aFilterMgr->CreateBelongToGenSurface();
3054           }
3055           tmpPred->SetShape( aThresholdID, aThresholdStr, aTypeOfElem );
3056           tmpPred->SetTolerance( aTolerance );
3057           aPredicate = tmpPred;
3058         }
3059         break;
3060       case SMESH::FT_LyingOnGeom:
3061         {
3062           SMESH::LyingOnGeom_ptr tmpPred = aFilterMgr->CreateLyingOnGeom();
3063           tmpPred->SetElementType( aTypeOfElem );
3064           tmpPred->SetShape( aThresholdID, aThresholdStr );
3065           tmpPred->SetTolerance( aTolerance );
3066           aPredicate = tmpPred;
3067         }
3068         break;
3069       case SMESH::FT_RangeOfIds:
3070         {
3071           SMESH::RangeOfIds_ptr tmpPred = aFilterMgr->CreateRangeOfIds();
3072           tmpPred->SetRangeStr( aThresholdStr );
3073           tmpPred->SetElementType( aTypeOfElem );
3074           aPredicate = tmpPred;
3075         }
3076         break;
3077       case SMESH::FT_BadOrientedVolume:
3078         {
3079           aPredicate = aFilterMgr->CreateBadOrientedVolume();
3080         }
3081         break;
3082       case SMESH::FT_BareBorderVolume:
3083         {
3084           aPredicate = aFilterMgr->CreateBareBorderVolume();
3085         }
3086         break;
3087       case SMESH::FT_BareBorderFace:
3088         {
3089           aPredicate = aFilterMgr->CreateBareBorderFace();
3090         }
3091         break;
3092       case SMESH::FT_OverConstrainedVolume:
3093         {
3094           aPredicate = aFilterMgr->CreateOverConstrainedVolume();
3095         }
3096         break;
3097       case SMESH::FT_OverConstrainedFace:
3098         {
3099           aPredicate = aFilterMgr->CreateOverConstrainedFace();
3100         }
3101         break;
3102       case SMESH::FT_LinearOrQuadratic:
3103         {
3104           SMESH::LinearOrQuadratic_ptr tmpPred = aFilterMgr->CreateLinearOrQuadratic();
3105           tmpPred->SetElementType( aTypeOfElem );
3106           aPredicate = tmpPred;
3107           break;
3108         }
3109       case SMESH::FT_GroupColor:
3110         {
3111           SMESH::GroupColor_ptr tmpPred = aFilterMgr->CreateGroupColor();
3112           tmpPred->SetElementType( aTypeOfElem );
3113           tmpPred->SetColorStr( aThresholdStr );
3114           aPredicate = tmpPred;
3115           break;
3116         }
3117       case SMESH::FT_ElemGeomType:
3118         {
3119           SMESH::ElemGeomType_ptr tmpPred = aFilterMgr->CreateElemGeomType();
3120           tmpPred->SetElementType( aTypeOfElem );
3121           tmpPred->SetGeometryType( (GeometryType)(int)(aThreshold + 0.5) );
3122           aPredicate = tmpPred;
3123           break;
3124         }
3125       case SMESH::FT_EntityType:
3126         {
3127           SMESH::ElemEntityType_ptr tmpPred = aFilterMgr->CreateElemEntityType();
3128           tmpPred->SetElementType( aTypeOfElem );
3129           tmpPred->SetEntityType( EntityType( (int (aThreshold + 0.5))));
3130           aPredicate = tmpPred;
3131           break;
3132         }
3133       case SMESH::FT_CoplanarFaces:
3134         {
3135           SMESH::CoplanarFaces_ptr tmpPred = aFilterMgr->CreateCoplanarFaces();
3136           tmpPred->SetFace( atol (aThresholdID ));
3137           tmpPred->SetTolerance( aTolerance );
3138           aPredicate = tmpPred;
3139           break;
3140         }
3141       case SMESH::FT_ConnectedElements:
3142         {
3143           SMESH::ConnectedElements_ptr tmpPred = aFilterMgr->CreateConnectedElements();
3144           if ( strlen( aThresholdID ) > 0 ) // shape ID
3145             tmpPred->SetThreshold( aThresholdID, SMESH::ConnectedElements::VERTEX );
3146           else if ( strlen( aThresholdStr ) > 0 ) // point coords
3147             tmpPred->SetThreshold( aThresholdStr, SMESH::ConnectedElements::POINT );
3148           else if ( aThreshold >= 1 )
3149             tmpPred->SetNode( (CORBA::Long) aThreshold ); // node ID
3150           tmpPred->SetElementType( aTypeOfElem );
3151           aPredicate = tmpPred;
3152           break;
3153         }
3154
3155       default:
3156         continue;
3157     }
3158
3159     // Comparator
3160     if ( !aFunctor->_is_nil() && aPredicate->_is_nil() )
3161     {
3162       SMESH::Comparator_ptr aComparator = SMESH::Comparator::_nil();
3163
3164       if ( aCompare == SMESH::FT_LessThan )
3165         aComparator = aFilterMgr->CreateLessThan();
3166       else if ( aCompare == SMESH::FT_MoreThan )
3167         aComparator = aFilterMgr->CreateMoreThan();
3168       else if ( aCompare == SMESH::FT_EqualTo )
3169         aComparator = aFilterMgr->CreateEqualTo();
3170       else
3171         continue;
3172
3173       aComparator->SetNumFunctor( aFunctor );
3174       aComparator->SetMargin( aThreshold );
3175
3176       if ( aCompare == FT_EqualTo )
3177       {
3178         SMESH::EqualTo_var anEqualTo = SMESH::EqualTo::_narrow( aComparator );
3179         anEqualTo->SetTolerance( aTolerance );
3180       }
3181
3182       aPredicate = aComparator;
3183
3184       aFunctor->SetPrecision( aPrecision );
3185     }
3186
3187     // Logical not
3188     if ( aUnary == FT_LogicalNOT )
3189     {
3190       SMESH::LogicalNOT_ptr aNotPred = aFilterMgr->CreateLogicalNOT();
3191       aNotPred->SetPredicate( aPredicate );
3192       aPredicate = aNotPred;
3193     }
3194
3195     // logical op
3196     aPredicates.push_back( aPredicate );
3197     aBinaries.push_back( aBinary );
3198     pd <<"aCriteria.append(aCriterion)";
3199
3200   } // end of for
3201   TPythonDump pd; pd<<this<<".SetCriteria(aCriteria)";
3202
3203   // CREATE ONE PREDICATE FROM PREVIOUSLY CREATED MAP
3204
3205   // combine all "AND" operations
3206
3207   std::list<SMESH::Predicate_ptr> aResList;
3208
3209   std::list<SMESH::Predicate_ptr>::iterator aPredIter;
3210   std::list<int>::iterator                  aBinaryIter;
3211
3212   SMESH::Predicate_ptr aPrevPredicate = SMESH::Predicate::_nil();
3213   int aPrevBinary = SMESH::FT_Undefined;
3214   if ( !aBinaries.empty() )
3215     aBinaries.back() = SMESH::FT_Undefined;
3216
3217   for ( aPredIter = aPredicates.begin(), aBinaryIter = aBinaries.begin();
3218         aPredIter != aPredicates.end() && aBinaryIter != aBinaries.end();
3219         ++aPredIter, ++aBinaryIter )
3220   {
3221     int aCurrBinary = *aBinaryIter;
3222
3223     SMESH::Predicate_ptr aCurrPred = SMESH::Predicate::_nil();
3224
3225     if ( aPrevBinary == SMESH::FT_LogicalAND )
3226     {
3227
3228       SMESH::LogicalBinary_ptr aBinaryPred = aFilterMgr->CreateLogicalAND();
3229       aBinaryPred->SetPredicate1( aPrevPredicate );
3230       aBinaryPred->SetPredicate2( *aPredIter );
3231       aCurrPred = aBinaryPred;
3232     }
3233     else
3234       aCurrPred = *aPredIter;
3235
3236     if ( aCurrBinary != SMESH::FT_LogicalAND )
3237       aResList.push_back( aCurrPred );
3238
3239     aPrevPredicate = aCurrPred;
3240     aPrevBinary = aCurrBinary;
3241   }
3242
3243   // combine all "OR" operations
3244
3245   SMESH::Predicate_ptr aResPredicate = SMESH::Predicate::_nil();
3246
3247   if ( aResList.size() == 1 )
3248     aResPredicate = *aResList.begin();
3249   else if ( aResList.size() > 1 )
3250   {
3251     std::list<SMESH::Predicate_ptr>::iterator anIter = aResList.begin();
3252     aResPredicate = *anIter;
3253     anIter++;
3254     for ( ; anIter != aResList.end(); ++anIter )
3255     {
3256       SMESH::LogicalBinary_ptr aBinaryPred = aFilterMgr->CreateLogicalOR();
3257       aBinaryPred->SetPredicate1( aResPredicate );
3258       aBinaryPred->SetPredicate2( *anIter );
3259       aResPredicate = aBinaryPred;
3260     }
3261   }
3262
3263   SetPredicate( aResPredicate );
3264   if ( !aResPredicate->_is_nil() )
3265     aResPredicate->UnRegister();
3266
3267   return !aResPredicate->_is_nil();
3268 }
3269
3270 //=======================================================================
3271 // name    : Filter_i::GetPredicate_i
3272 // Purpose : Get implementation of predicate
3273 //=======================================================================
3274 Predicate_i* Filter_i::GetPredicate_i()
3275 {
3276   return myPredicate;
3277 }
3278
3279 //=======================================================================
3280 // name    : Filter_i::GetPredicate
3281 // Purpose : Get predicate
3282 //=======================================================================
3283 Predicate_ptr Filter_i::GetPredicate()
3284 {
3285   if ( myPredicate == 0 )
3286     return SMESH::Predicate::_nil();
3287   else
3288   {
3289     SMESH::Predicate_var anObj = myPredicate->_this();
3290     // if ( SMESH::Functor_i* fun = SMESH::DownCast<SMESH::Functor_i*>( anObj ))
3291     //   TPythonDump() << fun << " = " << this << ".GetPredicate()";
3292     return anObj._retn();
3293   }
3294 }
3295
3296 //================================================================================
3297 /*!
3298  * \brief Find groups it depends on
3299  */
3300 //================================================================================
3301
3302 void Filter_i::FindBaseObjects()
3303 {
3304   // release current groups
3305   for ( size_t i = 0; i < myBaseGroups.size(); ++i )
3306     if ( myBaseGroups[i] )
3307     {
3308       myBaseGroups[i]->RemoveModifWaiter( this );
3309       myBaseGroups[i]->UnRegister();
3310     }
3311
3312   // remember new groups
3313   myBaseGroups.clear();
3314   if ( myPredicate )
3315   {
3316     std::vector<Predicate_i*> predicates;
3317     getPrediacates( myPredicate, predicates );
3318     for ( size_t i = 0; i < predicates.size(); ++i )
3319       if ( BelongToMeshGroup_i* bmg = dynamic_cast< BelongToMeshGroup_i* >( predicates[i] ))
3320       {
3321         SMESH::SMESH_GroupBase_var g = bmg->GetGroup();
3322         SMESH_GroupBase_i* g_i = SMESH::DownCast< SMESH_GroupBase_i*>( g );
3323         if ( g_i )
3324         {
3325           g_i->AddModifWaiter( this );
3326           g_i->Register();
3327           myBaseGroups.push_back( g_i );
3328         }
3329       }
3330   }
3331 }
3332
3333 //================================================================================
3334 /*!
3335  * \brief When notified on removal of myBaseGroups[i], remove a reference to a
3336  *        group from a predicate
3337  */
3338 //================================================================================
3339
3340 void Filter_i::OnBaseObjModified(NotifyerAndWaiter* group, bool removed)
3341 {
3342   if ( !removed )
3343     return; // a GroupOnFilter holding this filter is notified automatically
3344
3345   if ( myPredicate )
3346   {
3347     std::vector<Predicate_i*> predicates;
3348     getPrediacates( myPredicate, predicates );
3349     for ( size_t i = 0; i < predicates.size(); ++i )
3350       if ( BelongToMeshGroup_i* bmg = dynamic_cast< BelongToMeshGroup_i* >( predicates[i] ))
3351       {
3352         SMESH::SMESH_GroupBase_var g = bmg->GetGroup();
3353         SMESH_GroupBase_i* g_i = SMESH::DownCast< SMESH_GroupBase_i*>( g );
3354         if ( g_i == group )
3355         {
3356           bmg->SetGroup( SMESH::SMESH_GroupBase::_nil() );
3357           bmg->SetGroupID( "" );
3358         }
3359       }
3360   }
3361
3362   FindBaseObjects(); // release and update myBaseGroups;
3363 }
3364
3365 /*
3366                             FILTER LIBRARY
3367 */
3368
3369 #define ATTR_TYPE          "type"
3370 #define ATTR_COMPARE       "compare"
3371 #define ATTR_THRESHOLD     "threshold"
3372 #define ATTR_UNARY         "unary"
3373 #define ATTR_BINARY        "binary"
3374 #define ATTR_THRESHOLD_STR "threshold_str"
3375 #define ATTR_TOLERANCE     "tolerance"
3376 #define ATTR_ELEMENT_TYPE  "ElementType"
3377
3378 //=======================================================================
3379 // name    : toString
3380 // Purpose : Convert bool to LDOMString
3381 //=======================================================================
3382 static inline LDOMString toString( CORBA::Boolean val )
3383 {
3384   return val ? "logical not" : "";
3385 }
3386
3387 //=======================================================================
3388 // name    : toBool
3389 // Purpose : Convert LDOMString to bool
3390 //=======================================================================
3391 static inline bool toBool( const LDOMString& theStr )
3392 {
3393   return theStr.equals( "logical not" );
3394 }
3395
3396 //=======================================================================
3397 // name    : toString
3398 // Purpose : Convert double to LDOMString
3399 //=======================================================================
3400 static inline LDOMString toString( CORBA::Double val )
3401 {
3402   char a[ 255 ];
3403   sprintf( a, "%e", val );
3404   return LDOMString( a );
3405 }
3406
3407 //=======================================================================
3408 // name    : toDouble
3409 // Purpose : Convert LDOMString to double
3410 //=======================================================================
3411 static inline double toDouble( const LDOMString& theStr )
3412 {
3413   return atof( theStr.GetString() );
3414 }
3415
3416 //=======================================================================
3417 // name    : toString
3418 // Purpose : Convert functor type to LDOMString
3419 //=======================================================================
3420 static inline LDOMString toString( CORBA::Long theType )
3421 {
3422   switch ( theType )
3423   {
3424     case FT_AspectRatio           : return "Aspect ratio";
3425     case FT_Warping               : return "Warping";
3426     case FT_MinimumAngle          : return "Minimum angle";
3427     case FT_Taper                 : return "Taper";
3428     case FT_Skew                  : return "Skew";
3429     case FT_Area                  : return "Area";
3430     case FT_Volume3D              : return "Volume3D";
3431     case FT_MaxElementLength2D    : return "Max element length 2D";
3432     case FT_MaxElementLength3D    : return "Max element length 3D";
3433     case FT_BelongToMeshGroup     : return "Belong to Mesh Group";
3434     case FT_BelongToGeom          : return "Belong to Geom";
3435     case FT_BelongToPlane         : return "Belong to Plane";
3436     case FT_BelongToCylinder      : return "Belong to Cylinder";
3437     case FT_BelongToGenSurface    : return "Belong to Generic Surface";
3438     case FT_LyingOnGeom           : return "Lying on Geom";
3439     case FT_BadOrientedVolume     : return "Bad Oriented Volume";
3440     case FT_BareBorderVolume      : return "Volumes with bare border";
3441     case FT_BareBorderFace        : return "Faces with bare border";
3442     case FT_OverConstrainedVolume : return "Over-constrained Volumes";
3443     case FT_OverConstrainedFace   : return "Over-constrained Faces";
3444     case FT_RangeOfIds            : return "Range of IDs";
3445     case FT_FreeBorders           : return "Free borders";
3446     case FT_FreeEdges             : return "Free edges";
3447     case FT_FreeFaces             : return "Free faces";
3448     case FT_FreeNodes             : return "Free nodes";
3449     case FT_EqualNodes            : return "Equal nodes";
3450     case FT_EqualEdges            : return "Equal edges";
3451     case FT_EqualFaces            : return "Equal faces";
3452     case FT_EqualVolumes          : return "Equal volumes";
3453     case FT_MultiConnection       : return "Borders at multi-connections";
3454     case FT_MultiConnection2D     :return "Borders at multi-connections 2D";
3455     case FT_Length                : return "Length";
3456     case FT_Length2D              : return "Length 2D";
3457     case FT_LessThan              : return "Less than";
3458     case FT_MoreThan              : return "More than";
3459     case FT_EqualTo               : return "Equal to";
3460     case FT_LogicalNOT            : return "Not";
3461     case FT_LogicalAND            : return "And";
3462     case FT_LogicalOR             : return "Or";
3463     case FT_GroupColor            : return "Color of Group";
3464     case FT_LinearOrQuadratic     : return "Linear or Quadratic";
3465     case FT_ElemGeomType          : return "Element geomtry type";
3466     case FT_EntityType            : return "Entity type";
3467     case FT_Undefined             : return "";
3468     default                       : return "";
3469   }
3470 }
3471
3472 //=======================================================================
3473 // name    : toFunctorType
3474 // Purpose : Convert LDOMString to functor type
3475 //=======================================================================
3476 static inline SMESH::FunctorType toFunctorType( const LDOMString& theStr )
3477 {
3478   if      ( theStr.equals( "Aspect ratio"                 ) ) return FT_AspectRatio;
3479   else if ( theStr.equals( "Warping"                      ) ) return FT_Warping;
3480   else if ( theStr.equals( "Minimum angle"                ) ) return FT_MinimumAngle;
3481   else if ( theStr.equals( "Taper"                        ) ) return FT_Taper;
3482   else if ( theStr.equals( "Skew"                         ) ) return FT_Skew;
3483   else if ( theStr.equals( "Area"                         ) ) return FT_Area;
3484   else if ( theStr.equals( "Volume3D"                     ) ) return FT_Volume3D;
3485   else if ( theStr.equals( "Max element length 2D"        ) ) return FT_MaxElementLength2D;
3486   else if ( theStr.equals( "Max element length 3D"        ) ) return FT_MaxElementLength3D;
3487   else if ( theStr.equals( "Belong to Mesh Group"         ) ) return FT_BelongToMeshGroup;
3488   else if ( theStr.equals( "Belong to Geom"               ) ) return FT_BelongToGeom;
3489   else if ( theStr.equals( "Belong to Plane"              ) ) return FT_BelongToPlane;
3490   else if ( theStr.equals( "Belong to Cylinder"           ) ) return FT_BelongToCylinder;
3491   else if ( theStr.equals( "Belong to Generic Surface"    ) ) return FT_BelongToGenSurface;
3492   else if ( theStr.equals( "Lying on Geom"                ) ) return FT_LyingOnGeom;
3493   else if ( theStr.equals( "Free borders"                 ) ) return FT_FreeBorders;
3494   else if ( theStr.equals( "Free edges"                   ) ) return FT_FreeEdges;
3495   else if ( theStr.equals( "Free faces"                   ) ) return FT_FreeFaces;
3496   else if ( theStr.equals( "Free nodes"                   ) ) return FT_FreeNodes;
3497   else if ( theStr.equals( "Equal nodes"                  ) ) return FT_EqualNodes;
3498   else if ( theStr.equals( "Equal edges"                  ) ) return FT_EqualEdges;
3499   else if ( theStr.equals( "Equal faces"                  ) ) return FT_EqualFaces;
3500   else if ( theStr.equals( "Equal volumes"                ) ) return FT_EqualVolumes;
3501   else if ( theStr.equals( "Borders at multi-connections" ) ) return FT_MultiConnection;
3502   //  else if ( theStr.equals( "Borders at multi-connections 2D" ) ) return FT_MultiConnection2D;
3503   else if ( theStr.equals( "Length"                       ) ) return FT_Length;
3504   //  else if ( theStr.equals( "Length2D"                     ) ) return FT_Length2D;
3505   else if ( theStr.equals( "Range of IDs"                 ) ) return FT_RangeOfIds;
3506   else if ( theStr.equals( "Bad Oriented Volume"          ) ) return FT_BadOrientedVolume;
3507   else if ( theStr.equals( "Volumes with bare border"     ) ) return FT_BareBorderVolume;
3508   else if ( theStr.equals( "Faces with bare border"       ) ) return FT_BareBorderFace;
3509   else if ( theStr.equals( "Over-constrained Volumes"     ) ) return FT_OverConstrainedVolume;
3510   else if ( theStr.equals( "Over-constrained Faces"       ) ) return FT_OverConstrainedFace;
3511   else if ( theStr.equals( "Less than"                    ) ) return FT_LessThan;
3512   else if ( theStr.equals( "More than"                    ) ) return FT_MoreThan;
3513   else if ( theStr.equals( "Equal to"                     ) ) return FT_EqualTo;
3514   else if ( theStr.equals( "Not"                          ) ) return FT_LogicalNOT;
3515   else if ( theStr.equals( "And"                          ) ) return FT_LogicalAND;
3516   else if ( theStr.equals( "Or"                           ) ) return FT_LogicalOR;
3517   else if ( theStr.equals( "Color of Group"               ) ) return FT_GroupColor;
3518   else if ( theStr.equals( "Linear or Quadratic"          ) ) return FT_LinearOrQuadratic;
3519   else if ( theStr.equals( "Element geomtry type"         ) ) return FT_ElemGeomType;
3520   else if ( theStr.equals( "Entity type"                  ) ) return FT_EntityType;
3521   else if ( theStr.equals( ""                             ) ) return FT_Undefined;
3522   else  return FT_Undefined;
3523 }
3524
3525 //=======================================================================
3526 // name    : toFunctorType
3527 // Purpose : Convert LDOMString to value of ElementType enumeration
3528 //=======================================================================
3529 static inline SMESH::ElementType toElementType( const LDOMString& theStr )
3530 {
3531   if      ( theStr.equals( "NODE"   ) ) return SMESH::NODE;
3532   else if ( theStr.equals( "EDGE"   ) ) return SMESH::EDGE;
3533   else if ( theStr.equals( "FACE"   ) ) return SMESH::FACE;
3534   else if ( theStr.equals( "VOLUME" ) ) return SMESH::VOLUME;
3535   else                                  return SMESH::ALL;
3536 }
3537
3538 //=======================================================================
3539 // name    : toString
3540 // Purpose : Convert ElementType to string
3541 //=======================================================================
3542 static inline LDOMString toString( const SMESH::ElementType theType )
3543 {
3544   switch ( theType )
3545   {
3546     case SMESH::NODE   : return "NODE";
3547     case SMESH::EDGE   : return "EDGE";
3548     case SMESH::FACE   : return "FACE";
3549     case SMESH::VOLUME : return "VOLUME";
3550     case SMESH::ALL    : return "ALL";
3551     default            : return "";
3552   }
3553 }
3554
3555 //=======================================================================
3556 // name    : findFilter
3557 // Purpose : Find filter in document
3558 //=======================================================================
3559 static LDOM_Element findFilter( const char* theFilterName,
3560                                 const LDOM_Document& theDoc,
3561                                 LDOM_Node* theParent = 0 )
3562 {
3563   LDOM_Element aRootElement = theDoc.getDocumentElement();
3564   if ( aRootElement.isNull() || !aRootElement.hasChildNodes() )
3565     return LDOM_Element();
3566
3567   for ( LDOM_Node aTypeNode = aRootElement.getFirstChild();
3568         !aTypeNode.isNull(); aTypeNode = aTypeNode.getNextSibling() )
3569   {
3570     for ( LDOM_Node aFilter = aTypeNode.getFirstChild();
3571           !aFilter.isNull(); aFilter = aFilter.getNextSibling() )
3572     {
3573       LDOM_Element* anElem = ( LDOM_Element* )&aFilter;
3574       if ( anElem->getTagName().equals( LDOMString( "filter" ) ) &&
3575            anElem->getAttribute( "name" ).equals( LDOMString( theFilterName ) ) )
3576       {
3577         if ( theParent != 0  )
3578           *theParent = aTypeNode;
3579         return (LDOM_Element&)aFilter;
3580       }
3581     }
3582   }
3583   return LDOM_Element();
3584 }
3585
3586 //=======================================================================
3587 // name    : getSectionName
3588 // Purpose : Get name of section of filters
3589 //=======================================================================
3590 static const char* getSectionName( const ElementType theType )
3591 {
3592   switch ( theType )
3593   {
3594     case SMESH::NODE   : return "Filters for nodes";
3595     case SMESH::EDGE   : return "Filters for edges";
3596     case SMESH::FACE   : return "Filters for faces";
3597     case SMESH::VOLUME : return "Filters for volumes";
3598     case SMESH::ALL    : return "Filters for elements";
3599     default            : return "";
3600   }
3601 }
3602
3603 //=======================================================================
3604 // name    : getSection
3605 // Purpose : Create section for filters corresponding to the entity type
3606 //=======================================================================
3607 static LDOM_Node getSection( const ElementType theType,
3608                              LDOM_Document&    theDoc,
3609                              const bool        toCreate = false )
3610 {
3611   LDOM_Element aRootElement = theDoc.getDocumentElement();
3612   if ( aRootElement.isNull() )
3613     return LDOM_Node();
3614
3615   // Find section
3616   bool anExist = false;
3617   const char* aSectionName = getSectionName( theType );
3618   if ( strcmp( aSectionName, "" ) == 0 )
3619     return LDOM_Node();
3620
3621   LDOM_NodeList aSections = theDoc.getElementsByTagName( "section" );
3622   LDOM_Node aNode;
3623   for ( int i = 0, n = aSections.getLength(); i < n; i++ )
3624   {
3625     aNode = aSections.item( i );
3626     LDOM_Element& anItem = ( LDOM_Element& )aNode;
3627     if ( anItem.getAttribute( "name" ).equals( LDOMString( aSectionName ) ) )
3628     {
3629       anExist = true;
3630       break;
3631     }
3632   }
3633
3634   // Create new section if necessary
3635   if ( !anExist )
3636   {
3637     if ( toCreate )
3638     {
3639       LDOM_Element aNewItem = theDoc.createElement( "section" );
3640       aNewItem.setAttribute( "name", aSectionName );
3641       aRootElement.appendChild( aNewItem );
3642       return aNewItem;
3643     }
3644     else
3645       return LDOM_Node();
3646   }
3647   return
3648     aNode;
3649 }
3650
3651 //=======================================================================
3652 // name    : createFilterItem
3653 // Purpose : Create filter item or LDOM document
3654 //=======================================================================
3655 static LDOM_Element createFilterItem( const char*       theName,
3656                                       SMESH::Filter_ptr theFilter,
3657                                       LDOM_Document&    theDoc )
3658 {
3659   // create new filter in document
3660   LDOM_Element aFilterItem = theDoc.createElement( "filter" );
3661   aFilterItem.setAttribute( "name", theName );
3662
3663   // save filter criterions
3664   SMESH::Filter::Criteria_var aCriteria = new SMESH::Filter::Criteria;
3665
3666   if ( !theFilter->GetCriteria( aCriteria ) )
3667     return LDOM_Element();
3668
3669   for ( CORBA::ULong i = 0, n = aCriteria->length(); i < n; i++ )
3670   {
3671     LDOM_Element aCriterionItem = theDoc.createElement( "criterion" );
3672
3673     aCriterionItem.setAttribute( ATTR_TYPE         , toString(  aCriteria[ i ].Type) );
3674     aCriterionItem.setAttribute( ATTR_COMPARE      , toString(  aCriteria[ i ].Compare ) );
3675     aCriterionItem.setAttribute( ATTR_THRESHOLD    , toString(  aCriteria[ i ].Threshold ) );
3676     aCriterionItem.setAttribute( ATTR_UNARY        , toString(  aCriteria[ i ].UnaryOp ) );
3677     aCriterionItem.setAttribute( ATTR_BINARY       , toString(  aCriteria[ i ].BinaryOp ) );
3678
3679     aCriterionItem.setAttribute( ATTR_THRESHOLD_STR, (const char*)aCriteria[ i ].ThresholdStr );
3680     aCriterionItem.setAttribute( ATTR_TOLERANCE    , toString( aCriteria[ i ].Tolerance ) );
3681     aCriterionItem.setAttribute( ATTR_ELEMENT_TYPE ,
3682       toString( (SMESH::ElementType)aCriteria[ i ].TypeOfElement ) );
3683
3684     aFilterItem.appendChild( aCriterionItem );
3685   }
3686
3687   return aFilterItem;
3688 }
3689
3690 //=======================================================================
3691 // name    : FilterLibrary_i::FilterLibrary_i
3692 // Purpose : Constructor
3693 //=======================================================================
3694 FilterLibrary_i::FilterLibrary_i( const char* theFileName )
3695 {
3696   myFileName = CORBA::string_dup( theFileName );
3697   SMESH::FilterManager_i* aFilterMgr = new SMESH::FilterManager_i();
3698   myFilterMgr = aFilterMgr->_this();
3699
3700   LDOMParser aParser;
3701
3702   // Try to use existing library file
3703   bool anExists = false;
3704   if ( !aParser.parse( myFileName ) )
3705   {
3706     myDoc = aParser.getDocument();
3707     anExists = true;
3708   }
3709   // Create a new XML document if it doesn't exist
3710   else
3711     myDoc = LDOM_Document::createDocument( LDOMString() );
3712
3713   LDOM_Element aRootElement = myDoc.getDocumentElement();
3714   if ( aRootElement.isNull() )
3715   {
3716     // If the existing document is empty --> try to create a new one
3717     if ( anExists )
3718       myDoc = LDOM_Document::createDocument( LDOMString() );
3719   }
3720 }
3721
3722 //=======================================================================
3723 // name    : FilterLibrary_i::FilterLibrary_i
3724 // Purpose : Constructor
3725 //=======================================================================
3726 FilterLibrary_i::FilterLibrary_i()
3727 {
3728   myFileName = 0;
3729   SMESH::FilterManager_i* aFilter = new SMESH::FilterManager_i();
3730   myFilterMgr = aFilter->_this();
3731
3732   myDoc = LDOM_Document::createDocument( LDOMString() );
3733 }
3734
3735 FilterLibrary_i::~FilterLibrary_i()
3736 {
3737   CORBA::string_free( myFileName );
3738   //TPythonDump()<<this<<".UnRegister()";
3739 }
3740
3741 //=======================================================================
3742 // name    : FilterLibrary_i::Copy
3743 // Purpose : Create filter and initialize it with values from library
3744 //=======================================================================
3745 Filter_ptr FilterLibrary_i::Copy( const char* theFilterName )
3746 {
3747   Filter_ptr aRes = Filter::_nil();
3748   LDOM_Node aFilter = findFilter( theFilterName, myDoc );
3749
3750   if ( aFilter.isNull() )
3751     return aRes;
3752
3753   std::list<SMESH::Filter::Criterion> aCriteria;
3754
3755   for ( LDOM_Node aCritNode = aFilter.getFirstChild();
3756         !aCritNode.isNull() ; aCritNode = aCritNode.getNextSibling() )
3757   {
3758     LDOM_Element* aCrit = (LDOM_Element*)&aCritNode;
3759
3760     const char* aTypeStr      = aCrit->getAttribute( ATTR_TYPE          ).GetString();
3761     const char* aCompareStr   = aCrit->getAttribute( ATTR_COMPARE       ).GetString();
3762     const char* aUnaryStr     = aCrit->getAttribute( ATTR_UNARY         ).GetString();
3763     const char* aBinaryStr    = aCrit->getAttribute( ATTR_BINARY        ).GetString();
3764     const char* anElemTypeStr = aCrit->getAttribute( ATTR_ELEMENT_TYPE  ).GetString();
3765
3766     SMESH::Filter::Criterion aCriterion = createCriterion();
3767
3768     aCriterion.Type          = toFunctorType( aTypeStr );
3769     aCriterion.Compare       = toFunctorType( aCompareStr );
3770     aCriterion.UnaryOp       = toFunctorType( aUnaryStr );
3771     aCriterion.BinaryOp      = toFunctorType( aBinaryStr );
3772
3773     aCriterion.TypeOfElement = toElementType( anElemTypeStr );
3774
3775     LDOMString str = aCrit->getAttribute( ATTR_THRESHOLD );
3776     int val = 0;
3777     aCriterion.Threshold = str.Type() == LDOMBasicString::LDOM_Integer && str.GetInteger( val )
3778       ? val : atof( str.GetString() );
3779
3780     str = aCrit->getAttribute( ATTR_TOLERANCE );
3781     aCriterion.Tolerance = str.Type() == LDOMBasicString::LDOM_Integer && str.GetInteger( val )
3782       ? val : atof( str.GetString() );
3783
3784     str = aCrit->getAttribute( ATTR_THRESHOLD_STR );
3785     if ( str.Type() == LDOMBasicString::LDOM_Integer && str.GetInteger( val ) )
3786     {
3787       char a[ 255 ];
3788       sprintf( a, "%d", val );
3789       aCriterion.ThresholdStr = CORBA::string_dup( a );
3790     }
3791     else
3792       aCriterion.ThresholdStr = str.GetString();
3793
3794     aCriteria.push_back( aCriterion );
3795   }
3796
3797   SMESH::Filter::Criteria_var aCriteriaVar = new SMESH::Filter::Criteria;
3798   aCriteriaVar->length( aCriteria.size() );
3799
3800   CORBA::ULong i = 0;
3801   std::list<SMESH::Filter::Criterion>::iterator anIter = aCriteria.begin();
3802
3803   for( ; anIter != aCriteria.end(); ++anIter )
3804     aCriteriaVar[ i++ ] = *anIter;
3805
3806   aRes = myFilterMgr->CreateFilter();
3807   aRes->SetCriteria( aCriteriaVar.inout() );
3808
3809   TPythonDump()<<this<<".Copy('"<<theFilterName<<"')";
3810
3811   return aRes;
3812 }
3813
3814 //=======================================================================
3815 // name    : FilterLibrary_i::SetFileName
3816 // Purpose : Set file name for library
3817 //=======================================================================
3818 void FilterLibrary_i::SetFileName( const char* theFileName )
3819 {
3820   CORBA::string_free( myFileName );
3821   myFileName = CORBA::string_dup( theFileName );
3822   TPythonDump()<<this<<".SetFileName('"<<theFileName<<"')";
3823 }
3824
3825 //=======================================================================
3826 // name    : FilterLibrary_i::GetFileName
3827 // Purpose : Get file name of library
3828 //=======================================================================
3829 char* FilterLibrary_i::GetFileName()
3830 {
3831   return CORBA::string_dup( myFileName );
3832 }
3833
3834 //=======================================================================
3835 // name    : FilterLibrary_i::Add
3836 // Purpose : Add new filter to library
3837 //=======================================================================
3838 CORBA::Boolean FilterLibrary_i::Add( const char* theFilterName, Filter_ptr theFilter )
3839 {
3840   // if filter already in library or entry filter is null do nothing
3841   LDOM_Node aFilterNode = findFilter( theFilterName, myDoc );
3842   if ( !aFilterNode.isNull() || theFilter->_is_nil() )
3843     return false;
3844
3845   // get section corresponding to the filter type
3846   ElementType anEntType = theFilter->GetElementType();
3847
3848   LDOM_Node aSection = getSection( anEntType, myDoc, true );
3849   if ( aSection.isNull() )
3850     return false;
3851
3852   // create filter item
3853   LDOM_Element aFilterItem = createFilterItem( theFilterName, theFilter, myDoc );
3854   if ( aFilterItem.isNull() )
3855     return false;
3856   else
3857   {
3858     aSection.appendChild( aFilterItem );
3859     if(Filter_i* aFilter = DownCast<Filter_i*>(theFilter))
3860       TPythonDump()<<this<<".Add('"<<theFilterName<<"',"<<aFilter<<")";
3861     return true;
3862   }
3863 }
3864
3865 //=======================================================================
3866 // name    : FilterLibrary_i::Add
3867 // Purpose : Add new filter to library
3868 //=======================================================================
3869 CORBA::Boolean FilterLibrary_i::AddEmpty( const char* theFilterName, ElementType theType )
3870 {
3871   // if filter already in library or entry filter is null do nothing
3872   LDOM_Node aFilterNode = findFilter( theFilterName, myDoc );
3873   if ( !aFilterNode.isNull() )
3874     return false;
3875
3876   LDOM_Node aSection = getSection( theType, myDoc, true );
3877   if ( aSection.isNull() )
3878     return false;
3879
3880   // create filter item
3881   Filter_var aFilter = myFilterMgr->CreateFilter();
3882
3883   LDOM_Element aFilterItem = createFilterItem( theFilterName, aFilter, myDoc );
3884   if ( aFilterItem.isNull() )
3885     return false;
3886   else
3887   {
3888     aSection.appendChild( aFilterItem );
3889     TPythonDump()<<this<<".AddEmpty('"<<theFilterName<<"',"<<theType<<")";
3890     return true;
3891   }
3892 }
3893
3894 //=======================================================================
3895 // name    : FilterLibrary_i::Delete
3896 // Purpose : Delete filter from library
3897 //=======================================================================
3898 CORBA::Boolean FilterLibrary_i::Delete ( const char* theFilterName )
3899 {
3900   LDOM_Node aParentNode;
3901   LDOM_Node aFilterNode = findFilter( theFilterName, myDoc, &aParentNode );
3902   if ( aFilterNode.isNull() || aParentNode.isNull() )
3903     return false;
3904
3905   aParentNode.removeChild( aFilterNode );
3906   TPythonDump()<<this<<".Delete('"<<theFilterName<<"')";
3907   return true;
3908 }
3909
3910 //=======================================================================
3911 // name      : FilterLibrary_i::Replace
3912 // Purpose   : Replace existing filter with entry filter.
3913 // IMPORTANT : If filter does not exist it is not created
3914 //=======================================================================
3915 CORBA::Boolean FilterLibrary_i::Replace( const char* theFilterName,
3916                                          const char* theNewName,
3917                                          Filter_ptr  theFilter )
3918 {
3919   LDOM_Element aFilterItem = findFilter( theFilterName, myDoc );
3920   if ( aFilterItem.isNull() || theFilter->_is_nil() )
3921     return false;
3922
3923   LDOM_Element aNewItem = createFilterItem( theNewName, theFilter, myDoc );
3924   if ( aNewItem.isNull() )
3925     return false;
3926   else
3927   {
3928     aFilterItem.ReplaceElement( aNewItem );
3929     if(Filter_i* aFilter = DownCast<Filter_i*>(theFilter))
3930       TPythonDump()<<this<<".Replace('"<<theFilterName<<"','"<<theNewName<<"',"<<aFilter<<")";
3931     return true;
3932   }
3933 }
3934
3935 //=======================================================================
3936 // name    : FilterLibrary_i::Save
3937 // Purpose : Save library on disk
3938 //=======================================================================
3939 CORBA::Boolean FilterLibrary_i::Save()
3940 {
3941   if ( myFileName == 0 || strlen( myFileName ) == 0 )
3942     return false;
3943
3944 #if OCC_VERSION_MAJOR < 7
3945   FILE* aOutFile = fopen( myFileName, "wt" );
3946   if ( !aOutFile )
3947     return false;
3948
3949   LDOM_XmlWriter aWriter( aOutFile );
3950   aWriter.SetIndentation( 2 );
3951   aWriter << myDoc;
3952   fclose( aOutFile );
3953 #else
3954   std::filebuf fb;
3955   fb.open( myFileName, std::ios::out );
3956
3957   Standard_OStream os( &fb );
3958
3959   LDOM_XmlWriter aWriter;
3960   aWriter.SetIndentation( 2 );
3961   aWriter.Write( os, myDoc );
3962   fb.close();
3963 #endif
3964
3965   TPythonDump()<<this<<".Save()";
3966   return true;
3967 }
3968
3969 //=======================================================================
3970 // name    : FilterLibrary_i::SaveAs
3971 // Purpose : Save library on disk
3972 //=======================================================================
3973 CORBA::Boolean FilterLibrary_i::SaveAs( const char* aFileName )
3974 {
3975   myFileName = strdup ( aFileName );
3976   TPythonDump()<<this<<".SaveAs('"<<aFileName<<"')";
3977   return Save();
3978 }
3979
3980 //=======================================================================
3981 // name    : FilterLibrary_i::IsPresent
3982 // Purpose : Verify whether filter is in library
3983 //=======================================================================
3984 CORBA::Boolean FilterLibrary_i::IsPresent( const char* theFilterName )
3985 {
3986   return !findFilter( theFilterName, myDoc ).isNull();
3987 }
3988
3989 //=======================================================================
3990 // name    : FilterLibrary_i::NbFilters
3991 // Purpose : Return amount of filters in library
3992 //=======================================================================
3993 CORBA::Long FilterLibrary_i::NbFilters( ElementType theType )
3994 {
3995   string_array_var aNames = GetNames( theType );
3996   return aNames->length();
3997 }
3998
3999 //=======================================================================
4000 // name    : FilterLibrary_i::GetNames
4001 // Purpose : Get names of filters from library
4002 //=======================================================================
4003 string_array* FilterLibrary_i::GetNames( ElementType theType )
4004 {
4005   string_array_var anArray = new string_array;
4006   TColStd_SequenceOfHAsciiString aSeq;
4007
4008   LDOM_Node aSection = getSection( theType, myDoc, false );
4009
4010   if ( !aSection.isNull() )
4011   {
4012     for ( LDOM_Node aFilter = aSection.getFirstChild();
4013           !aFilter.isNull(); aFilter = aFilter.getNextSibling() )
4014     {
4015       LDOM_Element& anElem = ( LDOM_Element& )aFilter;
4016       aSeq.Append( new TCollection_HAsciiString(
4017          (Standard_CString)anElem.getAttribute( "name" ).GetString() ) );
4018     }
4019   }
4020
4021   anArray->length( aSeq.Length() );
4022   for ( int i = 1, n = aSeq.Length(); i <= n; i++ )
4023     anArray[ i - 1 ] = CORBA::string_dup( aSeq( i )->ToCString() );
4024
4025   return anArray._retn();
4026 }
4027
4028 //=======================================================================
4029 // name    : FilterLibrary_i::GetAllNames
4030 // Purpose : Get names of filters from library
4031 //=======================================================================
4032 string_array* FilterLibrary_i::GetAllNames()
4033 {
4034   string_array_var aResArray = new string_array;
4035   for ( int type = SMESH::ALL; type <= SMESH::VOLUME; type++ )
4036   {
4037     SMESH::string_array_var aNames = GetNames( (SMESH::ElementType)type );
4038
4039     int aPrevLength = aResArray->length();
4040     aResArray->length( aPrevLength + aNames->length() );
4041     for ( int i = 0, n = aNames->length(); i < n; i++ )
4042       aResArray[ aPrevLength + i ] = aNames[ i ];
4043   }
4044
4045   return aResArray._retn();
4046 }
4047
4048 //================================================================================
4049 /*!
4050  * \brief Return an array of strings corresponding to items of enum FunctorType
4051  */
4052 //================================================================================
4053
4054 static const char** getFunctNames()
4055 {
4056   static const char* functName[] = {
4057     // IT's necessary to update this array according to enum FunctorType (SMESH_Filter.idl)
4058     // The order is IMPORTANT !!!
4059     "FT_AspectRatio",
4060     "FT_AspectRatio3D",
4061     "FT_Warping",
4062     "FT_MinimumAngle",
4063     "FT_Taper",
4064     "FT_Skew",
4065     "FT_Area",
4066     "FT_Volume3D",
4067     "FT_MaxElementLength2D",
4068     "FT_MaxElementLength3D",
4069     "FT_FreeBorders",
4070     "FT_FreeEdges",
4071     "FT_FreeNodes",
4072     "FT_FreeFaces",
4073     "FT_EqualNodes",
4074     "FT_EqualEdges",
4075     "FT_EqualFaces",
4076     "FT_EqualVolumes",
4077     "FT_MultiConnection",
4078     "FT_MultiConnection2D",
4079     "FT_Length",
4080     "FT_Length2D",
4081     "FT_NodeConnectivityNumber",
4082     "FT_BelongToMeshGroup",
4083     "FT_BelongToGeom",
4084     "FT_BelongToPlane",
4085     "FT_BelongToCylinder",
4086     "FT_BelongToGenSurface",
4087     "FT_LyingOnGeom",
4088     "FT_RangeOfIds",
4089     "FT_BadOrientedVolume",
4090     "FT_BareBorderVolume",
4091     "FT_BareBorderFace",
4092     "FT_OverConstrainedVolume",
4093     "FT_OverConstrainedFace",
4094     "FT_LinearOrQuadratic",
4095     "FT_GroupColor",
4096     "FT_ElemGeomType",
4097     "FT_EntityType",
4098     "FT_CoplanarFaces",
4099     "FT_BallDiameter",
4100     "FT_ConnectedElements",
4101     "FT_LessThan",
4102     "FT_MoreThan",
4103     "FT_EqualTo",
4104     "FT_LogicalNOT",
4105     "FT_LogicalAND",
4106     "FT_LogicalOR",
4107     "FT_Undefined"};
4108
4109 #ifdef _DEBUG_
4110   // check if functName is complete, compilation failure means that enum FunctorType changed
4111   const int nbFunctors = sizeof(functName) / sizeof(const char*);
4112   int _assert[( nbFunctors == SMESH::FT_Undefined + 1 ) ? 2 : -1 ]; _assert[0]=_assert[1];
4113 #endif
4114
4115   return functName;
4116 }
4117
4118 //================================================================================
4119 /*!
4120  * \brief Return a string corresponding to an item of enum FunctorType
4121  */
4122 //================================================================================
4123
4124 const char* SMESH::FunctorTypeToString(SMESH::FunctorType ft)
4125 {
4126   if ( ft < 0 || ft > SMESH::FT_Undefined )
4127     return "FT_Undefined";
4128   return getFunctNames()[ ft ];
4129 }
4130
4131 //================================================================================
4132 /*!
4133  * \brief Converts a string to FunctorType. This is reverse of FunctorTypeToString()
4134  */
4135 //================================================================================
4136
4137 SMESH::FunctorType SMESH::StringToFunctorType(const char* str)
4138 {
4139   std::string name( str + 3 ); // skip "FT_"
4140   const char** functNames = getFunctNames();
4141   int ft = 0;
4142   for ( ; ft < SMESH::FT_Undefined; ++ft )
4143     if ( name == ( functNames[ft] + 3 ))
4144       break;
4145
4146   //ASSERT( strcmp( str, FunctorTypeToString( SMESH::FunctorType( ft ))) == 0 );
4147
4148   return SMESH::FunctorType( ft );
4149 }
4150
4151 //================================================================================
4152 /*!
4153  * \brief calls OnBaseObjModified(), if who != this, and myWaiters[i]->Modified(who)
4154  */
4155 //================================================================================
4156
4157 void NotifyerAndWaiter::Modified( bool removed, NotifyerAndWaiter* who )
4158 {
4159   if ( who != 0 && who != this )
4160     OnBaseObjModified( who, removed );
4161   else
4162     who = this;
4163
4164   std::list<NotifyerAndWaiter*> waiters = myWaiters; // myWaiters can be changed by Modified()
4165   std::list<NotifyerAndWaiter*>::iterator i = waiters.begin();
4166   for ( ; i != waiters.end(); ++i )
4167     (*i)->Modified( removed, who );
4168 }
4169
4170 //================================================================================
4171 /*!
4172  * \brief Stores an object to be notified on change of predicate
4173  */
4174 //================================================================================
4175
4176 void NotifyerAndWaiter::AddModifWaiter( NotifyerAndWaiter* waiter )
4177 {
4178   if ( waiter )
4179     myWaiters.push_back( waiter );
4180 }
4181
4182 //================================================================================
4183 /*!
4184  * \brief Removes an object to be notified on change of predicate
4185  */
4186 //================================================================================
4187
4188 void NotifyerAndWaiter::RemoveModifWaiter( NotifyerAndWaiter* waiter )
4189 {
4190   myWaiters.remove( waiter );
4191 }
4192
4193 //================================================================================
4194 /*!
4195  * \brief Checks if a waiter is among myWaiters, maybe nested
4196  */
4197 //================================================================================
4198
4199 bool NotifyerAndWaiter::ContainModifWaiter( NotifyerAndWaiter* waiter )
4200 {
4201   bool is = ( waiter == this );
4202
4203   std::list<NotifyerAndWaiter*>::iterator w = myWaiters.begin();
4204   for ( ; !is && w != myWaiters.end(); ++w )
4205     is = (*w)->ContainModifWaiter( waiter );
4206
4207   return is;
4208 }