src/SMESH_I/SMESH_Filter_i.cxx

   1 //  Copyright (C) 2007-2010  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 subshape 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 subshape 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<<".Destroy()";
 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 void NumericalFunctor_i::SetPrecision( CORBA::Long thePrecision )
 600 {
 601   myNumericalFunctorPtr->SetPrecision( thePrecision );
 602   TPythonDump()<<this<<".SetPrecision("<<thePrecision<<")";
 603 }
 604
 605 CORBA::Long NumericalFunctor_i::GetPrecision()
 606 {
 607  return myNumericalFunctorPtr->GetPrecision();
 608 }
 609
 610 Controls::NumericalFunctorPtr NumericalFunctor_i::GetNumericalFunctor()
 611 {
 612   return myNumericalFunctorPtr;
 613 }
 614
 615
 616 /*
 617   Class       : SMESH_MinimumAngle
 618   Description : Functor for calculation of minimum angle
 619 */
 620 MinimumAngle_i::MinimumAngle_i()
 621 {
 622   myNumericalFunctorPtr.reset( new Controls::MinimumAngle() );
 623   myFunctorPtr = myNumericalFunctorPtr;
 624 }
 625
 626 FunctorType MinimumAngle_i::GetFunctorType()
 627 {
 628   return SMESH::FT_MinimumAngle;
 629 }
 630
 631
 632 /*
 633   Class       : AspectRatio
 634   Description : Functor for calculating aspect ratio
 635 */
 636 AspectRatio_i::AspectRatio_i()
 637 {
 638   myNumericalFunctorPtr.reset( new Controls::AspectRatio() );
 639   myFunctorPtr = myNumericalFunctorPtr;
 640 }
 641
 642 FunctorType AspectRatio_i::GetFunctorType()
 643 {
 644   return SMESH::FT_AspectRatio;
 645 }
 646
 647
 648 /*
 649   Class       : AspectRatio3D
 650   Description : Functor for calculating aspect ratio 3D
 651 */
 652 AspectRatio3D_i::AspectRatio3D_i()
 653 {
 654   myNumericalFunctorPtr.reset( new Controls::AspectRatio3D() );
 655   myFunctorPtr = myNumericalFunctorPtr;
 656 }
 657
 658 FunctorType AspectRatio3D_i::GetFunctorType()
 659 {
 660   return SMESH::FT_AspectRatio3D;
 661 }
 662
 663
 664 /*
 665   Class       : Warping_i
 666   Description : Functor for calculating warping
 667 */
 668 Warping_i::Warping_i()
 669 {
 670   myNumericalFunctorPtr.reset( new Controls::Warping() );
 671   myFunctorPtr = myNumericalFunctorPtr;
 672 }
 673
 674 FunctorType Warping_i::GetFunctorType()
 675 {
 676   return SMESH::FT_Warping;
 677 }
 678
 679
 680 /*
 681   Class       : Taper_i
 682   Description : Functor for calculating taper
 683 */
 684 Taper_i::Taper_i()
 685 {
 686   myNumericalFunctorPtr.reset( new Controls::Taper() );
 687   myFunctorPtr = myNumericalFunctorPtr;
 688 }
 689
 690 FunctorType Taper_i::GetFunctorType()
 691 {
 692   return SMESH::FT_Taper;
 693 }
 694
 695
 696 /*
 697   Class       : Skew_i
 698   Description : Functor for calculating skew in degrees
 699 */
 700 Skew_i::Skew_i()
 701 {
 702   myNumericalFunctorPtr.reset( new Controls::Skew() );
 703   myFunctorPtr = myNumericalFunctorPtr;
 704 }
 705
 706 FunctorType Skew_i::GetFunctorType()
 707 {
 708   return SMESH::FT_Skew;
 709 }
 710
 711 /*
 712   Class       : Area_i
 713   Description : Functor for calculating area
 714 */
 715 Area_i::Area_i()
 716 {
 717   myNumericalFunctorPtr.reset( new Controls::Area() );
 718   myFunctorPtr = myNumericalFunctorPtr;
 719 }
 720
 721 FunctorType Area_i::GetFunctorType()
 722 {
 723   return SMESH::FT_Area;
 724 }
 725
 726 /*
 727   Class       : Volume3D_i
 728   Description : Functor for calculating volume of 3D element
 729 */
 730 Volume3D_i::Volume3D_i()
 731 {
 732   myNumericalFunctorPtr.reset( new Controls::Volume() );
 733   myFunctorPtr = myNumericalFunctorPtr;
 734 }
 735
 736 FunctorType Volume3D_i::GetFunctorType()
 737 {
 738   return SMESH::FT_Volume3D;
 739 }
 740
 741 /*
 742   Class       : Length_i
 743   Description : Functor for calculating length off edge
 744 */
 745 Length_i::Length_i()
 746 {
 747   myNumericalFunctorPtr.reset( new Controls::Length() );
 748   myFunctorPtr = myNumericalFunctorPtr;
 749 }
 750
 751 FunctorType Length_i::GetFunctorType()
 752 {
 753   return SMESH::FT_Length;
 754 }
 755
 756 /*
 757   Class       : Length2D_i
 758   Description : Functor for calculating length of edge
 759 */
 760 Length2D_i::Length2D_i()
 761 {
 762   myNumericalFunctorPtr.reset( new Controls::Length2D() );
 763   myFunctorPtr = myNumericalFunctorPtr;
 764 }
 765
 766 FunctorType Length2D_i::GetFunctorType()
 767 {
 768   return SMESH::FT_Length2D;
 769 }
 770
 771 SMESH::Length2D::Values* Length2D_i::GetValues()
 772 {
 773   INFOS("Length2D_i::GetValues");
 774   SMESH::Controls::Length2D::TValues aValues;
 775   myLength2DPtr->GetValues( aValues );
 776
 777   long i = 0, iEnd = aValues.size();
 778
 779   SMESH::Length2D::Values_var aResult = new SMESH::Length2D::Values(iEnd);
 780
 781   SMESH::Controls::Length2D::TValues::const_iterator anIter;
 782   for ( anIter = aValues.begin() ; anIter != aValues.end(); anIter++, i++ )
 783   {
 784     const SMESH::Controls::Length2D::Value&  aVal = *anIter;
 785     SMESH::Length2D::Value &aValue = aResult[ i ];
 786
 787     aValue.myLength = aVal.myLength;
 788     aValue.myPnt1 = aVal.myPntId[ 0 ];
 789     aValue.myPnt2 = aVal.myPntId[ 1 ];
 790   }
 791
 792   INFOS("Length2D_i::GetValuess~");
 793   return aResult._retn();
 794 }
 795
 796 /*
 797   Class       : MultiConnection_i
 798   Description : Functor for calculating number of faces conneted to the edge
 799 */
 800 MultiConnection_i::MultiConnection_i()
 801 {
 802   myNumericalFunctorPtr.reset( new Controls::MultiConnection() );
 803   myFunctorPtr = myNumericalFunctorPtr;
 804 }
 805
 806 FunctorType MultiConnection_i::GetFunctorType()
 807 {
 808   return SMESH::FT_MultiConnection;
 809 }
 810
 811 /*
 812   Class       : MultiConnection2D_i
 813   Description : Functor for calculating number of faces conneted to the edge
 814 */
 815 MultiConnection2D_i::MultiConnection2D_i()
 816 {
 817   myNumericalFunctorPtr.reset( new Controls::MultiConnection2D() );
 818   myFunctorPtr = myNumericalFunctorPtr;
 819 }
 820
 821 FunctorType MultiConnection2D_i::GetFunctorType()
 822 {
 823   return SMESH::FT_MultiConnection2D;
 824 }
 825
 826 SMESH::MultiConnection2D::Values* MultiConnection2D_i::GetValues()
 827 {
 828   INFOS("MultiConnection2D_i::GetValues");
 829   SMESH::Controls::MultiConnection2D::MValues aValues;
 830   myMulticonnection2DPtr->GetValues( aValues );
 831
 832   long i = 0, iEnd = aValues.size();
 833
 834   SMESH::MultiConnection2D::Values_var aResult = new SMESH::MultiConnection2D::Values(iEnd);
 835
 836   SMESH::Controls::MultiConnection2D::MValues::const_iterator anIter;
 837   for ( anIter = aValues.begin() ; anIter != aValues.end(); anIter++, i++ )
 838   {
 839     const SMESH::Controls::MultiConnection2D::Value&  aVal = (*anIter).first;
 840     SMESH::MultiConnection2D::Value &aValue = aResult[ i ];
 841
 842     aValue.myPnt1 = aVal.myPntId[ 0 ];
 843     aValue.myPnt2 = aVal.myPntId[ 1 ];
 844     aValue.myNbConnects = (*anIter).second;
 845   }
 846
 847   INFOS("Multiconnection2D_i::GetValuess~");
 848   return aResult._retn();
 849 }
 850
 851 /*
 852                             PREDICATES
 853 */
 854
 855
 856 /*
 857   Class       : Predicate_i
 858   Description : Base class for all predicates
 859 */
 860 CORBA::Boolean Predicate_i::IsSatisfy( CORBA::Long theId )
 861 {
 862   return myPredicatePtr->IsSatisfy( theId );
 863 }
 864
 865 Controls::PredicatePtr Predicate_i::GetPredicate()
 866 {
 867   return myPredicatePtr;
 868 }
 869
 870 /*
 871   Class       : BadOrientedVolume_i
 872   Description : Verify whether a mesh volume is incorrectly oriented from
 873                 the point of view of MED convention
 874 */
 875 BadOrientedVolume_i::BadOrientedVolume_i()
 876 {
 877   Controls::PredicatePtr control( new Controls::BadOrientedVolume() );
 878   myFunctorPtr = myPredicatePtr = control;
 879 };
 880
 881 FunctorType BadOrientedVolume_i::GetFunctorType()
 882 {
 883   return SMESH::FT_BadOrientedVolume;
 884 }
 885
 886 /*
 887   Class       : BelongToGeom_i
 888   Description : Predicate for selection on geometrical support
 889 */
 890 BelongToGeom_i::BelongToGeom_i()
 891 {
 892   myBelongToGeomPtr.reset( new Controls::BelongToGeom() );
 893   myFunctorPtr = myPredicatePtr = myBelongToGeomPtr;
 894   myShapeName = 0;
 895   myShapeID   = 0;
 896 }
 897
 898 BelongToGeom_i::~BelongToGeom_i()
 899 {
 900   delete myShapeName;
 901   delete myShapeID;
 902 }
 903
 904 void BelongToGeom_i::SetGeom( GEOM::GEOM_Object_ptr theGeom )
 905 {
 906   if ( theGeom->_is_nil() )
 907     return;
 908   SMESH_Gen_i* aSMESHGen = SMESH_Gen_i::GetSMESHGen();
 909   GEOM::GEOM_Gen_ptr aGEOMGen = SMESH_Gen_i::GetGeomEngine();
 910   TopoDS_Shape aLocShape = aSMESHGen->GetShapeReader()->GetShape( aGEOMGen, theGeom );
 911   myBelongToGeomPtr->SetGeom( aLocShape );
 912   TPythonDump()<<this<<".SetGeom("<<theGeom<<")";
 913 }
 914
 915 void BelongToGeom_i::SetGeom( const TopoDS_Shape& theShape )
 916 {
 917   myBelongToGeomPtr->SetGeom( theShape );
 918 }
 919
 920 void BelongToGeom_i::SetElementType(ElementType theType){
 921   myBelongToGeomPtr->SetType(SMDSAbs_ElementType(theType));
 922   TPythonDump()<<this<<".SetElementType("<<theType<<")";
 923 }
 924
 925 FunctorType BelongToGeom_i::GetFunctorType()
 926 {
 927   return SMESH::FT_BelongToGeom;
 928 }
 929
 930 void BelongToGeom_i::SetShapeName( const char* theName )
 931 {
 932   delete myShapeName;
 933   myShapeName = strdup( theName );
 934   myBelongToGeomPtr->SetGeom( getShapeByName( myShapeName ) );
 935   TPythonDump()<<this<<".SetShapeName('"<<theName<<"')";
 936 }
 937
 938 void BelongToGeom_i::SetShape( const char* theID, const char* theName )
 939 {
 940   delete myShapeName;
 941   myShapeName = strdup( theName );
 942   delete myShapeID;
 943   if ( theID )
 944     myShapeID = strdup( theID );
 945   else
 946     myShapeID = 0;
 947
 948   if ( myShapeID && strcmp(myShapeName, getShapeNameByID(myShapeID)) == 0 )
 949     myBelongToGeomPtr->SetGeom( getShapeByID(myShapeID) );
 950   else
 951     myBelongToGeomPtr->SetGeom( getShapeByName( myShapeName ) );
 952 }
 953
 954 char* BelongToGeom_i::GetShapeName()
 955 {
 956   return CORBA::string_dup( myShapeName );
 957 }
 958
 959 char* BelongToGeom_i::GetShapeID()
 960 {
 961   return CORBA::string_dup( myShapeID );
 962 }
 963
 964 void BelongToGeom_i::SetTolerance( CORBA::Double theToler )
 965 {
 966   myBelongToGeomPtr->SetTolerance( theToler );
 967   TPythonDump()<<this<<".SetTolerance("<<theToler<<")";
 968 }
 969
 970 CORBA::Double BelongToGeom_i::GetTolerance()
 971 {
 972   return myBelongToGeomPtr->GetTolerance();
 973 }
 974
 975 /*
 976   Class       : BelongToSurface_i
 977   Description : Predicate for selection on geometrical support
 978 */
 979 BelongToSurface_i::BelongToSurface_i( const Handle(Standard_Type)& theSurfaceType )
 980 {
 981   myElementsOnSurfacePtr.reset( new Controls::ElementsOnSurface() );
 982   myFunctorPtr = myPredicatePtr = myElementsOnSurfacePtr;
 983   myShapeName = 0;
 984   myShapeID   = 0;
 985   mySurfaceType = theSurfaceType;
 986 }
 987
 988 BelongToSurface_i::~BelongToSurface_i()
 989 {
 990   delete myShapeName;
 991   delete myShapeID;
 992 }
 993
 994 void BelongToSurface_i::SetSurface( GEOM::GEOM_Object_ptr theGeom, ElementType theType )
 995 {
 996   if ( theGeom->_is_nil() )
 997     return;
 998   SMESH_Gen_i* aSMESHGen = SMESH_Gen_i::GetSMESHGen();
 999   GEOM::GEOM_Gen_ptr aGEOMGen = SMESH_Gen_i::GetGeomEngine();
1000   TopoDS_Shape aLocShape = aSMESHGen->GetShapeReader()->GetShape( aGEOMGen, theGeom );
1001
1002   if ( aLocShape.ShapeType() == TopAbs_FACE )
1003   {
1004     Handle(Geom_Surface) aSurf = BRep_Tool::Surface( TopoDS::Face( aLocShape ) );
1005     if ( !aSurf.IsNull() && aSurf->DynamicType() == mySurfaceType )
1006     {
1007       myElementsOnSurfacePtr->SetSurface( aLocShape, (SMDSAbs_ElementType)theType );
1008       return;
1009     }
1010   }
1011
1012   myElementsOnSurfacePtr->SetSurface( TopoDS_Shape(), (SMDSAbs_ElementType)theType );
1013 }
1014
1015 void BelongToSurface_i::SetShapeName( const char* theName, ElementType theType )
1016 {
1017   delete myShapeName;
1018   myShapeName = strdup( theName );
1019   myElementsOnSurfacePtr->SetSurface( getShapeByName( myShapeName ), (SMDSAbs_ElementType)theType );
1020   TPythonDump()<<this<<".SetShapeName('"<<theName<<"',"<<theType<<")";
1021 }
1022
1023 void BelongToSurface_i::SetShape( const char* theID,  const char* theName, ElementType theType )
1024 {
1025   delete myShapeName;
1026   myShapeName = strdup( theName );
1027   delete myShapeID;
1028   if ( theID )
1029     myShapeID = strdup( theID );
1030   else
1031     myShapeID = 0;
1032
1033   if ( myShapeID && strcmp(myShapeName, getShapeNameByID(myShapeID)) == 0 )
1034     myElementsOnSurfacePtr->SetSurface( getShapeByID(myShapeID), (SMDSAbs_ElementType)theType );
1035   else
1036     myElementsOnSurfacePtr->SetSurface( getShapeByName( myShapeName ), (SMDSAbs_ElementType)theType );
1037 }
1038
1039 char* BelongToSurface_i::GetShapeName()
1040 {
1041   return CORBA::string_dup( myShapeName );
1042 }
1043
1044 char* BelongToSurface_i::GetShapeID()
1045 {
1046   return CORBA::string_dup( myShapeID );
1047 }
1048
1049 void BelongToSurface_i::SetTolerance( CORBA::Double theToler )
1050 {
1051   myElementsOnSurfacePtr->SetTolerance( theToler );
1052   TPythonDump()<<this<<".SetTolerance("<<theToler<<")";
1053 }
1054
1055 CORBA::Double BelongToSurface_i::GetTolerance()
1056 {
1057   return myElementsOnSurfacePtr->GetTolerance();
1058 }
1059
1060 void BelongToSurface_i::SetUseBoundaries( CORBA::Boolean theUseBndRestrictions )
1061 {
1062   myElementsOnSurfacePtr->SetUseBoundaries( theUseBndRestrictions );
1063   TPythonDump()<<this<<".SetUseBoundaries( " << theUseBndRestrictions << " )";
1064 }
1065
1066 CORBA::Boolean BelongToSurface_i::GetUseBoundaries()
1067 {
1068   return myElementsOnSurfacePtr->GetUseBoundaries();
1069 }
1070
1071
1072 /*
1073   Class       : BelongToPlane_i
1074   Description : Verify whether mesh element lie in pointed Geom planar object
1075 */
1076
1077 BelongToPlane_i::BelongToPlane_i()
1078 : BelongToSurface_i( STANDARD_TYPE( Geom_Plane ) )
1079 {
1080 }
1081
1082 void BelongToPlane_i::SetPlane( GEOM::GEOM_Object_ptr theGeom, ElementType theType )
1083 {
1084   BelongToSurface_i::SetSurface( theGeom, theType );
1085   TPythonDump()<<this<<".SetPlane("<<theGeom<<","<<theType<<")";
1086 }
1087
1088 FunctorType BelongToPlane_i::GetFunctorType()
1089 {
1090   return FT_BelongToPlane;
1091 }
1092
1093 /*
1094   Class       : BelongToCylinder_i
1095   Description : Verify whether mesh element lie in pointed Geom planar object
1096 */
1097
1098 BelongToCylinder_i::BelongToCylinder_i()
1099 : BelongToSurface_i( STANDARD_TYPE( Geom_CylindricalSurface ) )
1100 {
1101 }
1102
1103 void BelongToCylinder_i::SetCylinder( GEOM::GEOM_Object_ptr theGeom, ElementType theType )
1104 {
1105   BelongToSurface_i::SetSurface( theGeom, theType );
1106   TPythonDump()<<this<<".SetCylinder("<<theGeom<<","<<theType<<")";
1107 }
1108
1109 FunctorType BelongToCylinder_i::GetFunctorType()
1110 {
1111   return FT_BelongToCylinder;
1112 }
1113
1114 /*
1115   Class       : BelongToGenSurface_i
1116   Description : Verify whether mesh element lie in pointed Geom planar object
1117 */
1118
1119 BelongToGenSurface_i::BelongToGenSurface_i()
1120 : BelongToSurface_i( STANDARD_TYPE( Geom_CylindricalSurface ) )
1121 {
1122 }
1123
1124 void BelongToGenSurface_i::SetSurface( GEOM::GEOM_Object_ptr theGeom, ElementType theType )
1125 {
1126   if ( theGeom->_is_nil() )
1127     return;
1128   TopoDS_Shape aLocShape = SMESH_Gen_i::GetSMESHGen()->GeomObjectToShape( theGeom );
1129   if ( !aLocShape.IsNull() && aLocShape.ShapeType() != TopAbs_FACE )
1130     aLocShape.Nullify();
1131
1132   BelongToSurface_i::myElementsOnSurfacePtr->SetSurface( aLocShape, (SMDSAbs_ElementType)theType );
1133   TPythonDump()<<this<<".SetGenSurface("<<theGeom<<","<<theType<<")";
1134 }
1135
1136 FunctorType BelongToGenSurface_i::GetFunctorType()
1137 {
1138   return FT_BelongToGenSurface;
1139 }
1140
1141 /*
1142   Class       : LyingOnGeom_i
1143   Description : Predicate for selection on geometrical support
1144 */
1145 LyingOnGeom_i::LyingOnGeom_i()
1146 {
1147   myLyingOnGeomPtr.reset( new Controls::LyingOnGeom() );
1148   myFunctorPtr = myPredicatePtr = myLyingOnGeomPtr;
1149   myShapeName = 0;
1150   myShapeID = 0;
1151 }
1152
1153 LyingOnGeom_i::~LyingOnGeom_i()
1154 {
1155   delete myShapeName;
1156   delete myShapeID;
1157 }
1158
1159 void LyingOnGeom_i::SetGeom( GEOM::GEOM_Object_ptr theGeom )
1160 {
1161   if ( theGeom->_is_nil() )
1162     return;
1163   SMESH_Gen_i* aSMESHGen = SMESH_Gen_i::GetSMESHGen();
1164   GEOM::GEOM_Gen_ptr aGEOMGen = SMESH_Gen_i::GetGeomEngine();
1165   TopoDS_Shape aLocShape = aSMESHGen->GetShapeReader()->GetShape( aGEOMGen, theGeom );
1166   myLyingOnGeomPtr->SetGeom( aLocShape );
1167   TPythonDump()<<this<<".SetGeom("<<theGeom<<")";
1168 }
1169
1170 void LyingOnGeom_i::SetGeom( const TopoDS_Shape& theShape )
1171 {
1172   myLyingOnGeomPtr->SetGeom( theShape );
1173 }
1174
1175 void LyingOnGeom_i::SetElementType(ElementType theType){
1176   myLyingOnGeomPtr->SetType(SMDSAbs_ElementType(theType));
1177   TPythonDump()<<this<<".SetElementType("<<theType<<")";
1178 }
1179
1180 FunctorType LyingOnGeom_i::GetFunctorType()
1181 {
1182   return SMESH::FT_LyingOnGeom;
1183 }
1184
1185 void LyingOnGeom_i::SetShapeName( const char* theName )
1186 {
1187   delete myShapeName;
1188   myShapeName = strdup( theName );
1189   myLyingOnGeomPtr->SetGeom( getShapeByName( myShapeName ) );
1190   TPythonDump()<<this<<".SetShapeName('"<<theName<<"')";
1191 }
1192
1193 void LyingOnGeom_i::SetShape( const char* theID, const char* theName )
1194 {
1195   delete myShapeName;
1196   myShapeName = strdup( theName );
1197   delete myShapeID;
1198   if ( theID )
1199     myShapeID = strdup( theID );
1200   else
1201     myShapeID = 0;
1202
1203   if ( myShapeID && strcmp(myShapeName, getShapeNameByID(myShapeID)) == 0 )
1204     myLyingOnGeomPtr->SetGeom( getShapeByID(myShapeID) );
1205   else
1206     myLyingOnGeomPtr->SetGeom( getShapeByName( myShapeName ) );
1207 }
1208
1209 char* LyingOnGeom_i::GetShapeName()
1210 {
1211   return CORBA::string_dup( myShapeName );
1212 }
1213
1214 char* LyingOnGeom_i::GetShapeID()
1215 {
1216   return CORBA::string_dup( myShapeID );
1217 }
1218
1219 void LyingOnGeom_i::SetTolerance( CORBA::Double theToler )
1220 {
1221   myLyingOnGeomPtr->SetTolerance( theToler );
1222   TPythonDump()<<this<<".SetTolerance("<<theToler<<")";
1223 }
1224
1225 CORBA::Double LyingOnGeom_i::GetTolerance()
1226 {
1227   return myLyingOnGeomPtr->GetTolerance();
1228 }
1229
1230 /*
1231   Class       : FreeBorders_i
1232   Description : Predicate for free borders
1233 */
1234 FreeBorders_i::FreeBorders_i()
1235 {
1236   myPredicatePtr.reset(new Controls::FreeBorders());
1237   myFunctorPtr = myPredicatePtr;
1238 }
1239
1240 FunctorType FreeBorders_i::GetFunctorType()
1241 {
1242   return SMESH::FT_FreeBorders;
1243 }
1244
1245 /*
1246   Class       : FreeEdges_i
1247   Description : Predicate for free borders
1248 */
1249 FreeEdges_i::FreeEdges_i()
1250 : myFreeEdgesPtr( new Controls::FreeEdges() )
1251 {
1252   myFunctorPtr = myPredicatePtr = myFreeEdgesPtr;
1253 }
1254
1255 SMESH::FreeEdges::Borders* FreeEdges_i::GetBorders()
1256 {
1257   INFOS("FreeEdges_i::GetBorders");
1258   SMESH::Controls::FreeEdges::TBorders aBorders;
1259   myFreeEdgesPtr->GetBoreders( aBorders );
1260
1261   long i = 0, iEnd = aBorders.size();
1262
1263   SMESH::FreeEdges::Borders_var aResult = new SMESH::FreeEdges::Borders;
1264   aResult->length(iEnd);
1265
1266   SMESH::Controls::FreeEdges::TBorders::const_iterator anIter;
1267   for ( anIter = aBorders.begin() ; anIter != aBorders.end(); anIter++, i++ )
1268   {
1269     const SMESH::Controls::FreeEdges::Border&  aBord = *anIter;
1270     SMESH::FreeEdges::Border &aBorder = aResult[ i ];
1271
1272     aBorder.myElemId = aBord.myElemId;
1273     aBorder.myPnt1 = aBord.myPntId[ 0 ];
1274     aBorder.myPnt2 = aBord.myPntId[ 1 ];
1275   }
1276
1277   INFOS("FreeEdges_i::GetBorders~");
1278   return aResult._retn();
1279 }
1280
1281 FunctorType FreeEdges_i::GetFunctorType()
1282 {
1283   return SMESH::FT_FreeEdges;
1284 }
1285
1286 /*
1287   Class       : FreeFaces_i
1288   Description : Predicate for free faces
1289 */
1290 FreeFaces_i::FreeFaces_i()
1291 {
1292   myPredicatePtr.reset(new Controls::FreeFaces());
1293   myFunctorPtr = myPredicatePtr;
1294 }
1295
1296 FunctorType FreeFaces_i::GetFunctorType()
1297 {
1298   return SMESH::FT_FreeFaces;
1299 }
1300
1301 /*
1302   Class       : FreeNodes_i
1303   Description : Predicate for free nodes
1304 */
1305 FreeNodes_i::FreeNodes_i()
1306 {
1307   myPredicatePtr.reset(new Controls::FreeNodes());
1308   myFunctorPtr = myPredicatePtr;
1309 }
1310
1311 FunctorType FreeNodes_i::GetFunctorType()
1312 {
1313   return SMESH::FT_FreeNodes;
1314 }
1315
1316 /*
1317   Class       : RangeOfIds_i
1318   Description : Predicate for Range of Ids.
1319                 Range may be specified with two ways.
1320                 1. Using AddToRange method
1321                 2. With SetRangeStr method. Parameter of this method is a string
1322                    like as "1,2,3,50-60,63,67,70-"
1323 */
1324
1325 RangeOfIds_i::RangeOfIds_i()
1326 {
1327   myRangeOfIdsPtr.reset( new Controls::RangeOfIds() );
1328   myFunctorPtr = myPredicatePtr = myRangeOfIdsPtr;
1329 }
1330
1331 void RangeOfIds_i::SetRange( const SMESH::long_array& theIds )
1332 {
1333   CORBA::Long iEnd = theIds.length();
1334   for ( CORBA::Long i = 0; i < iEnd; i++ )
1335     myRangeOfIdsPtr->AddToRange( theIds[ i ] );
1336   TPythonDump()<<this<<".SetRange("<<theIds<<")";
1337 }
1338
1339 CORBA::Boolean RangeOfIds_i::SetRangeStr( const char* theRange )
1340 {
1341   TPythonDump()<<this<<".SetRangeStr('"<<theRange<<"')";
1342   return myRangeOfIdsPtr->SetRangeStr(
1343     TCollection_AsciiString( (Standard_CString)theRange ) );
1344 }
1345
1346 char* RangeOfIds_i::GetRangeStr()
1347 {
1348   TCollection_AsciiString aStr;
1349   myRangeOfIdsPtr->GetRangeStr( aStr );
1350   return CORBA::string_dup( aStr.ToCString() );
1351 }
1352
1353 void RangeOfIds_i::SetElementType( ElementType theType )
1354 {
1355   myRangeOfIdsPtr->SetType( SMDSAbs_ElementType( theType ) );
1356   TPythonDump()<<this<<".SetElementType("<<theType<<")";
1357 }
1358
1359 FunctorType RangeOfIds_i::GetFunctorType()
1360 {
1361   return SMESH::FT_RangeOfIds;
1362 }
1363
1364 /*
1365   Class       : LinearOrQuadratic_i
1366   Description : Predicate to verify whether a mesh element is linear
1367 */
1368 LinearOrQuadratic_i::LinearOrQuadratic_i()
1369 {
1370   myLinearOrQuadraticPtr.reset(new Controls::LinearOrQuadratic());
1371   myFunctorPtr = myPredicatePtr = myLinearOrQuadraticPtr;
1372 }
1373
1374 void LinearOrQuadratic_i::SetElementType(ElementType theType)
1375 {
1376   myLinearOrQuadraticPtr->SetType(SMDSAbs_ElementType(theType));
1377   TPythonDump()<<this<<".SetElementType("<<theType<<")";
1378 }
1379
1380 FunctorType LinearOrQuadratic_i::GetFunctorType()
1381 {
1382   return SMESH::FT_LinearOrQuadratic;
1383 }
1384
1385 /*
1386   Class       : GroupColor_i
1387   Description : Functor for check color of group to whic mesh element belongs to
1388 */
1389 GroupColor_i::GroupColor_i()
1390 {
1391   myGroupColorPtr.reset(new Controls::GroupColor());
1392   myFunctorPtr = myPredicatePtr = myGroupColorPtr;
1393 }
1394
1395 FunctorType GroupColor_i::GetFunctorType()
1396 {
1397   return SMESH::FT_GroupColor;
1398 }
1399
1400 void GroupColor_i::SetColorStr( const char* theColor )
1401 {
1402   myGroupColorPtr->SetColorStr(
1403     TCollection_AsciiString( (Standard_CString)theColor ) );
1404   TPythonDump()<<this<<".SetColorStr('"<<theColor<<"')";
1405 }
1406
1407 char* GroupColor_i::GetColorStr()
1408 {
1409   TCollection_AsciiString aStr;
1410   myGroupColorPtr->GetColorStr( aStr );
1411   return CORBA::string_dup( aStr.ToCString() );
1412 }
1413
1414 void GroupColor_i::SetElementType(ElementType theType)
1415 {
1416   myGroupColorPtr->SetType(SMDSAbs_ElementType(theType));
1417   TPythonDump()<<this<<".SetElementType("<<theType<<")";
1418 }
1419
1420 /*
1421   Class       : ElemGeomType_i
1422   Description : Predicate check is element has indicated geometry type
1423 */
1424 ElemGeomType_i::ElemGeomType_i()
1425 {
1426   myElemGeomTypePtr.reset(new Controls::ElemGeomType());
1427   myFunctorPtr = myPredicatePtr = myElemGeomTypePtr;
1428 }
1429
1430 void ElemGeomType_i::SetElementType(ElementType theType)
1431 {
1432   myElemGeomTypePtr->SetType(SMDSAbs_ElementType(theType));
1433   TPythonDump()<<this<<".SetElementType("<<theType<<")";
1434 }
1435
1436 void ElemGeomType_i::SetGeometryType(GeometryType theType)
1437 {
1438   myElemGeomTypePtr->SetGeomType(SMDSAbs_GeometryType(theType));
1439   TPythonDump()<<this<<".SetGeometryType("<<theType<<")";
1440 }
1441
1442 GeometryType ElemGeomType_i::GetGeometryType() const
1443 {
1444   return (GeometryType)myElemGeomTypePtr->GetGeomType();;
1445 }
1446
1447 FunctorType ElemGeomType_i::GetFunctorType()
1448 {
1449   return SMESH::FT_ElemGeomType;
1450 }
1451
1452 /*
1453   Class       : Comparator_i
1454   Description : Base class for comparators
1455 */
1456 Comparator_i::Comparator_i():
1457   myNumericalFunctor( NULL )
1458 {}
1459
1460 Comparator_i::~Comparator_i()
1461 {
1462   if ( myNumericalFunctor )
1463     myNumericalFunctor->Destroy();
1464 }
1465
1466 void Comparator_i::SetMargin( CORBA::Double theValue )
1467 {
1468   myComparatorPtr->SetMargin( theValue );
1469   TPythonDump()<<this<<".SetMargin("<<theValue<<")";
1470 }
1471
1472 CORBA::Double Comparator_i::GetMargin()
1473 {
1474   return myComparatorPtr->GetMargin();
1475 }
1476
1477 void Comparator_i::SetNumFunctor( NumericalFunctor_ptr theFunct )
1478 {
1479   if ( myNumericalFunctor )
1480     myNumericalFunctor->Destroy();
1481
1482   myNumericalFunctor = DownCast<NumericalFunctor_i*>(theFunct);
1483
1484   if ( myNumericalFunctor )
1485   {
1486     myComparatorPtr->SetNumFunctor( myNumericalFunctor->GetNumericalFunctor() );
1487     myNumericalFunctor->Register();
1488     TPythonDump()<<this<<".SetNumFunctor("<<myNumericalFunctor<<")";
1489   }
1490 }
1491
1492 Controls::ComparatorPtr Comparator_i::GetComparator()
1493 {
1494   return myComparatorPtr;
1495 }
1496
1497 NumericalFunctor_i* Comparator_i::GetNumFunctor_i()
1498 {
1499   return myNumericalFunctor;
1500 }
1501
1502
1503 /*
1504   Class       : LessThan_i
1505   Description : Comparator "<"
1506 */
1507 LessThan_i::LessThan_i()
1508 {
1509   myComparatorPtr.reset( new Controls::LessThan() );
1510   myFunctorPtr = myPredicatePtr = myComparatorPtr;
1511 }
1512
1513 FunctorType LessThan_i::GetFunctorType()
1514 {
1515   return SMESH::FT_LessThan;
1516 }
1517
1518
1519 /*
1520   Class       : MoreThan_i
1521   Description : Comparator ">"
1522 */
1523 MoreThan_i::MoreThan_i()
1524 {
1525   myComparatorPtr.reset( new Controls::MoreThan() );
1526   myFunctorPtr = myPredicatePtr = myComparatorPtr;
1527 }
1528
1529 FunctorType MoreThan_i::GetFunctorType()
1530 {
1531   return SMESH::FT_MoreThan;
1532 }
1533
1534
1535 /*
1536   Class       : EqualTo_i
1537   Description : Comparator "="
1538 */
1539 EqualTo_i::EqualTo_i()
1540 : myEqualToPtr( new Controls::EqualTo() )
1541 {
1542   myFunctorPtr = myPredicatePtr = myComparatorPtr = myEqualToPtr;
1543 }
1544
1545 void EqualTo_i::SetTolerance( CORBA::Double theToler )
1546 {
1547   myEqualToPtr->SetTolerance( theToler );
1548   TPythonDump()<<this<<".SetTolerance("<<theToler<<")";
1549 }
1550
1551 CORBA::Double EqualTo_i::GetTolerance()
1552 {
1553   return myEqualToPtr->GetTolerance();
1554 }
1555
1556 FunctorType EqualTo_i::GetFunctorType()
1557 {
1558   return SMESH::FT_EqualTo;
1559 }
1560
1561 /*
1562   Class       : LogicalNOT_i
1563   Description : Logical NOT predicate
1564 */
1565 LogicalNOT_i::LogicalNOT_i()
1566 : myPredicate( NULL ),
1567   myLogicalNOTPtr( new Controls::LogicalNOT() )
1568 {
1569   myFunctorPtr = myPredicatePtr = myLogicalNOTPtr;
1570 }
1571
1572 LogicalNOT_i::~LogicalNOT_i()
1573 {
1574   if ( myPredicate )
1575     myPredicate->Destroy();
1576 }
1577
1578 void LogicalNOT_i::SetPredicate( Predicate_ptr thePredicate )
1579 {
1580   if ( myPredicate )
1581     myPredicate->Destroy();
1582
1583   myPredicate = SMESH::GetPredicate(thePredicate);
1584
1585   if ( myPredicate ){
1586     myLogicalNOTPtr->SetPredicate(myPredicate->GetPredicate());
1587     myPredicate->Register();
1588     TPythonDump()<<this<<".SetPredicate("<<myPredicate<<")";
1589   }
1590 }
1591
1592 FunctorType LogicalNOT_i::GetFunctorType()
1593 {
1594   return SMESH::FT_LogicalNOT;
1595 }
1596
1597 Predicate_i* LogicalNOT_i::GetPredicate_i()
1598 {
1599   return myPredicate;
1600 }
1601
1602
1603 /*
1604   Class       : LogicalBinary_i
1605   Description : Base class for binary logical predicate
1606 */
1607 LogicalBinary_i::LogicalBinary_i()
1608 : myPredicate1( NULL ),
1609   myPredicate2( NULL )
1610 {}
1611
1612 LogicalBinary_i::~LogicalBinary_i()
1613 {
1614   if ( myPredicate1 )
1615     myPredicate1->Destroy();
1616
1617   if ( myPredicate2 )
1618     myPredicate2->Destroy();
1619 }
1620
1621 void LogicalBinary_i::SetMesh( SMESH_Mesh_ptr theMesh )
1622 {
1623   if ( myPredicate1 )
1624     myPredicate1->SetMesh( theMesh );
1625
1626   if ( myPredicate2 )
1627     myPredicate2->SetMesh( theMesh );
1628 }
1629
1630 void LogicalBinary_i::SetPredicate1( Predicate_ptr thePredicate )
1631 {
1632   if ( myPredicate1 )
1633     myPredicate1->Destroy();
1634
1635   myPredicate1 = SMESH::GetPredicate(thePredicate);
1636
1637   if ( myPredicate1 ){
1638     myLogicalBinaryPtr->SetPredicate1(myPredicate1->GetPredicate());
1639     myPredicate1->Register();
1640     TPythonDump()<<this<<".SetPredicate1("<<myPredicate1<<")";
1641   }
1642 }
1643
1644 void LogicalBinary_i::SetPredicate2( Predicate_ptr thePredicate )
1645 {
1646   if ( myPredicate2 )
1647     myPredicate2->Destroy();
1648
1649   myPredicate2 = SMESH::GetPredicate(thePredicate);
1650
1651   if ( myPredicate2 ){
1652     myLogicalBinaryPtr->SetPredicate2(myPredicate2->GetPredicate());
1653     myPredicate2->Register();
1654     TPythonDump()<<this<<".SetPredicate2("<<myPredicate2<<")";
1655   }
1656 }
1657
1658 Controls::LogicalBinaryPtr LogicalBinary_i::GetLogicalBinary()
1659 {
1660   return myLogicalBinaryPtr;
1661 }
1662
1663 Predicate_i* LogicalBinary_i::GetPredicate1_i()
1664 {
1665   return myPredicate1;
1666 }
1667 Predicate_i* LogicalBinary_i::GetPredicate2_i()
1668 {
1669   return myPredicate2;
1670 }
1671
1672
1673 /*
1674   Class       : LogicalAND_i
1675   Description : Logical AND
1676 */
1677 LogicalAND_i::LogicalAND_i()
1678 {
1679   myLogicalBinaryPtr.reset( new Controls::LogicalAND() );
1680   myFunctorPtr = myPredicatePtr = myLogicalBinaryPtr;
1681 }
1682
1683 FunctorType LogicalAND_i::GetFunctorType()
1684 {
1685   return SMESH::FT_LogicalAND;
1686 }
1687
1688
1689 /*
1690   Class       : LogicalOR_i
1691   Description : Logical OR
1692 */
1693 LogicalOR_i::LogicalOR_i()
1694 {
1695   myLogicalBinaryPtr.reset( new Controls::LogicalOR() );
1696   myFunctorPtr = myPredicatePtr = myLogicalBinaryPtr;
1697 }
1698
1699 FunctorType LogicalOR_i::GetFunctorType()
1700 {
1701   return SMESH::FT_LogicalOR;
1702 }
1703
1704
1705 /*
1706                             FILTER MANAGER
1707 */
1708
1709 FilterManager_i::FilterManager_i()
1710 : SALOME::GenericObj_i( SMESH_Gen_i::GetPOA() )
1711 {
1712   //Base class Salome_GenericObject do it inmplicitly by overriding PortableServer::POA_ptr _default_POA() method
1713   //PortableServer::ObjectId_var anObjectId =
1714   //  SMESH_Gen_i::GetPOA()->activate_object( this );
1715 }
1716
1717
1718 FilterManager_i::~FilterManager_i()
1719 {
1720   //TPythonDump()<<this<<".Destroy()";
1721 }
1722
1723
1724 MinimumAngle_ptr FilterManager_i::CreateMinimumAngle()
1725 {
1726   SMESH::MinimumAngle_i* aServant = new SMESH::MinimumAngle_i();
1727   SMESH::MinimumAngle_var anObj = aServant->_this();
1728   TPythonDump()<<aServant<<" = "<<this<<".CreateMinimumAngle()";
1729   return anObj._retn();
1730 }
1731
1732
1733 AspectRatio_ptr FilterManager_i::CreateAspectRatio()
1734 {
1735   SMESH::AspectRatio_i* aServant = new SMESH::AspectRatio_i();
1736   SMESH::AspectRatio_var anObj = aServant->_this();
1737   TPythonDump()<<aServant<<" = "<<this<<".CreateAspectRatio()";
1738   return anObj._retn();
1739 }
1740
1741
1742 AspectRatio3D_ptr FilterManager_i::CreateAspectRatio3D()
1743 {
1744   SMESH::AspectRatio3D_i* aServant = new SMESH::AspectRatio3D_i();
1745   SMESH::AspectRatio3D_var anObj = aServant->_this();
1746   TPythonDump()<<aServant<<" = "<<this<<".CreateAspectRatio3D()";
1747   return anObj._retn();
1748 }
1749
1750
1751 Warping_ptr FilterManager_i::CreateWarping()
1752 {
1753   SMESH::Warping_i* aServant = new SMESH::Warping_i();
1754   SMESH::Warping_var anObj = aServant->_this();
1755   TPythonDump()<<aServant<<" = "<<this<<".CreateWarping()";
1756   return anObj._retn();
1757 }
1758
1759
1760 Taper_ptr FilterManager_i::CreateTaper()
1761 {
1762   SMESH::Taper_i* aServant = new SMESH::Taper_i();
1763   SMESH::Taper_var anObj = aServant->_this();
1764   TPythonDump()<<aServant<<" = "<<this<<".CreateTaper()";
1765   return anObj._retn();
1766 }
1767
1768
1769 Skew_ptr FilterManager_i::CreateSkew()
1770 {
1771   SMESH::Skew_i* aServant = new SMESH::Skew_i();
1772   SMESH::Skew_var anObj = aServant->_this();
1773   TPythonDump()<<aServant<<" = "<<this<<".CreateSkew()";
1774   return anObj._retn();
1775 }
1776
1777
1778 Area_ptr FilterManager_i::CreateArea()
1779 {
1780   SMESH::Area_i* aServant = new SMESH::Area_i();
1781   SMESH::Area_var anObj = aServant->_this();
1782   TPythonDump()<<aServant<<" = "<<this<<".CreateArea()";
1783   return anObj._retn();
1784 }
1785
1786
1787 Volume3D_ptr FilterManager_i::CreateVolume3D()
1788 {
1789   SMESH::Volume3D_i* aServant = new SMESH::Volume3D_i();
1790   SMESH::Volume3D_var anObj = aServant->_this();
1791   TPythonDump()<<aServant<<" = "<<this<<".CreateVolume3D()";
1792   return anObj._retn();
1793 }
1794
1795
1796 Length_ptr FilterManager_i::CreateLength()
1797 {
1798   SMESH::Length_i* aServant = new SMESH::Length_i();
1799   SMESH::Length_var anObj = aServant->_this();
1800   TPythonDump()<<aServant<<" = "<<this<<".CreateLength()";
1801   return anObj._retn();
1802 }
1803
1804 Length2D_ptr FilterManager_i::CreateLength2D()
1805 {
1806   SMESH::Length2D_i* aServant = new SMESH::Length2D_i();
1807   SMESH::Length2D_var anObj = aServant->_this();
1808   TPythonDump()<<aServant<<" = "<<this<<".CreateLength2D()";
1809   return anObj._retn();
1810 }
1811
1812 MultiConnection_ptr FilterManager_i::CreateMultiConnection()
1813 {
1814   SMESH::MultiConnection_i* aServant = new SMESH::MultiConnection_i();
1815   SMESH::MultiConnection_var anObj = aServant->_this();
1816   TPythonDump()<<aServant<<" = "<<this<<".CreateMultiConnection()";
1817   return anObj._retn();
1818 }
1819
1820 MultiConnection2D_ptr FilterManager_i::CreateMultiConnection2D()
1821 {
1822   SMESH::MultiConnection2D_i* aServant = new SMESH::MultiConnection2D_i();
1823   SMESH::MultiConnection2D_var anObj = aServant->_this();
1824   TPythonDump()<<aServant<<" = "<<this<<".CreateMultiConnection2D()";
1825   return anObj._retn();
1826 }
1827
1828 BelongToGeom_ptr FilterManager_i::CreateBelongToGeom()
1829 {
1830   SMESH::BelongToGeom_i* aServant = new SMESH::BelongToGeom_i();
1831   SMESH::BelongToGeom_var anObj = aServant->_this();
1832   TPythonDump()<<aServant<<" = "<<this<<".CreateBelongToGeom()";
1833   return anObj._retn();
1834 }
1835
1836 BelongToPlane_ptr FilterManager_i::CreateBelongToPlane()
1837 {
1838   SMESH::BelongToPlane_i* aServant = new SMESH::BelongToPlane_i();
1839   SMESH::BelongToPlane_var anObj = aServant->_this();
1840   TPythonDump()<<aServant<<" = "<<this<<".CreateBelongToPlane()";
1841   return anObj._retn();
1842 }
1843
1844 BelongToCylinder_ptr FilterManager_i::CreateBelongToCylinder()
1845 {
1846   SMESH::BelongToCylinder_i* aServant = new SMESH::BelongToCylinder_i();
1847   SMESH::BelongToCylinder_var anObj = aServant->_this();
1848   TPythonDump()<<aServant<<" = "<<this<<".CreateBelongToCylinder()";
1849   return anObj._retn();
1850 }
1851
1852 BelongToGenSurface_ptr FilterManager_i::CreateBelongToGenSurface()
1853 {
1854   SMESH::BelongToGenSurface_i* aServant = new SMESH::BelongToGenSurface_i();
1855   SMESH::BelongToGenSurface_var anObj = aServant->_this();
1856   TPythonDump()<<aServant<<" = "<<this<<".CreateBelongToGenSurface()";
1857   return anObj._retn();
1858 }
1859
1860 LyingOnGeom_ptr FilterManager_i::CreateLyingOnGeom()
1861 {
1862   SMESH::LyingOnGeom_i* aServant = new SMESH::LyingOnGeom_i();
1863   SMESH::LyingOnGeom_var anObj = aServant->_this();
1864   TPythonDump()<<aServant<<" = "<<this<<".CreateLyingOnGeom()";
1865   return anObj._retn();
1866 }
1867
1868 FreeBorders_ptr FilterManager_i::CreateFreeBorders()
1869 {
1870   SMESH::FreeBorders_i* aServant = new SMESH::FreeBorders_i();
1871   SMESH::FreeBorders_var anObj = aServant->_this();
1872   TPythonDump()<<aServant<<" = "<<this<<".CreateFreeBorders()";
1873   return anObj._retn();
1874 }
1875
1876 FreeEdges_ptr FilterManager_i::CreateFreeEdges()
1877 {
1878   SMESH::FreeEdges_i* aServant = new SMESH::FreeEdges_i();
1879   SMESH::FreeEdges_var anObj = aServant->_this();
1880   TPythonDump()<<aServant<<" = "<<this<<".CreateFreeEdges()";
1881   return anObj._retn();
1882 }
1883
1884 FreeFaces_ptr FilterManager_i::CreateFreeFaces()
1885 {
1886   SMESH::FreeFaces_i* aServant = new SMESH::FreeFaces_i();
1887   SMESH::FreeFaces_var anObj = aServant->_this();
1888   TPythonDump()<<aServant<<" = "<<this<<".CreateFreeFaces()";
1889   return anObj._retn();
1890 }
1891
1892 FreeNodes_ptr FilterManager_i::CreateFreeNodes()
1893 {
1894   SMESH::FreeNodes_i* aServant = new SMESH::FreeNodes_i();
1895   SMESH::FreeNodes_var anObj = aServant->_this();
1896   TPythonDump()<<aServant<<" = "<<this<<".CreateFreeNodes()";
1897   return anObj._retn();
1898 }
1899
1900 RangeOfIds_ptr FilterManager_i::CreateRangeOfIds()
1901 {
1902   SMESH::RangeOfIds_i* aServant = new SMESH::RangeOfIds_i();
1903   SMESH::RangeOfIds_var anObj = aServant->_this();
1904   TPythonDump()<<aServant<<" = "<<this<<".CreateRangeOfIds()";
1905   return anObj._retn();
1906 }
1907
1908 BadOrientedVolume_ptr FilterManager_i::CreateBadOrientedVolume()
1909 {
1910   SMESH::BadOrientedVolume_i* aServant = new SMESH::BadOrientedVolume_i();
1911   SMESH::BadOrientedVolume_var anObj = aServant->_this();
1912   TPythonDump()<<aServant<<" = "<<this<<".CreateBadOrientedVolume()";
1913   return anObj._retn();
1914 }
1915
1916 LessThan_ptr FilterManager_i::CreateLessThan()
1917 {
1918   SMESH::LessThan_i* aServant = new SMESH::LessThan_i();
1919   SMESH::LessThan_var anObj = aServant->_this();
1920   TPythonDump()<<aServant<<" = "<<this<<".CreateLessThan()";
1921   return anObj._retn();
1922 }
1923
1924 MoreThan_ptr FilterManager_i::CreateMoreThan()
1925 {
1926   SMESH::MoreThan_i* aServant = new SMESH::MoreThan_i();
1927   SMESH::MoreThan_var anObj = aServant->_this();
1928   TPythonDump()<<aServant<<" = "<<this<<".CreateMoreThan()";
1929   return anObj._retn();
1930 }
1931
1932 EqualTo_ptr FilterManager_i::CreateEqualTo()
1933 {
1934   SMESH::EqualTo_i* aServant = new SMESH::EqualTo_i();
1935   SMESH::EqualTo_var anObj = aServant->_this();
1936   TPythonDump()<<aServant<<" = "<<this<<".CreateEqualTo()";
1937   return anObj._retn();
1938 }
1939
1940 LogicalNOT_ptr FilterManager_i::CreateLogicalNOT()
1941 {
1942   SMESH::LogicalNOT_i* aServant = new SMESH::LogicalNOT_i();
1943   SMESH::LogicalNOT_var anObj = aServant->_this();
1944   TPythonDump()<<aServant<<" = "<<this<<".CreateLogicalNOT()";
1945   return anObj._retn();
1946 }
1947
1948 LogicalAND_ptr FilterManager_i::CreateLogicalAND()
1949 {
1950   SMESH::LogicalAND_i* aServant = new SMESH::LogicalAND_i();
1951   SMESH::LogicalAND_var anObj = aServant->_this();
1952   TPythonDump()<<aServant<<" = "<<this<<".CreateLogicalAND()";
1953   return anObj._retn();
1954 }
1955
1956 LogicalOR_ptr FilterManager_i::CreateLogicalOR()
1957 {
1958   SMESH::LogicalOR_i* aServant = new SMESH::LogicalOR_i();
1959   SMESH::LogicalOR_var anObj = aServant->_this();
1960   TPythonDump()<<aServant<<" = "<<this<<".CreateLogicalOR()";
1961   return anObj._retn();
1962 }
1963
1964 LinearOrQuadratic_ptr FilterManager_i::CreateLinearOrQuadratic()
1965 {
1966   SMESH::LinearOrQuadratic_i* aServant = new SMESH::LinearOrQuadratic_i();
1967   SMESH::LinearOrQuadratic_var anObj = aServant->_this();
1968   TPythonDump()<<aServant<<" = "<<this<<".CreateLinearOrQuadratic()";
1969   return anObj._retn();
1970 }
1971
1972 GroupColor_ptr FilterManager_i::CreateGroupColor()
1973 {
1974   SMESH::GroupColor_i* aServant = new SMESH::GroupColor_i();
1975   SMESH::GroupColor_var anObj = aServant->_this();
1976   TPythonDump()<<aServant<<" = "<<this<<".CreateGroupColor()";
1977   return anObj._retn();
1978 }
1979
1980 ElemGeomType_ptr FilterManager_i::CreateElemGeomType()
1981 {
1982   SMESH::ElemGeomType_i* aServant = new SMESH::ElemGeomType_i();
1983   SMESH::ElemGeomType_var anObj = aServant->_this();
1984   TPythonDump()<<aServant<<" = "<<this<<".CreateElemGeomType()";
1985   return anObj._retn();
1986 }
1987
1988 Filter_ptr FilterManager_i::CreateFilter()
1989 {
1990   SMESH::Filter_i* aServant = new SMESH::Filter_i();
1991   SMESH::Filter_var anObj = aServant->_this();
1992   TPythonDump()<<aServant<<" = "<<this<<".CreateFilter()";
1993   return anObj._retn();
1994 }
1995
1996 FilterLibrary_ptr FilterManager_i::LoadLibrary( const char* aFileName )
1997 {
1998   SMESH::FilterLibrary_i* aServant = new SMESH::FilterLibrary_i( aFileName );
1999   SMESH::FilterLibrary_var anObj = aServant->_this();
2000   TPythonDump()<<aServant<<" = "<<this<<".LoadLibrary('"<<aFileName<<"')";
2001   return anObj._retn();
2002 }
2003
2004 FilterLibrary_ptr FilterManager_i::CreateLibrary()
2005 {
2006   SMESH::FilterLibrary_i* aServant = new SMESH::FilterLibrary_i();
2007   SMESH::FilterLibrary_var anObj = aServant->_this();
2008   TPythonDump()<<aServant<<" = "<<this<<".CreateLibrary()";
2009   return anObj._retn();
2010 }
2011
2012 CORBA::Boolean FilterManager_i::DeleteLibrary( const char* aFileName )
2013 {
2014   TPythonDump()<<this<<".DeleteLibrary("<<aFileName<<")";
2015   return remove( aFileName ) ? false : true;
2016 }
2017
2018 //=============================================================================
2019 /*!
2020  *  SMESH_Gen_i::CreateFilterManager
2021  *
2022  *  Create filter manager
2023  */
2024 //=============================================================================
2025
2026 SMESH::FilterManager_ptr SMESH_Gen_i::CreateFilterManager()
2027 {
2028   SMESH::FilterManager_i* aFilter = new SMESH::FilterManager_i();
2029   SMESH::FilterManager_var anObj = aFilter->_this();
2030   return anObj._retn();
2031 }
2032
2033
2034 /*
2035                               FILTER
2036 */
2037
2038 //=======================================================================
2039 // name    : Filter_i::Filter_i
2040 // Purpose : Constructor
2041 //=======================================================================
2042 Filter_i::Filter_i()
2043 : myPredicate( NULL )
2044 {}
2045
2046 //=======================================================================
2047 // name    : Filter_i::~Filter_i
2048 // Purpose : Destructor
2049 //=======================================================================
2050 Filter_i::~Filter_i()
2051 {
2052   if ( myPredicate )
2053     myPredicate->Destroy();
2054
2055   if(!CORBA::is_nil(myMesh))
2056     myMesh->Destroy();
2057
2058   //TPythonDump()<<this<<".Destroy()";
2059 }
2060
2061 //=======================================================================
2062 // name    : Filter_i::SetPredicate
2063 // Purpose : Set predicate
2064 //=======================================================================
2065 void Filter_i::SetPredicate( Predicate_ptr thePredicate )
2066 {
2067   if ( myPredicate )
2068     myPredicate->Destroy();
2069
2070   myPredicate = SMESH::GetPredicate(thePredicate);
2071
2072   if ( myPredicate )
2073   {
2074     myFilter.SetPredicate( myPredicate->GetPredicate() );
2075     myPredicate->Register();
2076     TPythonDump()<<this<<".SetPredicate("<<myPredicate<<")";
2077   }
2078 }
2079
2080 //=======================================================================
2081 // name    : Filter_i::GetElementType
2082 // Purpose : Get entity type
2083 //=======================================================================
2084 SMESH::ElementType Filter_i::GetElementType()
2085 {
2086   return myPredicate != 0 ? myPredicate->GetElementType() : SMESH::ALL;
2087 }
2088
2089 //=======================================================================
2090 // name    : Filter_i::SetMesh
2091 // Purpose : Set mesh
2092 //=======================================================================
2093 void
2094 Filter_i::
2095 SetMesh( SMESH_Mesh_ptr theMesh )
2096 {
2097   if(!CORBA::is_nil(theMesh))
2098     theMesh->Register();
2099
2100   if(!CORBA::is_nil(myMesh))
2101     myMesh->Destroy();
2102
2103   myMesh = SMESH_Mesh::_duplicate( theMesh );
2104   TPythonDump()<<this<<".SetMesh("<<theMesh<<")";
2105 }
2106
2107 SMESH::long_array*
2108 Filter_i::
2109 GetIDs()
2110 {
2111   return GetElementsId(myMesh);
2112 }
2113
2114 //=======================================================================
2115 // name    : Filter_i::GetElementsId
2116 // Purpose : Get ids of entities
2117 //=======================================================================
2118 void
2119 Filter_i::
2120 GetElementsId( Predicate_i* thePredicate,
2121                const SMDS_Mesh* theMesh,
2122                Controls::Filter::TIdSequence& theSequence )
2123 {
2124   if (thePredicate)
2125     Controls::Filter::GetElementsId(theMesh,thePredicate->GetPredicate(),theSequence);
2126 }
2127
2128 void
2129 Filter_i::
2130 GetElementsId( Predicate_i* thePredicate,
2131                SMESH_Mesh_ptr theMesh,
2132                Controls::Filter::TIdSequence& theSequence )
2133 {
2134   if (thePredicate)
2135     if(const SMDS_Mesh* aMesh = MeshPtr2SMDSMesh(theMesh))
2136       Controls::Filter::GetElementsId(aMesh,thePredicate->GetPredicate(),theSequence);
2137 }
2138
2139 SMESH::long_array*
2140 Filter_i::
2141 GetElementsId( SMESH_Mesh_ptr theMesh )
2142 {
2143   SMESH::long_array_var anArray = new SMESH::long_array;
2144   if(!CORBA::is_nil(theMesh) && myPredicate){
2145     Controls::Filter::TIdSequence aSequence;
2146     GetElementsId(myPredicate,theMesh,aSequence);
2147     long i = 0, iEnd = aSequence.size();
2148     anArray->length( iEnd );
2149     for ( ; i < iEnd; i++ )
2150       anArray[ i ] = aSequence[i];
2151   }
2152   return anArray._retn();
2153 }
2154
2155 template<class TElement, class TIterator, class TPredicate>
2156 static void collectMeshInfo(const TIterator& theItr,
2157                             TPredicate& thePred,
2158                             SMESH::long_array& theRes)
2159 {
2160   if (!theItr)
2161     return;
2162   while (theItr->more()) {
2163     const SMDS_MeshElement* anElem = theItr->next();
2164     if ( thePred->IsSatisfy( anElem->GetID() ) )
2165       theRes[ anElem->GetEntityType() ]++;
2166   }
2167 }
2168
2169 //=============================================================================
2170 /*!
2171  * \brief Returns statistic of mesh elements
2172  */
2173 //=============================================================================
2174 SMESH::long_array* ::Filter_i::GetMeshInfo()
2175 {
2176   SMESH::long_array_var aRes = new SMESH::long_array();
2177   aRes->length(SMESH::Entity_Last);
2178   for (int i = SMESH::Entity_Node; i < SMESH::Entity_Last; i++)
2179     aRes[i] = 0;
2180
2181   if(!CORBA::is_nil(myMesh) && myPredicate) {
2182     const SMDS_Mesh* aMesh = MeshPtr2SMDSMesh(myMesh);
2183     SMDS_ElemIteratorPtr it;
2184     switch( GetElementType() )
2185     {
2186   case SMDSAbs_Node:
2187     collectMeshInfo<const SMDS_MeshNode*>(aMesh->nodesIterator(),myPredicate,aRes);
2188     break;
2189   case SMDSAbs_Edge:
2190     collectMeshInfo<const SMDS_MeshElement*>(aMesh->edgesIterator(),myPredicate,aRes);
2191     break;
2192   case SMDSAbs_Face:
2193     collectMeshInfo<const SMDS_MeshElement*>(aMesh->facesIterator(),myPredicate,aRes);
2194     break;
2195   case SMDSAbs_Volume:
2196     collectMeshInfo<const SMDS_MeshElement*>(aMesh->volumesIterator(),myPredicate,aRes);
2197     break;
2198   case SMDSAbs_All:
2199   default:
2200     collectMeshInfo<const SMDS_MeshElement*>(aMesh->elementsIterator(),myPredicate,aRes);
2201     break;
2202     }
2203   }
2204
2205
2206   return aRes._retn();
2207 }
2208
2209 //================================================================================
2210 /*!
2211  * \brief Return GetElementType() within an array
2212  * Implement SMESH_IDSource interface
2213  */
2214 //================================================================================
2215
2216 SMESH::array_of_ElementType* Filter_i::GetTypes()
2217 {
2218   SMESH::array_of_ElementType_var types = new SMESH::array_of_ElementType;
2219   types->length( 1 );
2220   types[0] = GetElementType();
2221   return types._retn();
2222 }
2223
2224 //=======================================================================
2225 // name    : getCriteria
2226 // Purpose : Retrieve criterions from predicate
2227 //=======================================================================
2228 static inline bool getCriteria( Predicate_i*                thePred,
2229                                 SMESH::Filter::Criteria_out theCriteria )
2230 {
2231   int aFType = thePred->GetFunctorType();
2232
2233   switch ( aFType )
2234   {
2235   case FT_FreeBorders:
2236   case FT_FreeEdges:
2237   case FT_FreeFaces:
2238   case FT_LinearOrQuadratic:
2239   case FT_FreeNodes:
2240     {
2241       CORBA::ULong i = theCriteria->length();
2242       theCriteria->length( i + 1 );
2243
2244       theCriteria[ i ] = createCriterion();
2245
2246       theCriteria[ i ].Type = aFType;
2247       theCriteria[ i ].TypeOfElement = thePred->GetElementType();
2248       return true;
2249     }
2250   case FT_BelongToGeom:
2251     {
2252       BelongToGeom_i* aPred = dynamic_cast<BelongToGeom_i*>( thePred );
2253
2254       CORBA::ULong i = theCriteria->length();
2255       theCriteria->length( i + 1 );
2256
2257       theCriteria[ i ] = createCriterion();
2258
2259       theCriteria[ i ].Type          = FT_BelongToGeom;
2260       theCriteria[ i ].ThresholdStr  = aPred->GetShapeName();
2261       theCriteria[ i ].ThresholdID   = aPred->GetShapeID();
2262       theCriteria[ i ].TypeOfElement = aPred->GetElementType();
2263       theCriteria[ i ].Tolerance     = aPred->GetTolerance();
2264
2265       return true;
2266     }
2267   case FT_BelongToPlane:
2268   case FT_BelongToCylinder:
2269   case FT_BelongToGenSurface:
2270     {
2271       BelongToSurface_i* aPred = dynamic_cast<BelongToSurface_i*>( thePred );
2272
2273       CORBA::ULong i = theCriteria->length();
2274       theCriteria->length( i + 1 );
2275
2276       theCriteria[ i ] = createCriterion();
2277
2278       theCriteria[ i ].Type          = aFType;
2279       theCriteria[ i ].ThresholdStr  = aPred->GetShapeName();
2280       theCriteria[ i ].ThresholdID   = aPred->GetShapeID();
2281       theCriteria[ i ].TypeOfElement = aPred->GetElementType();
2282       theCriteria[ i ].Tolerance     = aPred->GetTolerance();
2283
2284       return true;
2285     }
2286    case FT_LyingOnGeom:
2287     {
2288       LyingOnGeom_i* aPred = dynamic_cast<LyingOnGeom_i*>( thePred );
2289
2290       CORBA::ULong i = theCriteria->length();
2291       theCriteria->length( i + 1 );
2292
2293       theCriteria[ i ] = createCriterion();
2294
2295       theCriteria[ i ].Type          = FT_LyingOnGeom;
2296       theCriteria[ i ].ThresholdStr  = aPred->GetShapeName();
2297       theCriteria[ i ].ThresholdID   = aPred->GetShapeID();
2298       theCriteria[ i ].TypeOfElement = aPred->GetElementType();
2299       theCriteria[ i ].Tolerance     = aPred->GetTolerance();
2300
2301       return true;
2302     }
2303   case FT_RangeOfIds:
2304     {
2305       RangeOfIds_i* aPred = dynamic_cast<RangeOfIds_i*>( thePred );
2306
2307       CORBA::ULong i = theCriteria->length();
2308       theCriteria->length( i + 1 );
2309
2310       theCriteria[ i ] = createCriterion();
2311
2312       theCriteria[ i ].Type          = FT_RangeOfIds;
2313       theCriteria[ i ].ThresholdStr  = aPred->GetRangeStr();
2314       theCriteria[ i ].TypeOfElement = aPred->GetElementType();
2315
2316       return true;
2317     }
2318   case FT_BadOrientedVolume:
2319     {
2320       BadOrientedVolume_i* aPred = dynamic_cast<BadOrientedVolume_i*>( thePred );
2321
2322       CORBA::ULong i = theCriteria->length();
2323       theCriteria->length( i + 1 );
2324
2325       theCriteria[ i ] = createCriterion();
2326
2327       theCriteria[ i ].Type          = FT_BadOrientedVolume;
2328       theCriteria[ i ].TypeOfElement = aPred->GetElementType();
2329
2330       return true;
2331     }
2332   case FT_LessThan:
2333   case FT_MoreThan:
2334   case FT_EqualTo:
2335     {
2336       Comparator_i* aCompar = dynamic_cast<Comparator_i*>( thePred );
2337
2338       CORBA::ULong i = theCriteria->length();
2339       theCriteria->length( i + 1 );
2340
2341       theCriteria[ i ] = createCriterion();
2342
2343       theCriteria[ i ].Type      = aCompar->GetNumFunctor_i()->GetFunctorType();
2344       theCriteria[ i ].Compare   = aFType;
2345       theCriteria[ i ].Threshold = aCompar->GetMargin();
2346       theCriteria[ i ].TypeOfElement = aCompar->GetElementType();
2347
2348       if ( aFType == FT_EqualTo )
2349       {
2350         EqualTo_i* aCompar = dynamic_cast<EqualTo_i*>( thePred );
2351         theCriteria[ i ].Tolerance = aCompar->GetTolerance();
2352       }
2353     }
2354     return true;
2355
2356   case FT_LogicalNOT:
2357     {
2358       Predicate_i* aPred = ( dynamic_cast<LogicalNOT_i*>( thePred ) )->GetPredicate_i();
2359       getCriteria( aPred, theCriteria );
2360       theCriteria[ theCriteria->length() - 1 ].UnaryOp = FT_LogicalNOT;
2361     }
2362     return true;
2363
2364   case FT_LogicalAND:
2365   case FT_LogicalOR:
2366     {
2367       Predicate_i* aPred1 = ( dynamic_cast<LogicalBinary_i*>( thePred ) )->GetPredicate1_i();
2368       Predicate_i* aPred2 = ( dynamic_cast<LogicalBinary_i*>( thePred ) )->GetPredicate2_i();
2369       if ( !getCriteria( aPred1, theCriteria ) )
2370         return false;
2371       theCriteria[ theCriteria->length() - 1 ].BinaryOp = aFType;
2372       return getCriteria( aPred2, theCriteria );
2373     }
2374   case FT_GroupColor:
2375     {
2376       CORBA::ULong i = theCriteria->length();
2377       theCriteria->length( i + 1 );
2378
2379       theCriteria[ i ] = createCriterion();
2380
2381       GroupColor_i* aPred = dynamic_cast<GroupColor_i*>( thePred );
2382       theCriteria[ i ].Type          = aFType;
2383       theCriteria[ i ].TypeOfElement = aPred->GetElementType();
2384       theCriteria[ i ].ThresholdStr  = aPred->GetColorStr();
2385
2386       return true;
2387     }
2388   case FT_ElemGeomType:
2389     {
2390       CORBA::ULong i = theCriteria->length();
2391       theCriteria->length( i + 1 );
2392
2393       theCriteria[ i ] = createCriterion();
2394
2395       ElemGeomType_i* aPred = dynamic_cast<ElemGeomType_i*>( thePred );
2396       theCriteria[ i ].Type          = aFType;
2397       theCriteria[ i ].TypeOfElement = aPred->GetElementType();
2398       theCriteria[ i ].Threshold     = (double)aPred->GetGeometryType();
2399       return true;
2400     }
2401
2402   case FT_Undefined:
2403     return false;
2404   default:
2405     return false;
2406   }
2407 }
2408
2409 //=======================================================================
2410 // name    : Filter_i::GetCriteria
2411 // Purpose : Retrieve criterions from predicate
2412 //=======================================================================
2413 CORBA::Boolean Filter_i::GetCriteria( SMESH::Filter::Criteria_out theCriteria )
2414 {
2415   theCriteria = new SMESH::Filter::Criteria;
2416   return myPredicate != 0 ? getCriteria( myPredicate, theCriteria ) : true;
2417 }
2418
2419 //=======================================================================
2420 // name    : Filter_i::SetCriteria
2421 // Purpose : Create new predicate and set criterions in it
2422 //=======================================================================
2423 CORBA::Boolean Filter_i::SetCriteria( const SMESH::Filter::Criteria& theCriteria )
2424 {
2425   if ( myPredicate != 0 )
2426     myPredicate->Destroy();
2427
2428   SMESH::FilterManager_i* aFilter = new SMESH::FilterManager_i();
2429   FilterManager_ptr aFilterMgr = aFilter->_this();
2430
2431   // CREATE two lists ( PREDICATES  and LOG OP )
2432
2433   // Criterion
2434   TPythonDump()<<"aCriteria = []";
2435   std::list<SMESH::Predicate_ptr> aPredicates;
2436   std::list<int>                  aBinaries;
2437   for ( int i = 0, n = theCriteria.length(); i < n; i++ )
2438   {
2439     int         aCriterion    = theCriteria[ i ].Type;
2440     int         aCompare      = theCriteria[ i ].Compare;
2441     double      aThreshold    = theCriteria[ i ].Threshold;
2442     const char* aThresholdStr = theCriteria[ i ].ThresholdStr;
2443     const char* aThresholdID  = theCriteria[ i ].ThresholdID;
2444     int         aUnary        = theCriteria[ i ].UnaryOp;
2445     int         aBinary       = theCriteria[ i ].BinaryOp;
2446     double      aTolerance    = theCriteria[ i ].Tolerance;
2447     ElementType aTypeOfElem   = theCriteria[ i ].TypeOfElement;
2448     long        aPrecision    = theCriteria[ i ].Precision;
2449
2450     {
2451       TPythonDump pd;
2452       pd << "aCriterion = SMESH.Filter.Criterion(" << aCriterion << "," << aCompare
2453          << "," << aThreshold << ",'" << aThresholdStr;
2454       if (aThresholdID && strlen(aThresholdID))
2455         //pd << "',salome.ObjectToID(" << aThresholdID
2456         pd << "','" << aThresholdID
2457            << "'," << aUnary << "," << aBinary << "," << aTolerance
2458            << "," << aTypeOfElem << "," << aPrecision << ")";
2459       else
2460         pd << "',''," << aUnary << "," << aBinary << "," << aTolerance
2461            << "," << aTypeOfElem << "," << aPrecision << ")";
2462     }
2463
2464     SMESH::Predicate_ptr aPredicate = SMESH::Predicate::_nil();
2465     SMESH::NumericalFunctor_ptr aFunctor = SMESH::NumericalFunctor::_nil();
2466
2467     switch ( aCriterion )
2468     {
2469       // Functors
2470
2471       case SMESH::FT_MultiConnection:
2472         aFunctor = aFilterMgr->CreateMultiConnection();
2473         break;
2474       case SMESH::FT_MultiConnection2D:
2475         aFunctor = aFilterMgr->CreateMultiConnection2D();
2476         break;
2477       case SMESH::FT_Length:
2478         aFunctor = aFilterMgr->CreateLength();
2479         break;
2480       case SMESH::FT_Length2D:
2481         aFunctor = aFilterMgr->CreateLength2D();
2482         break;
2483       case SMESH::FT_AspectRatio:
2484         aFunctor = aFilterMgr->CreateAspectRatio();
2485         break;
2486       case SMESH::FT_AspectRatio3D:
2487         aFunctor = aFilterMgr->CreateAspectRatio3D();
2488         break;
2489       case SMESH::FT_Warping:
2490         aFunctor = aFilterMgr->CreateWarping();
2491         break;
2492       case SMESH::FT_MinimumAngle:
2493         aFunctor = aFilterMgr->CreateMinimumAngle();
2494         break;
2495       case SMESH::FT_Taper:
2496         aFunctor = aFilterMgr->CreateTaper();
2497         break;
2498       case SMESH::FT_Skew:
2499         aFunctor = aFilterMgr->CreateSkew();
2500         break;
2501       case SMESH::FT_Area:
2502         aFunctor = aFilterMgr->CreateArea();
2503         break;
2504       case SMESH::FT_Volume3D:
2505         aFunctor = aFilterMgr->CreateVolume3D();
2506         break;
2507
2508       // Predicates
2509
2510       case SMESH::FT_FreeBorders:
2511         aPredicate = aFilterMgr->CreateFreeBorders();
2512         break;
2513       case SMESH::FT_FreeEdges:
2514         aPredicate = aFilterMgr->CreateFreeEdges();
2515         break;
2516       case SMESH::FT_FreeFaces:
2517         aPredicate = aFilterMgr->CreateFreeFaces();
2518         break;
2519       case SMESH::FT_FreeNodes:
2520         aPredicate = aFilterMgr->CreateFreeNodes();
2521         break;
2522       case SMESH::FT_BelongToGeom:
2523         {
2524           SMESH::BelongToGeom_ptr tmpPred = aFilterMgr->CreateBelongToGeom();
2525           tmpPred->SetElementType( aTypeOfElem );
2526           tmpPred->SetShape( aThresholdID, aThresholdStr );
2527           tmpPred->SetTolerance( aTolerance );
2528           aPredicate = tmpPred;
2529         }
2530         break;
2531       case SMESH::FT_BelongToPlane:
2532       case SMESH::FT_BelongToCylinder:
2533       case SMESH::FT_BelongToGenSurface:
2534         {
2535           SMESH::BelongToSurface_ptr tmpPred;
2536           switch ( aCriterion ) {
2537           case SMESH::FT_BelongToPlane:
2538             tmpPred = aFilterMgr->CreateBelongToPlane(); break;
2539           case SMESH::FT_BelongToCylinder:
2540             tmpPred = aFilterMgr->CreateBelongToCylinder(); break;
2541           default:
2542             tmpPred = aFilterMgr->CreateBelongToGenSurface();
2543           }
2544           tmpPred->SetShape( aThresholdID, aThresholdStr, aTypeOfElem );
2545           tmpPred->SetTolerance( aTolerance );
2546           aPredicate = tmpPred;
2547         }
2548         break;
2549       case SMESH::FT_LyingOnGeom:
2550         {
2551           SMESH::LyingOnGeom_ptr tmpPred = aFilterMgr->CreateLyingOnGeom();
2552           tmpPred->SetElementType( aTypeOfElem );
2553           tmpPred->SetShape( aThresholdID, aThresholdStr );
2554           tmpPred->SetTolerance( aTolerance );
2555           aPredicate = tmpPred;
2556         }
2557         break;
2558       case SMESH::FT_RangeOfIds:
2559         {
2560           SMESH::RangeOfIds_ptr tmpPred = aFilterMgr->CreateRangeOfIds();
2561           tmpPred->SetRangeStr( aThresholdStr );
2562           tmpPred->SetElementType( aTypeOfElem );
2563           aPredicate = tmpPred;
2564         }
2565         break;
2566       case SMESH::FT_BadOrientedVolume:
2567         {
2568           aPredicate = aFilterMgr->CreateBadOrientedVolume();
2569         }
2570         break;
2571       case SMESH::FT_LinearOrQuadratic:
2572         {
2573           SMESH::LinearOrQuadratic_ptr tmpPred = aFilterMgr->CreateLinearOrQuadratic();
2574           tmpPred->SetElementType( aTypeOfElem );
2575           aPredicate = tmpPred;
2576           break;
2577         }
2578       case SMESH::FT_GroupColor:
2579         {
2580           SMESH::GroupColor_ptr tmpPred = aFilterMgr->CreateGroupColor();
2581           tmpPred->SetElementType( aTypeOfElem );
2582           tmpPred->SetColorStr( aThresholdStr );
2583           aPredicate = tmpPred;
2584           break;
2585         }
2586       case SMESH::FT_ElemGeomType:
2587         {
2588           SMESH::ElemGeomType_ptr tmpPred = aFilterMgr->CreateElemGeomType();
2589           tmpPred->SetElementType( aTypeOfElem );
2590           tmpPred->SetGeometryType( (GeometryType)(int)(aThreshold + 0.5) );
2591           aPredicate = tmpPred;
2592           break;
2593         }
2594
2595       default:
2596         continue;
2597     }
2598
2599     // Comparator
2600     if ( !aFunctor->_is_nil() && aPredicate->_is_nil() )
2601     {
2602       SMESH::Comparator_ptr aComparator = SMESH::Comparator::_nil();
2603
2604       if ( aCompare == SMESH::FT_LessThan )
2605         aComparator = aFilterMgr->CreateLessThan();
2606       else if ( aCompare == SMESH::FT_MoreThan )
2607         aComparator = aFilterMgr->CreateMoreThan();
2608       else if ( aCompare == SMESH::FT_EqualTo )
2609         aComparator = aFilterMgr->CreateEqualTo();
2610       else
2611         continue;
2612
2613       aComparator->SetNumFunctor( aFunctor );
2614       aComparator->SetMargin( aThreshold );
2615
2616       if ( aCompare == FT_EqualTo )
2617       {
2618         SMESH::EqualTo_var anEqualTo = SMESH::EqualTo::_narrow( aComparator );
2619         anEqualTo->SetTolerance( aTolerance );
2620       }
2621
2622       aPredicate = aComparator;
2623
2624       aFunctor->SetPrecision( aPrecision );
2625     }
2626
2627     // Logical not
2628     if ( aUnary == FT_LogicalNOT )
2629     {
2630       SMESH::LogicalNOT_ptr aNotPred = aFilterMgr->CreateLogicalNOT();
2631       aNotPred->SetPredicate( aPredicate );
2632       aPredicate = aNotPred;
2633     }
2634
2635     // logical op
2636     aPredicates.push_back( aPredicate );
2637     aBinaries.push_back( aBinary );
2638     TPythonDump()<<"aCriteria.append(aCriterion)";
2639
2640   } // end of for
2641   TPythonDump()<<this<<".SetCriteria(aCriteria)";
2642
2643   // CREATE ONE PREDICATE FROM PREVIOUSLY CREATED MAP
2644
2645   // combine all "AND" operations
2646
2647   std::list<SMESH::Predicate_ptr> aResList;
2648
2649   std::list<SMESH::Predicate_ptr>::iterator aPredIter;
2650   std::list<int>::iterator                  aBinaryIter;
2651
2652   SMESH::Predicate_ptr aPrevPredicate = SMESH::Predicate::_nil();
2653   int aPrevBinary = SMESH::FT_Undefined;
2654
2655   for ( aPredIter = aPredicates.begin(), aBinaryIter = aBinaries.begin();
2656         aPredIter != aPredicates.end() && aBinaryIter != aBinaries.end();
2657         ++aPredIter, ++aBinaryIter )
2658   {
2659     int aCurrBinary = *aBinaryIter;
2660
2661     SMESH::Predicate_ptr aCurrPred = SMESH::Predicate::_nil();
2662
2663     if ( aPrevBinary == SMESH::FT_LogicalAND )
2664     {
2665
2666       SMESH::LogicalBinary_ptr aBinaryPred = aFilterMgr->CreateLogicalAND();
2667       aBinaryPred->SetPredicate1( aPrevPredicate );
2668       aBinaryPred->SetPredicate2( *aPredIter );
2669       aCurrPred = aBinaryPred;
2670     }
2671     else
2672       aCurrPred = *aPredIter;
2673
2674     if ( aCurrBinary != SMESH::FT_LogicalAND )
2675       aResList.push_back( aCurrPred );
2676
2677     aPrevPredicate = aCurrPred;
2678     aPrevBinary = aCurrBinary;
2679   }
2680
2681   // combine all "OR" operations
2682
2683   SMESH::Predicate_ptr aResPredicate = SMESH::Predicate::_nil();
2684
2685   if ( aResList.size() == 1 )
2686     aResPredicate = *aResList.begin();
2687   else if ( aResList.size() > 1 )
2688   {
2689     std::list<SMESH::Predicate_ptr>::iterator anIter = aResList.begin();
2690     aResPredicate = *anIter;
2691     anIter++;
2692     for ( ; anIter != aResList.end(); ++anIter )
2693     {
2694       SMESH::LogicalBinary_ptr aBinaryPred = aFilterMgr->CreateLogicalOR();
2695       aBinaryPred->SetPredicate1( aResPredicate );
2696       aBinaryPred->SetPredicate2( *anIter );
2697       aResPredicate = aBinaryPred;
2698     }
2699   }
2700
2701   SetPredicate( aResPredicate );
2702
2703   return !aResPredicate->_is_nil();
2704 }
2705
2706 //=======================================================================
2707 // name    : Filter_i::GetPredicate_i
2708 // Purpose : Get implementation of predicate
2709 //=======================================================================
2710 Predicate_i* Filter_i::GetPredicate_i()
2711 {
2712   return myPredicate;
2713 }
2714
2715 //=======================================================================
2716 // name    : Filter_i::GetPredicate
2717 // Purpose : Get predicate
2718 //=======================================================================
2719 Predicate_ptr Filter_i::GetPredicate()
2720 {
2721   if ( myPredicate == 0 )
2722     return SMESH::Predicate::_nil();
2723   else
2724   {
2725     SMESH::Predicate_var anObj = myPredicate->_this();
2726     return anObj._retn();
2727   }
2728 }
2729
2730 /*
2731                             FILTER LIBRARY
2732 */
2733
2734 #define ATTR_TYPE          "type"
2735 #define ATTR_COMPARE       "compare"
2736 #define ATTR_THRESHOLD     "threshold"
2737 #define ATTR_UNARY         "unary"
2738 #define ATTR_BINARY        "binary"
2739 #define ATTR_THRESHOLD_STR "threshold_str"
2740 #define ATTR_TOLERANCE     "tolerance"
2741 #define ATTR_ELEMENT_TYPE  "ElementType"
2742
2743 //=======================================================================
2744 // name    : toString
2745 // Purpose : Convert bool to LDOMString
2746 //=======================================================================
2747 static inline LDOMString toString( CORBA::Boolean val )
2748 {
2749   return val ? "logical not" : "";
2750 }
2751
2752 //=======================================================================
2753 // name    : toBool
2754 // Purpose : Convert LDOMString to bool
2755 //=======================================================================
2756 static inline bool toBool( const LDOMString& theStr )
2757 {
2758   return theStr.equals( "logical not" );
2759 }
2760
2761 //=======================================================================
2762 // name    : toString
2763 // Purpose : Convert double to LDOMString
2764 //=======================================================================
2765 static inline LDOMString toString( CORBA::Double val )
2766 {
2767   char a[ 255 ];
2768   sprintf( a, "%e", val );
2769   return LDOMString( a );
2770 }
2771
2772 //=======================================================================
2773 // name    : toDouble
2774 // Purpose : Convert LDOMString to double
2775 //=======================================================================
2776 static inline double toDouble( const LDOMString& theStr )
2777 {
2778   return atof( theStr.GetString() );
2779 }
2780
2781 //=======================================================================
2782 // name    : toString
2783 // Purpose : Convert functor type to LDOMString
2784 //=======================================================================
2785 static inline LDOMString toString( CORBA::Long theType )
2786 {
2787   switch ( theType )
2788   {
2789     case FT_AspectRatio     : return "Aspect ratio";
2790     case FT_Warping         : return "Warping";
2791     case FT_MinimumAngle    : return "Minimum angle";
2792     case FT_Taper           : return "Taper";
2793     case FT_Skew            : return "Skew";
2794     case FT_Area            : return "Area";
2795     case FT_Volume3D        : return "Volume3D";
2796     case FT_BelongToGeom    : return "Belong to Geom";
2797     case FT_BelongToPlane   : return "Belong to Plane";
2798     case FT_BelongToCylinder: return "Belong to Cylinder";
2799     case FT_BelongToGenSurface: return "Belong to Generic Surface";
2800     case FT_LyingOnGeom     : return "Lying on Geom";
2801     case FT_BadOrientedVolume: return "Bad Oriented Volume";
2802     case FT_RangeOfIds      : return "Range of IDs";
2803     case FT_FreeBorders     : return "Free borders";
2804     case FT_FreeEdges       : return "Free edges";
2805     case FT_FreeFaces       : return "Free faces";
2806     case FT_FreeNodes       : return "Free nodes";
2807     case FT_MultiConnection : return "Borders at multi-connections";
2808     case FT_MultiConnection2D: return "Borders at multi-connections 2D";
2809     case FT_Length          : return "Length";
2810     case FT_Length2D        : return "Length 2D";
2811     case FT_LessThan        : return "Less than";
2812     case FT_MoreThan        : return "More than";
2813     case FT_EqualTo         : return "Equal to";
2814     case FT_LogicalNOT      : return "Not";
2815     case FT_LogicalAND      : return "And";
2816     case FT_LogicalOR       : return "Or";
2817     case FT_GroupColor      : return "Color of Group";
2818     case FT_LinearOrQuadratic : return "Linear or Quadratic";
2819     case FT_ElemGeomType    : return "Element geomtry type";
2820     case FT_Undefined       : return "";
2821     default                 : return "";
2822   }
2823 }
2824
2825 //=======================================================================
2826 // name    : toFunctorType
2827 // Purpose : Convert LDOMString to functor type
2828 //=======================================================================
2829 static inline SMESH::FunctorType toFunctorType( const LDOMString& theStr )
2830 {
2831   if      ( theStr.equals( "Aspect ratio"                 ) ) return FT_AspectRatio;
2832   else if ( theStr.equals( "Warping"                      ) ) return FT_Warping;
2833   else if ( theStr.equals( "Minimum angle"                ) ) return FT_MinimumAngle;
2834   else if ( theStr.equals( "Taper"                        ) ) return FT_Taper;
2835   else if ( theStr.equals( "Skew"                         ) ) return FT_Skew;
2836   else if ( theStr.equals( "Area"                         ) ) return FT_Area;
2837   else if ( theStr.equals( "Volume3D"                     ) ) return FT_Volume3D;
2838   else if ( theStr.equals( "Belong to Geom"               ) ) return FT_BelongToGeom;
2839   else if ( theStr.equals( "Belong to Plane"              ) ) return FT_BelongToPlane;
2840   else if ( theStr.equals( "Belong to Cylinder"           ) ) return FT_BelongToCylinder;
2841   else if ( theStr.equals( "Belong to Generic Surface"    ) ) return FT_BelongToGenSurface;
2842   else if ( theStr.equals( "Lying on Geom"                ) ) return FT_LyingOnGeom;
2843   else if ( theStr.equals( "Free borders"                 ) ) return FT_FreeBorders;
2844   else if ( theStr.equals( "Free edges"                   ) ) return FT_FreeEdges;
2845   else if ( theStr.equals( "Free faces"                   ) ) return FT_FreeFaces;
2846   else if ( theStr.equals( "Free nodes"                   ) ) return FT_FreeNodes;
2847   else if ( theStr.equals( "Borders at multi-connections" ) ) return FT_MultiConnection;
2848   //  else if ( theStr.equals( "Borders at multi-connections 2D" ) ) return FT_MultiConnection2D;
2849   else if ( theStr.equals( "Length"                       ) ) return FT_Length;
2850   //  else if ( theStr.equals( "Length2D"                     ) ) return FT_Length2D;
2851   else if ( theStr.equals( "Range of IDs"                 ) ) return FT_RangeOfIds;
2852   else if ( theStr.equals( "Bad Oriented Volume"          ) ) return FT_BadOrientedVolume;
2853   else if ( theStr.equals( "Less than"                    ) ) return FT_LessThan;
2854   else if ( theStr.equals( "More than"                    ) ) return FT_MoreThan;
2855   else if ( theStr.equals( "Equal to"                     ) ) return FT_EqualTo;
2856   else if ( theStr.equals( "Not"                          ) ) return FT_LogicalNOT;
2857   else if ( theStr.equals( "And"                          ) ) return FT_LogicalAND;
2858   else if ( theStr.equals( "Or"                           ) ) return FT_LogicalOR;
2859   else if ( theStr.equals( "Color of Group"               ) ) return FT_GroupColor;
2860   else if ( theStr.equals( "Linear or Quadratic"          ) ) return FT_LinearOrQuadratic;
2861   else if ( theStr.equals( "Element geomtry type"         ) ) return FT_ElemGeomType;
2862   else if ( theStr.equals( ""                             ) ) return FT_Undefined;
2863   else  return FT_Undefined;
2864 }
2865
2866 //=======================================================================
2867 // name    : toFunctorType
2868 // Purpose : Convert LDOMString to value of ElementType enumeration
2869 //=======================================================================
2870 static inline SMESH::ElementType toElementType( const LDOMString& theStr )
2871 {
2872   if      ( theStr.equals( "NODE"   ) ) return SMESH::NODE;
2873   else if ( theStr.equals( "EDGE"   ) ) return SMESH::EDGE;
2874   else if ( theStr.equals( "FACE"   ) ) return SMESH::FACE;
2875   else if ( theStr.equals( "VOLUME" ) ) return SMESH::VOLUME;
2876   else                                  return SMESH::ALL;
2877 }
2878
2879 //=======================================================================
2880 // name    : toString
2881 // Purpose : Convert ElementType to string
2882 //=======================================================================
2883 static inline LDOMString toString( const SMESH::ElementType theType )
2884 {
2885   switch ( theType )
2886   {
2887     case SMESH::NODE   : return "NODE";
2888     case SMESH::EDGE   : return "EDGE";
2889     case SMESH::FACE   : return "FACE";
2890     case SMESH::VOLUME : return "VOLUME";
2891     case SMESH::ALL    : return "ALL";
2892     default            : return "";
2893   }
2894 }
2895
2896 //=======================================================================
2897 // name    : findFilter
2898 // Purpose : Find filter in document
2899 //=======================================================================
2900 static LDOM_Element findFilter( const char* theFilterName,
2901                                 const LDOM_Document& theDoc,
2902                                 LDOM_Node* theParent = 0 )
2903 {
2904   LDOM_Element aRootElement = theDoc.getDocumentElement();
2905   if ( aRootElement.isNull() || !aRootElement.hasChildNodes() )
2906     return LDOM_Element();
2907
2908   for ( LDOM_Node aTypeNode = aRootElement.getFirstChild();
2909         !aTypeNode.isNull(); aTypeNode = aTypeNode.getNextSibling() )
2910   {
2911     for ( LDOM_Node aFilter = aTypeNode.getFirstChild();
2912           !aFilter.isNull(); aFilter = aFilter.getNextSibling() )
2913     {
2914       LDOM_Element* anElem = ( LDOM_Element* )&aFilter;
2915       if ( anElem->getTagName().equals( LDOMString( "filter" ) ) &&
2916            anElem->getAttribute( "name" ).equals( LDOMString( theFilterName ) ) )
2917       {
2918         if ( theParent != 0  )
2919           *theParent = aTypeNode;
2920         return (LDOM_Element&)aFilter;
2921       }
2922     }
2923   }
2924   return LDOM_Element();
2925 }
2926
2927 //=======================================================================
2928 // name    : getSectionName
2929 // Purpose : Get name of section of filters
2930 //=======================================================================
2931 static const char* getSectionName( const ElementType theType )
2932 {
2933   switch ( theType )
2934   {
2935     case SMESH::NODE   : return "Filters for nodes";
2936     case SMESH::EDGE   : return "Filters for edges";
2937     case SMESH::FACE   : return "Filters for faces";
2938     case SMESH::VOLUME : return "Filters for volumes";
2939     case SMESH::ALL    : return "Filters for elements";
2940     default            : return "";
2941   }
2942 }
2943
2944 //=======================================================================
2945 // name    : getSection
2946 // Purpose : Create section for filters corresponding to the entity type
2947 //=======================================================================
2948 static LDOM_Node getSection( const ElementType theType,
2949                              LDOM_Document&    theDoc,
2950                              const bool        toCreate = false )
2951 {
2952   LDOM_Element aRootElement = theDoc.getDocumentElement();
2953   if ( aRootElement.isNull() )
2954     return LDOM_Node();
2955
2956   // Find section
2957   bool anExist = false;
2958   const char* aSectionName = getSectionName( theType );
2959   if ( strcmp( aSectionName, "" ) == 0 )
2960     return LDOM_Node();
2961
2962   LDOM_NodeList aSections = theDoc.getElementsByTagName( "section" );
2963   LDOM_Node aNode;
2964   for ( int i = 0, n = aSections.getLength(); i < n; i++ )
2965   {
2966     aNode = aSections.item( i );
2967     LDOM_Element& anItem = ( LDOM_Element& )aNode;
2968     if ( anItem.getAttribute( "name" ).equals( LDOMString( aSectionName ) ) )
2969     {
2970       anExist = true;
2971       break;
2972     }
2973   }
2974
2975   // Create new section if necessary
2976   if ( !anExist )
2977   {
2978     if ( toCreate )
2979     {
2980       LDOM_Element aNewItem = theDoc.createElement( "section" );
2981       aNewItem.setAttribute( "name", aSectionName );
2982       aRootElement.appendChild( aNewItem );
2983       return aNewItem;
2984     }
2985     else
2986       return LDOM_Node();
2987   }
2988   return
2989     aNode;
2990 }
2991
2992 //=======================================================================
2993 // name    : createFilterItem
2994 // Purpose : Create filter item or LDOM document
2995 //=======================================================================
2996 static LDOM_Element createFilterItem( const char*       theName,
2997                                       SMESH::Filter_ptr theFilter,
2998                                       LDOM_Document&    theDoc )
2999 {
3000   // create new filter in document
3001   LDOM_Element aFilterItem = theDoc.createElement( "filter" );
3002   aFilterItem.setAttribute( "name", theName );
3003
3004   // save filter criterions
3005   SMESH::Filter::Criteria_var aCriteria = new SMESH::Filter::Criteria;
3006
3007   if ( !theFilter->GetCriteria( aCriteria ) )
3008     return LDOM_Element();
3009
3010   for ( CORBA::ULong i = 0, n = aCriteria->length(); i < n; i++ )
3011   {
3012     LDOM_Element aCriterionItem = theDoc.createElement( "criterion" );
3013
3014     aCriterionItem.setAttribute( ATTR_TYPE         , toString(  aCriteria[ i ].Type) );
3015     aCriterionItem.setAttribute( ATTR_COMPARE      , toString(  aCriteria[ i ].Compare ) );
3016     aCriterionItem.setAttribute( ATTR_THRESHOLD    , toString(  aCriteria[ i ].Threshold ) );
3017     aCriterionItem.setAttribute( ATTR_UNARY        , toString(  aCriteria[ i ].UnaryOp ) );
3018     aCriterionItem.setAttribute( ATTR_BINARY       , toString(  aCriteria[ i ].BinaryOp ) );
3019
3020     aCriterionItem.setAttribute( ATTR_THRESHOLD_STR, (const char*)aCriteria[ i ].ThresholdStr );
3021     aCriterionItem.setAttribute( ATTR_TOLERANCE    , toString( aCriteria[ i ].Tolerance ) );
3022     aCriterionItem.setAttribute( ATTR_ELEMENT_TYPE ,
3023       toString( (SMESH::ElementType)aCriteria[ i ].TypeOfElement ) );
3024
3025     aFilterItem.appendChild( aCriterionItem );
3026   }
3027
3028   return aFilterItem;
3029 }
3030
3031 //=======================================================================
3032 // name    : FilterLibrary_i::FilterLibrary_i
3033 // Purpose : Constructor
3034 //=======================================================================
3035 FilterLibrary_i::FilterLibrary_i( const char* theFileName )
3036 {
3037   myFileName = strdup( theFileName );
3038   SMESH::FilterManager_i* aFilterMgr = new SMESH::FilterManager_i();
3039   myFilterMgr = aFilterMgr->_this();
3040
3041   LDOMParser aParser;
3042
3043   // Try to use existing library file
3044   bool anExists = false;
3045   if ( !aParser.parse( myFileName ) )
3046   {
3047     myDoc = aParser.getDocument();
3048     anExists = true;
3049   }
3050   // Create a new XML document if it doesn't exist
3051   else
3052     myDoc = LDOM_Document::createDocument( LDOMString() );
3053
3054   LDOM_Element aRootElement = myDoc.getDocumentElement();
3055   if ( aRootElement.isNull() )
3056   {
3057     // If the existing document is empty --> try to create a new one
3058     if ( anExists )
3059       myDoc = LDOM_Document::createDocument( LDOMString() );
3060   }
3061 }
3062
3063 //=======================================================================
3064 // name    : FilterLibrary_i::FilterLibrary_i
3065 // Purpose : Constructor
3066 //=======================================================================
3067 FilterLibrary_i::FilterLibrary_i()
3068 {
3069   myFileName = 0;
3070   SMESH::FilterManager_i* aFilter = new SMESH::FilterManager_i();
3071   myFilterMgr = aFilter->_this();
3072
3073   myDoc = LDOM_Document::createDocument( LDOMString() );
3074 }
3075
3076 FilterLibrary_i::~FilterLibrary_i()
3077 {
3078   delete myFileName;
3079   //TPythonDump()<<this<<".Destroy()";
3080 }
3081
3082 //=======================================================================
3083 // name    : FilterLibrary_i::Copy
3084 // Purpose : Create filter and initialize it with values from library
3085 //=======================================================================
3086 Filter_ptr FilterLibrary_i::Copy( const char* theFilterName )
3087 {
3088   Filter_ptr aRes = Filter::_nil();
3089   LDOM_Node aFilter = findFilter( theFilterName, myDoc );
3090
3091   if ( aFilter.isNull() )
3092     return aRes;
3093
3094   std::list<SMESH::Filter::Criterion> aCriteria;
3095
3096   for ( LDOM_Node aCritNode = aFilter.getFirstChild();
3097         !aCritNode.isNull() ; aCritNode = aCritNode.getNextSibling() )
3098   {
3099     LDOM_Element* aCrit = (LDOM_Element*)&aCritNode;
3100
3101     const char* aTypeStr      = aCrit->getAttribute( ATTR_TYPE          ).GetString();
3102     const char* aCompareStr   = aCrit->getAttribute( ATTR_COMPARE       ).GetString();
3103     const char* aUnaryStr     = aCrit->getAttribute( ATTR_UNARY         ).GetString();
3104     const char* aBinaryStr    = aCrit->getAttribute( ATTR_BINARY        ).GetString();
3105     const char* anElemTypeStr = aCrit->getAttribute( ATTR_ELEMENT_TYPE  ).GetString();
3106
3107     SMESH::Filter::Criterion aCriterion = createCriterion();
3108
3109     aCriterion.Type          = toFunctorType( aTypeStr );
3110     aCriterion.Compare       = toFunctorType( aCompareStr );
3111     aCriterion.UnaryOp       = toFunctorType( aUnaryStr );
3112     aCriterion.BinaryOp      = toFunctorType( aBinaryStr );
3113
3114     aCriterion.TypeOfElement = toElementType( anElemTypeStr );
3115
3116     LDOMString str = aCrit->getAttribute( ATTR_THRESHOLD );
3117     int val = 0;
3118     aCriterion.Threshold = str.Type() == LDOMBasicString::LDOM_Integer && str.GetInteger( val )
3119       ? val : atof( str.GetString() );
3120
3121     str = aCrit->getAttribute( ATTR_TOLERANCE );
3122     aCriterion.Tolerance = str.Type() == LDOMBasicString::LDOM_Integer && str.GetInteger( val )
3123       ? val : atof( str.GetString() );
3124
3125     str = aCrit->getAttribute( ATTR_THRESHOLD_STR );
3126     if ( str.Type() == LDOMBasicString::LDOM_Integer && str.GetInteger( val ) )
3127     {
3128       char a[ 255 ];
3129       sprintf( a, "%d", val );
3130       aCriterion.ThresholdStr = strdup( a );
3131     }
3132     else
3133       aCriterion.ThresholdStr = str.GetString();
3134
3135     aCriteria.push_back( aCriterion );
3136   }
3137
3138   SMESH::Filter::Criteria_var aCriteriaVar = new SMESH::Filter::Criteria;
3139   aCriteriaVar->length( aCriteria.size() );
3140
3141   CORBA::ULong i = 0;
3142   std::list<SMESH::Filter::Criterion>::iterator anIter = aCriteria.begin();
3143
3144   for( ; anIter != aCriteria.end(); ++anIter )
3145     aCriteriaVar[ i++ ] = *anIter;
3146
3147   aRes = myFilterMgr->CreateFilter();
3148   aRes->SetCriteria( aCriteriaVar.inout() );
3149
3150   TPythonDump()<<this<<".Copy('"<<theFilterName<<"')";
3151
3152   return aRes;
3153 }
3154
3155 //=======================================================================
3156 // name    : FilterLibrary_i::SetFileName
3157 // Purpose : Set file name for library
3158 //=======================================================================
3159 void FilterLibrary_i::SetFileName( const char* theFileName )
3160 {
3161   delete myFileName;
3162   myFileName = strdup( theFileName );
3163   TPythonDump()<<this<<".SetFileName('"<<theFileName<<"')";
3164 }
3165
3166 //=======================================================================
3167 // name    : FilterLibrary_i::GetFileName
3168 // Purpose : Get file name of library
3169 //=======================================================================
3170 char* FilterLibrary_i::GetFileName()
3171 {
3172   return CORBA::string_dup( myFileName );
3173 }
3174
3175 //=======================================================================
3176 // name    : FilterLibrary_i::Add
3177 // Purpose : Add new filter to library
3178 //=======================================================================
3179 CORBA::Boolean FilterLibrary_i::Add( const char* theFilterName, Filter_ptr theFilter )
3180 {
3181   // if filter already in library or entry filter is null do nothing
3182   LDOM_Node aFilterNode = findFilter( theFilterName, myDoc );
3183   if ( !aFilterNode.isNull() || theFilter->_is_nil() )
3184     return false;
3185
3186   // get section corresponding to the filter type
3187   ElementType anEntType = theFilter->GetElementType();
3188
3189   LDOM_Node aSection = getSection( anEntType, myDoc, true );
3190   if ( aSection.isNull() )
3191     return false;
3192
3193   // create filter item
3194   LDOM_Element aFilterItem = createFilterItem( theFilterName, theFilter, myDoc );
3195   if ( aFilterItem.isNull() )
3196     return false;
3197   else
3198   {
3199     aSection.appendChild( aFilterItem );
3200     if(Filter_i* aFilter = DownCast<Filter_i*>(theFilter))
3201       TPythonDump()<<this<<".Add('"<<theFilterName<<"',"<<aFilter<<")";
3202     return true;
3203   }
3204 }
3205
3206 //=======================================================================
3207 // name    : FilterLibrary_i::Add
3208 // Purpose : Add new filter to library
3209 //=======================================================================
3210 CORBA::Boolean FilterLibrary_i::AddEmpty( const char* theFilterName, ElementType theType )
3211 {
3212   // if filter already in library or entry filter is null do nothing
3213   LDOM_Node aFilterNode = findFilter( theFilterName, myDoc );
3214   if ( !aFilterNode.isNull() )
3215     return false;
3216
3217   LDOM_Node aSection = getSection( theType, myDoc, true );
3218   if ( aSection.isNull() )
3219     return false;
3220
3221   // create filter item
3222   Filter_var aFilter = myFilterMgr->CreateFilter();
3223
3224   LDOM_Element aFilterItem = createFilterItem( theFilterName, aFilter, myDoc );
3225   if ( aFilterItem.isNull() )
3226     return false;
3227   else
3228   {
3229     aSection.appendChild( aFilterItem );
3230     TPythonDump()<<this<<".AddEmpty('"<<theFilterName<<"',"<<theType<<")";
3231     return true;
3232   }
3233 }
3234
3235 //=======================================================================
3236 // name    : FilterLibrary_i::Delete
3237 // Purpose : Delete filter from library
3238 //=======================================================================
3239 CORBA::Boolean FilterLibrary_i::Delete ( const char* theFilterName )
3240 {
3241   LDOM_Node aParentNode;
3242   LDOM_Node aFilterNode = findFilter( theFilterName, myDoc, &aParentNode );
3243   if ( aFilterNode.isNull() || aParentNode.isNull() )
3244     return false;
3245
3246   aParentNode.removeChild( aFilterNode );
3247   TPythonDump()<<this<<".Delete('"<<theFilterName<<"')";
3248   return true;
3249 }
3250
3251 //=======================================================================
3252 // name      : FilterLibrary_i::Replace
3253 // Purpose   : Replace existing filter with entry filter.
3254 // IMPORTANT : If filter does not exist it is not created
3255 //=======================================================================
3256 CORBA::Boolean FilterLibrary_i::Replace( const char* theFilterName,
3257                                          const char* theNewName,
3258                                          Filter_ptr  theFilter )
3259 {
3260   LDOM_Element aFilterItem = findFilter( theFilterName, myDoc );
3261   if ( aFilterItem.isNull() || theFilter->_is_nil() )
3262     return false;
3263
3264   LDOM_Element aNewItem = createFilterItem( theNewName, theFilter, myDoc );
3265   if ( aNewItem.isNull() )
3266     return false;
3267   else
3268   {
3269     aFilterItem.ReplaceElement( aNewItem );
3270     if(Filter_i* aFilter = DownCast<Filter_i*>(theFilter))
3271       TPythonDump()<<this<<".Replace('"<<theFilterName<<"',"<<theNewName<<"',"<<aFilter<<")";
3272     return true;
3273   }
3274 }
3275
3276 //=======================================================================
3277 // name    : FilterLibrary_i::Save
3278 // Purpose : Save library on disk
3279 //=======================================================================
3280 CORBA::Boolean FilterLibrary_i::Save()
3281 {
3282   if ( myFileName == 0 || strlen( myFileName ) == 0 )
3283     return false;
3284
3285   FILE* aOutFile = fopen( myFileName, "wt" );
3286   if ( !aOutFile )
3287     return false;
3288
3289   LDOM_XmlWriter aWriter( aOutFile );
3290   aWriter.SetIndentation( 2 );
3291   aWriter << myDoc;
3292   fclose( aOutFile );
3293
3294   TPythonDump()<<this<<".Save()";
3295   return true;
3296 }
3297
3298 //=======================================================================
3299 // name    : FilterLibrary_i::SaveAs
3300 // Purpose : Save library on disk
3301 //=======================================================================
3302 CORBA::Boolean FilterLibrary_i::SaveAs( const char* aFileName )
3303 {
3304   myFileName = strdup ( aFileName );
3305   TPythonDump()<<this<<".SaveAs('"<<aFileName<<"')";
3306   return Save();
3307 }
3308
3309 //=======================================================================
3310 // name    : FilterLibrary_i::IsPresent
3311 // Purpose : Verify whether filter is in library
3312 //=======================================================================
3313 CORBA::Boolean FilterLibrary_i::IsPresent( const char* theFilterName )
3314 {
3315   return !findFilter( theFilterName, myDoc ).isNull();
3316 }
3317
3318 //=======================================================================
3319 // name    : FilterLibrary_i::NbFilters
3320 // Purpose : Return amount of filters in library
3321 //=======================================================================
3322 CORBA::Long FilterLibrary_i::NbFilters( ElementType theType )
3323 {
3324   string_array_var aNames = GetNames( theType );
3325   return aNames->length();
3326 }
3327
3328 //=======================================================================
3329 // name    : FilterLibrary_i::GetNames
3330 // Purpose : Get names of filters from library
3331 //=======================================================================
3332 string_array* FilterLibrary_i::GetNames( ElementType theType )
3333 {
3334   string_array_var anArray = new string_array;
3335   TColStd_SequenceOfHAsciiString aSeq;
3336
3337   LDOM_Node aSection = getSection( theType, myDoc, false );
3338
3339   if ( !aSection.isNull() )
3340   {
3341     for ( LDOM_Node aFilter = aSection.getFirstChild();
3342           !aFilter.isNull(); aFilter = aFilter.getNextSibling() )
3343     {
3344       LDOM_Element& anElem = ( LDOM_Element& )aFilter;
3345       aSeq.Append( new TCollection_HAsciiString(
3346          (Standard_CString)anElem.getAttribute( "name" ).GetString() ) );
3347     }
3348   }
3349
3350   anArray->length( aSeq.Length() );
3351   for ( int i = 1, n = aSeq.Length(); i <= n; i++ )
3352     anArray[ i - 1 ] = CORBA::string_dup( aSeq( i )->ToCString() );
3353
3354   return anArray._retn();
3355 }
3356
3357 //=======================================================================
3358 // name    : FilterLibrary_i::GetAllNames
3359 // Purpose : Get names of filters from library
3360 //=======================================================================
3361 string_array* FilterLibrary_i::GetAllNames()
3362 {
3363   string_array_var aResArray = new string_array;
3364   for ( int type = SMESH::ALL; type <= SMESH::VOLUME; type++ )
3365   {
3366     SMESH::string_array_var aNames = GetNames( (SMESH::ElementType)type );
3367
3368     int aPrevLength = aResArray->length();
3369     aResArray->length( aPrevLength + aNames->length() );
3370     for ( int i = 0, n = aNames->length(); i < n; i++ )
3371       aResArray[ aPrevLength + i ] = aNames[ i ];
3372   }
3373
3374   return aResArray._retn();
3375 }