1 // Copyright (C) 2007-2015 CEA/DEN, EDF R&D, OPEN CASCADE
3 // Copyright (C) 2003-2007 OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
4 // CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
6 // This library is free software; you can redistribute it and/or
7 // modify it under the terms of the GNU Lesser General Public
8 // License as published by the Free Software Foundation; either
9 // version 2.1 of the License, or (at your option) any later version.
11 // This library is distributed in the hope that it will be useful,
12 // but WITHOUT ANY WARRANTY; without even the implied warranty of
13 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
14 // Lesser General Public License for more details.
16 // You should have received a copy of the GNU Lesser General Public
17 // License along with this library; if not, write to the Free Software
18 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
20 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
23 // SMESH SMESH : implementaion of SMESH idl descriptions
24 // File : StdMeshers_NumberOfSegments.cxx
25 // Moved here from SMESH_NumberOfSegments.cxx
26 // Author : Paul RASCLE, EDF
29 #include "StdMeshers_NumberOfSegments.hxx"
31 #include "StdMeshers_Distribution.hxx"
32 #include "SMESHDS_SubMesh.hxx"
33 #include "SMESH_Mesh.hxx"
34 #include "SMESH_Comment.hxx"
36 #include <ExprIntrp_GenExp.hxx>
37 #include <Expr_Array1OfNamedUnknown.hxx>
38 #include <Expr_NamedUnknown.hxx>
39 #include <TColStd_Array1OfReal.hxx>
40 #include <TCollection_AsciiString.hxx>
42 #include <TopTools_IndexedMapOfShape.hxx>
44 #if (OCC_VERSION_MAJOR << 16 | OCC_VERSION_MINOR << 8 | OCC_VERSION_MAINTENANCE) > 0x060100
48 #include <Standard_Failure.hxx>
51 #include <Standard_ErrorHandler.hxx>
54 #include <Basics_Utils.hxx>
56 using namespace StdMeshers;
59 const double PRECISION = 1e-7;
61 //=============================================================================
65 //=============================================================================
67 StdMeshers_NumberOfSegments::StdMeshers_NumberOfSegments(int hypId,
70 : SMESH_Hypothesis(hypId, studyId, gen),
71 _numberOfSegments(15),//issue 19923
72 _distrType(DT_Regular),
74 _convMode(1) //cut negative by default
76 _name = "NumberOfSegments";
80 //=============================================================================
84 //=============================================================================
86 StdMeshers_NumberOfSegments::~StdMeshers_NumberOfSegments()
90 //=============================================================================
94 //=============================================================================
96 StdMeshers_NumberOfSegments::BuildDistributionExpr( const char* expr,int nbSeg,int conv )
97 throw ( SALOME_Exception )
99 if( !buildDistribution( TCollection_AsciiString( ( Standard_CString )expr ), conv, 0.0, 1.0, nbSeg, _distr, 1E-4 ) )
104 const vector<double>&
105 StdMeshers_NumberOfSegments::BuildDistributionTab( const vector<double>& tab,
108 throw ( SALOME_Exception )
110 if( !buildDistribution( tab, conv, 0.0, 1.0, nbSeg, _distr, 1E-4 ) )
115 //=============================================================================
119 //=============================================================================
121 void StdMeshers_NumberOfSegments::SetNumberOfSegments(int segmentsNumber)
122 throw(SALOME_Exception)
124 int oldNumberOfSegments = _numberOfSegments;
125 if (segmentsNumber <= 0)
126 throw SALOME_Exception(LOCALIZED("number of segments must be positive"));
127 _numberOfSegments = segmentsNumber;
129 if (oldNumberOfSegments != _numberOfSegments)
130 NotifySubMeshesHypothesisModification();
133 //=============================================================================
137 //=============================================================================
139 int StdMeshers_NumberOfSegments::GetNumberOfSegments() const
141 return _numberOfSegments;
144 //================================================================================
148 //================================================================================
150 void StdMeshers_NumberOfSegments::SetDistrType(DistrType typ)
151 throw(SALOME_Exception)
153 if (typ < DT_Regular || typ > DT_ExprFunc)
154 throw SALOME_Exception(LOCALIZED("distribution type is out of range"));
156 if (typ != _distrType)
159 NotifySubMeshesHypothesisModification();
163 //================================================================================
167 //================================================================================
169 StdMeshers_NumberOfSegments::DistrType StdMeshers_NumberOfSegments::GetDistrType() const
174 //================================================================================
178 //================================================================================
180 void StdMeshers_NumberOfSegments::SetScaleFactor(double scaleFactor)
181 throw(SALOME_Exception)
183 if (scaleFactor < PRECISION)
184 throw SALOME_Exception(LOCALIZED("scale factor must be positive"));
186 if (_distrType != DT_Scale)
187 _distrType = DT_Scale;
189 // commented by mpa for IPAL 52986
190 // if ( fabs(scaleFactor - 1.0) < PRECISION )
191 // _distrType = DT_Regular;
193 if ( fabs(_scaleFactor - scaleFactor) > PRECISION )
195 _scaleFactor = scaleFactor;
196 NotifySubMeshesHypothesisModification();
200 //================================================================================
204 //================================================================================
206 double StdMeshers_NumberOfSegments::GetScaleFactor() const
207 throw(SALOME_Exception)
209 if (_distrType != DT_Scale)
210 throw SALOME_Exception(LOCALIZED("not a scale distribution"));
214 //================================================================================
218 //================================================================================
220 void StdMeshers_NumberOfSegments::SetTableFunction(const vector<double>& table)
221 throw(SALOME_Exception)
223 if (_distrType != DT_TabFunc)
224 _distrType = DT_TabFunc;
225 //throw SALOME_Exception(LOCALIZED("not a table function distribution"));
226 if ( (table.size() % 2) != 0 )
227 throw SALOME_Exception(LOCALIZED("odd size of vector of table function"));
229 double prev = -PRECISION;
230 bool isSame = table.size() == _table.size();
233 for ( size_t i = 0; i < table.size() / 2; i++ )
235 double par = table[i*2];
236 double val = table[i*2+1];
243 val = pow( 10.0, val );
245 catch(Standard_Failure) {
246 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
247 throw SALOME_Exception( LOCALIZED( "invalid value"));
251 else if( _convMode==1 && val<0.0 )
254 if ( par < 0 || par > 1)
255 throw SALOME_Exception(LOCALIZED("parameter of table function is out of range [0,1]"));
256 if ( fabs(par-prev) < PRECISION )
257 throw SALOME_Exception(LOCALIZED("two parameters are the same"));
259 throw SALOME_Exception(LOCALIZED("value of table function is not positive"));
260 if( val > PRECISION )
264 double oldpar = _table[i*2];
265 double oldval = _table[i*2+1];
266 if ( fabs(par - oldpar) > PRECISION || fabs(val - oldval) > PRECISION )
273 throw SALOME_Exception(LOCALIZED("value of table function is not positive"));
275 if ( pos && !isSame )
278 NotifySubMeshesHypothesisModification();
282 //================================================================================
286 //================================================================================
288 const vector<double>& StdMeshers_NumberOfSegments::GetTableFunction() const
289 throw(SALOME_Exception)
291 if (_distrType != DT_TabFunc)
292 throw SALOME_Exception(LOCALIZED("not a table function distribution"));
296 //================================================================================
297 /*! check if only 't' is unknown variable in expression
299 //================================================================================
300 bool isCorrectArg( const Handle( Expr_GeneralExpression )& expr )
302 Handle( Expr_NamedUnknown ) sub = Handle( Expr_NamedUnknown )::DownCast( expr );
304 return sub->GetName()=="t";
307 for( int i=1, n=expr->NbSubExpressions(); i<=n && res; i++ )
309 Handle( Expr_GeneralExpression ) sub = expr->SubExpression( i );
310 Handle( Expr_NamedUnknown ) name = Handle( Expr_NamedUnknown )::DownCast( sub );
313 if( name->GetName()!="t" )
317 res = isCorrectArg( sub );
322 //================================================================================
323 /*! this function parses the expression 'str' in order to check if syntax is correct
324 * ( result in 'syntax' ) and if only 't' is unknown variable in expression ( result in 'args' )
326 //================================================================================
327 bool process( const TCollection_AsciiString& str, int convMode,
328 bool& syntax, bool& args,
329 bool& non_neg, bool& non_zero,
330 bool& singulars, double& sing_point )
332 Kernel_Utils::Localizer loc;
334 bool parsed_ok = true;
335 Handle( ExprIntrp_GenExp ) myExpr;
340 myExpr = ExprIntrp_GenExp::Create();
341 myExpr->Process( str.ToCString() );
342 } catch(Standard_Failure) {
343 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
349 if( parsed_ok && myExpr->IsDone() )
352 args = isCorrectArg( myExpr->Expression() );
355 bool res = parsed_ok && syntax && args;
365 FunctionExpr f( str.ToCString(), convMode );
367 for( int i=0; i<=max; i++ )
369 double t = double(i)/double(max), val;
370 if( !f.value( t, val ) )
386 return res && non_neg && non_zero && ( !singulars );
389 //================================================================================
393 //================================================================================
395 void StdMeshers_NumberOfSegments::SetExpressionFunction(const char* expr)
396 throw(SALOME_Exception)
398 if (_distrType != DT_ExprFunc)
399 _distrType = DT_ExprFunc;
401 string func = CheckExpressionFunction( expr, _convMode );
405 NotifySubMeshesHypothesisModification();
409 //=======================================================================
410 //function : CheckExpressionFunction
411 //purpose : Checks validity of the expression of the function f(t), e.g. "sin(t)".
412 // In case of validity returns a cleaned expression
413 //=======================================================================
416 StdMeshers_NumberOfSegments::CheckExpressionFunction( const std::string& expr,
418 throw (SALOME_Exception)
420 // remove white spaces
421 TCollection_AsciiString str((Standard_CString)expr.c_str());
427 bool syntax, args, non_neg, singulars, non_zero;
429 bool res = process( str, convMode, syntax, args, non_neg, non_zero, singulars, sing_point );
433 throw SALOME_Exception(SMESH_Comment("invalid expression syntax: ") << str );
435 throw SALOME_Exception(LOCALIZED("only 't' may be used as function argument"));
437 throw SALOME_Exception(LOCALIZED("only non-negative function can be used"));
441 sprintf( buf, "Function has singular point in %.3f", sing_point );
442 throw SALOME_Exception( buf );
445 throw SALOME_Exception(LOCALIZED("f(t)=0 cannot be used"));
448 return str.ToCString();
451 //================================================================================
455 //================================================================================
457 const char* StdMeshers_NumberOfSegments::GetExpressionFunction() const
458 throw(SALOME_Exception)
460 if (_distrType != DT_ExprFunc)
461 throw SALOME_Exception(LOCALIZED("not an expression function distribution"));
462 return _func.c_str();
465 //================================================================================
469 //================================================================================
471 void StdMeshers_NumberOfSegments::SetConversionMode( int conv )
472 throw(SALOME_Exception)
474 // if (_distrType != DT_TabFunc && _distrType != DT_ExprFunc)
475 // throw SALOME_Exception(LOCALIZED("not a functional distribution"));
477 if( conv != _convMode )
480 NotifySubMeshesHypothesisModification();
484 //================================================================================
488 //================================================================================
490 int StdMeshers_NumberOfSegments::ConversionMode() const
491 throw(SALOME_Exception)
493 // if (_distrType != DT_TabFunc && _distrType != DT_ExprFunc)
494 // throw SALOME_Exception(LOCALIZED("not a functional distribution"));
498 //=============================================================================
502 //=============================================================================
504 ostream & StdMeshers_NumberOfSegments::SaveTo(ostream & save)
506 int listSize = _edgeIDs.size();
507 save << _numberOfSegments << " " << (int)_distrType;
511 save << " " << _scaleFactor;
514 save << " " << _table.size();
515 for ( size_t i = 0; i < _table.size(); i++ )
516 save << " " << _table[i];
519 save << " " << _func;
526 if (_distrType == DT_TabFunc || _distrType == DT_ExprFunc)
527 save << " " << _convMode;
529 if ( _distrType != DT_Regular && listSize > 0 ) {
530 save << " " << listSize;
531 for ( int i = 0; i < listSize; i++ )
532 save << " " << _edgeIDs[i];
533 save << " " << _objEntry;
539 //=============================================================================
543 //=============================================================================
545 istream & StdMeshers_NumberOfSegments::LoadFrom(istream & load)
550 // read number of segments
553 _numberOfSegments = a;
555 load.clear(ios::badbit | load.rdstate());
557 // read second stored value. It can be two variants here:
558 // 1. If the hypothesis is stored in old format (nb.segments and scale factor),
559 // we wait here the scale factor, which is double.
560 // 2. If the hypothesis is stored in new format
561 // (nb.segments, distr.type, some other params.),
562 // we wait here the ditribution type, which is integer
564 isOK = (load >> scale_factor);
565 a = (int)scale_factor;
567 // try to interprete ditribution type,
568 // supposing that this hypothesis was written in the new format
571 if (a < DT_Regular || a > DT_ExprFunc)
572 _distrType = DT_Regular;
574 _distrType = (DistrType) a;
577 load.clear(ios::badbit | load.rdstate());
579 // parameters of distribution
590 load.clear(ios::badbit | load.rdstate());
591 // this can mean, that the hypothesis is stored in old format
592 _distrType = DT_Regular;
593 _scaleFactor = scale_factor;
602 _table.resize(a, 0.);
603 for ( size_t i=0; i < _table.size(); i++ )
609 load.clear(ios::badbit | load.rdstate());
614 load.clear(ios::badbit | load.rdstate());
615 // this can mean, that the hypothesis is stored in old format
616 _distrType = DT_Regular;
617 _scaleFactor = scale_factor;
624 isOK = (load >> str);
629 load.clear(ios::badbit | load.rdstate());
630 // this can mean, that the hypothesis is stored in old format
631 _distrType = DT_Regular;
632 _scaleFactor = scale_factor;
641 if (_distrType == DT_TabFunc || _distrType == DT_ExprFunc)
647 load.clear(ios::badbit | load.rdstate());
650 // load reversed edges IDs
652 isOK = (load >> intVal);
653 if ( isOK && _distrType != DT_Regular && intVal > 0 ) {
654 _edgeIDs.reserve( intVal );
655 for ( size_t i = 0; i < _edgeIDs.capacity() && isOK; i++) {
656 isOK = (load >> intVal);
657 if ( isOK ) _edgeIDs.push_back( intVal );
659 isOK = (load >> _objEntry);
665 //=============================================================================
669 //=============================================================================
671 ostream & operator <<(ostream & save, StdMeshers_NumberOfSegments & hyp)
673 return hyp.SaveTo( save );
676 //=============================================================================
680 //=============================================================================
682 istream & operator >>(istream & load, StdMeshers_NumberOfSegments & hyp)
684 return hyp.LoadFrom( load );
687 //================================================================================
689 * \brief Initialize number of segments by the mesh built on the geometry
690 * \param theMesh - the built mesh
691 * \param theShape - the geometry of interest
692 * \retval bool - true if parameter values have been successfully defined
694 //================================================================================
696 bool StdMeshers_NumberOfSegments::SetParametersByMesh(const SMESH_Mesh* theMesh,
697 const TopoDS_Shape& theShape)
699 if ( !theMesh || theShape.IsNull() )
702 _numberOfSegments = 0;
703 _distrType = DT_Regular;
706 TopTools_IndexedMapOfShape edgeMap;
707 TopExp::MapShapes( theShape, TopAbs_EDGE, edgeMap );
708 SMESHDS_Mesh* aMeshDS = const_cast< SMESH_Mesh* >( theMesh )->GetMeshDS();
709 for ( int i = 1; i <= edgeMap.Extent(); ++i )
711 // get current segment length
712 SMESHDS_SubMesh * eSubMesh = aMeshDS->MeshElements( edgeMap( i ));
713 if ( eSubMesh && eSubMesh->NbElements())
714 _numberOfSegments += eSubMesh->NbElements();
719 _numberOfSegments /= nbEdges;
721 if (_numberOfSegments == 0) _numberOfSegments = 1;
725 //================================================================================
727 * \brief Initialize my parameter values by default parameters.
728 * \retval bool - true if parameter values have been successfully defined
730 //================================================================================
732 bool StdMeshers_NumberOfSegments::SetParametersByDefaults(const TDefaults& dflts,
733 const SMESH_Mesh* /*theMesh*/)
735 return (_numberOfSegments = dflts._nbSegments );
738 //=============================================================================
742 //=============================================================================
744 void StdMeshers_NumberOfSegments::SetReversedEdges( std::vector<int>& ids )
746 if ( ids != _edgeIDs ) {
749 NotifySubMeshesHypothesisModification();