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