src/SMESH_I/SMESH_Filter_i.cxx

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