X-Git-Url: http://git.salome-platform.org/gitweb/?a=blobdiff_plain;f=src%2FSMESH_I%2FSMESH_Filter_i.cxx;h=ae97ce1c664fd6f444e9127a111055f543ae2a02;hb=HEAD;hp=724d602dd7cfc3b2864f22f8012b8bf4369f1ecc;hpb=c3bf92bd87b770fd81631a3853f7f5bb1ac6a4e8;p=modules%2Fsmesh.git diff --git a/src/SMESH_I/SMESH_Filter_i.cxx b/src/SMESH_I/SMESH_Filter_i.cxx index 724d602dd..5448c6699 100644 --- a/src/SMESH_I/SMESH_Filter_i.cxx +++ b/src/SMESH_I/SMESH_Filter_i.cxx @@ -1,914 +1,4342 @@ -// SMESH SMESH_I : idl implementation based on 'SMESH' unit's calsses +// Copyright (C) 2007-2024 CEA, EDF, OPEN CASCADE // -// Copyright (C) 2003 OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN, -// CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS -// -// This library is free software; you can redistribute it and/or -// modify it under the terms of the GNU Lesser General Public -// License as published by the Free Software Foundation; either -// version 2.1 of the License. -// -// This library is distributed in the hope that it will be useful, -// but WITHOUT ANY WARRANTY; without even the implied warranty of -// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU -// Lesser General Public License for more details. -// -// You should have received a copy of the GNU Lesser General Public -// License along with this library; if not, write to the Free Software -// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA -// -// See http://www.opencascade.org/SALOME/ or email : webmaster.salome@opencascade.org +// Copyright (C) 2003-2007 OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN, +// CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS // +// This library is free software; you can redistribute it and/or +// modify it under the terms of the GNU Lesser General Public +// License as published by the Free Software Foundation; either +// version 2.1 of the License, or (at your option) any later version. // +// This library is distributed in the hope that it will be useful, +// but WITHOUT ANY WARRANTY; without even the implied warranty of +// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU +// Lesser General Public License for more details. // +// You should have received a copy of the GNU Lesser General Public +// License along with this library; if not, write to the Free Software +// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA +// +// See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com +// + +// SMESH SMESH_I : idl implementation based on 'SMESH' unit's classes // File : SMESH_Filter_i.cxx // Author : Alexey Petrov, OCC // Module : SMESH - #include "SMESH_Filter_i.hxx" -#include "SMDS_Iterator.hxx" +#include "SMDS_ElemIterator.hxx" +#include "SMDS_Mesh.hxx" #include "SMDS_MeshElement.hxx" #include "SMDS_MeshNode.hxx" -#include "SMDSAbs_ElementType.hxx" -#include "SMESH_Gen_i.hxx" +#include "SMESHDS_GroupBase.hxx" #include "SMESHDS_Mesh.hxx" - -#include -#include -#include +#include "SMESH_Gen_i.hxx" +#include "SMESH_Group_i.hxx" +#include "SMESH_PythonDump.hxx" + +#include +#include + +#include +#include +#include +#include +#include +#include +#include +#include +#include #include -#include -#include -#include #include +#include +#include +#include +#include +#include +#include +#include +#include + +//using namespace SMESH; +//using namespace SMESH::Controls; + + +namespace SMESH +{ + Predicate_i* + GetPredicate( Predicate_ptr thePredicate ) + { + return DownCast(thePredicate); + } +} /* - AUXILIARY METHODS + AUXILIARY METHODS */ -static inline double getAngle( const gp_XYZ& P1, const gp_XYZ& P2, const gp_XYZ& P3 ) +inline +const SMDS_Mesh* +MeshPtr2SMDSMesh( SMESH::SMESH_Mesh_ptr theMesh ) { - return gp_Vec( P1 - P2 ).Angle( gp_Vec( P3 - P2 ) ); + SMESH_Mesh_i* anImplPtr = SMESH::DownCast(theMesh); + return anImplPtr ? anImplPtr->GetImpl().GetMeshDS() : 0; } -static inline double getArea( const gp_XYZ& P1, const gp_XYZ& P2, const gp_XYZ& P3 ) +inline +SMESH::long_array* +toArray( const TColStd_ListOfInteger& aList ) { - gp_Vec aVec1( P2 - P1 ); - gp_Vec aVec2( P3 - P1 ); - return ( aVec1 ^ aVec2 ).Magnitude() * 0.5; + SMESH::long_array_var anArray = new SMESH::long_array; + anArray->length( aList.Extent() ); + TColStd_ListIteratorOfListOfInteger anIter( aList ); + int i = 0; + for( ; anIter.More(); anIter.Next() ) + anArray[ i++ ] = anIter.Value(); + + return anArray._retn(); } -static inline double getArea( const gp_Pnt& P1, const gp_Pnt& P2, const gp_Pnt& P3 ) +inline +SMESH::double_array* +toArray( const TColStd_ListOfReal& aList ) { - return getArea( P1.XYZ(), P2.XYZ(), P3.XYZ() ); + SMESH::double_array_var anArray = new SMESH::double_array; + anArray->length( aList.Extent() ); + TColStd_ListIteratorOfListOfReal anIter( aList ); + int i = 0; + for( ; anIter.More(); anIter.Next() ) + anArray[ i++ ] = anIter.Value(); + + return anArray._retn(); } -static inline double getDistance( const gp_XYZ& P1, const gp_XYZ& P2 ) +static SMESH::Filter::Criterion createCriterion() { - double aDist = gp_Pnt( P1 ).Distance( gp_Pnt( P2 ) ); - return aDist; + SMESH::Filter::Criterion aCriterion; + + aCriterion.Type = SMESH::FT_Undefined; + aCriterion.Compare = SMESH::FT_Undefined; + aCriterion.Threshold = 0; + aCriterion.UnaryOp = SMESH::FT_Undefined; + aCriterion.BinaryOp = SMESH::FT_Undefined; + aCriterion.ThresholdStr = ""; + aCriterion.ThresholdID = ""; + aCriterion.Tolerance = Precision::Confusion(); + aCriterion.TypeOfElement = SMESH::ALL; + aCriterion.Precision = -1; + + return aCriterion; } -static int getNbMultiConnection( SMESHDS_Mesh* theMesh, const int theId ) +static TopoDS_Shape getShapeByName( const char* theName ) { - if ( theMesh == 0 ) - return 0; - - const SMDS_MeshElement* anEdge = theMesh->FindElement( theId ); - if ( anEdge == 0 || anEdge->GetType() != SMDSAbs_Edge || anEdge->NbNodes() != 2 ) - return 0; - - TColStd_MapOfInteger aMap; - - int aResult = 0; - SMDS_ElemIteratorPtr anIter = anEdge->nodesIterator(); - if ( anIter != 0 ) + if ( theName != 0 ) { - while( anIter->more() ) + SMESH_Gen_i* aSMESHGen = SMESH_Gen_i::GetSMESHGen(); + SALOMEDS::Study::ListOfSObject_var aList = SMESH_Gen_i::GetSMESHGen()->getStudyServant()->FindObjectByName( theName, "GEOM" ); + if ( aList->length() == 0 ) + aList = SMESH_Gen_i::GetSMESHGen()->getStudyServant()->FindObjectByName( theName, "SHAPERSTUDY" ); + if ( aList->length() > 0 ) { - const SMDS_MeshNode* aNode = (SMDS_MeshNode*)anIter->next(); - if ( aNode == 0 ) - return 0; - SMDS_ElemIteratorPtr anElemIter = aNode->GetInverseElementIterator(); - while( anElemIter->more() ) - { - const SMDS_MeshElement* anElem = anElemIter->next(); - if ( anElem != 0 && anElem->GetType() != SMDSAbs_Edge ) - { - int anId = anElem->GetID(); - - if ( anIter->more() ) // i.e. first node - aMap.Add( anId ); - else if ( aMap.Contains( anId ) ) - aResult++; - } - } -// delete anElemIter; + CORBA::Object_var anObj = aList[ 0 ]->GetObject(); + GEOM::GEOM_Object_var aGeomObj = GEOM::GEOM_Object::_narrow( anObj ); + TopoDS_Shape shape = aSMESHGen->GeomObjectToShape( aGeomObj ); + SALOME::UnRegister( aList ); // UnRegister() objects in aList + return shape; } -// delete anIter; } + return TopoDS_Shape(); +} - return aResult; +static TopoDS_Shape getShapeByID (const char* theID) +{ + if ( theID && strlen( theID ) > 0 ) { + SMESH_Gen_i* aSMESHGen = SMESH_Gen_i::GetSMESHGen(); + SALOMEDS::SObject_wrap aSObj = SMESH_Gen_i::GetSMESHGen()->getStudyServant()->FindObjectID(theID); + if ( !aSObj->_is_nil() ) { + CORBA::Object_var obj = aSObj->GetObject(); + GEOM::GEOM_Object_var aGeomObj = GEOM::GEOM_Object::_narrow(obj); + return aSMESHGen->GeomObjectToShape( aGeomObj ); + } + } + return TopoDS_Shape(); } -using namespace std; -using namespace SMESH; +// static std::string getShapeNameByID (const char* theID) +// { +// if ( theID && strlen( theID ) > 0 ) { +// SALOMEDS::SObject_wrap aSObj = SMESH_Gen_i::getStudyServant()->FindObjectID(theID); +// if ( !aSObj->_is_nil() ) { +// CORBA::String_var name = aSObj->GetName(); +// return name.in(); +// } +// } +// return ""; +// } /* FUNCTORS */ -/* - Class : NumericalFunctor_i - Description : Base class for numerical functors -*/ +//============================================================================= +/*! + * SMESH_Gen_i::CreateFilterManager + * + * Create filter manager + */ + //============================================================================= -NumericalFunctor_i::NumericalFunctor_i() -: SALOME::GenericObj_i( SMESH_Gen_i::GetPOA() ) +SMESH::FilterManager_ptr SMESH_Gen_i::CreateFilterManager() { - myMesh = 0; - SMESH_Gen_i::GetPOA()->activate_object( this ); + SMESH::FilterManager_i* aFilter = new SMESH::FilterManager_i(); + SMESH::FilterManager_var anObj = aFilter->_this(); + return anObj._retn(); } -void NumericalFunctor_i::SetMesh( SMESH_Mesh_ptr theMesh ) -{ - SMESH_Mesh_i* anImplPtr = - dynamic_cast( SMESH_Gen_i::GetServant( theMesh ).in() ); - myMesh = anImplPtr ? anImplPtr->GetImpl().GetMeshDS() : 0; -} +namespace SMESH { + /* + Class : Functor_i + Description : An abstract class for all functors + */ + Functor_i::Functor_i(): + SALOME::GenericObj_i( SMESH_Gen_i::GetPOA() ) + { + //Base class Salome_GenericObject do it inmplicitly by overriding PortableServer::POA_ptr _default_POA() method + //PortableServer::ObjectId_var anObjectId = + // SMESH_Gen_i::GetPOA()->activate_object( this ); + } -bool NumericalFunctor_i::getPoints( const int theId, - TColgp_SequenceOfXYZ& theRes ) const -{ - theRes.Clear(); + Functor_i::~Functor_i() + { + //TPythonDump()<SetMesh( MeshPtr2SMDSMesh( theMesh ) ); + TPythonDump()<GetType(); + } - if ( myMesh == 0 ) - return false; - // Get nodes of the face - const SMDS_MeshElement* anElem = myMesh->FindElement( theId ); - if ( anElem == 0 || anElem->GetType() != GetType() ) - return false; + /* + Class : NumericalFunctor_i + Description : Base class for numerical functors + */ + CORBA::Double NumericalFunctor_i::GetValue( SMESH::smIdType theId ) + { + return myNumericalFunctorPtr->GetValue( theId ); + } - int nbNodes = anElem->NbNodes(); + CORBA::Boolean NumericalFunctor_i::IsApplicable( SMESH::smIdType theElementId ) + { + return myNumericalFunctorPtr->IsApplicable( theElementId ); + } - SMDS_ElemIteratorPtr anIter = anElem->nodesIterator(); - if ( anIter != 0 ) + SMESH::Histogram* NumericalFunctor_i::GetHistogram(CORBA::Short nbIntervals, CORBA::Boolean isLogarithmic) { - while( anIter->more() ) + std::vector nbEvents; + std::vector funValues; + std::vector<::smIdType> elements; + myNumericalFunctorPtr->GetHistogram(nbIntervals, nbEvents, funValues ,elements, 0, isLogarithmic); + + SMESH::Histogram_var histogram = new SMESH::Histogram; + + nbIntervals = CORBA::Short( Min( int( nbEvents.size()), + int( funValues.size() - 1 ))); + if ( nbIntervals > 0 ) { - const SMDS_MeshNode* aNode = (SMDS_MeshNode*)anIter->next(); - if ( aNode != 0 ) - theRes.Append( gp_XYZ( aNode->X(), aNode->Y(), aNode->Z() ) ); + histogram->length( nbIntervals ); + for ( int i = 0; i < nbIntervals; ++i ) + { + HistogramRectangle& rect = histogram[i]; + rect.nbEvents = nbEvents[i]; + rect.min = funValues[i]; + rect.max = funValues[i+1]; + } } - -// delete anIter; + return histogram._retn(); } - return true; -} + SMESH::Histogram* NumericalFunctor_i::GetLocalHistogram(CORBA::Short nbIntervals, + CORBA::Boolean isLogarithmic, + SMESH::SMESH_IDSource_ptr object) + { + SMESH::Histogram_var histogram = new SMESH::Histogram; + std::vector nbEvents; + std::vector funValues; + std::vector<::smIdType> elements; -/* - Class : SMESH_MinimumAngleFunct - Description : Functor for calculation of minimum angle -*/ + SMDS_ElemIteratorPtr elemIt; + if ( SMESH::DownCast< SMESH_GroupOnFilter_i* >( object ) || + SMESH::DownCast< SMESH::Filter_i* >( object )) + { + elemIt = SMESH_Mesh_i::GetElements( object, GetElementType() ); + if ( !elemIt ) return histogram._retn(); + } + else + { + SMESH::SMESH_Mesh_var mesh = object->GetMesh(); + SMESH::smIdType_array_var objNbElems = object->GetNbElementsByType(); + SMESH::smIdType_array_var meshNbElems = mesh-> GetNbElementsByType(); + if ( meshNbElems[ GetElementType() ] != + objNbElems [ GetElementType() ] ) + { + elements.reserve( objNbElems[ GetElementType() ]); + elemIt = SMESH_Mesh_i::GetElements( object, GetElementType() ); + } + } + if ( elemIt ) + { + while ( elemIt->more() ) + elements.push_back( elemIt->next()->GetID() ); + if ( elements.empty() ) return histogram._retn(); + } -CORBA::Double MinimumAngle_i::GetValue( CORBA::Long theId ) -{ - TColgp_SequenceOfXYZ P; - if ( !getPoints( theId, P ) || P.Length() != 3 && P.Length() != 4 ) - return 0; + myNumericalFunctorPtr->GetHistogram(nbIntervals,nbEvents,funValues, elements, 0, isLogarithmic); - double aMin; + nbIntervals = CORBA::Short( Min( int( nbEvents.size()), + int( funValues.size() - 1 ))); + if ( nbIntervals > 0 ) + { + histogram->length( nbIntervals ); + for ( int i = 0; i < nbIntervals; ++i ) + { + HistogramRectangle& rect = histogram[i]; + rect.nbEvents = nbEvents[i]; + rect.min = funValues[i]; + rect.max = funValues[i+1]; + } + } + return histogram._retn(); + } - if ( P.Length() == 3 ) + void NumericalFunctor_i::SetPrecision( CORBA::Long thePrecision ) { - double A0 = getAngle( P( 3 ), P( 1 ), P( 2 ) ); - double A1 = getAngle( P( 1 ), P( 2 ), P( 3 ) ); - double A2 = getAngle( P( 2 ), P( 3 ), P( 1 ) ); + myNumericalFunctorPtr->SetPrecision( thePrecision ); + TPythonDump()<GetPrecision(); } - else + + Controls::NumericalFunctorPtr NumericalFunctor_i::GetNumericalFunctor() { - double A0 = getAngle( P( 4 ), P( 1 ), P( 2 ) ); - double A1 = getAngle( P( 1 ), P( 2 ), P( 3 ) ); - double A2 = getAngle( P( 2 ), P( 3 ), P( 4 ) ); - double A3 = getAngle( P( 3 ), P( 4 ), P( 1 ) ); - - aMin = Min( Min( A0, A1 ), Min( A2, A3 ) ); + return myNumericalFunctorPtr; } - - return aMin * 180 / PI; -} -int MinimumAngle_i::GetType() const -{ - return SMDSAbs_Face; -} -/* - Class : AspectRatio_i - Description : Functor for calculating aspect ratio -*/ + /* + Class : SMESH_MinimumAngle + Description : Functor for calculation of minimum angle + */ + MinimumAngle_i::MinimumAngle_i() + { + myNumericalFunctorPtr.reset( new Controls::MinimumAngle() ); + myFunctorPtr = myNumericalFunctorPtr; + } -CORBA::Double AspectRatio_i::GetValue( CORBA::Long theId ) -{ - TColgp_SequenceOfXYZ P; - if ( !getPoints( theId, P ) || P.Length() != 3 && P.Length() != 4 ) - return 0; + FunctorType MinimumAngle_i::GetFunctorType() + { + return SMESH::FT_MinimumAngle; + } - int nbNodes = P.Length(); - // Compute lengths of the sides + /* + Class : AspectRatio + Description : Functor for calculating aspect ratio + */ + AspectRatio_i::AspectRatio_i() + { + myNumericalFunctorPtr.reset( new Controls::AspectRatio() ); + myFunctorPtr = myNumericalFunctorPtr; + } - double aLen[ nbNodes ]; - for ( int i = 0; i < nbNodes - 1; i++ ) - aLen[ i ] = getDistance( P( i + 1 ), P( i + 2 ) ); - aLen[ nbNodes - 1 ] = getDistance( P( 1 ), P( nbNodes ) ); + FunctorType AspectRatio_i::GetFunctorType() + { + return SMESH::FT_AspectRatio; + } - // Compute aspect ratio - if ( nbNodes == 3 ) + /* + Class : AspectRatio3D + Description : Functor for calculating aspect ratio 3D + */ + AspectRatio3D_i::AspectRatio3D_i() { - double aMaxLen = Max( aLen[ 0 ], Max( aLen[ 1 ], aLen[ 2 ] ) ); - double anArea = getArea( P( 1 ), P( 2 ), P( 3 ) ); - static double aCoef = sqrt( 3. ) / 4; - - return anArea != 0 ? aCoef * aMaxLen * aMaxLen / anArea : 0; + myNumericalFunctorPtr.reset( new Controls::AspectRatio3D() ); + myFunctorPtr = myNumericalFunctorPtr; } - else + + FunctorType AspectRatio3D_i::GetFunctorType() { - double aMaxLen = Max( Max( aLen[ 0 ], aLen[ 1 ] ), Max( aLen[ 2 ], aLen[ 3 ] ) ); - double aMinLen = Min( Min( aLen[ 0 ], aLen[ 1 ] ), Min( aLen[ 2 ], aLen[ 3 ] ) ); - - return aMinLen != 0 ? aMaxLen / aMinLen : 0; + return SMESH::FT_AspectRatio3D; } -} -int AspectRatio_i::GetType() const -{ - return SMDSAbs_Face; -} -/* - Class : Warping_i - Description : Functor for calculating warping -*/ + /* + Class : Warping_i + Description : Functor for calculating warping + */ + Warping_i::Warping_i() + { + myNumericalFunctorPtr.reset( new Controls::Warping() ); + myFunctorPtr = myNumericalFunctorPtr; + } -CORBA::Double Warping_i::GetValue( CORBA::Long theId ) -{ - TColgp_SequenceOfXYZ P; - if ( !getPoints( theId, P ) || P.Length() != 4 ) - return 0; + FunctorType Warping_i::GetFunctorType() + { + return SMESH::FT_Warping; + } - gp_XYZ G = ( P( 1 ) + P( 2 ) + P( 3 ) + P( 4 ) ) / 4; + /* + Class : Warping3D_i + Description : Functor for calculating warping + */ + Warping3D_i::Warping3D_i() + { + myNumericalFunctorPtr.reset(new Controls::Warping3D()); + myFunctorPtr = myNumericalFunctorPtr; + } - double A1 = ComputeA( P( 1 ), P( 2 ), P( 3 ), G ); - double A2 = ComputeA( P( 2 ), P( 3 ), P( 4 ), G ); - double A3 = ComputeA( P( 3 ), P( 4 ), P( 1 ), G ); - double A4 = ComputeA( P( 4 ), P( 1 ), P( 2 ), G ); + FunctorType Warping3D_i::GetFunctorType() + { + return SMESH::FT_Warping3D; + } - return Max( Max( A1, A2 ), Max( A3, A4 ) ); -} + /* + Class : Taper_i + Description : Functor for calculating taper + */ + Taper_i::Taper_i() + { + myNumericalFunctorPtr.reset( new Controls::Taper() ); + myFunctorPtr = myNumericalFunctorPtr; + } -double Warping_i::ComputeA( const gp_XYZ& thePnt1, - const gp_XYZ& thePnt2, - const gp_XYZ& thePnt3, - const gp_XYZ& theG ) const -{ - double aLen1 = gp_Pnt( thePnt1 ).Distance( gp_Pnt( thePnt2 ) ); - double aLen2 = gp_Pnt( thePnt2 ).Distance( gp_Pnt( thePnt3 ) ); - double L = Min( aLen1, aLen2 ) * 0.5; + FunctorType Taper_i::GetFunctorType() + { + return SMESH::FT_Taper; + } - gp_XYZ GI = ( thePnt2 - thePnt1 ) / 2. - theG; - gp_XYZ GJ = ( thePnt3 - thePnt2 ) / 2. - theG; - gp_XYZ N = GI.Crossed( GJ ); - N.Normalize(); + /* + Class : Skew_i + Description : Functor for calculating skew in degrees + */ + Skew_i::Skew_i() + { + myNumericalFunctorPtr.reset( new Controls::Skew() ); + myFunctorPtr = myNumericalFunctorPtr; + } - double H = gp_Vec( thePnt2 - theG ).Dot( gp_Vec( N ) ); - return asin( fabs( H / L ) ) * 180 / PI; -} + FunctorType Skew_i::GetFunctorType() + { + return SMESH::FT_Skew; + } -int Warping_i::GetType() const -{ - return SMDSAbs_Face; -} + /* + Class : Area_i + Description : Functor for calculating area + */ + Area_i::Area_i() + { + myNumericalFunctorPtr.reset( new Controls::Area() ); + myFunctorPtr = myNumericalFunctorPtr; + } -/* - Class : Taper_i - Description : Functor for calculating taper -*/ + FunctorType Area_i::GetFunctorType() + { + return SMESH::FT_Area; + } -CORBA::Double Taper_i::GetValue( CORBA::Long theId ) -{ - TColgp_SequenceOfXYZ P; - if ( !getPoints( theId, P ) || P.Length() != 4 ) - return 0; + /* + Class : Volume3D_i + Description : Functor for calculating volume of 3D element + */ + Volume3D_i::Volume3D_i() + { + myNumericalFunctorPtr.reset( new Controls::Volume() ); + myFunctorPtr = myNumericalFunctorPtr; + } - // Compute taper - double J1 = getArea( P( 4 ), P( 1 ), P( 2 ) ) / 2; - double J2 = getArea( P( 3 ), P( 1 ), P( 2 ) ) / 2; - double J3 = getArea( P( 2 ), P( 3 ), P( 4 ) ) / 2; - double J4 = getArea( P( 3 ), P( 4 ), P( 1 ) ) / 2; + FunctorType Volume3D_i::GetFunctorType() + { + return SMESH::FT_Volume3D; + } - double JA = 0.25 * ( J1 + J2 + J3 + J4 ); + /* + Class : ScaledJacobian_i + Description : Functor for calculating volume of 3D element + */ + ScaledJacobian_i::ScaledJacobian_i() + { + myNumericalFunctorPtr.reset( new Controls::ScaledJacobian() ); + myFunctorPtr = myNumericalFunctorPtr; + } - double T1 = fabs( ( J1 - JA ) / JA ); - double T2 = fabs( ( J2 - JA ) / JA ); - double T3 = fabs( ( J3 - JA ) / JA ); - double T4 = fabs( ( J4 - JA ) / JA ); + FunctorType ScaledJacobian_i::GetFunctorType() + { + return SMESH::FT_ScaledJacobian; + } - return Max( Max( T1, T2 ), Max( T3, T4 ) ); -} + /* + Class : MaxElementLength2D_i + Description : Functor for calculating maximum length of 2D element + */ + MaxElementLength2D_i::MaxElementLength2D_i() + { + myNumericalFunctorPtr.reset( new Controls::MaxElementLength2D() ); + myFunctorPtr = myNumericalFunctorPtr; + } -int Taper_i::GetType() const -{ - return SMDSAbs_Face; -} + FunctorType MaxElementLength2D_i::GetFunctorType() + { + return SMESH::FT_MaxElementLength2D; + } -/* - Class : Skew_i - Description : Functor for calculating skew in degrees -*/ + /* + Class : MaxElementLength3D_i + Description : Functor for calculating maximum length of 3D element + */ + MaxElementLength3D_i::MaxElementLength3D_i() + { + myNumericalFunctorPtr.reset( new Controls::MaxElementLength3D() ); + myFunctorPtr = myNumericalFunctorPtr; + } -static inline double skewAngle( const gp_XYZ& p1, const gp_XYZ& p2, const gp_XYZ& p3 ) -{ - gp_XYZ p12 = ( p2 + p1 ) / 2; - gp_XYZ p23 = ( p3 + p2 ) / 2; - gp_XYZ p31 = ( p3 + p1 ) / 2; + FunctorType MaxElementLength3D_i::GetFunctorType() + { + return SMESH::FT_MaxElementLength3D; + } - return gp_Vec( p31 - p2 ).Angle( p12 - p23 ); -} + /* + Class : Length_i + Description : Functor for calculating length off edge + */ + Length_i::Length_i() + { + myNumericalFunctorPtr.reset( new Controls::Length() ); + myFunctorPtr = myNumericalFunctorPtr; + } -CORBA::Double Skew_i::GetValue( CORBA::Long theId ) -{ - TColgp_SequenceOfXYZ P; - if ( !getPoints( theId, P ) || P.Length() != 3 && P.Length() != 4 ) - return 0; + FunctorType Length_i::GetFunctorType() + { + return SMESH::FT_Length; + } - // Compute skew - static double PI2 = PI / 2; - if ( P.Length() == 3 ) + /* + Class : Length2D_i + Description : Functor for calculating length of edge + */ + Length2D_i::Length2D_i() { - double A0 = fabs( PI2 - skewAngle( P( 3 ), P( 1 ), P( 2 ) ) ); - double A1 = fabs( PI2 - skewAngle( P( 1 ), P( 2 ), P( 3 ) ) ); - double A2 = fabs( PI2 - skewAngle( P( 2 ), P( 3 ), P( 1 ) ) ); + myNumericalFunctorPtr.reset( new Controls::Length2D() ); + myFunctorPtr = myNumericalFunctorPtr; + } - return Max( A0, Max( A1, A2 ) ) * 180 / PI; + FunctorType Length2D_i::GetFunctorType() + { + return SMESH::FT_Length2D; } - else + + SMESH::Length2D::Values* Length2D_i::GetValues() { - gp_XYZ p12 = ( P( 1 ) + P( 2 ) ) / 2; - gp_XYZ p23 = ( P( 2 ) + P( 3 ) ) / 2; - gp_XYZ p34 = ( P( 3 ) + P( 4 ) ) / 2; - gp_XYZ p41 = ( P( 4 ) + P( 1 ) ) / 2; - - double A = fabs( PI2 - gp_Vec( p34 - p12 ).Angle( p23 - p41 ) ); + SMESH::Controls::Length2D::TValues aValues; + (dynamic_cast(myFunctorPtr.get()))->GetValues( aValues ); - return A * 180 / PI; - } -} + long i = 0, iEnd = aValues.size(); -int Skew_i::GetType() const -{ - return SMDSAbs_Face; -} + SMESH::Length2D::Values_var aResult = new SMESH::Length2D::Values(iEnd); + aResult->length(iEnd); -/* - Class : Area_i - Description : Functor for calculating area -*/ + SMESH::Controls::Length2D::TValues::const_iterator anIter; + for ( anIter = aValues.begin() ; anIter != aValues.end(); anIter++, i++ ) + { + const SMESH::Controls::Length2D::Value& aVal = *anIter; + SMESH::Length2D::Value &aValue = aResult[ i ]; -CORBA::Double Area_i::GetValue( CORBA::Long theId ) -{ - TColgp_SequenceOfXYZ P; - if ( !getPoints( theId, P ) || P.Length() != 3 && P.Length() != 4 ) - return 0; - - if ( P.Length() == 3 ) - return getArea( P( 1 ), P( 2 ), P( 3 ) ); - else - return getArea( P( 1 ), P( 2 ), P( 3 ) ) + getArea( P( 1 ), P( 3 ), P( 4 ) ); -} + aValue.myLength = aVal.myLength; + aValue.myPnt1 = aVal.myPntId[ 0 ]; + aValue.myPnt2 = aVal.myPntId[ 1 ]; + } -int Area_i::GetType() const -{ - return SMDSAbs_Face; -} + return aResult._retn(); + } -/* - Class : Length_i - Description : Functor for calculating length off edge -*/ -CORBA::Double Length_i::GetValue( CORBA::Long theId ) -{ - TColgp_SequenceOfXYZ P; - return getPoints( theId, P ) && P.Length() == 2 ? getDistance( P( 1 ), P( 2 ) ) : 0; -} + /* + Class : Length3D_i + Description : Functor for calculating length of edge + */ + Length3D_i::Length3D_i() + { + myNumericalFunctorPtr.reset( new Controls::Length3D() ); + myFunctorPtr = myNumericalFunctorPtr; + } -int Length_i::GetType() const -{ - return SMDSAbs_Edge; -} + FunctorType Length3D_i::GetFunctorType() + { + return SMESH::FT_Length3D; + } -/* - Class : MultiConnection_i - Description : Functor for calculating number of faces conneted to the edge -*/ + // SMESH::Length3D::Values* Length3D_i::GetValues() + // { + // SMESH::Controls::Length3D::TValues aValues; + // (dynamic_cast(myFunctorPtr.get()))->GetValues( aValues ); -CORBA::Double MultiConnection_i::GetValue( CORBA::Long theId ) -{ - return getNbMultiConnection( myMesh, theId ); -} + // long i = 0, iEnd = aValues.size(); -int MultiConnection_i::GetType() const -{ - return SMDSAbs_Edge; -} + // SMESH::Length3D::Values_var aResult = new SMESH::Length3D::Values(iEnd); + // aResult->length(iEnd); -/* - PREDICATES -*/ + // SMESH::Controls::Length3D::TValues::const_iterator anIter; + // for ( anIter = aValues.begin() ; anIter != aValues.end(); anIter++, i++ ) + // { + // const SMESH::Controls::Length3D::Value& aVal = *anIter; + // SMESH::Length3D::Value &aValue = aResult[ i ]; -/* - Class : Predicate_i - Description : Base class for all predicates -*/ -Predicate_i::Predicate_i() -: SALOME::GenericObj_i( SMESH_Gen_i::GetPOA() ) -{ - SMESH_Gen_i::GetPOA()->activate_object( this ); -} + // aValue.myLength = aVal.myLength; + // aValue.myPnt1 = aVal.myPntId[ 0 ]; + // aValue.myPnt2 = aVal.myPntId[ 1 ]; + // } + // return aResult._retn(); + // } -/* - Class : FreeBorders_i - Description : Predicate for free borders -*/ + /* + Class : Deflection2D_i + Description : Functor for calculating distance between a face and geometry + */ + Deflection2D_i::Deflection2D_i() + { + myNumericalFunctorPtr.reset( new Controls::Deflection2D() ); + myFunctorPtr = myNumericalFunctorPtr; + } -FreeBorders_i::FreeBorders_i() -{ - myMesh = 0; -} + FunctorType Deflection2D_i::GetFunctorType() + { + return SMESH::FT_Deflection2D; + } -void FreeBorders_i::SetMesh( SMESH_Mesh_ptr theMesh ) -{ - SMESH_Mesh_i* anImplPtr = - dynamic_cast( SMESH_Gen_i::GetServant( theMesh ).in() ); - myMesh = anImplPtr ? anImplPtr->GetImpl().GetMeshDS() : 0; -} + /* + Class : MultiConnection_i + Description : Functor for calculating number of faces conneted to the edge + */ + MultiConnection_i::MultiConnection_i() + { + myNumericalFunctorPtr.reset( new Controls::MultiConnection() ); + myFunctorPtr = myNumericalFunctorPtr; + } -CORBA::Boolean FreeBorders_i::IsSatisfy( CORBA::Long theId ) -{ - return getNbMultiConnection( myMesh, theId ) == 1; -} + FunctorType MultiConnection_i::GetFunctorType() + { + return SMESH::FT_MultiConnection; + } -int FreeBorders_i::GetType() const -{ - return SMDSAbs_Edge; -} + /* + Class : BallDiameter_i + Description : Functor returning diameter of a ball element + */ + BallDiameter_i::BallDiameter_i() + { + myNumericalFunctorPtr.reset( new Controls::BallDiameter() ); + myFunctorPtr = myNumericalFunctorPtr; + } -/* - Class : Comparator_i - Description : Base class for comparators -*/ + FunctorType BallDiameter_i::GetFunctorType() + { + return SMESH::FT_BallDiameter; + } -Comparator_i::Comparator_i() -{ - myMargin = 0; - myFunctor = 0; -} + /* + Class : NodeConnectivityNumber_i + Description : Functor returning diameter of a ball element + */ + NodeConnectivityNumber_i::NodeConnectivityNumber_i() + { + myNumericalFunctorPtr.reset( new Controls::NodeConnectivityNumber() ); + myFunctorPtr = myNumericalFunctorPtr; + } -Comparator_i::~Comparator_i() -{ - if ( myFunctor != 0 ) - myFunctor->Destroy(); -} + FunctorType NodeConnectivityNumber_i::GetFunctorType() + { + return SMESH::FT_NodeConnectivityNumber; + } -void Comparator_i::SetMesh( SMESH_Mesh_ptr theMesh ) -{ - if ( myFunctor != 0 ) - myFunctor->SetMesh( theMesh ); -} + /* + Class : MultiConnection2D_i + Description : Functor for calculating number of faces conneted to the edge + */ + MultiConnection2D_i::MultiConnection2D_i() + { + myNumericalFunctorPtr.reset( new Controls::MultiConnection2D() ); + myFunctorPtr = myNumericalFunctorPtr; + } -void Comparator_i::SetMargin( CORBA::Double theValue ) -{ - myMargin = theValue; -} + FunctorType MultiConnection2D_i::GetFunctorType() + { + return SMESH::FT_MultiConnection2D; + } -void Comparator_i::SetNumFunctor( NumericalFunctor_ptr theFunct ) -{ - if ( myFunctor != 0 ) - myFunctor->Destroy(); + SMESH::MultiConnection2D::Values* MultiConnection2D_i::GetValues() + { + SMESH::Controls::MultiConnection2D::MValues aValues; + (dynamic_cast(myFunctorPtr.get()))->GetValues( aValues ); + + long i = 0, iEnd = aValues.size(); - myFunctor = dynamic_cast( SMESH_Gen_i::GetServant( theFunct ).in() ); + SMESH::MultiConnection2D::Values_var aResult = new SMESH::MultiConnection2D::Values(iEnd); + aResult->length(iEnd); - if ( myFunctor != 0 ) - myFunctor->Register(); -} + SMESH::Controls::MultiConnection2D::MValues::const_iterator anIter; + for ( anIter = aValues.begin() ; anIter != aValues.end(); anIter++, i++ ) + { + const SMESH::Controls::MultiConnection2D::Value& aVal = (*anIter).first; + SMESH::MultiConnection2D::Value &aValue = aResult[ i ]; -int Comparator_i::GetType() const -{ - return myFunctor != 0 ? myFunctor->GetType() : SMDSAbs_All; -} + aValue.myPnt1 = aVal.myPntId[ 0 ]; + aValue.myPnt2 = aVal.myPntId[ 1 ]; + aValue.myNbConnects = (*anIter).second; + } -/* - Class : LessThan_i - Description : Comparator "<" -*/ + return aResult._retn(); + } + + /* + PREDICATES + */ + + + /* + Class : Predicate_i + Description : Base class for all predicates + */ + CORBA::Boolean Predicate_i::IsSatisfy( CORBA::Long theId ) + { + return myPredicatePtr->IsSatisfy( theId ); + } + + CORBA::Long Predicate_i::NbSatisfying( SMESH::SMESH_IDSource_ptr obj ) + { + SMESH::SMESH_Mesh_var meshVar = obj->GetMesh(); + const SMDS_Mesh* meshDS = MeshPtr2SMDSMesh( meshVar ); + if ( !meshDS ) + return 0; + myPredicatePtr->SetMesh( meshDS ); + + SMDSAbs_ElementType elemType = SMDSAbs_ElementType( GetElementType() ); + + int nb = 0; + SMDS_ElemIteratorPtr elemIt = + SMESH::DownCast( meshVar )->GetElements( obj, GetElementType() ); + if ( elemIt ) + while ( elemIt->more() ) + { + const SMDS_MeshElement* e = elemIt->next(); + if ( e && e->GetType() == elemType ) + nb += myPredicatePtr->IsSatisfy( e->GetID() ); + } + return nb; + } + + Controls::PredicatePtr Predicate_i::GetPredicate() + { + return myPredicatePtr; + } + + /* + Class : BadOrientedVolume_i + Description : Verify whether a mesh volume is incorrectly oriented from + the point of view of MED convention + */ + BadOrientedVolume_i::BadOrientedVolume_i() + { + Controls::PredicatePtr control( new Controls::BadOrientedVolume() ); + myFunctorPtr = myPredicatePtr = control; + } + + FunctorType BadOrientedVolume_i::GetFunctorType() + { + return SMESH::FT_BadOrientedVolume; + } + + /* + Class : BareBorderVolume_i + Description : Verify whether a mesh volume has a free facet without a face on it + */ + BareBorderVolume_i::BareBorderVolume_i() + { + Controls::PredicatePtr control( new Controls::BareBorderVolume() ); + myFunctorPtr = myPredicatePtr = control; + } + + FunctorType BareBorderVolume_i::GetFunctorType() + { + return SMESH::FT_BareBorderVolume; + } + + /* + Class : BareBorderFace_i + Description : Verify whether a mesh face has a free border without an edge on it + */ + BareBorderFace_i::BareBorderFace_i() + { + Controls::PredicatePtr control( new Controls::BareBorderFace() ); + myFunctorPtr = myPredicatePtr = control; + } + + FunctorType BareBorderFace_i::GetFunctorType() + { + return SMESH::FT_BareBorderFace; + } + + /* + Class : OverConstrainedVolume_i + Description : Verify whether a mesh volume has only one facet shared with other volumes + */ + OverConstrainedVolume_i::OverConstrainedVolume_i() + { + Controls::PredicatePtr control( new Controls::OverConstrainedVolume() ); + myFunctorPtr = myPredicatePtr = control; + } + + FunctorType OverConstrainedVolume_i::GetFunctorType() + { + return SMESH::FT_OverConstrainedVolume; + } + + /* + Class : OverConstrainedFace_i + Description : Verify whether a mesh face has only one border shared with other faces + */ + OverConstrainedFace_i::OverConstrainedFace_i() + { + Controls::PredicatePtr control( new Controls::OverConstrainedFace() ); + myFunctorPtr = myPredicatePtr = control; + } + + FunctorType OverConstrainedFace_i::GetFunctorType() + { + return SMESH::FT_OverConstrainedFace; + } + + /* + Class : BelongToMeshGroup_i + Description : Verify whether a mesh element is included into a mesh group + */ + BelongToMeshGroup_i::BelongToMeshGroup_i() + { + myBelongToMeshGroup = Controls::BelongToMeshGroupPtr( new Controls::BelongToMeshGroup() ); + myFunctorPtr = myPredicatePtr = myBelongToMeshGroup; + } + + BelongToMeshGroup_i::~BelongToMeshGroup_i() + { + SetGroup( SMESH::SMESH_GroupBase::_nil() ); + } + + void BelongToMeshGroup_i::SetGroup( SMESH::SMESH_GroupBase_ptr theGroup ) + { + if ( myGroup->_is_equivalent( theGroup )) + return; + + if ( ! myGroup->_is_nil() ) + myGroup->UnRegister(); + + myGroup = SMESH_GroupBase::_duplicate( theGroup ); + + myBelongToMeshGroup->SetGroup( 0 ); + if ( SMESH_GroupBase_i* gr_i = SMESH::DownCast< SMESH_GroupBase_i* >( myGroup )) + { + myBelongToMeshGroup->SetGroup( gr_i->GetGroupDS() ); + myGroup->Register(); + } + } + + void BelongToMeshGroup_i::SetGroupID( const char* theID ) // IOR or StoreName + { + myID = theID; + if ( strncmp( "IOR:", myID.c_str(), 4 ) == 0 ) // transient mode, no GUI + { + CORBA::Object_var obj = SMESH_Gen_i::GetORB()->string_to_object( myID.c_str() ); + SetGroup( SMESH::SMESH_GroupBase::_narrow( obj )); + } + else if ( strncmp( "0:", myID.c_str(), 2 ) == 0 ) // transient mode + GUI + { + SALOMEDS::SObject_wrap aSObj = SMESH_Gen_i::GetSMESHGen()->getStudyServant()->FindObjectID( myID.c_str() ); + if ( !aSObj->_is_nil() ) { + CORBA::Object_var obj = aSObj->GetObject(); + SetGroup( SMESH::SMESH_GroupBase::_narrow( obj )); + } + } + else if ( !myID.empty() ) // persistent mode + { + myBelongToMeshGroup->SetStoreName( myID ); + } + } + + std::string BelongToMeshGroup_i::GetGroupID() + { + if ( myGroup->_is_nil() ) + SMESH::SMESH_GroupBase_var( GetGroup() ); // decref the returned pointer + + if ( !myGroup->_is_nil() ) + myID = SMESH_Gen_i::GetORB()->object_to_string( myGroup ); + + return myID; + } + + SMESH::SMESH_GroupBase_ptr BelongToMeshGroup_i::GetGroup() + { + if ( myGroup->_is_nil() && myBelongToMeshGroup->GetGroup() ) + { + // search for a group in a current study + SMESH_Gen_i* aSMESHGen = SMESH_Gen_i::GetSMESHGen(); + if ( StudyContext* sc = aSMESHGen->GetStudyContext() ) + { + int id = 1; + std::string ior; + while (true) + { + ior = sc->getIORbyId( id++ ); + if ( ior.empty() ) break; + CORBA::Object_var obj = aSMESHGen->GetORB()->string_to_object( ior.c_str() ); + if ( SMESH_GroupBase_i* g_i = SMESH::DownCast( obj )) + if ( g_i->GetGroupDS() == myBelongToMeshGroup->GetGroup() ) + { + SetGroup( g_i->_this() ); + break; + } + } + } + } + return SMESH::SMESH_GroupBase::_duplicate( myGroup ); + } + + FunctorType BelongToMeshGroup_i::GetFunctorType() + { + return SMESH::FT_BelongToMeshGroup; + } + + /* + Class : BelongToGeom_i + Description : Predicate for selection on geometrical support + */ + BelongToGeom_i::BelongToGeom_i() + { + myBelongToGeomPtr.reset( new Controls::BelongToGeom() ); + myFunctorPtr = myPredicatePtr = myBelongToGeomPtr; + myShapeName = 0; + myShapeID = 0; + } + + BelongToGeom_i::~BelongToGeom_i() + { + CORBA::string_free( myShapeName ); + CORBA::string_free( myShapeID ); + } + + void BelongToGeom_i::SetGeom( GEOM::GEOM_Object_ptr theGeom ) + { + if ( theGeom->_is_nil() ) + return; + TopoDS_Shape aLocShape = SMESH_Gen_i::GetSMESHGen()->GeomObjectToShape( theGeom ); + myBelongToGeomPtr->SetGeom( aLocShape ); + TPythonDump()<SetGeom( theShape ); + } + + void BelongToGeom_i::SetElementType(ElementType theType) + { + myBelongToGeomPtr->SetType(SMDSAbs_ElementType(theType)); + TPythonDump()<SetGeom( getShapeByName( myShapeName ) ); + TPythonDump()<SetGeom( S ); + } + + char* BelongToGeom_i::GetShapeName() + { + return CORBA::string_dup( myShapeName ); + } + + char* BelongToGeom_i::GetShapeID() + { + return CORBA::string_dup( myShapeID ); + } + + void BelongToGeom_i::SetTolerance( CORBA::Double theToler ) + { + myBelongToGeomPtr->SetTolerance( theToler ); + TPythonDump()<GetTolerance(); + } + + /* + Class : BelongToSurface_i + Description : Predicate for selection on geometrical support + */ + BelongToSurface_i::BelongToSurface_i( const Handle(Standard_Type)& theSurfaceType ) + { + myElementsOnSurfacePtr.reset( new Controls::ElementsOnSurface() ); + myFunctorPtr = myPredicatePtr = myElementsOnSurfacePtr; + myShapeName = 0; + myShapeID = 0; + mySurfaceType = theSurfaceType; + } + + BelongToSurface_i::~BelongToSurface_i() + { + CORBA::string_free( myShapeName ); + CORBA::string_free( myShapeID ); + } + + void BelongToSurface_i::SetSurface( GEOM::GEOM_Object_ptr theGeom, ElementType theType ) + { + if ( theGeom->_is_nil() ) + return; + TopoDS_Shape aLocShape = SMESH_Gen_i::GetSMESHGen()->GeomObjectToShape( theGeom ); + + if ( aLocShape.ShapeType() == TopAbs_FACE ) + { + Handle(Geom_Surface) aSurf = BRep_Tool::Surface( TopoDS::Face( aLocShape ) ); + if ( !aSurf.IsNull() && aSurf->DynamicType() == mySurfaceType ) + { + myElementsOnSurfacePtr->SetSurface( aLocShape, (SMDSAbs_ElementType)theType ); + return; + } + } + + myElementsOnSurfacePtr->SetSurface( TopoDS_Shape(), (SMDSAbs_ElementType)theType ); + } + + void BelongToSurface_i::SetShapeName( const char* theName, ElementType theType ) + { + CORBA::string_free( myShapeName ); + myShapeName = CORBA::string_dup( theName ); + myElementsOnSurfacePtr->SetSurface( getShapeByName( myShapeName ), (SMDSAbs_ElementType)theType ); + TPythonDump()<SetSurface( S, (SMDSAbs_ElementType)theType ); + } + + char* BelongToSurface_i::GetShapeName() + { + return CORBA::string_dup( myShapeName ); + } + + char* BelongToSurface_i::GetShapeID() + { + return CORBA::string_dup( myShapeID ); + } + + void BelongToSurface_i::SetTolerance( CORBA::Double theToler ) + { + myElementsOnSurfacePtr->SetTolerance( theToler ); + TPythonDump()<GetTolerance(); + } + + void BelongToSurface_i::SetUseBoundaries( CORBA::Boolean theUseBndRestrictions ) + { + myElementsOnSurfacePtr->SetUseBoundaries( theUseBndRestrictions ); + TPythonDump()<GetUseBoundaries(); + } + + + /* + Class : BelongToPlane_i + Description : Verify whether mesh element lie in pointed Geom planar object + */ + + BelongToPlane_i::BelongToPlane_i() + : BelongToSurface_i( STANDARD_TYPE( Geom_Plane ) ) + { + } + + void BelongToPlane_i::SetPlane( GEOM::GEOM_Object_ptr theGeom, ElementType theType ) + { + BelongToSurface_i::SetSurface( theGeom, theType ); + TPythonDump()<_is_nil() ) + return; + TopoDS_Shape aLocShape = SMESH_Gen_i::GetSMESHGen()->GeomObjectToShape( theGeom ); + if ( !aLocShape.IsNull() && aLocShape.ShapeType() != TopAbs_FACE ) + aLocShape.Nullify(); + + BelongToSurface_i::myElementsOnSurfacePtr->SetSurface( aLocShape, (SMDSAbs_ElementType)theType ); + TPythonDump()<_is_nil() ) + return; + TopoDS_Shape aLocShape = SMESH_Gen_i::GetSMESHGen()->GeomObjectToShape( theGeom ); + myLyingOnGeomPtr->SetGeom( aLocShape ); + TPythonDump()<SetGeom( theShape ); + } + + void LyingOnGeom_i::SetElementType(ElementType theType){ + myLyingOnGeomPtr->SetType(SMDSAbs_ElementType(theType)); + TPythonDump()<SetGeom( getShapeByName( myShapeName ) ); + TPythonDump()<SetGeom( S ); + } + + char* LyingOnGeom_i::GetShapeName() + { + return CORBA::string_dup( myShapeName ); + } + + char* LyingOnGeom_i::GetShapeID() + { + return CORBA::string_dup( myShapeID ); + } + + void LyingOnGeom_i::SetTolerance( CORBA::Double theToler ) + { + myLyingOnGeomPtr->SetTolerance( theToler ); + TPythonDump()<GetTolerance(); + } + + /* + Class : FreeBorders_i + Description : Predicate for free borders + */ + FreeBorders_i::FreeBorders_i() + { + myPredicatePtr.reset(new Controls::FreeBorders()); + myFunctorPtr = myPredicatePtr; + } + + FunctorType FreeBorders_i::GetFunctorType() + { + return SMESH::FT_FreeBorders; + } + + /* + Class : FreeEdges_i + Description : Predicate for free borders + */ + FreeEdges_i::FreeEdges_i() + : myFreeEdgesPtr( new Controls::FreeEdges() ) + { + myFunctorPtr = myPredicatePtr = myFreeEdgesPtr; + } + + SMESH::FreeEdges::Borders* FreeEdges_i::GetBorders() + { + SMESH::Controls::FreeEdges::TBorders aBorders; + myFreeEdgesPtr->GetBoreders( aBorders ); + + long i = 0, iEnd = aBorders.size(); + + SMESH::FreeEdges::Borders_var aResult = new SMESH::FreeEdges::Borders; + aResult->length(iEnd); + + SMESH::Controls::FreeEdges::TBorders::const_iterator anIter; + for ( anIter = aBorders.begin() ; anIter != aBorders.end(); anIter++, i++ ) + { + const SMESH::Controls::FreeEdges::Border& aBord = *anIter; + SMESH::FreeEdges::Border &aBorder = aResult[ i ]; + + aBorder.myElemId = aBord.myElemId; + aBorder.myPnt1 = aBord.myPntId[ 0 ]; + aBorder.myPnt2 = aBord.myPntId[ 1 ]; + } + return aResult._retn(); + } + + FunctorType FreeEdges_i::GetFunctorType() + { + return SMESH::FT_FreeEdges; + } + + /* + Class : FreeFaces_i + Description : Predicate for free faces + */ + FreeFaces_i::FreeFaces_i() + { + myPredicatePtr.reset(new Controls::FreeFaces()); + myFunctorPtr = myPredicatePtr; + } + + FunctorType FreeFaces_i::GetFunctorType() + { + return SMESH::FT_FreeFaces; + } + + /* + Class : FreeNodes_i + Description : Predicate for free nodes + */ + FreeNodes_i::FreeNodes_i() + { + myPredicatePtr.reset(new Controls::FreeNodes()); + myFunctorPtr = myPredicatePtr; + } + + FunctorType FreeNodes_i::GetFunctorType() + { + return SMESH::FT_FreeNodes; + } + + /* + Class : EqualNodes_i + Description : Predicate for Equal nodes + */ + EqualNodes_i::EqualNodes_i() + { + myCoincidentNodesPtr.reset(new Controls::CoincidentNodes()); + myFunctorPtr = myPredicatePtr = myCoincidentNodesPtr; + } + + FunctorType EqualNodes_i::GetFunctorType() + { + return SMESH::FT_EqualNodes; + } + + void EqualNodes_i::SetTolerance( double tol ) + { + myCoincidentNodesPtr->SetTolerance( tol ); + } + + double EqualNodes_i::GetTolerance() + { + return myCoincidentNodesPtr->GetTolerance(); + } + + /* + Class : EqualEdges_i + Description : Predicate for Equal Edges + */ + EqualEdges_i::EqualEdges_i() + { + myPredicatePtr.reset(new Controls::CoincidentElements1D()); + myFunctorPtr = myPredicatePtr; + } + + FunctorType EqualEdges_i::GetFunctorType() + { + return SMESH::FT_EqualEdges; + } + + /* + Class : EqualFaces_i + Description : Predicate for Equal Faces + */ + EqualFaces_i::EqualFaces_i() + { + myPredicatePtr.reset(new Controls::CoincidentElements2D()); + myFunctorPtr = myPredicatePtr; + } + + FunctorType EqualFaces_i::GetFunctorType() + { + return SMESH::FT_EqualFaces; + } + + /* + Class : EqualVolumes_i + Description : Predicate for Equal Volumes + */ + EqualVolumes_i::EqualVolumes_i() + { + myPredicatePtr.reset(new Controls::CoincidentElements3D()); + myFunctorPtr = myPredicatePtr; + } + + FunctorType EqualVolumes_i::GetFunctorType() + { + return SMESH::FT_EqualVolumes; + } + + + /* + Class : RangeOfIds_i + Description : Predicate for Range of Ids. + Range may be specified with two ways. + 1. Using AddToRange method + 2. With SetRangeStr method. Parameter of this method is a string + like as "1,2,3,50-60,63,67,70-" + */ + + RangeOfIds_i::RangeOfIds_i() + { + myRangeOfIdsPtr.reset( new Controls::RangeOfIds() ); + myFunctorPtr = myPredicatePtr = myRangeOfIdsPtr; + } + + void RangeOfIds_i::SetRange( const SMESH::smIdType_array& theIds ) + { + SMESH::smIdType iEnd = theIds.length(); + for ( SMESH::smIdType i = 0; i < iEnd; i++ ) + myRangeOfIdsPtr->AddToRange( theIds[ i ] ); + TPythonDump()<SetRangeStr( + TCollection_AsciiString( (Standard_CString)theRange ) ); + } + + char* RangeOfIds_i::GetRangeStr() + { + TCollection_AsciiString aStr; + myRangeOfIdsPtr->GetRangeStr( aStr ); + return CORBA::string_dup( aStr.ToCString() ); + } + + void RangeOfIds_i::SetElementType( ElementType theType ) + { + myRangeOfIdsPtr->SetType( SMDSAbs_ElementType( theType ) ); + TPythonDump()<SetType(SMDSAbs_ElementType(theType)); + TPythonDump()<SetColorStr( + TCollection_AsciiString( (Standard_CString)theColor ) ); + TPythonDump()<GetColorStr( aStr ); + return CORBA::string_dup( aStr.ToCString() ); + } + + void GroupColor_i::SetElementType(ElementType theType) + { + myGroupColorPtr->SetType(SMDSAbs_ElementType(theType)); + TPythonDump()<SetType(SMDSAbs_ElementType(theType)); + TPythonDump()<SetGeomType(SMDSAbs_GeometryType(theType)); + TPythonDump()<GetGeomType(); + } + + FunctorType ElemGeomType_i::GetFunctorType() + { + return SMESH::FT_ElemGeomType; + } + + /* + Class : ElemEntityType_i + Description : Predicate check is element has indicated entity type + */ + ElemEntityType_i::ElemEntityType_i() + { + myElemEntityTypePtr.reset(new Controls::ElemEntityType()); + myFunctorPtr = myPredicatePtr = myElemEntityTypePtr; + } + + void ElemEntityType_i::SetElementType(ElementType theType) + { + myElemEntityTypePtr->SetType(SMDSAbs_ElementType(theType)); + TPythonDump()<SetElemEntityType(SMDSAbs_EntityType (theEntityType)); + TPythonDump()<GetElemEntityType(); + } + + FunctorType ElemEntityType_i::GetFunctorType() + { + return SMESH::FT_EntityType; + } + + /* + Class : CoplanarFaces_i + Description : Returns true if a mesh face is a coplanar neighbour to a given one + */ + CoplanarFaces_i::CoplanarFaces_i() + { + myCoplanarFacesPtr.reset(new Controls::CoplanarFaces()); + myFunctorPtr = myPredicatePtr = myCoplanarFacesPtr; + } + + void CoplanarFaces_i::SetFace ( CORBA::Long theFaceID ) + { + myCoplanarFacesPtr->SetFace(theFaceID); + TPythonDump()<SetTolerance(theToler); + TPythonDump()<GetFace(); + } + + char* CoplanarFaces_i::GetFaceAsString () const + { + TCollection_AsciiString str(Standard_Integer(myCoplanarFacesPtr->GetFace())); + return CORBA::string_dup( str.ToCString() ); + } + + CORBA::Double CoplanarFaces_i::GetTolerance() const + { + return myCoplanarFacesPtr->GetTolerance(); + } + + FunctorType CoplanarFaces_i::GetFunctorType() + { + return SMESH::FT_CoplanarFaces; + } + + /* + * Class : ConnectedElements_i + * Description : Returns true if an element is connected via other elements to the element + * located at a given point. + */ + ConnectedElements_i::ConnectedElements_i() + { + myConnectedElementsPtr.reset(new Controls::ConnectedElements()); + myFunctorPtr = myPredicatePtr = myConnectedElementsPtr; + } + + FunctorType ConnectedElements_i::GetFunctorType() + { + return FT_ConnectedElements; + } + + void ConnectedElements_i::SetElementType( ElementType theType ) + { + myConnectedElementsPtr->SetType( SMDSAbs_ElementType( theType )); + TPythonDump() << this << ".SetElementType( " << theType << " )"; + } + + void ConnectedElements_i::SetPoint( CORBA::Double x, CORBA::Double y, CORBA::Double z ) + { + myConnectedElementsPtr->SetPoint( x,y,z ); + myVertexID.clear(); + TPythonDump() << this << ".SetPoint( " << x << ", " << y << ", " << z << " )"; + } + + void ConnectedElements_i::SetVertex( GEOM::GEOM_Object_ptr vertex ) + { + TopoDS_Shape shape = SMESH_Gen_i::GetSMESHGen()->GeomObjectToShape( vertex ); + if ( shape.IsNull() ) + THROW_SALOME_CORBA_EXCEPTION( "ConnectedElements_i::SetVertex(): NULL Vertex", + SALOME::BAD_PARAM ); + + TopExp_Explorer v( shape, TopAbs_VERTEX ); + if ( !v.More() ) + THROW_SALOME_CORBA_EXCEPTION( "ConnectedElements_i::SetVertex(): empty vertex", + SALOME::BAD_PARAM ); + + gp_Pnt p = BRep_Tool::Pnt( TopoDS::Vertex( v.Current() )); + myConnectedElementsPtr->SetPoint( p.X(), p.Y(), p.Z() ); + // + CORBA::String_var id = vertex->GetStudyEntry(); + myVertexID = id.in(); + + TPythonDump() << this << ".SetVertex( " << vertex << " )"; + } + + void ConnectedElements_i::SetNode ( SMESH::smIdType nodeID ) + { + if ( nodeID < 1 ) + THROW_SALOME_CORBA_EXCEPTION( "ConnectedElements_i::SetNode(): nodeID must be > 0", + SALOME::BAD_PARAM ); + + myConnectedElementsPtr->SetNode( nodeID ); + myVertexID.clear(); + TPythonDump() << this << ".SetNode( " << nodeID << " )"; + } + + /*! + * \brief This is a comfort method for Filter dialog + */ + void ConnectedElements_i::SetThreshold ( const char* threshold, + SMESH::ConnectedElements::ThresholdType type ) + { + if ( !threshold ) + THROW_SALOME_CORBA_EXCEPTION( "ConnectedElements_i::SetThreshold(): NULL threshold", + SALOME::BAD_PARAM ); + switch ( type ) + { + case SMESH::ConnectedElements::POINT: // read 3 node coordinates /////////////////// + { + std::vector< double > xyz; + char* endptr; + do + { + // skip a separator + while ( *threshold && + *threshold != '+' && + *threshold != '-' && + !isdigit( *threshold )) + ++threshold; + if ( !*threshold ) + break; + // read a coordinate + xyz.push_back( strtod( threshold, &endptr )); + if ( threshold == endptr ) + { + xyz.resize( xyz.size() - 1 ); + break; + } + threshold = endptr; + } + while ( xyz.size() < 3 ); + + if ( xyz.size() < 3 ) + THROW_SALOME_CORBA_EXCEPTION + ( "ConnectedElements_i::SetThreshold(): invalid point coordinates", SALOME::BAD_PARAM ); + + SetPoint( xyz[0], xyz[1], xyz[2] ); + break; + } + case SMESH::ConnectedElements::VERTEX: // get a VERTEX by its entry ///////////////// + { + SALOMEDS::SObject_wrap sobj = SMESH_Gen_i::GetSMESHGen()->getStudyServant()->FindObjectID( threshold ); + if ( sobj->_is_nil() ) + THROW_SALOME_CORBA_EXCEPTION + ( "ConnectedElements_i::SetThreshold(): invalid vertex study entry", SALOME::BAD_PARAM ); + CORBA::Object_var obj = sobj->GetObject(); + GEOM::GEOM_Object_var vertex = GEOM::GEOM_Object::_narrow( obj ); + if ( vertex->_is_nil() ) + THROW_SALOME_CORBA_EXCEPTION + ( "ConnectedElements_i::SetThreshold(): no GEOM_Object in SObject", SALOME::BAD_PARAM ); + SetVertex( vertex ); + break; + } + case SMESH::ConnectedElements::NODE: // read a node ID //////////////////////////// + { + char* endptr; + int id = strtol( threshold, &endptr, 10 ); + if ( threshold == endptr ) + THROW_SALOME_CORBA_EXCEPTION + ( "ConnectedElements_i::SetThreshold(): invalid node ID", SALOME::BAD_PARAM ); + SetNode( id ); + break; + } + default: + THROW_SALOME_CORBA_EXCEPTION + ( "ConnectedElements_i::SetThreshold(): invalid ThresholdType", SALOME::BAD_PARAM ); + } + } + + char* ConnectedElements_i::GetThreshold ( SMESH::ConnectedElements::ThresholdType& type ) + { + std::string threshold; + if ( !myVertexID.empty() ) + { + threshold = myVertexID; + type = SMESH::ConnectedElements::VERTEX; + } + else + { + std::vector xyz = myConnectedElementsPtr->GetPoint(); + if ( xyz.size() == 3 ) + { + threshold = SMESH_Comment( xyz[0] ) << "; " << xyz[1] << "; " << xyz[2]; + type = SMESH::ConnectedElements::POINT; + } + else + { + threshold = SMESH_Comment( myConnectedElementsPtr->GetNode() ); + type = SMESH::ConnectedElements::NODE; + } + } + return CORBA::string_dup( threshold.c_str() ); + } + + /* + Class : Comparator_i + Description : Base class for comparators + */ + Comparator_i::Comparator_i(): + myNumericalFunctor( NULL ) + {} + + Comparator_i::~Comparator_i() + { + if ( myNumericalFunctor ) + myNumericalFunctor->UnRegister(); + } + + void Comparator_i::SetMargin( CORBA::Double theValue ) + { + myComparatorPtr->SetMargin( theValue ); + TPythonDump()<GetMargin(); + } + + void Comparator_i::SetNumFunctor( NumericalFunctor_ptr theFunct ) + { + if ( myNumericalFunctor ) + myNumericalFunctor->UnRegister(); + + myNumericalFunctor = DownCast(theFunct); + + if ( myNumericalFunctor ) + { + myComparatorPtr->SetNumFunctor( myNumericalFunctor->GetNumericalFunctor() ); + myNumericalFunctor->Register(); + TPythonDump()<" + */ + MoreThan_i::MoreThan_i() + { + myComparatorPtr.reset( new Controls::MoreThan() ); + myFunctorPtr = myPredicatePtr = myComparatorPtr; + } + + FunctorType MoreThan_i::GetFunctorType() + { + return SMESH::FT_MoreThan; + } + + + /* + Class : EqualTo_i + Description : Comparator "=" + */ + EqualTo_i::EqualTo_i() + : myEqualToPtr( new Controls::EqualTo() ) + { + myFunctorPtr = myPredicatePtr = myComparatorPtr = myEqualToPtr; + } + + void EqualTo_i::SetTolerance( CORBA::Double theToler ) + { + myEqualToPtr->SetTolerance( theToler ); + TPythonDump()<GetTolerance(); + } + + FunctorType EqualTo_i::GetFunctorType() + { + return SMESH::FT_EqualTo; + } + + /* + Class : LogicalNOT_i + Description : Logical NOT predicate + */ + LogicalNOT_i::LogicalNOT_i(): + myLogicalNOTPtr( new Controls::LogicalNOT() ), + myPredicate( NULL ) + { + myFunctorPtr = myPredicatePtr = myLogicalNOTPtr; + } + + LogicalNOT_i::~LogicalNOT_i() + { + if ( myPredicate ) + myPredicate->UnRegister(); + } + + void LogicalNOT_i::SetPredicate( Predicate_ptr thePredicate ) + { + if ( myPredicate ) + myPredicate->UnRegister(); + + myPredicate = SMESH::GetPredicate(thePredicate); + + if ( myPredicate ){ + myLogicalNOTPtr->SetPredicate(myPredicate->GetPredicate()); + myPredicate->Register(); + TPythonDump()<UnRegister(); + + if ( myPredicate2 ) + myPredicate2->UnRegister(); + } + + void LogicalBinary_i::SetMesh( SMESH_Mesh_ptr theMesh ) + { + if ( myPredicate1 ) + myPredicate1->SetMesh( theMesh ); + + if ( myPredicate2 ) + myPredicate2->SetMesh( theMesh ); + } + + void LogicalBinary_i::SetPredicate1( Predicate_ptr thePredicate ) + { + if ( myPredicate1 ) + myPredicate1->UnRegister(); + + myPredicate1 = SMESH::GetPredicate(thePredicate); + + if ( myPredicate1 ){ + myLogicalBinaryPtr->SetPredicate1(myPredicate1->GetPredicate()); + myPredicate1->Register(); + TPythonDump()<UnRegister(); + + myPredicate2 = SMESH::GetPredicate(thePredicate); + + if ( myPredicate2 ){ + myLogicalBinaryPtr->SetPredicate2(myPredicate2->GetPredicate()); + myPredicate2->Register(); + TPythonDump()<activate_object( this ); + } + + + FilterManager_i::~FilterManager_i() + { + //TPythonDump()<_this(); + TPythonDump()<_this(); + TPythonDump()<_this(); + TPythonDump()<_this(); + TPythonDump()<_this(); + TPythonDump() << aServant << " = " << this << ".CreateWarping3D()"; + return anObj._retn(); + } + + Taper_ptr FilterManager_i::CreateTaper() + { + SMESH::Taper_i* aServant = new SMESH::Taper_i(); + SMESH::Taper_var anObj = aServant->_this(); + TPythonDump()<_this(); + TPythonDump()<_this(); + TPythonDump()<_this(); + TPythonDump()<_this(); + TPythonDump()<_this(); + TPythonDump()<_this(); + TPythonDump()<_this(); + TPythonDump()<_this(); + TPythonDump()<_this(); + TPythonDump()<_this(); + TPythonDump()<_this(); + TPythonDump()<_this(); + TPythonDump()<_this(); + TPythonDump()<_this(); + TPythonDump()<_this(); + TPythonDump()<_this(); + TPythonDump()<_this(); + TPythonDump()<_this(); + TPythonDump()<_this(); + TPythonDump()<_this(); + TPythonDump()<_this(); + TPythonDump()<_this(); + TPythonDump()<_this(); + TPythonDump()<_this(); + TPythonDump()<_this(); + TPythonDump()<_this(); + TPythonDump()<_this(); + TPythonDump()<_this(); + TPythonDump()<_this(); + TPythonDump()<_this(); + TPythonDump()<_this(); + TPythonDump()<_this(); + TPythonDump()<_this(); + TPythonDump()<_this(); + TPythonDump()<_this(); + TPythonDump()<_this(); + TPythonDump()<_this(); + TPythonDump()<_this(); + TPythonDump()<_this(); + TPythonDump()<_this(); + TPythonDump()<_this(); + TPythonDump()<_this(); + TPythonDump()<_this(); + TPythonDump()<_this(); + TPythonDump()<_this(); + TPythonDump()<_this(); + TPythonDump()<_this(); + TPythonDump()<_this(); + TPythonDump()<_this(); + TPythonDump()<UnRegister(); + + if(!CORBA::is_nil(myMesh)) + myMesh->UnRegister(); + + myPredicate = 0; + FindBaseObjects(); + } + + //======================================================================= + // name : Filter_i::SetPredicate + // Purpose : Set predicate + //======================================================================= + void Filter_i::SetPredicate( Predicate_ptr thePredicate ) + { + if ( myPredicate ) + myPredicate->UnRegister(); + + myPredicate = SMESH::GetPredicate(thePredicate); + + if ( myPredicate ) + { + myFilter.SetPredicate( myPredicate->GetPredicate() ); + myPredicate->Register(); + if ( const SMDS_Mesh* aMesh = MeshPtr2SMDSMesh(myMesh)) + myPredicate->GetPredicate()->SetMesh( aMesh ); + TPythonDump()<GetElementType() : SMESH::ALL; + } + + //======================================================================= + // name : Filter_i::SetMesh + // Purpose : Set mesh + //======================================================================= + void + Filter_i:: + SetMesh( SMESH_Mesh_ptr theMesh ) + { + if(!CORBA::is_nil(theMesh)) + theMesh->Register(); + + if(!CORBA::is_nil(myMesh)) + myMesh->UnRegister(); + + myMesh = SMESH_Mesh::_duplicate( theMesh ); + TPythonDump()<GetPredicate()->SetMesh( aMesh ); + } + + SMESH::smIdType_array* + Filter_i:: + GetIDs() + { + return GetElementsId(myMesh); + } + + //======================================================================= + // name : Filter_i::GetElementsId + // Purpose : Get ids of entities + //======================================================================= + void + Filter_i:: + GetElementsId( Predicate_i* thePredicate, + const SMDS_Mesh* theMesh, + Controls::Filter::TIdSequence& theSequence ) + { + if (thePredicate) + Controls::Filter::GetElementsId(theMesh,thePredicate->GetPredicate(),theSequence); + } + + void + Filter_i:: + GetElementsId( Predicate_i* thePredicate, + SMESH_Mesh_ptr theMesh, + Controls::Filter::TIdSequence& theSequence ) + { + if (thePredicate) + if(const SMDS_Mesh* aMesh = MeshPtr2SMDSMesh(theMesh)) + Controls::Filter::GetElementsId(aMesh,thePredicate->GetPredicate(),theSequence); + } + + SMESH::smIdType_array* + Filter_i:: + GetElementsId( SMESH_Mesh_ptr theMesh ) + { + SMESH::smIdType_array_var anArray = new SMESH::smIdType_array; + if(!CORBA::is_nil(theMesh) && myPredicate){ + theMesh->Load(); + Controls::Filter::TIdSequence aSequence; + GetElementsId(myPredicate,theMesh,aSequence); + long i = 0, iEnd = aSequence.size(); + anArray->length( iEnd ); + for ( ; i < iEnd; i++ ) + anArray[ i ] = aSequence[i]; + } + return anArray._retn(); + } + + SMESH::smIdType_array* + Filter_i:: + GetElementsIdFromParts( const ListOfIDSources& theParts ) + { + SMESH::smIdType_array_var array = new SMESH::smIdType_array; + if ( theParts.length() > 0 && myPredicate ) + { + SMESH_Mesh_ptr mesh = theParts[0]->GetMesh(); + mesh->Load(); + const SMDS_Mesh* meshDS = MeshPtr2SMDSMesh( mesh ); + Controls::Filter::TIdSequence totalSequence; + for ( CORBA::ULong i = 0; i < theParts.length(); ++i ) + { + if ( SMESH::Filter_i* filter = SMESH::DownCast( theParts[i] )) + filter->SetMesh( mesh ); + SMDS_ElemIteratorPtr iter = SMESH_Mesh_i::GetElements( theParts[i], GetElementType() ); + if ( iter && meshDS ) + { + Controls::Filter::TIdSequence sequence; + Controls::Filter::GetElementsId( meshDS, myPredicate->GetPredicate(), sequence, iter ); + totalSequence.insert( totalSequence.end(), sequence.begin(), sequence.end() ); + } + } + array->length( totalSequence.size() ); + for ( size_t i = 0; i < totalSequence.size(); ++i ) + array[ i ] = totalSequence[ i ]; + } + return array._retn(); + } + + //============================================================================= + /*! + * \brief Returns number of mesh elements per each \a EntityType + */ + //============================================================================= + + SMESH::smIdType_array* Filter_i::GetMeshInfo() + { + SMESH::smIdType_array_var aRes = new SMESH::smIdType_array(); + aRes->length(SMESH::Entity_Last); + for (int i = 0; i < SMESH::Entity_Last; i++) + aRes[i] = 0; + + if ( !CORBA::is_nil(myMesh) && myPredicate ) + { + const SMDS_Mesh* aMesh = MeshPtr2SMDSMesh(myMesh); + SMDS_ElemIteratorPtr it = aMesh->elementsIterator( SMDSAbs_ElementType( GetElementType() )); + while ( it->more() ) + { + const SMDS_MeshElement* anElem = it->next(); + if ( myPredicate->IsSatisfy( anElem->GetID() ) ) + aRes[ anElem->GetEntityType() ]++; + } + } + + return aRes._retn(); + } + + //============================================================================= + /*! + * \brief Returns number of mesh elements of each \a ElementType + */ + //============================================================================= + + SMESH::smIdType_array* Filter_i::GetNbElementsByType() + { + SMESH::smIdType_array_var aRes = new SMESH::smIdType_array(); + aRes->length(SMESH::NB_ELEMENT_TYPES); + for (int i = 0; i < SMESH::NB_ELEMENT_TYPES; i++) + aRes[i] = 0; + + if ( !CORBA::is_nil(myMesh) && myPredicate ) { + const SMDS_Mesh* aMesh = MeshPtr2SMDSMesh(myMesh); + SMDS_ElemIteratorPtr it = aMesh->elementsIterator( SMDSAbs_ElementType( GetElementType() )); + SMESH::smIdType& nbElems = aRes[ GetElementType() ]; + while ( it->more() ) + { + const SMDS_MeshElement* anElem = it->next(); + if ( myPredicate->IsSatisfy( anElem->GetID() ) ) + nbElems++; + } + } + + return aRes._retn(); + } + + + //================================================================================ + /*! + * \brief Return GetElementType() within an array + * Implement SMESH_IDSource interface + */ + //================================================================================ + + SMESH::array_of_ElementType* Filter_i::GetTypes() + { + SMESH::array_of_ElementType_var types = new SMESH::array_of_ElementType; + + // check if any element passes through the filter + if ( !CORBA::is_nil(myMesh) && myPredicate ) + { + const SMDS_Mesh* aMesh = MeshPtr2SMDSMesh(myMesh); + SMDS_ElemIteratorPtr it = aMesh->elementsIterator( SMDSAbs_ElementType( GetElementType() )); + bool satisfies = false; + while ( !satisfies && it->more() ) + satisfies = myPredicate->IsSatisfy( it->next()->GetID() ); + if ( satisfies ) { + types->length( 1 ); + types[0] = GetElementType(); + } + } + return types._retn(); + } + + //======================================================================= + //function : GetMesh + //purpose : Returns mesh + //======================================================================= + + SMESH::SMESH_Mesh_ptr Filter_i::GetMesh() + { + return SMESH_Mesh::_duplicate( myMesh ); + } + + //======================================================================= + //function : GetVtkUgStream + //purpose : Return data vtk unstructured grid (not implemented) + //======================================================================= + + SALOMEDS::TMPFile* Filter_i::GetVtkUgStream() + { + SALOMEDS::TMPFile_var SeqFile; + return SeqFile._retn(); + } + //======================================================================= + // name : getCriteria + // Purpose : Retrieve criterions from predicate + //======================================================================= + static inline void getPrediacates( Predicate_i* thePred, + std::vector & thePredVec ) + { + const int aFType = thePred->GetFunctorType(); + + switch ( aFType ) + { + case FT_LogicalNOT: + { + Predicate_i* aPred = ( dynamic_cast( thePred ) )->GetPredicate_i(); + getPrediacates( aPred, thePredVec ); + break; + } + case FT_LogicalAND: + case FT_LogicalOR: + { + Predicate_i* aPred1 = ( dynamic_cast( thePred ) )->GetPredicate1_i(); + Predicate_i* aPred2 = ( dynamic_cast( thePred ) )->GetPredicate2_i(); + getPrediacates( aPred1, thePredVec ); + getPrediacates( aPred2, thePredVec ); + break; + } + default:; + } + thePredVec.push_back( thePred ); + } + + //======================================================================= + // name : getCriteria + // Purpose : Retrieve criterions from predicate + //======================================================================= + static inline bool getCriteria( Predicate_i* thePred, + SMESH::Filter::Criteria_out theCriteria ) + { + const int aFType = thePred->GetFunctorType(); + + switch ( aFType ) + { + case FT_LogicalNOT: + { + Predicate_i* aPred = ( dynamic_cast( thePred ) )->GetPredicate_i(); + getCriteria( aPred, theCriteria ); + theCriteria[ theCriteria->length() - 1 ].UnaryOp = FT_LogicalNOT; + } + return true; + + case FT_LogicalAND: + case FT_LogicalOR: + { + Predicate_i* aPred1 = ( dynamic_cast( thePred ) )->GetPredicate1_i(); + Predicate_i* aPred2 = ( dynamic_cast( thePred ) )->GetPredicate2_i(); + if ( !getCriteria( aPred1, theCriteria ) ) + return false; + theCriteria[ theCriteria->length() - 1 ].BinaryOp = aFType; + return getCriteria( aPred2, theCriteria ); + } + case FT_Undefined: + return false; + } + + // resize theCriteria + CORBA::ULong i = theCriteria->length(); + theCriteria->length( i + 1 ); + theCriteria[ i ] = createCriterion(); + + // set members of the added Criterion + + theCriteria[ i ].Type = aFType; + theCriteria[ i ].TypeOfElement = thePred->GetElementType(); + + switch ( aFType ) + { + case FT_FreeBorders: + case FT_FreeEdges: + case FT_FreeFaces: + case FT_LinearOrQuadratic: + case FT_FreeNodes: + case FT_EqualEdges: + case FT_EqualFaces: + case FT_EqualVolumes: + case FT_BadOrientedVolume: + case FT_BareBorderVolume: + case FT_BareBorderFace: + case FT_OverConstrainedVolume: + case FT_OverConstrainedFace: + { + return true; + } + case FT_BelongToMeshGroup: + { + BelongToMeshGroup_i* aPred = dynamic_cast( thePred ); + SMESH::SMESH_GroupBase_var grp = aPred->GetGroup(); + if ( !grp->_is_nil() ) + { + theCriteria[ i ].ThresholdStr = grp->GetName(); + theCriteria[ i ].ThresholdID = aPred->GetGroupID().c_str(); + } + return true; + } + case FT_BelongToGeom: + { + BelongToGeom_i* aPred = dynamic_cast( thePred ); + theCriteria[ i ].ThresholdStr = aPred->GetShapeName(); + theCriteria[ i ].ThresholdID = aPred->GetShapeID(); + theCriteria[ i ].Tolerance = aPred->GetTolerance(); + return true; + } + case FT_BelongToPlane: + case FT_BelongToCylinder: + case FT_BelongToGenSurface: + { + BelongToSurface_i* aPred = dynamic_cast( thePred ); + theCriteria[ i ].ThresholdStr = aPred->GetShapeName(); + theCriteria[ i ].ThresholdID = aPred->GetShapeID(); + theCriteria[ i ].Tolerance = aPred->GetTolerance(); + return true; + } + case FT_LyingOnGeom: + { + LyingOnGeom_i* aPred = dynamic_cast( thePred ); + theCriteria[ i ].ThresholdStr = aPred->GetShapeName(); + theCriteria[ i ].ThresholdID = aPred->GetShapeID(); + theCriteria[ i ].Tolerance = aPred->GetTolerance(); + return true; + } + case FT_CoplanarFaces: + { + CoplanarFaces_i* aPred = dynamic_cast( thePred ); + theCriteria[ i ].ThresholdID = aPred->GetFaceAsString(); + theCriteria[ i ].Tolerance = aPred->GetTolerance(); + return true; + } + case FT_ConnectedElements: + { + ConnectedElements_i* aPred = dynamic_cast( thePred ); + SMESH::ConnectedElements::ThresholdType type; + CORBA::String_var threshold = aPred->GetThreshold( type ); + switch ( type ) { + case SMESH::ConnectedElements::POINT: + theCriteria[ i ].ThresholdStr = threshold; break; + case SMESH::ConnectedElements::VERTEX: + theCriteria[ i ].ThresholdID = threshold; break; + case SMESH::ConnectedElements::NODE: + theCriteria[ i ].Threshold = atof( threshold.in() ); break; + default:; + } + return true; + } + case FT_EqualNodes: + { + EqualNodes_i* aPred = dynamic_cast( thePred ); + theCriteria[ i ].Tolerance = aPred->GetTolerance(); + return true; + } + case FT_RangeOfIds: + { + RangeOfIds_i* aPred = dynamic_cast( thePred ); + theCriteria[ i ].ThresholdStr = aPred->GetRangeStr(); + return true; + } + case FT_LessThan: + case FT_MoreThan: + case FT_EqualTo: + { + Comparator_i* aCompar = dynamic_cast( thePred ); + theCriteria[ i ].Type = aCompar->GetNumFunctor_i()->GetFunctorType(); + theCriteria[ i ].Compare = aFType; + theCriteria[ i ].Threshold = aCompar->GetMargin(); + if ( aFType == FT_EqualTo ) + { + EqualTo_i* aCompar = dynamic_cast( thePred ); + theCriteria[ i ].Tolerance = aCompar->GetTolerance(); + } + return true; + } + case FT_GroupColor: + { + GroupColor_i* aPred = dynamic_cast( thePred ); + theCriteria[ i ].ThresholdStr = aPred->GetColorStr(); + return true; + } + case FT_ElemGeomType: + { + ElemGeomType_i* aPred = dynamic_cast( thePred ); + theCriteria[ i ].Threshold = (double)aPred->GetGeometryType(); + return true; + } + case FT_EntityType: + { + ElemEntityType_i* aPred = dynamic_cast( thePred ); + theCriteria[ i ].Threshold = (double)aPred->GetEntityType(); + return true; + } + default: + return false; + } + } + + //======================================================================= + // name : Filter_i::GetCriteria + // Purpose : Retrieve criterions from predicate + //======================================================================= + CORBA::Boolean Filter_i::GetCriteria( SMESH::Filter::Criteria_out theCriteria ) + { + theCriteria = new SMESH::Filter::Criteria; + return myPredicate != 0 ? getCriteria( myPredicate, theCriteria ) : true; + } + + //======================================================================= + // name : Filter_i::SetCriteria + // Purpose : Create new predicate and set criterions in it + //======================================================================= + CORBA::Boolean Filter_i::SetCriteria( const SMESH::Filter::Criteria& theCriteria ) + { + SetPredicate( SMESH::Predicate::_nil() ); + + SMESH::FilterManager_i* aFilter = new SMESH::FilterManager_i(); + FilterManager_ptr aFilterMgr = aFilter->_this(); + + // CREATE two lists ( PREDICATES and LOG OP ) + + // Criterion + TPythonDump()<<"aCriteria = []"; + std::list aPredicates; + std::list aBinaries; + for ( int i = 0, n = theCriteria.length(); i < n; i++ ) + { + int aCriterion = theCriteria[ i ].Type; + int aCompare = theCriteria[ i ].Compare; + double aThreshold = theCriteria[ i ].Threshold; + const char* aThresholdStr = theCriteria[ i ].ThresholdStr; + const char* aThresholdID = theCriteria[ i ].ThresholdID; + int aUnary = theCriteria[ i ].UnaryOp; + int aBinary = theCriteria[ i ].BinaryOp; + double aTolerance = theCriteria[ i ].Tolerance; + ElementType aTypeOfElem = theCriteria[ i ].TypeOfElement; + long aPrecision = theCriteria[ i ].Precision; + + { + TPythonDump pd; + pd << "aCriterion = SMESH.Filter.Criterion(" + << aCriterion << ", " + << aCompare << ", " + << aThreshold << ", '" + << aThresholdStr << "', '" + << aThresholdID << "', " + << aUnary << ", " + << aBinary << ", " + << aTolerance << ", " + << aTypeOfElem << ", " + << aPrecision << ")"; + } + TPythonDump pd; + + SMESH::Predicate_ptr aPredicate = SMESH::Predicate::_nil(); + SMESH::NumericalFunctor_ptr aFunctor = SMESH::NumericalFunctor::_nil(); + + switch ( aCriterion ) + { + // Functors + + case SMESH::FT_MultiConnection: + aFunctor = aFilterMgr->CreateMultiConnection(); + break; + case SMESH::FT_MultiConnection2D: + aFunctor = aFilterMgr->CreateMultiConnection2D(); + break; + case SMESH::FT_Length: + aFunctor = aFilterMgr->CreateLength(); + break; + case SMESH::FT_Length2D: + aFunctor = aFilterMgr->CreateLength2D(); + break; + case SMESH::FT_Length3D: + aFunctor = aFilterMgr->CreateLength3D(); + break; + case SMESH::FT_Deflection2D: + aFunctor = aFilterMgr->CreateDeflection2D(); + break; + case SMESH::FT_AspectRatio: + aFunctor = aFilterMgr->CreateAspectRatio(); + break; + case SMESH::FT_AspectRatio3D: + aFunctor = aFilterMgr->CreateAspectRatio3D(); + break; + case SMESH::FT_Warping: + aFunctor = aFilterMgr->CreateWarping(); + break; + case SMESH::FT_Warping3D: + aFunctor = aFilterMgr->CreateWarping3D(); + break; + case SMESH::FT_MinimumAngle: + aFunctor = aFilterMgr->CreateMinimumAngle(); + break; + case SMESH::FT_Taper: + aFunctor = aFilterMgr->CreateTaper(); + break; + case SMESH::FT_Skew: + aFunctor = aFilterMgr->CreateSkew(); + break; + case SMESH::FT_Area: + aFunctor = aFilterMgr->CreateArea(); + break; + case SMESH::FT_Volume3D: + aFunctor = aFilterMgr->CreateVolume3D(); + break; + case SMESH::FT_MaxElementLength2D: + aFunctor = aFilterMgr->CreateMaxElementLength2D(); + break; + case SMESH::FT_MaxElementLength3D: + aFunctor = aFilterMgr->CreateMaxElementLength3D(); + break; + case SMESH::FT_BallDiameter: + aFunctor = aFilterMgr->CreateBallDiameter(); + break; + case SMESH::FT_NodeConnectivityNumber: + aFunctor = aFilterMgr->CreateNodeConnectivityNumber(); + break; + case SMESH::FT_ScaledJacobian: + aFunctor = aFilterMgr->CreateScaledJacobian(); + break; + + // Predicates + + case SMESH::FT_FreeBorders: + aPredicate = aFilterMgr->CreateFreeBorders(); + break; + case SMESH::FT_FreeEdges: + aPredicate = aFilterMgr->CreateFreeEdges(); + break; + case SMESH::FT_FreeFaces: + aPredicate = aFilterMgr->CreateFreeFaces(); + break; + case SMESH::FT_FreeNodes: + aPredicate = aFilterMgr->CreateFreeNodes(); + break; + case SMESH::FT_EqualNodes: + { + SMESH::EqualNodes_ptr pred = aFilterMgr->CreateEqualNodes(); + pred->SetTolerance( aTolerance ); + aPredicate = pred; + break; + } + case SMESH::FT_EqualEdges: + aPredicate = aFilterMgr->CreateEqualEdges(); + break; + case SMESH::FT_EqualFaces: + aPredicate = aFilterMgr->CreateEqualFaces(); + break; + case SMESH::FT_EqualVolumes: + aPredicate = aFilterMgr->CreateEqualVolumes(); + break; + case SMESH::FT_BelongToMeshGroup: + { + SMESH::BelongToMeshGroup_ptr tmpPred = aFilterMgr->CreateBelongToMeshGroup(); + tmpPred->SetGroupID( aThresholdID ); + aPredicate = tmpPred; + } + break; + case SMESH::FT_BelongToGeom: + { + SMESH::BelongToGeom_ptr tmpPred = aFilterMgr->CreateBelongToGeom(); + tmpPred->SetElementType( aTypeOfElem ); + tmpPred->SetShape( aThresholdID, aThresholdStr ); + tmpPred->SetTolerance( aTolerance ); + aPredicate = tmpPred; + } + break; + case SMESH::FT_BelongToPlane: + case SMESH::FT_BelongToCylinder: + case SMESH::FT_BelongToGenSurface: + { + SMESH::BelongToSurface_ptr tmpPred; + switch ( aCriterion ) { + case SMESH::FT_BelongToPlane: + tmpPred = aFilterMgr->CreateBelongToPlane(); break; + case SMESH::FT_BelongToCylinder: + tmpPred = aFilterMgr->CreateBelongToCylinder(); break; + default: + tmpPred = aFilterMgr->CreateBelongToGenSurface(); + } + tmpPred->SetShape( aThresholdID, aThresholdStr, aTypeOfElem ); + tmpPred->SetTolerance( aTolerance ); + aPredicate = tmpPred; + } + break; + case SMESH::FT_LyingOnGeom: + { + SMESH::LyingOnGeom_ptr tmpPred = aFilterMgr->CreateLyingOnGeom(); + tmpPred->SetElementType( aTypeOfElem ); + tmpPred->SetShape( aThresholdID, aThresholdStr ); + tmpPred->SetTolerance( aTolerance ); + aPredicate = tmpPred; + } + break; + case SMESH::FT_RangeOfIds: + { + SMESH::RangeOfIds_ptr tmpPred = aFilterMgr->CreateRangeOfIds(); + tmpPred->SetRangeStr( aThresholdStr ); + tmpPred->SetElementType( aTypeOfElem ); + aPredicate = tmpPred; + } + break; + case SMESH::FT_BadOrientedVolume: + { + aPredicate = aFilterMgr->CreateBadOrientedVolume(); + } + break; + case SMESH::FT_BareBorderVolume: + { + aPredicate = aFilterMgr->CreateBareBorderVolume(); + } + break; + case SMESH::FT_BareBorderFace: + { + aPredicate = aFilterMgr->CreateBareBorderFace(); + } + break; + case SMESH::FT_OverConstrainedVolume: + { + aPredicate = aFilterMgr->CreateOverConstrainedVolume(); + } + break; + case SMESH::FT_OverConstrainedFace: + { + aPredicate = aFilterMgr->CreateOverConstrainedFace(); + } + break; + case SMESH::FT_LinearOrQuadratic: + { + SMESH::LinearOrQuadratic_ptr tmpPred = aFilterMgr->CreateLinearOrQuadratic(); + tmpPred->SetElementType( aTypeOfElem ); + aPredicate = tmpPred; + break; + } + case SMESH::FT_GroupColor: + { + SMESH::GroupColor_ptr tmpPred = aFilterMgr->CreateGroupColor(); + tmpPred->SetElementType( aTypeOfElem ); + tmpPred->SetColorStr( aThresholdStr ); + aPredicate = tmpPred; + break; + } + case SMESH::FT_ElemGeomType: + { + SMESH::ElemGeomType_ptr tmpPred = aFilterMgr->CreateElemGeomType(); + tmpPred->SetElementType( aTypeOfElem ); + tmpPred->SetGeometryType( (GeometryType)(int)(aThreshold + 0.5) ); + aPredicate = tmpPred; + break; + } + case SMESH::FT_EntityType: + { + SMESH::ElemEntityType_ptr tmpPred = aFilterMgr->CreateElemEntityType(); + tmpPred->SetElementType( aTypeOfElem ); + tmpPred->SetEntityType( EntityType( (int (aThreshold + 0.5)))); + aPredicate = tmpPred; + break; + } + case SMESH::FT_CoplanarFaces: + { + SMESH::CoplanarFaces_ptr tmpPred = aFilterMgr->CreateCoplanarFaces(); + tmpPred->SetFace( atol (aThresholdID )); + tmpPred->SetTolerance( aTolerance ); + aPredicate = tmpPred; + break; + } + case SMESH::FT_ConnectedElements: + { + SMESH::ConnectedElements_ptr tmpPred = aFilterMgr->CreateConnectedElements(); + if ( strlen( aThresholdID ) > 0 ) // shape ID + tmpPred->SetThreshold( aThresholdID, SMESH::ConnectedElements::VERTEX ); + else if ( strlen( aThresholdStr ) > 0 ) // point coords + tmpPred->SetThreshold( aThresholdStr, SMESH::ConnectedElements::POINT ); + else if ( aThreshold >= 1 ) + tmpPred->SetNode( (CORBA::Long) aThreshold ); // node ID + tmpPred->SetElementType( aTypeOfElem ); + aPredicate = tmpPred; + break; + } + + default: + continue; + } + + // Comparator + if ( !aFunctor->_is_nil() && aPredicate->_is_nil() ) + { + SMESH::Comparator_ptr aComparator = SMESH::Comparator::_nil(); + + if ( aCompare == SMESH::FT_LessThan ) + aComparator = aFilterMgr->CreateLessThan(); + else if ( aCompare == SMESH::FT_MoreThan ) + aComparator = aFilterMgr->CreateMoreThan(); + else if ( aCompare == SMESH::FT_EqualTo ) + aComparator = aFilterMgr->CreateEqualTo(); + else + continue; + + aComparator->SetNumFunctor( aFunctor ); + aComparator->SetMargin( aThreshold ); + + if ( aCompare == FT_EqualTo ) + { + SMESH::EqualTo_var anEqualTo = SMESH::EqualTo::_narrow( aComparator ); + anEqualTo->SetTolerance( aTolerance ); + } + + aPredicate = aComparator; + + aFunctor->SetPrecision( aPrecision ); + } + + // Logical not + if ( aUnary == FT_LogicalNOT ) + { + SMESH::LogicalNOT_ptr aNotPred = aFilterMgr->CreateLogicalNOT(); + aNotPred->SetPredicate( aPredicate ); + aPredicate = aNotPred; + } + + // logical op + aPredicates.push_back( aPredicate ); + aBinaries.push_back( aBinary ); + pd <<"aCriteria.append(aCriterion)"; + + } // end of for + TPythonDump pd; pd< aResList; + + std::list::iterator aPredIter; + std::list::iterator aBinaryIter; + + SMESH::Predicate_ptr aPrevPredicate = SMESH::Predicate::_nil(); + int aPrevBinary = SMESH::FT_Undefined; + if ( !aBinaries.empty() ) + aBinaries.back() = SMESH::FT_Undefined; + + for ( aPredIter = aPredicates.begin(), aBinaryIter = aBinaries.begin(); + aPredIter != aPredicates.end() && aBinaryIter != aBinaries.end(); + ++aPredIter, ++aBinaryIter ) + { + int aCurrBinary = *aBinaryIter; + + SMESH::Predicate_ptr aCurrPred = SMESH::Predicate::_nil(); + + if ( aPrevBinary == SMESH::FT_LogicalAND ) + { + + SMESH::LogicalBinary_ptr aBinaryPred = aFilterMgr->CreateLogicalAND(); + aBinaryPred->SetPredicate1( aPrevPredicate ); + aBinaryPred->SetPredicate2( *aPredIter ); + aCurrPred = aBinaryPred; + } + else + aCurrPred = *aPredIter; + + if ( aCurrBinary != SMESH::FT_LogicalAND ) + aResList.push_back( aCurrPred ); + + aPrevPredicate = aCurrPred; + aPrevBinary = aCurrBinary; + } + + // combine all "OR" operations + + SMESH::Predicate_ptr aResPredicate = SMESH::Predicate::_nil(); + + if ( aResList.size() == 1 ) + aResPredicate = *aResList.begin(); + else if ( aResList.size() > 1 ) + { + std::list::iterator anIter = aResList.begin(); + aResPredicate = *anIter; + anIter++; + for ( ; anIter != aResList.end(); ++anIter ) + { + SMESH::LogicalBinary_ptr aBinaryPred = aFilterMgr->CreateLogicalOR(); + aBinaryPred->SetPredicate1( aResPredicate ); + aBinaryPred->SetPredicate2( *anIter ); + aResPredicate = aBinaryPred; + } + } + + SetPredicate( aResPredicate ); + if ( !aResPredicate->_is_nil() ) + aResPredicate->UnRegister(); + + return !aResPredicate->_is_nil(); + } + + //======================================================================= + // name : Filter_i::GetPredicate_i + // Purpose : Get implementation of predicate + //======================================================================= + Predicate_i* Filter_i::GetPredicate_i() + { + return myPredicate; + } + + //======================================================================= + // name : Filter_i::GetPredicate + // Purpose : Get predicate + //======================================================================= + Predicate_ptr Filter_i::GetPredicate() + { + if ( myPredicate == 0 ) + return SMESH::Predicate::_nil(); + else + { + SMESH::Predicate_var anObj = myPredicate->_this(); + // if ( SMESH::Functor_i* fun = SMESH::DownCast( anObj )) + // TPythonDump() << fun << " = " << this << ".GetPredicate()"; + return anObj._retn(); + } + } + + //================================================================================ + /*! + * \brief Find groups it depends on + */ + //================================================================================ + + void Filter_i::FindBaseObjects() + { + // release current groups + for ( size_t i = 0; i < myBaseGroups.size(); ++i ) + if ( myBaseGroups[i] ) + { + myBaseGroups[i]->RemoveModifWaiter( this ); + myBaseGroups[i]->UnRegister(); + } + + // remember new groups + myBaseGroups.clear(); + if ( myPredicate ) + { + std::vector predicates; + getPrediacates( myPredicate, predicates ); + for ( size_t i = 0; i < predicates.size(); ++i ) + if ( BelongToMeshGroup_i* bmg = dynamic_cast< BelongToMeshGroup_i* >( predicates[i] )) + { + SMESH::SMESH_GroupBase_var g = bmg->GetGroup(); + SMESH_GroupBase_i* g_i = SMESH::DownCast< SMESH_GroupBase_i*>( g ); + if ( g_i ) + { + g_i->AddModifWaiter( this ); + g_i->Register(); + myBaseGroups.push_back( g_i ); + } + } + } + } + + //================================================================================ + /*! + * \brief When notified on removal of myBaseGroups[i], remove a reference to a + * group from a predicate + */ + //================================================================================ + + void Filter_i::OnBaseObjModified(NotifyerAndWaiter* group, bool removed) + { + if ( !removed ) + return; // a GroupOnFilter holding this filter is notified automatically -CORBA::Boolean LessThan_i::IsSatisfy( CORBA::Long theId ) -{ - return myFunctor != 0 && myFunctor->GetValue( theId ) < myMargin; -} + if ( myPredicate ) + { + std::vector predicates; + getPrediacates( myPredicate, predicates ); + for ( size_t i = 0; i < predicates.size(); ++i ) + if ( BelongToMeshGroup_i* bmg = dynamic_cast< BelongToMeshGroup_i* >( predicates[i] )) + { + SMESH::SMESH_GroupBase_var g = bmg->GetGroup(); + SMESH_GroupBase_i* g_i = SMESH::DownCast< SMESH_GroupBase_i*>( g ); + if ( g_i == group ) + { + bmg->SetGroup( SMESH::SMESH_GroupBase::_nil() ); + bmg->SetGroupID( "" ); + } + } + } -/* - Class : MoreThan_i - Description : Comparator ">" -*/ + FindBaseObjects(); // release and update myBaseGroups; + } -CORBA::Boolean MoreThan_i::IsSatisfy( CORBA::Long theId ) -{ - return myFunctor != 0 && myFunctor->GetValue( theId ) > myMargin; -} + /* + FILTER LIBRARY + */ + + #define ATTR_TYPE "type" + #define ATTR_COMPARE "compare" + #define ATTR_THRESHOLD "threshold" + #define ATTR_UNARY "unary" + #define ATTR_BINARY "binary" + #define ATTR_THRESHOLD_STR "threshold_str" + #define ATTR_TOLERANCE "tolerance" + #define ATTR_ELEMENT_TYPE "ElementType" + + //======================================================================= + // name : toString + // Purpose : Convert bool to LDOMString + //======================================================================= + static inline LDOMString toString( CORBA::Boolean val ) + { + return val ? "logical not" : ""; + } -/* - Class : EqualTo_i - Description : Comparator "=" -*/ -EqualTo_i::EqualTo_i() -{ - myToler = Precision::Confusion(); -} + //======================================================================= + // name : toBool + // Purpose : Convert LDOMString to bool + //======================================================================= + static inline bool toBool( const LDOMString& theStr ) + { + return theStr.equals( "logical not" ); + } -CORBA::Boolean EqualTo_i::IsSatisfy( CORBA::Long theId ) -{ - return myFunctor != 0 && fabs( myFunctor->GetValue( theId ) - myMargin ) < myToler; -} + //======================================================================= + // name : toString + // Purpose : Convert double to LDOMString + //======================================================================= + static inline LDOMString toString( CORBA::Double val ) + { + char a[ 255 ]; + sprintf( a, "%e", val ); + return LDOMString( a ); + } -void EqualTo_i::SetTolerance( CORBA::Double theToler ) -{ - myToler = theToler; -} + //======================================================================= + // name : toDouble + // Purpose : Convert LDOMString to double + //======================================================================= + static inline double toDouble( const LDOMString& theStr ) + { + return atof( theStr.GetString() ); + } + //======================================================================= + // name : toString + // Purpose : Convert functor type to LDOMString + //======================================================================= + static inline LDOMString toString( CORBA::Long theType ) + { + switch ( theType ) + { + case FT_AspectRatio : return "Aspect ratio"; + case FT_Warping : return "Warping"; + case FT_Warping3D : return "Warping 3D"; + case FT_MinimumAngle : return "Minimum angle"; + case FT_Taper : return "Taper"; + case FT_Skew : return "Skew"; + case FT_Area : return "Area"; + case FT_Volume3D : return "Volume3D"; + case FT_ScaledJacobian : return "ScaledJacobian"; + case FT_MaxElementLength2D : return "Max element length 2D"; + case FT_MaxElementLength3D : return "Max element length 3D"; + case FT_BelongToMeshGroup : return "Belong to Mesh Group"; + case FT_BelongToGeom : return "Belong to Geom"; + case FT_BelongToPlane : return "Belong to Plane"; + case FT_BelongToCylinder : return "Belong to Cylinder"; + case FT_BelongToGenSurface : return "Belong to Generic Surface"; + case FT_LyingOnGeom : return "Lying on Geom"; + case FT_BadOrientedVolume : return "Bad Oriented Volume"; + case FT_BareBorderVolume : return "Volumes with bare border"; + case FT_BareBorderFace : return "Faces with bare border"; + case FT_OverConstrainedVolume : return "Over-constrained Volumes"; + case FT_OverConstrainedFace : return "Over-constrained Faces"; + case FT_RangeOfIds : return "Range of IDs"; + case FT_FreeBorders : return "Free borders"; + case FT_FreeEdges : return "Free edges"; + case FT_FreeFaces : return "Free faces"; + case FT_FreeNodes : return "Free nodes"; + case FT_EqualNodes : return "Equal nodes"; + case FT_EqualEdges : return "Equal edges"; + case FT_EqualFaces : return "Equal faces"; + case FT_EqualVolumes : return "Equal volumes"; + case FT_MultiConnection : return "Borders at multi-connections"; + case FT_MultiConnection2D : return "Borders at multi-connections 2D"; + case FT_Length : return "Length"; + case FT_Length2D : return "Length 2D"; + case FT_Length3D : return "Length 3D"; + case FT_Deflection2D : return "Deflection 2D"; + case FT_LessThan : return "Less than"; + case FT_MoreThan : return "More than"; + case FT_EqualTo : return "Equal to"; + case FT_LogicalNOT : return "Not"; + case FT_LogicalAND : return "And"; + case FT_LogicalOR : return "Or"; + case FT_GroupColor : return "Color of Group"; + case FT_LinearOrQuadratic : return "Linear or Quadratic"; + case FT_ElemGeomType : return "Element geometry type"; + case FT_EntityType : return "Entity type"; + case FT_Undefined : return ""; + default : return ""; + } + } -/* - Class : LogicalNOT_i - Description : Logical NOT predicate -*/ + //======================================================================= + // name : toFunctorType + // Purpose : Convert LDOMString to functor type + //======================================================================= + static inline SMESH::FunctorType toFunctorType( const LDOMString& theStr ) + { + if ( theStr.equals( "Aspect ratio" ) ) return FT_AspectRatio; + else if ( theStr.equals( "Warping" ) ) return FT_Warping; + else if ( theStr.equals( "Warping 3D" ) ) return FT_Warping3D; + else if ( theStr.equals( "Minimum angle" ) ) return FT_MinimumAngle; + else if ( theStr.equals( "Taper" ) ) return FT_Taper; + else if ( theStr.equals( "Skew" ) ) return FT_Skew; + else if ( theStr.equals( "Area" ) ) return FT_Area; + else if ( theStr.equals( "Volume3D" ) ) return FT_Volume3D; + else if ( theStr.equals( "ScaledJacobian" ) ) return FT_ScaledJacobian; + else if ( theStr.equals( "Max element length 2D" ) ) return FT_MaxElementLength2D; + else if ( theStr.equals( "Max element length 3D" ) ) return FT_MaxElementLength3D; + else if ( theStr.equals( "Belong to Mesh Group" ) ) return FT_BelongToMeshGroup; + else if ( theStr.equals( "Belong to Geom" ) ) return FT_BelongToGeom; + else if ( theStr.equals( "Belong to Plane" ) ) return FT_BelongToPlane; + else if ( theStr.equals( "Belong to Cylinder" ) ) return FT_BelongToCylinder; + else if ( theStr.equals( "Belong to Generic Surface" ) ) return FT_BelongToGenSurface; + else if ( theStr.equals( "Lying on Geom" ) ) return FT_LyingOnGeom; + else if ( theStr.equals( "Free borders" ) ) return FT_FreeBorders; + else if ( theStr.equals( "Free edges" ) ) return FT_FreeEdges; + else if ( theStr.equals( "Free faces" ) ) return FT_FreeFaces; + else if ( theStr.equals( "Free nodes" ) ) return FT_FreeNodes; + else if ( theStr.equals( "Equal nodes" ) ) return FT_EqualNodes; + else if ( theStr.equals( "Equal edges" ) ) return FT_EqualEdges; + else if ( theStr.equals( "Equal faces" ) ) return FT_EqualFaces; + else if ( theStr.equals( "Equal volumes" ) ) return FT_EqualVolumes; + else if ( theStr.equals( "Borders at multi-connections" ) ) return FT_MultiConnection; + // else if ( theStr.equals( "Borders at multi-connections 2D" ) ) return FT_MultiConnection2D; + else if ( theStr.equals( "Length" ) ) return FT_Length; + // else if ( theStr.equals( "Length2D" ) ) return FT_Length2D; + // else if ( theStr.equals( "Length3D" ) ) return FT_Length3D; + else if ( theStr.equals( "Deflection" ) ) return FT_Deflection2D; + else if ( theStr.equals( "Range of IDs" ) ) return FT_RangeOfIds; + else if ( theStr.equals( "Bad Oriented Volume" ) ) return FT_BadOrientedVolume; + else if ( theStr.equals( "Volumes with bare border" ) ) return FT_BareBorderVolume; + else if ( theStr.equals( "Faces with bare border" ) ) return FT_BareBorderFace; + else if ( theStr.equals( "Over-constrained Volumes" ) ) return FT_OverConstrainedVolume; + else if ( theStr.equals( "Over-constrained Faces" ) ) return FT_OverConstrainedFace; + else if ( theStr.equals( "Less than" ) ) return FT_LessThan; + else if ( theStr.equals( "More than" ) ) return FT_MoreThan; + else if ( theStr.equals( "Equal to" ) ) return FT_EqualTo; + else if ( theStr.equals( "Not" ) ) return FT_LogicalNOT; + else if ( theStr.equals( "And" ) ) return FT_LogicalAND; + else if ( theStr.equals( "Or" ) ) return FT_LogicalOR; + else if ( theStr.equals( "Color of Group" ) ) return FT_GroupColor; + else if ( theStr.equals( "Linear or Quadratic" ) ) return FT_LinearOrQuadratic; + else if ( theStr.equals( "Element geometry type" ) ) return FT_ElemGeomType; + else if ( theStr.equals( "Entity type" ) ) return FT_EntityType; + else if ( theStr.equals( "" ) ) return FT_Undefined; + else return FT_Undefined; + } -LogicalNOT_i::LogicalNOT_i() -{ - myPredicate = 0; -} + //======================================================================= + // name : toFunctorType + // Purpose : Convert LDOMString to value of ElementType enumeration + //======================================================================= + static inline SMESH::ElementType toElementType( const LDOMString& theStr ) + { + if ( theStr.equals( "NODE" ) ) return SMESH::NODE; + else if ( theStr.equals( "EDGE" ) ) return SMESH::EDGE; + else if ( theStr.equals( "FACE" ) ) return SMESH::FACE; + else if ( theStr.equals( "VOLUME" ) ) return SMESH::VOLUME; + else return SMESH::ALL; + } -LogicalNOT_i::~LogicalNOT_i() -{ - if ( myPredicate ) - myPredicate->Destroy(); -} + //======================================================================= + // name : toString + // Purpose : Convert ElementType to string + //======================================================================= + static inline LDOMString toString( const SMESH::ElementType theType ) + { + switch ( theType ) + { + case SMESH::NODE : return "NODE"; + case SMESH::EDGE : return "EDGE"; + case SMESH::FACE : return "FACE"; + case SMESH::VOLUME : return "VOLUME"; + case SMESH::ALL : return "ALL"; + default : return ""; + } + } -CORBA::Boolean LogicalNOT_i::IsSatisfy( CORBA::Long theId ) -{ - return myPredicate !=0 && !myPredicate->IsSatisfy( theId ); -} + //======================================================================= + // name : findFilter + // Purpose : Find filter in document + //======================================================================= + static LDOM_Element findFilter( const char* theFilterName, + const LDOM_Document& theDoc, + LDOM_Node* theParent = 0 ) + { + LDOM_Element aRootElement = theDoc.getDocumentElement(); + if ( aRootElement.isNull() || !aRootElement.hasChildNodes() ) + return LDOM_Element(); -void LogicalNOT_i::SetMesh( SMESH_Mesh_ptr theMesh ) -{ - if ( myPredicate != 0 ) - myPredicate->SetMesh( theMesh ); -} + for ( LDOM_Node aTypeNode = aRootElement.getFirstChild(); + !aTypeNode.isNull(); aTypeNode = aTypeNode.getNextSibling() ) + { + for ( LDOM_Node aFilter = aTypeNode.getFirstChild(); + !aFilter.isNull(); aFilter = aFilter.getNextSibling() ) + { + LDOM_Element* anElem = ( LDOM_Element* )&aFilter; + if ( anElem->getTagName().equals( LDOMString( "filter" ) ) && + anElem->getAttribute( "name" ).equals( LDOMString( theFilterName ) ) ) + { + if ( theParent != 0 ) + *theParent = aTypeNode; + return (LDOM_Element&)aFilter; + } + } + } + return LDOM_Element(); + } -void LogicalNOT_i::SetPredicate( Predicate_ptr thePred ) -{ - if ( myPredicate != 0 ) - myPredicate->Destroy(); + //======================================================================= + // name : getSectionName + // Purpose : Get name of section of filters + //======================================================================= + static const char* getSectionName( const ElementType theType ) + { + switch ( theType ) + { + case SMESH::NODE : return "Filters for nodes"; + case SMESH::EDGE : return "Filters for edges"; + case SMESH::FACE : return "Filters for faces"; + case SMESH::VOLUME : return "Filters for volumes"; + case SMESH::ALL : return "Filters for elements"; + default : return ""; + } + } - myPredicate = dynamic_cast( SMESH_Gen_i::GetServant( thePred ).in() ); + //======================================================================= + // name : getSection + // Purpose : Create section for filters corresponding to the entity type + //======================================================================= + static LDOM_Node getSection( const ElementType theType, + LDOM_Document& theDoc, + const bool toCreate = false ) + { + LDOM_Element aRootElement = theDoc.getDocumentElement(); + if ( aRootElement.isNull() ) + return LDOM_Node(); + + // Find section + bool anExist = false; + const char* aSectionName = getSectionName( theType ); + if ( strcmp( aSectionName, "" ) == 0 ) + return LDOM_Node(); + + LDOM_NodeList aSections = theDoc.getElementsByTagName( "section" ); + LDOM_Node aNode; + for ( int i = 0, n = aSections.getLength(); i < n; i++ ) + { + aNode = aSections.item( i ); + LDOM_Element& anItem = ( LDOM_Element& )aNode; + if ( anItem.getAttribute( "name" ).equals( LDOMString( aSectionName ) ) ) + { + anExist = true; + break; + } + } - if ( myPredicate != 0 ) - myPredicate->Register(); -} + // Create new section if necessary + if ( !anExist ) + { + if ( toCreate ) + { + LDOM_Element aNewItem = theDoc.createElement( "section" ); + aNewItem.setAttribute( "name", aSectionName ); + aRootElement.appendChild( aNewItem ); + return aNewItem; + } + else + return LDOM_Node(); + } + return + aNode; + } -int LogicalNOT_i::GetType() const -{ - return myPredicate != 0 ? myPredicate->GetType() : SMDSAbs_All; -} + //======================================================================= + // name : createFilterItem + // Purpose : Create filter item or LDOM document + //======================================================================= + static LDOM_Element createFilterItem( const char* theName, + SMESH::Filter_ptr theFilter, + LDOM_Document& theDoc ) + { + // create new filter in document + LDOM_Element aFilterItem = theDoc.createElement( "filter" ); + aFilterItem.setAttribute( "name", theName ); + // save filter criterions + SMESH::Filter::Criteria_var aCriteria = new SMESH::Filter::Criteria; -/* - Class : LogicalBinary_i - Description : Base class for binary logical predicate -*/ + if ( !theFilter->GetCriteria( aCriteria ) ) + return LDOM_Element(); -LogicalBinary_i::LogicalBinary_i() -{ - myPredicate1 = 0; - myPredicate2 = 0; -} -LogicalBinary_i::~LogicalBinary_i() -{ - if ( myPredicate1 != 0 ) - myPredicate1->Destroy(); + for ( CORBA::ULong i = 0, n = aCriteria->length(); i < n; i++ ) + { + LDOM_Element aCriterionItem = theDoc.createElement( "criterion" ); + + aCriterionItem.setAttribute( ATTR_TYPE , toString( aCriteria[ i ].Type) ); + aCriterionItem.setAttribute( ATTR_COMPARE , toString( aCriteria[ i ].Compare ) ); + aCriterionItem.setAttribute( ATTR_THRESHOLD , toString( aCriteria[ i ].Threshold ) ); + aCriterionItem.setAttribute( ATTR_UNARY , toString( aCriteria[ i ].UnaryOp ) ); + aCriterionItem.setAttribute( ATTR_BINARY , toString( aCriteria[ i ].BinaryOp ) ); + + aCriterionItem.setAttribute( ATTR_THRESHOLD_STR, (const char*)aCriteria[ i ].ThresholdStr ); + aCriterionItem.setAttribute( ATTR_TOLERANCE , toString( aCriteria[ i ].Tolerance ) ); + aCriterionItem.setAttribute( ATTR_ELEMENT_TYPE , + toString( (SMESH::ElementType)aCriteria[ i ].TypeOfElement ) ); + + aFilterItem.appendChild( aCriterionItem ); + } - if ( myPredicate2 != 0 ) - myPredicate2->Destroy(); -} + return aFilterItem; + } -void LogicalBinary_i::SetMesh( SMESH_Mesh_ptr theMesh ) -{ - if ( myPredicate1 != 0 ) - myPredicate1->SetMesh( theMesh ); + //======================================================================= + // name : FilterLibrary_i::FilterLibrary_i + // Purpose : Constructor + //======================================================================= + FilterLibrary_i::FilterLibrary_i( const char* theFileName ) + { + myFileName = CORBA::string_dup( theFileName ); + SMESH::FilterManager_i* aFilterMgr = new SMESH::FilterManager_i(); + myFilterMgr = aFilterMgr->_this(); - if ( myPredicate2 != 0 ) - myPredicate2->SetMesh( theMesh ); -} + LDOMParser aParser; -void LogicalBinary_i::SetPredicate1( Predicate_ptr thePredicate ) -{ - if ( myPredicate1 != 0 ) - myPredicate1->Destroy(); + // Try to use existing library file + bool anExists = false; + if ( !aParser.parse( myFileName ) ) + { + myDoc = aParser.getDocument(); + anExists = true; + } + // Create a new XML document if it doesn't exist + else + myDoc = LDOM_Document::createDocument( LDOMString() ); - myPredicate1 = dynamic_cast( SMESH_Gen_i::GetServant( thePredicate ).in() ); + LDOM_Element aRootElement = myDoc.getDocumentElement(); + if ( aRootElement.isNull() ) + { + // If the existing document is empty --> try to create a new one + if ( anExists ) + myDoc = LDOM_Document::createDocument( LDOMString() ); + } + } - if ( myPredicate1 != 0 ) - myPredicate1->Register(); -} + //======================================================================= + // name : FilterLibrary_i::FilterLibrary_i + // Purpose : Constructor + //======================================================================= + FilterLibrary_i::FilterLibrary_i() + { + myFileName = 0; + SMESH::FilterManager_i* aFilter = new SMESH::FilterManager_i(); + myFilterMgr = aFilter->_this(); -void LogicalBinary_i::SetPredicate2( Predicate_ptr thePredicate ) -{ - if ( myPredicate2 != 0 ) - myPredicate2->Destroy(); + myDoc = LDOM_Document::createDocument( LDOMString() ); + } - myPredicate2 = dynamic_cast( SMESH_Gen_i::GetServant( thePredicate ).in() ); + FilterLibrary_i::~FilterLibrary_i() + { + CORBA::string_free( myFileName ); + //TPythonDump()<Register(); -} + //======================================================================= + // name : FilterLibrary_i::Copy + // Purpose : Create filter and initialize it with values from library + //======================================================================= + Filter_ptr FilterLibrary_i::Copy( const char* theFilterName ) + { + Filter_ptr aRes = Filter::_nil(); + LDOM_Node aFilter = findFilter( theFilterName, myDoc ); -int LogicalBinary_i::GetType() const -{ - if ( myPredicate1 == 0 || myPredicate2 == 0 ) - return SMDSAbs_All; + if ( aFilter.isNull() ) + return aRes; - int aType1 = myPredicate1->GetType(); - int aType2 = myPredicate2->GetType(); + std::list aCriteria; - return aType1 == aType2 ? aType1 : SMDSAbs_All; -} + for ( LDOM_Node aCritNode = aFilter.getFirstChild(); + !aCritNode.isNull() ; aCritNode = aCritNode.getNextSibling() ) + { + LDOM_Element* aCrit = (LDOM_Element*)&aCritNode; -/* - Class : LogicalAND_i - Description : Logical AND -*/ + const char* aTypeStr = aCrit->getAttribute( ATTR_TYPE ).GetString(); + const char* aCompareStr = aCrit->getAttribute( ATTR_COMPARE ).GetString(); + const char* aUnaryStr = aCrit->getAttribute( ATTR_UNARY ).GetString(); + const char* aBinaryStr = aCrit->getAttribute( ATTR_BINARY ).GetString(); + const char* anElemTypeStr = aCrit->getAttribute( ATTR_ELEMENT_TYPE ).GetString(); -CORBA::Boolean LogicalAND_i::IsSatisfy( CORBA::Long theId ) -{ - return myPredicate1 != 0 && - myPredicate2 != 0 && - myPredicate1->IsSatisfy( theId ) && myPredicate2->IsSatisfy( theId );; -} + SMESH::Filter::Criterion aCriterion = createCriterion(); -/* - Class : LogicalOR_i - Description : Logical OR -*/ + aCriterion.Type = toFunctorType( aTypeStr ); + aCriterion.Compare = toFunctorType( aCompareStr ); + aCriterion.UnaryOp = toFunctorType( aUnaryStr ); + aCriterion.BinaryOp = toFunctorType( aBinaryStr ); -CORBA::Boolean LogicalOR_i::IsSatisfy( CORBA::Long theId ) -{ - return myPredicate1 != 0 && - myPredicate2 != 0 && - myPredicate1->IsSatisfy( theId ) || myPredicate2->IsSatisfy( theId ); -} + aCriterion.TypeOfElement = toElementType( anElemTypeStr ); + LDOMString str = aCrit->getAttribute( ATTR_THRESHOLD ); + int val = 0; + aCriterion.Threshold = str.Type() == LDOMBasicString::LDOM_Integer && str.GetInteger( val ) + ? val : atof( str.GetString() ); -/* - FILTER -*/ + str = aCrit->getAttribute( ATTR_TOLERANCE ); + aCriterion.Tolerance = str.Type() == LDOMBasicString::LDOM_Integer && str.GetInteger( val ) + ? val : atof( str.GetString() ); -Filter_i::Filter_i() -{ - myPredicate = 0; -} + str = aCrit->getAttribute( ATTR_THRESHOLD_STR ); + if ( str.Type() == LDOMBasicString::LDOM_Integer && str.GetInteger( val ) ) + { + char a[ 255 ]; + sprintf( a, "%d", val ); + aCriterion.ThresholdStr = CORBA::string_dup( a ); + } + else + aCriterion.ThresholdStr = str.GetString(); -Filter_i::~Filter_i() -{ - if ( myPredicate != 0 ) - myPredicate->Destroy(); -} + aCriteria.push_back( aCriterion ); + } -void Filter_i::SetPredicate( Predicate_ptr thePredicate ) -{ - if ( myPredicate != 0 ) - myPredicate->Destroy(); + SMESH::Filter::Criteria_var aCriteriaVar = new SMESH::Filter::Criteria; + aCriteriaVar->length( aCriteria.size() ); - myPredicate = dynamic_cast( SMESH_Gen_i::GetServant( thePredicate ).in() ); + CORBA::ULong i = 0; + std::list::iterator anIter = aCriteria.begin(); - if ( myPredicate != 0 ) - myPredicate->Register(); -} + for( ; anIter != aCriteria.end(); ++anIter ) + aCriteriaVar[ i++ ] = *anIter; -void Filter_i::SetMesh( SMESH_Mesh_ptr theMesh ) -{ - if ( myPredicate != 0 ) - myPredicate->SetMesh( theMesh ); -} + aRes = myFilterMgr->CreateFilter(); + aRes->SetCriteria( aCriteriaVar.inout() ); -SMESH::long_array* Filter_i::GetElementsId( SMESH_Mesh_ptr theMesh ) -{ + TPythonDump()<( SMESH_Gen_i::GetServant( theMesh ).in() ); + //======================================================================= + // name : FilterLibrary_i::SetFileName + // Purpose : Set file name for library + //======================================================================= + void FilterLibrary_i::SetFileName( const char* theFileName ) + { + CORBA::string_free( myFileName ); + myFileName = CORBA::string_dup( theFileName ); + TPythonDump()<GetImpl().GetMeshDS(); + // if filter already in library or entry filter is null do nothing + LDOM_Node aFilterNode = findFilter( theFilterName, myDoc ); + if ( !aFilterNode.isNull() || theFilter->_is_nil() ) + return false; + + // get section corresponding to the filter type + ElementType anEntType = theFilter->GetElementType(); + + LDOM_Node aSection = getSection( anEntType, myDoc, true ); + if ( aSection.isNull() ) + return false; + + // create filter item + LDOM_Element aFilterItem = createFilterItem( theFilterName, theFilter, myDoc ); + if ( aFilterItem.isNull() ) + return false; + else + { + aSection.appendChild( aFilterItem ); + if(Filter_i* aFilter = DownCast(theFilter)) + TPythonDump()<CreateFilter(); + + LDOM_Element aFilterItem = createFilterItem( theFilterName, aFilter, myDoc ); + if ( aFilterItem.isNull() ) + return false; + else { - int aType = myPredicate->GetType(); - - if ( aType == SMDSAbs_Edge ) - { - SMDS_EdgeIteratorPtr anIter = aMesh->edgesIterator(); - if ( anIter != 0 ) - { - while( anIter->more() ) - { - const SMDS_MeshElement* anElem = anIter->next(); - if ( myPredicate->IsSatisfy( anElem->GetID() ) ) - aList.Append( anElem->GetID() ); - } - } -// delete anIter; - } - else if ( aType == SMDSAbs_Face ) - { - SMDS_FaceIteratorPtr anIter = aMesh->facesIterator(); - if ( anIter != 0 ) - { - while( anIter->more() ) - { - const SMDS_MeshElement* anElem = anIter->next(); - if ( myPredicate->IsSatisfy( anElem->GetID() ) ) - aList.Append( anElem->GetID() ); - } - } -// delete anIter; - } + aSection.appendChild( aFilterItem ); + TPythonDump()<length( aList.Extent() ); - TColStd_ListIteratorOfListOfInteger anIter( aList ); - int i = 0; - for( ; anIter.More(); anIter.Next() ) - anArray[ i++ ] = anIter.Value(); - return anArray._retn(); -} + //======================================================================= + // name : FilterLibrary_i::Delete + // Purpose : Delete filter from library + //======================================================================= + CORBA::Boolean FilterLibrary_i::Delete ( const char* theFilterName ) + { + LDOM_Node aParentNode; + LDOM_Node aFilterNode = findFilter( theFilterName, myDoc, &aParentNode ); + if ( aFilterNode.isNull() || aParentNode.isNull() ) + return false; + + aParentNode.removeChild( aFilterNode ); + TPythonDump()<_is_nil() ) + return false; + + LDOM_Element aNewItem = createFilterItem( theNewName, theFilter, myDoc ); + if ( aNewItem.isNull() ) + return false; + else + { + aFilterItem.ReplaceElement( aNewItem ); + if(Filter_i* aFilter = DownCast(theFilter)) + TPythonDump()<activate_object( this ); -} + //======================================================================= + // name : FilterLibrary_i::Save + // Purpose : Save library on disk + //======================================================================= + CORBA::Boolean FilterLibrary_i::Save() + { + if ( myFileName == 0 || strlen( myFileName ) == 0 ) + return false; -MinimumAngle_ptr FilterManager_i::CreateMinimumAngle() -{ - SMESH::MinimumAngle_i* aServant = new SMESH::MinimumAngle_i(); - SMESH::MinimumAngle_var anObj = aServant->_this(); - return anObj._retn(); -} + std::filebuf fb; + fb.open( myFileName, std::ios::out ); + Standard_OStream os( &fb ); -AspectRatio_ptr FilterManager_i::CreateAspectRatio() -{ - SMESH::AspectRatio_i* aServant = new SMESH::AspectRatio_i(); - SMESH::AspectRatio_var anObj = aServant->_this(); - return anObj._retn(); -} + LDOM_XmlWriter aWriter; + aWriter.SetIndentation( 2 ); + aWriter.Write( os, myDoc ); + fb.close(); + TPythonDump()<_this(); - return anObj._retn(); -} + //======================================================================= + // name : FilterLibrary_i::SaveAs + // Purpose : Save library on disk + //======================================================================= + CORBA::Boolean FilterLibrary_i::SaveAs( const char* aFileName ) + { + myFileName = strdup ( aFileName ); + TPythonDump()<_this(); - return anObj._retn(); -} + //======================================================================= + // name : FilterLibrary_i::NbFilters + // Purpose : Return amount of filters in library + //======================================================================= + CORBA::Long FilterLibrary_i::NbFilters( ElementType theType ) + { + string_array_var aNames = GetNames( theType ); + return aNames->length(); + } + //======================================================================= + // name : FilterLibrary_i::GetNames + // Purpose : Get names of filters from library + //======================================================================= + string_array* FilterLibrary_i::GetNames( ElementType theType ) + { + string_array_var anArray = new string_array; + TColStd_SequenceOfHAsciiString aSeq; -Skew_ptr FilterManager_i::CreateSkew() -{ - SMESH::Skew_i* aServant = new SMESH::Skew_i(); - SMESH::Skew_var anObj = aServant->_this(); - return anObj._retn(); -} + LDOM_Node aSection = getSection( theType, myDoc, false ); + if ( !aSection.isNull() ) + { + for ( LDOM_Node aFilter = aSection.getFirstChild(); + !aFilter.isNull(); aFilter = aFilter.getNextSibling() ) + { + LDOM_Element& anElem = ( LDOM_Element& )aFilter; + aSeq.Append( new TCollection_HAsciiString( + (Standard_CString)anElem.getAttribute( "name" ).GetString() ) ); + } + } -Area_ptr FilterManager_i::CreateArea() -{ - SMESH::Area_i* aServant = new SMESH::Area_i(); - SMESH::Area_var anObj = aServant->_this(); - return anObj._retn(); -} + anArray->length( aSeq.Length() ); + for ( int i = 1, n = aSeq.Length(); i <= n; i++ ) + anArray[ i - 1 ] = CORBA::string_dup( aSeq( i )->ToCString() ); + return anArray._retn(); + } -Length_ptr FilterManager_i::CreateLength() -{ - SMESH::Length_i* aServant = new SMESH::Length_i(); - SMESH::Length_var anObj = aServant->_this(); - return anObj._retn(); -} + //======================================================================= + // name : FilterLibrary_i::GetAllNames + // Purpose : Get names of filters from library + //======================================================================= + string_array* FilterLibrary_i::GetAllNames() + { + string_array_var aResArray = new string_array; + for ( int type = SMESH::ALL; type <= SMESH::VOLUME; type++ ) + { + SMESH::string_array_var aNames = GetNames( (SMESH::ElementType)type ); + int aPrevLength = aResArray->length(); + aResArray->length( aPrevLength + aNames->length() ); + for ( int i = 0, n = aNames->length(); i < n; i++ ) + aResArray[ aPrevLength + i ] = aNames[ i ]; + } -MultiConnection_ptr FilterManager_i::CreateMultiConnection() -{ - SMESH::MultiConnection_i* aServant = new SMESH::MultiConnection_i(); - SMESH::MultiConnection_var anObj = aServant->_this(); - return anObj._retn(); -} + return aResArray._retn(); + } + //================================================================================ + /*! + * \brief Return an array of strings corresponding to items of enum FunctorType + */ + //================================================================================ -FreeBorders_ptr FilterManager_i::CreateFreeBorders() -{ - SMESH::FreeBorders_i* aServant = new SMESH::FreeBorders_i(); - SMESH::FreeBorders_var anObj = aServant->_this(); - return anObj._retn(); -} + static const char** getFunctNames() + { + static const char* functName[] = { + // IT's necessary to update this array according to enum FunctorType (SMESH_Filter.idl) + // The order is IMPORTANT !!! + "FT_AspectRatio", + "FT_AspectRatio3D", + "FT_Warping", + "FT_Warping3D", + "FT_MinimumAngle", + "FT_Taper", + "FT_Skew", + "FT_Area", + "FT_Volume3D", + "FT_ScaledJacobian", + "FT_MaxElementLength2D", + "FT_MaxElementLength3D", + "FT_FreeBorders", + "FT_FreeEdges", + "FT_FreeNodes", + "FT_FreeFaces", + "FT_EqualNodes", + "FT_EqualEdges", + "FT_EqualFaces", + "FT_EqualVolumes", + "FT_MultiConnection", + "FT_MultiConnection2D", + "FT_Length", + "FT_Length2D", + "FT_Length3D", + "FT_Deflection2D", + "FT_NodeConnectivityNumber", + "FT_BelongToMeshGroup", + "FT_BelongToGeom", + "FT_BelongToPlane", + "FT_BelongToCylinder", + "FT_BelongToGenSurface", + "FT_LyingOnGeom", + "FT_RangeOfIds", + "FT_BadOrientedVolume", + "FT_BareBorderVolume", + "FT_BareBorderFace", + "FT_OverConstrainedVolume", + "FT_OverConstrainedFace", + "FT_LinearOrQuadratic", + "FT_GroupColor", + "FT_ElemGeomType", + "FT_EntityType", + "FT_CoplanarFaces", + "FT_BallDiameter", + "FT_ConnectedElements", + "FT_LessThan", + "FT_MoreThan", + "FT_EqualTo", + "FT_LogicalNOT", + "FT_LogicalAND", + "FT_LogicalOR", + "FT_Undefined"}; + + // check if functName is complete, compilation failure means that enum FunctorType changed + static_assert( sizeof(functName) / sizeof(const char*) == SMESH::FT_Undefined + 1, + "Update names of FunctorType's!!!" ); + + return functName; + } -LessThan_ptr FilterManager_i::CreateLessThan() -{ - SMESH::LessThan_i* aServant = new SMESH::LessThan_i(); - SMESH::LessThan_var anObj = aServant->_this(); - return anObj._retn(); -} + //================================================================================ + /*! + * \brief Return a string corresponding to an item of enum FunctorType + */ + //================================================================================ + + const char* FunctorTypeToString(FunctorType ft) + { + if ( ft < 0 || ft > SMESH::FT_Undefined ) + return "FT_Undefined"; + return getFunctNames()[ ft ]; + } + //================================================================================ + /*! + * \brief Converts a string to FunctorType. This is reverse of FunctorTypeToString() + */ + //================================================================================ -MoreThan_ptr FilterManager_i::CreateMoreThan() -{ - SMESH::MoreThan_i* aServant = new SMESH::MoreThan_i(); - SMESH::MoreThan_var anObj = aServant->_this(); - return anObj._retn(); -} + FunctorType StringToFunctorType(const char* str) + { + std::string name( str + 3 ); // skip "FT_" + const char** functNames = getFunctNames(); + int ft = 0; + for ( ; ft < SMESH::FT_Undefined; ++ft ) + if ( name == ( functNames[ft] + 3 )) + break; -EqualTo_ptr FilterManager_i::CreateEqualTo() -{ - SMESH::EqualTo_i* aServant = new SMESH::EqualTo_i(); - SMESH::EqualTo_var anObj = aServant->_this(); - return anObj._retn(); -} + //ASSERT( strcmp( str, FunctorTypeToString( SMESH::FunctorType( ft ))) == 0 ); + return SMESH::FunctorType( ft ); + } -LogicalNOT_ptr FilterManager_i::CreateLogicalNOT() -{ - SMESH::LogicalNOT_i* aServant = new SMESH::LogicalNOT_i(); - SMESH::LogicalNOT_var anObj = aServant->_this(); - return anObj._retn(); -} + //================================================================================ + /*! + * \brief calls OnBaseObjModified(), if who != this, and myWaiters[i]->Modified(who) + */ + //================================================================================ + void NotifyerAndWaiter::Modified( bool removed, NotifyerAndWaiter* who ) + { + if ( who != 0 && who != this ) + OnBaseObjModified( who, removed ); + else + who = this; + + std::list waiters = myWaiters; // myWaiters can be changed by Modified() + std::list::iterator i = waiters.begin(); + for ( ; i != waiters.end(); ++i ) + (*i)->Modified( removed, who ); + } -LogicalAND_ptr FilterManager_i::CreateLogicalAND() -{ - SMESH::LogicalAND_i* aServant = new SMESH::LogicalAND_i(); - SMESH::LogicalAND_var anObj = aServant->_this(); - return anObj._retn(); -} + //================================================================================ + /*! + * \brief Stores an object to be notified on change of predicate + */ + //================================================================================ + void NotifyerAndWaiter::AddModifWaiter( NotifyerAndWaiter* waiter ) + { + if ( waiter ) + myWaiters.push_back( waiter ); + } -LogicalOR_ptr FilterManager_i::CreateLogicalOR() -{ - SMESH::LogicalOR_i* aServant = new SMESH::LogicalOR_i(); - SMESH::LogicalOR_var anObj = aServant->_this(); - return anObj._retn(); -} + //================================================================================ + /*! + * \brief Removes an object to be notified on change of predicate + */ + //================================================================================ -Filter_ptr FilterManager_i::CreateFilter() -{ - SMESH::Filter_i* aServant = new SMESH::Filter_i(); - SMESH::Filter_var anObj = aServant->_this(); - return anObj._retn(); + void NotifyerAndWaiter::RemoveModifWaiter( NotifyerAndWaiter* waiter ) + { + myWaiters.remove( waiter ); + } + + //================================================================================ + /*! + * \brief Checks if a waiter is among myWaiters, maybe nested + */ + //================================================================================ + + bool NotifyerAndWaiter::ContainModifWaiter( NotifyerAndWaiter* waiter ) + { + bool is = ( waiter == this ); + + std::list::iterator w = myWaiters.begin(); + for ( ; !is && w != myWaiters.end(); ++w ) + is = (*w)->ContainModifWaiter( waiter ); + + return is; + } }