-// NETGENPlugin : C++ implementation
+// Copyright (C) 2007-2020 CEA/DEN, EDF R&D, OPEN CASCADE
+//
+// Copyright (C) 2003-2007 OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
+// CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
+//
+// This library is free software; you can redistribute it and/or
+// modify it under the terms of the GNU Lesser General Public
+// License as published by the Free Software Foundation; either
+// version 2.1 of the License, or (at your option) any later version.
//
-// Copyright (C) 2006 OPEN CASCADE, CEA/DEN, EDF R&D
-//
-// This library is free software; you can redistribute it and/or
-// modify it under the terms of the GNU Lesser General Public
-// License as published by the Free Software Foundation; either
-// version 2.1 of the License.
-//
-// This library is distributed in the hope that it will be useful,
-// but WITHOUT ANY WARRANTY; without even the implied warranty of
-// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
-// Lesser General Public License for more details.
-//
-// You should have received a copy of the GNU Lesser General Public
-// License along with this library; if not, write to the Free Software
-// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
-//
-// See http://www.opencascade.org/SALOME/ or email : webmaster.salome@opencascade.org
+// This library is distributed in the hope that it will be useful,
+// but WITHOUT ANY WARRANTY; without even the implied warranty of
+// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
+// Lesser General Public License for more details.
//
+// You should have received a copy of the GNU Lesser General Public
+// License along with this library; if not, write to the Free Software
+// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
//
+// See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
+//
+
+// NETGENPlugin : C++ implementation
// File : NETGENPlugin_NETGEN_2D3D.cxx
// Author : Michael Sazonov (OCN)
// Date : 20/03/2006
// Project : SALOME
-// $Header$
//=============================================================================
-using namespace std;
-
+//
#include "NETGENPlugin_NETGEN_2D3D.hxx"
#include "NETGENPlugin_Hypothesis.hxx"
+#include "NETGENPlugin_SimpleHypothesis_3D.hxx"
#include "NETGENPlugin_Mesher.hxx"
+#include <SMESHDS_Mesh.hxx>
+#include <SMESH_ControlsDef.hxx>
#include <SMESH_Gen.hxx>
#include <SMESH_Mesh.hxx>
-#include <SMESH_ControlsDef.hxx>
-#include <SMESHDS_Mesh.hxx>
+#include <StdMeshers_ViscousLayers.hxx>
+
#include <utilities.h>
#include <list>
+namespace nglib {
+#include <nglib.h>
+}
+#include <meshing.hpp>
+
+using namespace std;
+
//=============================================================================
/*!
*
*/
//=============================================================================
-NETGENPlugin_NETGEN_2D3D::NETGENPlugin_NETGEN_2D3D(int hypId, int studyId,
+NETGENPlugin_NETGEN_2D3D::NETGENPlugin_NETGEN_2D3D(int hypId,
SMESH_Gen* gen)
- : SMESH_3D_Algo(hypId, studyId, gen)
+ : SMESH_3D_Algo(hypId, gen)
{
- MESSAGE("NETGENPlugin_NETGEN_2D3D::NETGENPlugin_NETGEN_2D3D");
+ //MESSAGE("NETGENPlugin_NETGEN_2D3D::NETGENPlugin_NETGEN_2D3D");
_name = "NETGEN_2D3D";
_shapeType = (1 << TopAbs_SHELL) | (1 << TopAbs_SOLID);// 1 bit /shape type
_compatibleHypothesis.push_back("NETGEN_Parameters");
- _requireDescretBoundary = false;
+ _compatibleHypothesis.push_back("NETGEN_SimpleParameters_3D");
+ _compatibleHypothesis.push_back( StdMeshers_ViscousLayers::GetHypType() );
+ _requireDiscreteBoundary = false;
_onlyUnaryInput = false;
_hypothesis = NULL;
+ _supportSubmeshes = true;
}
//=============================================================================
/*!
- *
+ *
*/
//=============================================================================
NETGENPlugin_NETGEN_2D3D::~NETGENPlugin_NETGEN_2D3D()
{
- MESSAGE("NETGENPlugin_NETGEN_2D3D::~NETGENPlugin_NETGEN_2D3D");
+ //MESSAGE("NETGENPlugin_NETGEN_2D3D::~NETGENPlugin_NETGEN_2D3D");
}
//=============================================================================
/*!
- *
+ *
*/
//=============================================================================
-bool NETGENPlugin_NETGEN_2D3D::CheckHypothesis
- (SMESH_Mesh& aMesh,
- const TopoDS_Shape& aShape,
- SMESH_Hypothesis::Hypothesis_Status& aStatus)
+bool NETGENPlugin_NETGEN_2D3D::CheckHypothesis (SMESH_Mesh& aMesh,
+ const TopoDS_Shape& aShape,
+ Hypothesis_Status& aStatus)
{
- MESSAGE("NETGENPlugin_NETGEN_2D3D::CheckHypothesis");
+ _hypothesis = NULL;
+ _viscousLayersHyp = NULL;
+ _mesher = NULL;
- _hypothesis = NULL;
-
- list<const SMESHDS_Hypothesis*>::const_iterator itl;
- const SMESHDS_Hypothesis* theHyp;
-
- const list<const SMESHDS_Hypothesis*>& hyps = GetUsedHypothesis(aMesh, aShape);
- int nbHyp = hyps.size();
- if (!nbHyp)
+ const list<const SMESHDS_Hypothesis*>& hyps = GetUsedHypothesis(aMesh, aShape, /*noAux=*/false);
+ if ( hyps.empty() )
{
aStatus = SMESH_Hypothesis::HYP_OK;
return true; // can work with no hypothesis
}
- itl = hyps.begin();
- theHyp = (*itl); // use only the first hypothesis
-
- string hypName = theHyp->GetName();
-
- if (hypName == "NETGEN_Parameters")
+ list<const SMESHDS_Hypothesis*>::const_iterator h = hyps.begin();
+ for ( ; h != hyps.end(); ++h )
{
- _hypothesis = static_cast<const NETGENPlugin_Hypothesis*> (theHyp);
- ASSERT(_hypothesis);
- aStatus = SMESH_Hypothesis::HYP_OK;
+ const SMESHDS_Hypothesis* aHyp = *h;
+ std::string hypName = aHyp->GetName();
+
+ if ( std::find( _compatibleHypothesis.begin(), _compatibleHypothesis.end(),
+ hypName ) != _compatibleHypothesis.end() )
+ {
+ if ( hypName == StdMeshers_ViscousLayers::GetHypType() )
+ _viscousLayersHyp = dynamic_cast<const StdMeshers_ViscousLayers*>( aHyp );
+ else
+ _hypothesis = aHyp;
+ aStatus = SMESH_Hypothesis::HYP_OK;
+ }
+ else
+ {
+ aStatus = SMESH_Hypothesis::HYP_INCOMPATIBLE;
+ break;
+ }
}
- else
- aStatus = SMESH_Hypothesis::HYP_INCOMPATIBLE;
return aStatus == SMESH_Hypothesis::HYP_OK;
}
//=============================================================================
/*!
- *Here we are going to use the NETGEN mesher
+ * Here we are going to use the NETGEN mesher
*/
//=============================================================================
bool NETGENPlugin_NETGEN_2D3D::Compute(SMESH_Mesh& aMesh,
const TopoDS_Shape& aShape)
{
- SMESHDS_Mesh* meshDS = aMesh.GetMeshDS();
+ netgen::multithread.terminate = 0;
- NETGENPlugin_Mesher mesher(meshDS, aShape, true);
- mesher.SetParameters(_hypothesis);
+ NETGENPlugin_Mesher mesher(&aMesh, aShape, true);
+ mesher.SetParameters(dynamic_cast<const NETGENPlugin_Hypothesis*>(_hypothesis));
+ mesher.SetParameters(dynamic_cast<const NETGENPlugin_SimpleHypothesis_3D*>(_hypothesis));
+ mesher.SetParameters(_viscousLayersHyp);
+ mesher.SetSelfPointer( &_mesher );
return mesher.Compute();
}
//=============================================================================
/*!
- *
+ *
*/
//=============================================================================
-ostream & NETGENPlugin_NETGEN_2D3D::SaveTo(ostream & save)
+void NETGENPlugin_NETGEN_2D3D::CancelCompute()
{
- return save;
+ SMESH_Algo::CancelCompute();
+ netgen::multithread.terminate = 1;
}
-//=============================================================================
+//================================================================================
/*!
- *
+ * \brief Return progress of Compute() [0.,1]
*/
-//=============================================================================
+//================================================================================
-istream & NETGENPlugin_NETGEN_2D3D::LoadFrom(istream & load)
+double NETGENPlugin_NETGEN_2D3D::GetProgress() const
{
- return load;
-}
+ double & progress = (double &)_progress;
+ if ( _mesher )
+ progress = _mesher->GetProgress(this, &_progressTic, &_progress);
+ else if ( _progress > 0.001 )
+ progress = 0.99;
-//=============================================================================
-/*!
- *
- */
-//=============================================================================
-
-ostream & operator << (ostream & save, NETGENPlugin_NETGEN_2D3D & hyp)
-{
- return hyp.SaveTo( save );
+ return _progress;
}
//=============================================================================
/*!
- *
+ *
*/
//=============================================================================
-istream & operator >> (istream & load, NETGENPlugin_NETGEN_2D3D & hyp)
+bool NETGENPlugin_NETGEN_2D3D::Evaluate(SMESH_Mesh& aMesh,
+ const TopoDS_Shape& aShape,
+ MapShapeNbElems& aResMap)
{
- return hyp.LoadFrom( load );
+ NETGENPlugin_Mesher mesher(&aMesh, aShape, true);
+ mesher.SetParameters(dynamic_cast<const NETGENPlugin_Hypothesis*>(_hypothesis));
+ mesher.SetParameters(dynamic_cast<const NETGENPlugin_SimpleHypothesis_2D*>(_hypothesis));
+ return mesher.Evaluate(aResMap);
}
