1 // Copyright (C) 2007-2023 CEA/DEN, EDF R&D, OPEN CASCADE
3 // This library is free software; you can redistribute it and/or
4 // modify it under the terms of the GNU Lesser General Public
5 // License as published by the Free Software Foundation; either
6 // version 2.1 of the License, or (at your option) any later version.
8 // This library is distributed in the hope that it will be useful,
9 // but WITHOUT ANY WARRANTY; without even the implied warranty of
10 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
11 // Lesser General Public License for more details.
13 // You should have received a copy of the GNU Lesser General Public
14 // License along with this library; if not, write to the Free Software
15 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
17 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
20 // NETGENPlugin : C++ implementation
21 // File : NETGENPlugin_SimpleHypothesis_2D.cxx
22 // Author : Edward AGAPOV
24 //=============================================================================
26 #include "NETGENPlugin_SimpleHypothesis_2D.hxx"
27 #include "NETGENPlugin_Hypothesis.hxx"
29 #include <SMESHDS_SubMesh.hxx>
30 #include <SMESH_ControlsDef.hxx>
31 #include <SMESH_Mesh.hxx>
32 #include <SMESH_subMesh.hxx>
34 #include <TopExp_Explorer.hxx>
36 #include <utilities.h>
40 //=============================================================================
44 //=============================================================================
45 NETGENPlugin_SimpleHypothesis_2D::NETGENPlugin_SimpleHypothesis_2D (int hypId,
47 : SMESH_Hypothesis(hypId, gen),
48 _nbSegments ((int)NETGENPlugin_Hypothesis::GetDefaultNbSegPerEdge()),
53 _name = "NETGEN_SimpleParameters_2D";
57 //=============================================================================
61 //=============================================================================
62 void NETGENPlugin_SimpleHypothesis_2D::SetNumberOfSegments(int nb)
65 throw SALOME_Exception("Number of segments must be positive");
66 if (nb != _nbSegments)
69 if ( _nbSegments ) _segmentLength = 0.;
70 NotifySubMeshesHypothesisModification();
74 //=============================================================================
78 //=============================================================================
79 void NETGENPlugin_SimpleHypothesis_2D::SetLocalLength(double segmentLength)
81 if ( segmentLength < DBL_MIN )
82 throw SALOME_Exception("segment length must be more than zero");
83 if (segmentLength != _segmentLength)
85 _segmentLength = segmentLength;
86 if ( _segmentLength > DBL_MIN ) _nbSegments = 0;
87 NotifySubMeshesHypothesisModification();
91 //=============================================================================
95 //=============================================================================
96 void NETGENPlugin_SimpleHypothesis_2D::LengthFromEdges()
101 NotifySubMeshesHypothesisModification();
105 //=============================================================================
109 //=============================================================================
110 void NETGENPlugin_SimpleHypothesis_2D::SetMaxElementArea(double area)
112 if ( area < DBL_MIN )
117 NotifySubMeshesHypothesisModification();
121 //=======================================================================
122 //function : SetAllowQuadrangles
123 //purpose : Enables/disables generation of quadrangular faces
124 //=======================================================================
126 void NETGENPlugin_SimpleHypothesis_2D::SetAllowQuadrangles(bool toAllow)
128 if ( _allowQuad != toAllow )
130 _allowQuad = toAllow;
131 NotifySubMeshesHypothesisModification();
135 //=======================================================================
136 //function : GetAllowQuadrangles
137 //purpose : Returns true if generation of quadrangular faces is enabled
138 //=======================================================================
140 bool NETGENPlugin_SimpleHypothesis_2D::GetAllowQuadrangles() const
145 //=============================================================================
149 //=============================================================================
150 ostream & NETGENPlugin_SimpleHypothesis_2D::SaveTo(ostream & save)
152 save << _nbSegments << " " << _segmentLength << " " << _area << " " << _allowQuad;
157 //=============================================================================
161 //=============================================================================
162 istream & NETGENPlugin_SimpleHypothesis_2D::LoadFrom(istream & load)
167 isOK = static_cast<bool>(load >> val);
169 _nbSegments = (int) val;
171 load.clear(ios::badbit | load.rdstate());
173 isOK = static_cast<bool>(load >> val);
175 _segmentLength = val;
177 load.clear(ios::badbit | load.rdstate());
179 isOK = static_cast<bool>(load >> val);
183 load.clear(ios::badbit | load.rdstate());
190 //================================================================================
192 * \brief Does nothing
193 * \param theMesh - the built mesh
194 * \param theShape - the geometry of interest
195 * \retval bool - always false
197 //================================================================================
198 bool NETGENPlugin_SimpleHypothesis_2D::SetParametersByMesh(const SMESH_Mesh* theMesh,
199 const TopoDS_Shape& theShape)
201 // Find out nb of segments.
202 smIdType nbSeg = 0, nbEdges = 0;
203 TopExp_Explorer exp( theShape, TopAbs_EDGE );
204 for ( ; exp.More(); exp.Next() ) {
205 SMESH_subMesh* sm = theMesh->GetSubMeshContaining( exp.Current() );
206 if ( sm && !sm->IsEmpty() ) {
207 nbSeg += sm->GetSubMeshDS()->NbElements();
212 _nbSegments = nbSeg / nbEdges;
214 // Find out max face area
216 SMESH::Controls::Area areaControl;
217 SMESH::Controls::TSequenceOfXYZ nodesCoords;
218 const int nbFacesToCheck = 100;
219 for ( exp.Init( theShape, TopAbs_FACE ); exp.More(); exp.Next() ) {
220 SMESH_subMesh* sm = theMesh->GetSubMeshContaining( exp.Current() );
221 if ( sm && !sm->IsEmpty() ) {
222 SMDS_ElemIteratorPtr fIt = sm->GetSubMeshDS()->GetElements();
223 int nbCheckedFaces = 0;
224 while ( fIt->more() && nbCheckedFaces++ < nbFacesToCheck ) {
225 const SMDS_MeshElement* elem = fIt->next();
226 areaControl.GetPoints( elem, nodesCoords );
227 _area = max( _area, areaControl.GetValue( nodesCoords ));
234 //================================================================================
236 * \brief Initialize my parameter values by default parameters.
237 * \retval bool - true if parameter values have been successfully defined
239 //================================================================================
241 bool NETGENPlugin_SimpleHypothesis_2D::SetParametersByDefaults(const TDefaults& dflts,
242 const SMESH_Mesh* /*theMesh*/)
244 _nbSegments = dflts._nbSegments;
245 _segmentLength = dflts._elemLength;
246 return _nbSegments && _segmentLength > 0;