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