1 // SMESH StdMeshers_Deflection1D : implementaion of SMESH idl descriptions
3 // Copyright (C) 2003 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.
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
24 // File : StdMeshers_Deflection1D.cxx
30 #include "StdMeshers_Deflection1D.hxx"
31 #include "utilities.h"
33 #include "SMESH_Mesh.hxx"
34 #include "SMESH_Algo.hxx"
36 #include <BRep_Tool.hxx>
37 #include <GeomAdaptor_Curve.hxx>
38 #include <Geom_Curve.hxx>
40 #include <TopLoc_Location.hxx>
41 #include <TopTools_IndexedMapOfShape.hxx>
43 #include <TopoDS_Edge.hxx>
47 //=============================================================================
51 //=============================================================================
53 StdMeshers_Deflection1D::StdMeshers_Deflection1D(int hypId,
56 :SMESH_Hypothesis(hypId, studyId, gen)
59 _name = "Deflection1D";
60 _param_algo_dim = 1; // is used by SMESH_Regular_1D
63 //=============================================================================
67 //=============================================================================
69 StdMeshers_Deflection1D::~StdMeshers_Deflection1D()
73 //=============================================================================
77 //=============================================================================
79 void StdMeshers_Deflection1D::SetDeflection(double value)
80 throw(SALOME_Exception)
82 if (_value != value) {
84 throw SALOME_Exception(LOCALIZED("Value must be positive"));
86 NotifySubMeshesHypothesisModification();
92 //=============================================================================
96 //=============================================================================
98 double StdMeshers_Deflection1D::GetDeflection() const
103 //=============================================================================
107 //=============================================================================
109 ostream & StdMeshers_Deflection1D::SaveTo(ostream & save)
115 //=============================================================================
119 //=============================================================================
121 istream & StdMeshers_Deflection1D::LoadFrom(istream & load)
123 bool isOK = (load >> _value);
125 load.clear(ios::badbit | load.rdstate());
129 //=============================================================================
133 //=============================================================================
135 ostream & operator <<(ostream & save, StdMeshers_Deflection1D & hyp)
137 return hyp.SaveTo( save );
140 //=============================================================================
144 //=============================================================================
146 istream & operator >>(istream & load, StdMeshers_Deflection1D & hyp)
148 return hyp.LoadFrom( load );
150 //================================================================================
152 * \brief Evaluate curve deflection between two points
153 * \param theCurve - the curve
154 * \param theU1 - the parameter of the first point
155 * \param theU2 - the parameter of the second point
156 * \retval double - deflection value
158 //================================================================================
160 static double deflection(const GeomAdaptor_Curve & theCurve,
164 if ( theCurve.GetType() == GeomAbs_Line )
166 // line between theU1 and theU2
167 gp_Pnt p1 = theCurve.Value( theU1 ), p2 = theCurve.Value( theU2 );
168 gp_Lin segment( p1, gp_Vec( p1, p2 ));
170 // evaluate square distance of theCurve from the segment
171 Standard_Real dist2 = 0;
173 const double step = ( theU2 - theU1 ) / nbPnt;
174 while (( theU1 += step ) < theU2 )
175 dist2 = Max( dist2, segment.SquareDistance( theCurve.Value( theU1 )));
177 return sqrt( dist2 );
180 //================================================================================
182 * \brief Initialize deflection value by the mesh built on the geometry
183 * \param theMesh - the built mesh
184 * \param theShape - the geometry of interest
185 * \retval bool - true if parameter values have been successfully defined
187 //================================================================================
189 bool StdMeshers_Deflection1D::SetParametersByMesh(const SMESH_Mesh* theMesh,
190 const TopoDS_Shape& theShape)
192 if ( !theMesh || theShape.IsNull() )
197 Standard_Real UMin, UMax;
201 TopTools_IndexedMapOfShape edgeMap;
202 TopExp::MapShapes( theShape, TopAbs_EDGE, edgeMap );
204 for ( int iE = 1; iE <= edgeMap.Extent(); ++iE )
206 const TopoDS_Edge& edge = TopoDS::Edge( edgeMap( iE ));
207 Handle(Geom_Curve) C = BRep_Tool::Curve( edge, L, UMin, UMax );
208 GeomAdaptor_Curve AdaptCurve(C);
209 if ( AdaptCurve.GetType() != GeomAbs_Line )
211 vector< double > params;
212 SMESHDS_Mesh* aMeshDS = const_cast< SMESH_Mesh* >( theMesh )->GetMeshDS();
213 if ( SMESH_Algo::GetNodeParamOnEdge( aMeshDS, edge, params ))
216 for ( int i = 1; i < params.size(); ++i )
217 _value = Max( _value, deflection( AdaptCurve, params[ i-1 ], params[ i ]));