]> SALOME platform Git repositories - modules/smesh.git/blob - src/StdMeshers/StdMeshers_LocalLength.cxx
Salome HOME
0020983: EDF 1569 SMESH: Numbering Display Elements with several groups selected...
[modules/smesh.git] / src / StdMeshers / StdMeshers_LocalLength.cxx
1 //  Copyright (C) 2007-2010  CEA/DEN, EDF R&D, OPEN CASCADE
2 //
3 //  Copyright (C) 2003-2007  OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
4 //  CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
5 //
6 //  This library is free software; you can redistribute it and/or
7 //  modify it under the terms of the GNU Lesser General Public
8 //  License as published by the Free Software Foundation; either
9 //  version 2.1 of the License.
10 //
11 //  This library is distributed in the hope that it will be useful,
12 //  but WITHOUT ANY WARRANTY; without even the implied warranty of
13 //  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
14 //  Lesser General Public License for more details.
15 //
16 //  You should have received a copy of the GNU Lesser General Public
17 //  License along with this library; if not, write to the Free Software
18 //  Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
19 //
20 //  See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
21 //
22
23 //  SMESH SMESH : implementaion of SMESH idl descriptions
24 //  File   : StdMeshers_LocalLength.cxx
25 //           Moved here from SMESH_LocalLength.cxx
26 //  Author : Paul RASCLE, EDF
27 //  Module : SMESH
28 //
29 #include "StdMeshers_LocalLength.hxx"
30
31 #include "SMESH_Mesh.hxx"
32 #include "SMESH_Algo.hxx"
33
34 #include "utilities.h"
35
36 #include <BRep_Tool.hxx>
37 #include <GCPnts_AbscissaPoint.hxx>
38 #include <GeomAdaptor_Curve.hxx>
39 #include <Geom_Curve.hxx>
40 #include <TopExp.hxx>
41 #include <TopLoc_Location.hxx>
42 #include <TopTools_IndexedMapOfShape.hxx>
43 #include <TopoDS.hxx>
44 #include <TopoDS_Edge.hxx>
45 #include <Precision.hxx>
46
47 using namespace std;
48
49 //=============================================================================
50 /*!
51  *  
52  */
53 //=============================================================================
54
55 StdMeshers_LocalLength::StdMeshers_LocalLength(int hypId, int studyId, SMESH_Gen * gen)
56   :SMESH_Hypothesis(hypId, studyId, gen)
57 {
58   _length = 1.;
59   _precision = Precision::Confusion();
60   _name = "LocalLength";
61   _param_algo_dim = 1; // is used by SMESH_Regular_1D
62 }
63
64 //=============================================================================
65 /*!
66  *  
67  */
68 //=============================================================================
69
70 StdMeshers_LocalLength::~StdMeshers_LocalLength()
71 {
72 }
73
74 //=============================================================================
75 /*!
76  *  
77  */
78 //=============================================================================
79
80 void StdMeshers_LocalLength::SetLength(double length) throw(SALOME_Exception)
81 {
82   double oldLength = _length;
83   if (length <= 0)
84     throw SALOME_Exception(LOCALIZED("length must be positive"));
85   _length = length;
86   const double precision = 1e-7;
87   if (fabs(oldLength - _length) > precision)
88     NotifySubMeshesHypothesisModification();
89 }
90
91 //=============================================================================
92 /*!
93  *  
94  */
95 //=============================================================================
96
97 double StdMeshers_LocalLength::GetLength() const
98 {
99   return _length;
100 }
101
102 //=============================================================================
103 /*!
104  *  
105  */
106 //=============================================================================
107 void StdMeshers_LocalLength::SetPrecision (double thePrecision) throw(SALOME_Exception)
108 {
109   double oldPrecision = _precision;
110   if (_precision < 0)
111     throw SALOME_Exception(LOCALIZED("precision cannot be negative"));
112   _precision = thePrecision;
113   const double precision = 1e-8;
114   if (fabs(oldPrecision - _precision) > precision)
115     NotifySubMeshesHypothesisModification();
116 }
117
118 //=============================================================================
119 /*!
120  *  
121  */
122 //=============================================================================
123 double StdMeshers_LocalLength::GetPrecision() const
124 {
125   return _precision;
126 }
127
128 //=============================================================================
129 /*!
130  *  
131  */
132 //=============================================================================
133
134 ostream & StdMeshers_LocalLength::SaveTo(ostream & save)
135 {
136   save << this->_length << " " << this->_precision;
137   return save;
138 }
139
140 //=============================================================================
141 /*!
142  *  
143  */
144 //=============================================================================
145
146 istream & StdMeshers_LocalLength::LoadFrom(istream & load)
147 {
148   bool isOK = true;
149   double a;
150
151   isOK = (load >> a);
152   if (isOK)
153     this->_length = a;
154   else
155     load.clear(ios::badbit | load.rdstate());
156
157   isOK = (load >> a);
158   if (isOK)
159     this->_precision = a;
160   else
161   {
162     load.clear(ios::badbit | load.rdstate());
163     // old format, without precision
164     _precision = 0.;
165   }
166
167   return load;
168 }
169
170 //=============================================================================
171 /*!
172  *  
173  */
174 //=============================================================================
175
176 ostream & operator <<(ostream & save, StdMeshers_LocalLength & hyp)
177 {
178   return hyp.SaveTo( save );
179 }
180
181 //=============================================================================
182 /*!
183  *  
184  */
185 //=============================================================================
186
187 istream & operator >>(istream & load, StdMeshers_LocalLength & hyp)
188 {
189   return hyp.LoadFrom( load );
190 }
191
192 //================================================================================
193 /*!
194  * \brief Initialize segment length by the mesh built on the geometry
195  * \param theMesh - the built mesh
196  * \param theShape - the geometry of interest
197  * \retval bool - true if parameter values have been successfully defined
198  */
199 //================================================================================
200
201 bool StdMeshers_LocalLength::SetParametersByMesh(const SMESH_Mesh*   theMesh,
202                                                  const TopoDS_Shape& theShape)
203 {
204   if ( !theMesh || theShape.IsNull() )
205     return false;
206
207   _length = 0.;
208
209   Standard_Real UMin, UMax;
210   TopLoc_Location L;
211
212   int nbEdges = 0;
213   TopTools_IndexedMapOfShape edgeMap;
214   TopExp::MapShapes( theShape, TopAbs_EDGE, edgeMap );
215   for ( int iE = 1; iE <= edgeMap.Extent(); ++iE )
216   {
217     const TopoDS_Edge& edge = TopoDS::Edge( edgeMap( iE ));
218     Handle(Geom_Curve) C = BRep_Tool::Curve( edge, L, UMin, UMax );
219     GeomAdaptor_Curve AdaptCurve(C);
220
221     vector< double > params;
222     SMESHDS_Mesh* aMeshDS = const_cast< SMESH_Mesh* >( theMesh )->GetMeshDS();
223     if ( SMESH_Algo::GetNodeParamOnEdge( aMeshDS, edge, params ))
224     {
225       for ( int i = 1; i < params.size(); ++i )
226         _length += GCPnts_AbscissaPoint::Length( AdaptCurve, params[ i-1 ], params[ i ]);
227       nbEdges += params.size() - 1;
228     }
229   }
230   if ( nbEdges )
231     _length /= nbEdges;
232
233   _precision = Precision::Confusion();
234
235   return nbEdges;
236 }
237 //================================================================================
238 /*!
239  * \brief Initialize my parameter values by default parameters.
240  *  \retval bool - true if parameter values have been successfully defined
241  */
242 //================================================================================
243
244 bool StdMeshers_LocalLength::SetParametersByDefaults(const TDefaults&  dflts,
245                                                      const SMESH_Mesh* /*theMesh*/)
246 {
247   return ( _length = dflts._elemLength );
248 }
249