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