src/SMESH_I/SMESH_Filter_i.cxx

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