src/SMESH_I/SMESH_Filter_i.cxx

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