src/SMESH_I/SMESH_Filter_i.cxx

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