1 // Copyright (C) 2007-2016 CEA/DEN, EDF R&D, OPEN CASCADE
3 // Copyright (C) 2003-2007 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, or (at your option) any later version.
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
23 // NETGENPlugin : C++ implementation
24 // File : NETGENPlugin_NETGEN_2D.cxx
25 // Author : Michael Sazonov (OCN)
28 //=============================================================================
30 #include "NETGENPlugin_NETGEN_2D.hxx"
31 #include "NETGENPlugin_Hypothesis_2D.hxx"
32 #include "NETGENPlugin_SimpleHypothesis_2D.hxx"
33 #include "NETGENPlugin_Mesher.hxx"
35 #include <SMESHDS_Mesh.hxx>
36 #include <SMESH_ControlsDef.hxx>
37 #include <SMESH_Gen.hxx>
38 #include <SMESH_Mesh.hxx>
39 #include <StdMeshers_ViscousLayers2D.hxx>
40 #include <utilities.h>
47 #include <meshing.hpp>
51 //=============================================================================
55 //=============================================================================
57 NETGENPlugin_NETGEN_2D::NETGENPlugin_NETGEN_2D(int hypId, int studyId,
59 : SMESH_2D_Algo(hypId, studyId, gen)
62 _shapeType = (1 << TopAbs_FACE); // 1 bit /shape type
63 _compatibleHypothesis.push_back("NETGEN_Parameters_2D");
64 _compatibleHypothesis.push_back("NETGEN_SimpleParameters_2D");
65 _compatibleHypothesis.push_back( StdMeshers_ViscousLayers2D::GetHypType() );
66 _requireDiscreteBoundary = false;
67 _onlyUnaryInput = false;
69 _supportSubmeshes = true;
72 //=============================================================================
76 //=============================================================================
78 NETGENPlugin_NETGEN_2D::~NETGENPlugin_NETGEN_2D()
82 //=============================================================================
86 //=============================================================================
88 bool NETGENPlugin_NETGEN_2D::CheckHypothesis (SMESH_Mesh& aMesh,
89 const TopoDS_Shape& aShape,
90 Hypothesis_Status& aStatus)
93 _isViscousLayers2D = false;
96 // can work with no hypothesis
97 aStatus = SMESH_Hypothesis::HYP_OK;
99 const list<const SMESHDS_Hypothesis*>& hyps = GetUsedHypothesis(aMesh, aShape, /*skipAux=*/false);
100 list<const SMESHDS_Hypothesis*>::const_iterator h = hyps.begin();
101 for ( ; h != hyps.end(); ++h )
103 const SMESHDS_Hypothesis* theHyp = *h;
104 string hypName = theHyp->GetName();
105 if ( hypName == StdMeshers_ViscousLayers2D::GetHypType() )
106 _isViscousLayers2D = true;
107 else if ( _hypothesis )
108 aStatus = SMESH_Hypothesis::HYP_INCOMPATIBLE;
110 _hypothesis = theHyp;
113 if ( aStatus == HYP_OK && _isViscousLayers2D )
114 error( StdMeshers_ViscousLayers2D::CheckHypothesis( aMesh, aShape, aStatus ));
116 return aStatus == SMESH_Hypothesis::HYP_OK;
119 //=============================================================================
123 //=============================================================================
125 bool NETGENPlugin_NETGEN_2D::Compute(SMESH_Mesh& aMesh,
126 const TopoDS_Shape& aShape)
128 netgen::multithread.terminate = 0;
130 NETGENPlugin_Mesher mesher(&aMesh, aShape, /*is3D = */false);
131 mesher.SetParameters(dynamic_cast<const NETGENPlugin_Hypothesis*>(_hypothesis));
132 mesher.SetParameters(dynamic_cast<const NETGENPlugin_SimpleHypothesis_2D*>(_hypothesis));
133 mesher.SetViscousLayers2DAssigned( _isViscousLayers2D );
134 mesher.SetSelfPointer( &_mesher );
135 return mesher.Compute();
138 //=============================================================================
140 * Terminate Compute()
142 //=============================================================================
144 void NETGENPlugin_NETGEN_2D::CancelCompute()
146 SMESH_Algo::CancelCompute();
147 netgen::multithread.terminate = 1;
150 //================================================================================
152 * \brief Return progress of Compute() [0.,1]
154 //================================================================================
156 double NETGENPlugin_NETGEN_2D::GetProgress() const
158 return _mesher ? _mesher->GetProgress(this, &_progressTic, &_progress) : 0;
161 //=============================================================================
165 //=============================================================================
167 bool NETGENPlugin_NETGEN_2D::Evaluate(SMESH_Mesh& aMesh,
168 const TopoDS_Shape& aShape,
169 MapShapeNbElems& aResMap)
172 NETGENPlugin_Mesher mesher(&aMesh, aShape, false);
173 mesher.SetParameters(dynamic_cast<const NETGENPlugin_Hypothesis*>(_hypothesis));
174 mesher.SetParameters(dynamic_cast<const NETGENPlugin_SimpleHypothesis_2D*>(_hypothesis));
175 return mesher.Evaluate(aResMap);