]> SALOME platform Git repositories - modules/smesh.git/blob - src/StdMeshers/StdMeshers_MaxElementArea.cxx
Salome HOME
Merge from branch BR_Dev_For_4_0 (from tag mergeto_BR_QT4_Dev_12Feb08)
[modules/smesh.git] / src / StdMeshers / StdMeshers_MaxElementArea.cxx
1 //  SMESH SMESH : implementaion of SMESH idl descriptions
2 //
3 //  Copyright (C) 2003  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 //
24 //  File   : StdMeshers_MaxElementArea.cxx
25 //           Moved here from SMESH_MaxElementArea.cxx
26 //  Author : Paul RASCLE, EDF
27 //  Module : SMESH
28 //  $Header$
29
30 #include "StdMeshers_MaxElementArea.hxx"
31
32 #include "SMESH_ControlsDef.hxx"
33 #include "SMDS_MeshElement.hxx"
34 #include "SMESHDS_SubMesh.hxx"
35 #include "SMESH_Mesh.hxx"
36
37 #include <TopExp.hxx>
38 #include <TopTools_IndexedMapOfShape.hxx>
39
40 #include "utilities.h"
41
42 using namespace std;
43
44 //=============================================================================
45 /*!
46  *  
47  */
48 //=============================================================================
49
50 StdMeshers_MaxElementArea::StdMeshers_MaxElementArea(int hypId, int studyId, SMESH_Gen* gen)
51   : SMESH_Hypothesis(hypId, studyId, gen)
52 {
53   _maxArea =1.;
54   _name = "MaxElementArea";
55   _param_algo_dim = 2; 
56 }
57
58 //=============================================================================
59 /*!
60  *  
61  */
62 //=============================================================================
63
64 StdMeshers_MaxElementArea::~StdMeshers_MaxElementArea()
65 {
66 }
67
68 //=============================================================================
69 /*!
70  *  
71  */
72 //=============================================================================
73
74 void StdMeshers_MaxElementArea::SetMaxArea(double maxArea)
75   throw (SALOME_Exception)
76 {
77   double oldArea = _maxArea;
78   if (maxArea <= 0) 
79     throw SALOME_Exception(LOCALIZED("maxArea must be positive"));
80   _maxArea = maxArea;
81   if (_maxArea != oldArea)
82     NotifySubMeshesHypothesisModification();
83 }
84
85 //=============================================================================
86 /*!
87  *  
88  */
89 //=============================================================================
90
91 double StdMeshers_MaxElementArea::GetMaxArea() const
92 {
93   return _maxArea;
94 }
95
96 //=============================================================================
97 /*!
98  *  
99  */
100 //=============================================================================
101
102 ostream & StdMeshers_MaxElementArea::SaveTo(ostream & save)
103 {
104   save << this->_maxArea;
105   return save;
106 }
107
108 //=============================================================================
109 /*!
110  *  
111  */
112 //=============================================================================
113
114 istream & StdMeshers_MaxElementArea::LoadFrom(istream & load)
115 {
116   bool isOK = true;
117   double a;
118   isOK = (load >> a);
119   if (isOK) 
120     this->_maxArea = a;
121   else 
122     load.clear(ios::badbit | load.rdstate());
123   return load;
124 }
125
126 //=============================================================================
127 /*!
128  *  
129  */
130 //=============================================================================
131
132 ostream & operator << (ostream & save, StdMeshers_MaxElementArea & hyp)
133 {
134   return hyp.SaveTo( save );
135 }
136
137 //=============================================================================
138 /*!
139  *  
140  */
141 //=============================================================================
142
143 istream & operator >> (istream & load, StdMeshers_MaxElementArea & hyp)
144 {
145   return hyp.LoadFrom( load );
146 }
147
148 //================================================================================
149 /*!
150  * \brief Initialize maximal area by the mesh built on the geometry
151  * \param theMesh - the built mesh
152  * \param theShape - the geometry of interest
153  * \retval bool - true if parameter values have been successfully defined
154  */
155 //================================================================================
156
157 bool StdMeshers_MaxElementArea::SetParametersByMesh(const SMESH_Mesh*   theMesh,
158                                                     const TopoDS_Shape& theShape)
159 {
160   if ( !theMesh || theShape.IsNull() )
161     return false;
162
163   _maxArea = 0;
164
165   SMESH::Controls::Area areaControl;
166   SMESH::Controls::TSequenceOfXYZ nodesCoords;
167
168   SMESHDS_Mesh* aMeshDS = const_cast< SMESH_Mesh* >( theMesh )->GetMeshDS();
169
170   TopTools_IndexedMapOfShape faceMap;
171   TopExp::MapShapes( theShape, TopAbs_FACE, faceMap );
172   for ( int iF = 1; iF <= faceMap.Extent(); ++iF )
173   {
174     SMESHDS_SubMesh * subMesh = aMeshDS->MeshElements( faceMap( iF ));
175     if ( !subMesh )
176       return false;
177     SMDS_ElemIteratorPtr fIt = subMesh->GetElements();
178     while ( fIt->more() )
179     {
180       const SMDS_MeshElement* elem = fIt->next();
181       if ( elem->GetType() == SMDSAbs_Face ) {
182         areaControl.GetPoints( elem, nodesCoords );
183         _maxArea = max( _maxArea, areaControl.GetValue( nodesCoords ));
184       }
185     }
186   }
187   return _maxArea > 0;
188 }