src/SMESH_I/SMESH_Filter_i.cxx

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