194 lines
6.2 KiB
C++
194 lines
6.2 KiB
C++
// Copyright (C) 2007-2019 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.
|
|
//
|
|
// 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
|
|
//=============================================================================
|
|
//
|
|
#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 <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,
|
|
SMESH_Gen* gen)
|
|
: SMESH_3D_Algo(hypId, gen)
|
|
{
|
|
//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");
|
|
_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");
|
|
}
|
|
|
|
//=============================================================================
|
|
/*!
|
|
*
|
|
*/
|
|
//=============================================================================
|
|
|
|
bool NETGENPlugin_NETGEN_2D3D::CheckHypothesis (SMESH_Mesh& aMesh,
|
|
const TopoDS_Shape& aShape,
|
|
Hypothesis_Status& aStatus)
|
|
{
|
|
_hypothesis = NULL;
|
|
_viscousLayersHyp = NULL;
|
|
_mesher = NULL;
|
|
|
|
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
|
|
}
|
|
|
|
list<const SMESHDS_Hypothesis*>::const_iterator h = hyps.begin();
|
|
for ( ; h != hyps.end(); ++h )
|
|
{
|
|
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;
|
|
}
|
|
}
|
|
|
|
return aStatus == SMESH_Hypothesis::HYP_OK;
|
|
}
|
|
|
|
//=============================================================================
|
|
/*!
|
|
* Here we are going to use the NETGEN mesher
|
|
*/
|
|
//=============================================================================
|
|
|
|
bool NETGENPlugin_NETGEN_2D3D::Compute(SMESH_Mesh& aMesh,
|
|
const TopoDS_Shape& aShape)
|
|
{
|
|
netgen::multithread.terminate = 0;
|
|
|
|
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();
|
|
}
|
|
|
|
//=============================================================================
|
|
/*!
|
|
*
|
|
*/
|
|
//=============================================================================
|
|
|
|
void NETGENPlugin_NETGEN_2D3D::CancelCompute()
|
|
{
|
|
SMESH_Algo::CancelCompute();
|
|
netgen::multithread.terminate = 1;
|
|
}
|
|
|
|
//================================================================================
|
|
/*!
|
|
* \brief Return progress of Compute() [0.,1]
|
|
*/
|
|
//================================================================================
|
|
|
|
double NETGENPlugin_NETGEN_2D3D::GetProgress() const
|
|
{
|
|
double & progress = (double &)_progress;
|
|
if ( _mesher )
|
|
progress = _mesher->GetProgress(this, &_progressTic, &_progress);
|
|
else if ( _progress > 0.001 )
|
|
progress = 0.99;
|
|
|
|
return _progress;
|
|
}
|
|
|
|
//=============================================================================
|
|
/*!
|
|
*
|
|
*/
|
|
//=============================================================================
|
|
|
|
bool NETGENPlugin_NETGEN_2D3D::Evaluate(SMESH_Mesh& aMesh,
|
|
const TopoDS_Shape& aShape,
|
|
MapShapeNbElems& aResMap)
|
|
{
|
|
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);
|
|
}
|