src/SMESH_I/SMESH_Filter_i.cxx

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