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