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