]> SALOME platform Git repositories - modules/smesh.git/blob - src/StdMeshers/StdMeshers_MaxElementVolume.cxx
Salome HOME
0020963: EDF 1528 SMESH: How to combine Netgen1D2D with sizemap with Quadrangle(Mappi...
[modules/smesh.git] / src / StdMeshers / StdMeshers_MaxElementVolume.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_MaxElementVolume.cxx
25 //           Moved here from SMESH_MaxElementVolume.cxx
26 //  Author : Paul RASCLE, EDF
27 //  Module : SMESH
28 //
29 #include "StdMeshers_MaxElementVolume.hxx"
30
31 #include "SMDS_MeshElement.hxx"
32 #include "SMESHDS_SubMesh.hxx"
33 #include "SMESH_ControlsDef.hxx"
34 #include "SMESH_Mesh.hxx"
35
36 #include "utilities.h"
37
38 #include <TopExp.hxx>
39 #include <TopExp_Explorer.hxx>
40 #include <TopTools_IndexedMapOfShape.hxx>
41
42 using namespace std;
43
44 //=============================================================================
45 /*!
46  *  
47  */
48 //=============================================================================
49
50 StdMeshers_MaxElementVolume::StdMeshers_MaxElementVolume(int hypId, int studyId, SMESH_Gen* gen)
51   : SMESH_Hypothesis(hypId, studyId, gen)
52 {
53   _maxVolume = 1.;
54   _name = "MaxElementVolume";
55   _param_algo_dim = 3;
56 }
57
58 //=============================================================================
59 /*!
60  *  
61  */
62 //=============================================================================
63
64 StdMeshers_MaxElementVolume::~StdMeshers_MaxElementVolume()
65 {
66   MESSAGE("StdMeshers_MaxElementVolume::~StdMeshers_MaxElementVolume");
67 }
68
69 //=============================================================================
70 /*!
71  *  
72  */
73 //=============================================================================
74
75 void StdMeshers_MaxElementVolume::SetMaxVolume(double maxVolume)
76   throw (SALOME_Exception)
77 {
78   double oldVolume = _maxVolume;
79   if (maxVolume <= 0) 
80     throw SALOME_Exception(LOCALIZED("maxVolume must be positive"));
81   _maxVolume = maxVolume;
82   if (_maxVolume != oldVolume)
83     NotifySubMeshesHypothesisModification();
84 }
85
86 //=============================================================================
87 /*!
88  *  
89  */
90 //=============================================================================
91
92 double StdMeshers_MaxElementVolume::GetMaxVolume() const
93 {
94   return _maxVolume;
95 }
96
97 //=============================================================================
98 /*!
99  *  
100  */
101 //=============================================================================
102
103 ostream & StdMeshers_MaxElementVolume::SaveTo(ostream & save)
104 {
105   save << this->_maxVolume;
106   return save;
107 }
108
109 //=============================================================================
110 /*!
111  *  
112  */
113 //=============================================================================
114
115 istream & StdMeshers_MaxElementVolume::LoadFrom(istream & load)
116 {
117   bool isOK = true;
118   double a;
119   isOK = (load >> a);
120   if (isOK)
121     this->_maxVolume = a;
122   else 
123     load.clear(ios::badbit | load.rdstate());
124   return load;
125 }
126
127 //=============================================================================
128 /*!
129  *  
130  */
131 //=============================================================================
132
133 ostream & operator << (ostream & save, StdMeshers_MaxElementVolume & hyp)
134 {
135   return hyp.SaveTo( save );
136 }
137
138 //=============================================================================
139 /*!
140  *  
141  */
142 //=============================================================================
143
144 istream & operator >> (istream & load, StdMeshers_MaxElementVolume & hyp)
145 {
146   return hyp.LoadFrom( load );
147 }
148
149
150 //================================================================================
151 /*!
152  * \brief Initialize maximal area by the mesh built on the geometry
153  * \param theMesh - the built mesh
154  * \param theShape - the geometry of interest
155  * \retval bool - true if parameter values have been successfully defined
156  */
157 //================================================================================
158
159 bool StdMeshers_MaxElementVolume::SetParametersByMesh(const SMESH_Mesh*   theMesh,
160                                                       const TopoDS_Shape& theShape)
161 {
162   if ( !theMesh || theShape.IsNull() )
163     return false;
164
165   _maxVolume = 0;
166
167   SMESH::Controls::Volume volumeControl;
168
169   TopTools_IndexedMapOfShape volMap;
170   TopExp::MapShapes( theShape, TopAbs_SOLID, volMap );
171   if ( volMap.IsEmpty() )
172     TopExp::MapShapes( theShape, TopAbs_SHELL, volMap );
173   if ( volMap.IsEmpty() )
174     return false;
175
176   SMESHDS_Mesh* aMeshDS = const_cast< SMESH_Mesh* >( theMesh )->GetMeshDS();
177
178   for ( int iV = 1; iV <= volMap.Extent(); ++iV )
179   {
180     const TopoDS_Shape& S = volMap( iV );
181     SMESHDS_SubMesh * subMesh = aMeshDS->MeshElements( S );
182     if ( !subMesh && S.ShapeType() == TopAbs_SOLID ) {
183       TopExp_Explorer shellExp( S, TopAbs_SHELL );
184       if ( shellExp.More() )
185         subMesh = aMeshDS->MeshElements( shellExp.Current() );
186     }
187     if ( !subMesh) 
188       return false;
189     SMDS_ElemIteratorPtr vIt = subMesh->GetElements();
190     while ( vIt->more() )
191     {
192       const SMDS_MeshElement* elem = vIt->next();
193       if ( elem->GetType() == SMDSAbs_Volume ) {
194         _maxVolume = max( _maxVolume, volumeControl.GetValue( elem->GetID() ));
195       }
196     }
197   }
198   return _maxVolume > 0;
199 }
200 //================================================================================
201 /*!
202  * \brief Initialize my parameter values by default parameters.
203  *  \retval bool - true if parameter values have been successfully defined
204  */
205 //================================================================================
206
207 bool StdMeshers_MaxElementVolume::SetParametersByDefaults(const TDefaults&  dflts,
208                                                           const SMESH_Mesh* /*theMesh*/)
209 {
210   return ( _maxVolume = dflts._elemLength*dflts._elemLength*dflts._elemLength );
211 }
212