ForkMaster/netgenplugin
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Merge from V6_4_BR 05/12/2011

Browse Source
This commit is contained in:
vsr 2011-12-05 11:51:56 +00:00
parent 6923ab131d
commit d736f1af2c
6 changed files with 120 additions and 66 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

9
clean_configure
View File

@ -24,10 +24,5 @@ find . -name "*.pyc" -print -exec rm {} \;
#exit #exit
# ==================== ON SORT AVANT # ==================== ON SORT AVANT
find bin -name Makefile.in | xargs rm -f find . -name Makefile.in | xargs rm -f
find doc -name Makefile.in | xargs rm -f ( cd adm_local/unix/config_files && rm -f config.* depcomp install-sh ltmain.sh missing py-compile )
find idl -name Makefile.in | xargs rm -f
find resources -name Makefile.in | xargs rm -f
find salome_adm -name Makefile.in | xargs rm -f
find src -name Makefile.in | xargs rm -f
rm -f Makefile.in

2
configure.ac
View File

@ -29,7 +29,7 @@
# Reorganization for usage of autotools # Reorganization for usage of autotools
# Created from configure.in.base # Created from configure.in.base
# #
AC_INIT([Salome2 Project NETGENPLUGIN module], [6.3.1], [webmaster.salome@opencascade.com], [SalomeNETGENPLUGIN]) AC_INIT([Salome2 Project NETGENPLUGIN module], [6.4.0], [webmaster.salome@opencascade.com], [SalomeNETGENPLUGIN])
AC_CONFIG_AUX_DIR(adm_local/unix/config_files) AC_CONFIG_AUX_DIR(adm_local/unix/config_files)
AC_CANONICAL_HOST AC_CANONICAL_HOST
AC_CANONICAL_TARGET AC_CANONICAL_TARGET

63
src/NETGEN/netgen49ForSalome.patch
View File

@ -333,11 +333,24 @@ diff -Naur netgen-4.9.13_orig/libsrc/occ/occgenmesh.cpp netgen-4.9.13_new/libsrc
if (geom.vmap.FindIndex(TopExp::FirstVertex (edge)) == if (geom.vmap.FindIndex(TopExp::FirstVertex (edge)) ==
geom.vmap.FindIndex(TopExp::LastVertex (edge))) geom.vmap.FindIndex(TopExp::LastVertex (edge)))
@@ -481,13 +499,45 @@ @@ -479,15 +497,64 @@
}
else
{ {
Point<3> fp = occ2ng (BRep_Tool::Pnt (TopExp::FirstVertex (edge))); - Point<3> fp = occ2ng (BRep_Tool::Pnt (TopExp::FirstVertex (edge)));
Point<3> lp = occ2ng (BRep_Tool::Pnt (TopExp::LastVertex (edge))); - Point<3> lp = occ2ng (BRep_Tool::Pnt (TopExp::LastVertex (edge)));
+ TopoDS_Iterator vIt( edge, false );
+ TopoDS_Vertex v1 = TopoDS::Vertex( vIt.Value() ); vIt.Next();
+ TopoDS_Vertex v2 = TopoDS::Vertex( vIt.Value() );
+ if ( v1.Orientation() == TopAbs_REVERSED )
+ std::swap( v1, v2 );
+ const bool isClosedEdge = v1.IsSame( v2 );
+
+ Point<3> fp = occ2ng (BRep_Tool::Pnt (v1));
+ Point<3> lp = occ2ng (BRep_Tool::Pnt (v2));
+ double tol2 = std::min( eps*eps, 1e-6 * Dist2( fp, lp )); + double tol2 = std::min( eps*eps, 1e-6 * Dist2( fp, lp ));
+ if ( isClosedEdge )
+ tol2 = BRep_Tool::Tolerance( v1 ) * BRep_Tool::Tolerance( v1 );
pnums[0] = -1; pnums[0] = -1;
pnums.Last() = -1; pnums.Last() = -1;
@ -348,40 +361,48 @@ diff -Naur netgen-4.9.13_orig/libsrc/occ/occgenmesh.cpp netgen-4.9.13_new/libsrc
+ if (Dist2 (mesh[pi], fp) < tol2) pnums[0] = pi; + if (Dist2 (mesh[pi], fp) < tol2) pnums[0] = pi;
+ if (Dist2 (mesh[pi], lp) < tol2) pnums.Last() = pi; + if (Dist2 (mesh[pi], lp) < tol2) pnums.Last() = pi;
+ } + }
+ if ( pnums[0] == pnums.Last() ) + if (( isClosedEdge && pnums[0] != pnums.Last() ) ||
+ pnums[0] = -1; + ( !isClosedEdge && pnums[0] == pnums.Last() ))
+ pnums[0] = pnums.Last() = -1;
+ if ( pnums[0] == -1 || pnums.Last() == -1 ) + if ( pnums[0] == -1 || pnums.Last() == -1 )
+ { + {
+ // take into account a possible large gap between a vertex and an edge curve + // take into account a possible large gap between a vertex and an edge curve
+ // and a large vertex tolerance covering the whole edge + // end and a large vertex tolerance covering the whole edge
+ if ( pnums[0] == -1 ) + if ( pnums[0] == -1 )
+ { + {
+ double tol = BRep_Tool::Tolerance( TopExp::FirstVertex (edge)); + double tol = BRep_Tool::Tolerance( v1 );
+ for (PointIndex pi = 1; pi < first_ep; pi++) + for (PointIndex pi = 1; pi < first_ep; pi++)
+ if (pi != pnums.Last() && Dist2 (mesh[pi], fp) < 2*tol*tol) + if (pi != pnums.Last() && Dist2 (mesh[pi], fp) < 2*tol*tol)
+ pnums[0] = pi; + pnums[0] = pi;
+ +
+ if ( pnums[0] == -1 ) + if ( pnums[0] == -1 )
+ pnums[0] = geom.vmap.FindIndex (TopExp::FirstVertex (edge)); + pnums[0] = first_ep-1- nvertices + geom.vmap.FindIndex ( v1 );
+ } + }
+ if ( pnums.Last() == -1 ) + if ( isClosedEdge )
+ { + {
+ double tol = BRep_Tool::Tolerance( TopExp::LastVertex (edge)); + pnums.Last() = pnums[0];
+ for (PointIndex pi = 1; pi < first_ep; pi++)
+ if (pi != pnums[0] && Dist2 (mesh[pi], lp) < 2*tol*tol)
+ pnums.Last() = pi;
+
+ if ( pnums.Last() == -1 )
+ pnums.Last() = geom.vmap.FindIndex (TopExp::LastVertex (edge));
+ } + }
+ else
+ {
+ if ( pnums.Last() == -1 )
+ {
+ double tol = BRep_Tool::Tolerance( v2 );
+ for (PointIndex pi = 1; pi < first_ep; pi++)
+ if (pi != pnums[0] && Dist2 (mesh[pi], lp) < 2*tol*tol)
+ pnums.Last() = pi;
+ +
+ if ( Dist2( fp, mesh[PointIndex(pnums[0])]) > + if ( pnums.Last() == -1 )
+ Dist2( lp, mesh[PointIndex(pnums.Last())])) + pnums.Last() = first_ep-1-nvertices + geom.vmap.FindIndex ( v2 );
+ }
+
+ if ( Dist2( fp, mesh[PointIndex(pnums[0])]) >
+ Dist2( lp, mesh[PointIndex(pnums.Last())]))
+ std::swap( pnums[0], pnums.Last() ); + std::swap( pnums[0], pnums.Last() );
+ }
} }
} }
@@ -633,7 +683,8 @@ @@ -633,7 +700,8 @@
} }
(*testout) << "mesh face " << k << endl; (*testout) << "mesh face " << k << endl;
@ -391,7 +412,7 @@ diff -Naur netgen-4.9.13_orig/libsrc/occ/occgenmesh.cpp netgen-4.9.13_new/libsrc
geom.facemeshstatus[k-1] = -1; geom.facemeshstatus[k-1] = -1;
@@ -901,7 +952,8 @@ @@ -901,7 +969,8 @@
// if (k != 36) continue; // if (k != 36) continue;
// (*testout) << "optimize face " << k << endl; // (*testout) << "optimize face " << k << endl;
@ -401,7 +422,7 @@ diff -Naur netgen-4.9.13_orig/libsrc/occ/occgenmesh.cpp netgen-4.9.13_new/libsrc
FaceDescriptor & fd = mesh.GetFaceDescriptor(k); FaceDescriptor & fd = mesh.GetFaceDescriptor(k);
@@ -1456,3 +1508,4 @@ @@ -1456,3 +1525,4 @@
} }
#endif #endif

98
src/NETGENPlugin/NETGENPlugin_Mesher.cxx
View File

@ -119,6 +119,7 @@ NETGENPlugin_Mesher::NETGENPlugin_Mesher (SMESH_Mesh* mesh,
_shape (aShape), _shape (aShape),
_isVolume(isVolume), _isVolume(isVolume),
_optimize(true), _optimize(true),
_fineness(NETGENPlugin_Hypothesis::GetDefaultFineness()),
_simpleHyp(NULL) _simpleHyp(NULL)
{ {
defaultParameters(); defaultParameters();
@ -153,6 +154,7 @@ void NETGENPlugin_Mesher::defaultParameters()
mparams.quad = 0; mparams.quad = 0;
else else
mparams.quad = NETGENPlugin_Hypothesis_2D::GetDefaultQuadAllowed() ? 1 : 0; mparams.quad = NETGENPlugin_Hypothesis_2D::GetDefaultQuadAllowed() ? 1 : 0;
_fineness = NETGENPlugin_Hypothesis::GetDefaultFineness();
} }
//============================================================================= //=============================================================================
@ -196,6 +198,8 @@ void NETGENPlugin_Mesher::SetParameters(const NETGENPlugin_Hypothesis* hyp)
// Initialize global NETGEN parameters: // Initialize global NETGEN parameters:
// maximal mesh segment size // maximal mesh segment size
mparams.maxh = hyp->GetMaxSize(); mparams.maxh = hyp->GetMaxSize();
// maximal mesh element linear size
mparams.minh = hyp->GetMinSize();
// minimal number of segments per edge // minimal number of segments per edge
mparams.segmentsperedge = hyp->GetNbSegPerEdge(); mparams.segmentsperedge = hyp->GetNbSegPerEdge();
// rate of growth of size between elements // rate of growth of size between elements
@ -210,6 +214,7 @@ void NETGENPlugin_Mesher::SetParameters(const NETGENPlugin_Hypothesis* hyp)
mparams.quad = static_cast<const NETGENPlugin_Hypothesis_2D*> mparams.quad = static_cast<const NETGENPlugin_Hypothesis_2D*>
(hyp)->GetQuadAllowed() ? 1 : 0; (hyp)->GetQuadAllowed() ? 1 : 0;
_optimize = hyp->GetOptimize(); _optimize = hyp->GetOptimize();
_fineness = hyp->GetFineness();
_simpleHyp = NULL; _simpleHyp = NULL;
SMESH_Gen_i* smeshGen_i = SMESH_Gen_i::GetSMESHGen(); SMESH_Gen_i* smeshGen_i = SMESH_Gen_i::GetSMESHGen();
@ -413,30 +418,33 @@ namespace
void updateTriangulation( const TopoDS_Shape& shape ) void updateTriangulation( const TopoDS_Shape& shape )
{ {
static set< Poly_Triangulation* > updated; // static set< Poly_Triangulation* > updated;
TopLoc_Location loc; // TopLoc_Location loc;
TopExp_Explorer fExp( shape, TopAbs_FACE ); // TopExp_Explorer fExp( shape, TopAbs_FACE );
for ( ; fExp.More(); fExp.Next() ) // for ( ; fExp.More(); fExp.Next() )
{ // {
Handle(Poly_Triangulation) triangulation = // Handle(Poly_Triangulation) triangulation =
BRep_Tool::Triangulation ( TopoDS::Face( fExp.Current() ), loc); // BRep_Tool::Triangulation ( TopoDS::Face( fExp.Current() ), loc);
if ( triangulation.IsNull() || // if ( triangulation.IsNull() ||
updated.insert( triangulation.operator->() ).second ) // updated.insert( triangulation.operator->() ).second )
{ // {
BRepTools::Clean (shape); // BRepTools::Clean (shape);
try { try {
#if (OCC_VERSION_MAJOR << 16 | OCC_VERSION_MINOR << 8 | OCC_VERSION_MAINTENANCE) > 0x060100 #if (OCC_VERSION_MAJOR << 16 | OCC_VERSION_MINOR << 8 | OCC_VERSION_MAINTENANCE) > 0x060100
OCC_CATCH_SIGNALS; OCC_CATCH_SIGNALS;
#endif #endif
BRepMesh_IncrementalMesh e(shape, 0.01, true); BRepMesh_IncrementalMesh e(shape, 0.01, true);
} }
catch (Standard_Failure) catch (Standard_Failure)
{ {
updated.erase( triangulation.operator->() );
} }
} // updated.erase( triangulation.operator->() );
} // triangulation = BRep_Tool::Triangulation ( TopoDS::Face( fExp.Current() ), loc);
// updated.insert( triangulation.operator->() );
// }
// }
} }
} }
@ -549,6 +557,7 @@ double NETGENPlugin_Mesher::GetDefaultMinSize(const TopoDS_Shape& geom,
Handle(Poly_Triangulation) triangulation = Handle(Poly_Triangulation) triangulation =
BRep_Tool::Triangulation ( TopoDS::Face( fExp.Current() ), loc); BRep_Tool::Triangulation ( TopoDS::Face( fExp.Current() ), loc);
if ( triangulation.IsNull() ) continue; if ( triangulation.IsNull() ) continue;
const double fTol = BRep_Tool::Tolerance( TopoDS::Face( fExp.Current() ));
const TColgp_Array1OfPnt& points = triangulation->Nodes(); const TColgp_Array1OfPnt& points = triangulation->Nodes();
const Poly_Array1OfTriangle& trias = triangulation->Triangles(); const Poly_Array1OfTriangle& trias = triangulation->Triangles();
for ( int iT = trias.Lower(); iT <= trias.Upper(); ++iT ) for ( int iT = trias.Lower(); iT <= trias.Upper(); ++iT )
@ -557,7 +566,7 @@ double NETGENPlugin_Mesher::GetDefaultMinSize(const TopoDS_Shape& geom,
for ( int j = 0; j < 3; ++j ) for ( int j = 0; j < 3; ++j )
{ {
double dist2 = points(*pi[j]).SquareDistance( points( *pi[j+1] )); double dist2 = points(*pi[j]).SquareDistance( points( *pi[j+1] ));
if ( dist2 < minh ) if ( dist2 < minh && fTol*fTol < dist2 )
minh = dist2; minh = dist2;
bb.Add( points(*pi[j])); bb.Add( points(*pi[j]));
} }
@ -579,6 +588,23 @@ double NETGENPlugin_Mesher::GetDefaultMinSize(const TopoDS_Shape& geom,
return minh; return minh;
} }
//================================================================================
/*!
* \brief Restrict size of elements at a given point
*/
//================================================================================
void NETGENPlugin_Mesher::RestrictLocalSize(netgen::Mesh& ngMesh, const gp_XYZ& p, const double size)
{
if ( netgen::mparam.minh > size )
{
ngMesh.SetMinimalH( size );
netgen::mparam.minh = size;
}
netgen::Point3d pi(p.X(), p.Y(), p.Z());
ngMesh.RestrictLocalH( pi, size );
}
//================================================================================ //================================================================================
/*! /*!
* \brief fill ngMesh with nodes and elements of computed submeshes * \brief fill ngMesh with nodes and elements of computed submeshes
@ -591,9 +617,8 @@ bool NETGENPlugin_Mesher::fillNgMesh(const netgen::OCCGeometry& occgeom,
const list< SMESH_subMesh* > & meshedSM) const list< SMESH_subMesh* > & meshedSM)
{ {
TNode2IdMap nodeNgIdMap; TNode2IdMap nodeNgIdMap;
if ( !nodeVec.empty() ) for ( int i = 1; i < nodeVec.size(); ++i )
for ( int i = 1; i < nodeVec.size(); ++i ) nodeNgIdMap.insert( make_pair( nodeVec[i], i ));
nodeNgIdMap.insert( make_pair( nodeVec[i], i ));
TopTools_MapOfShape visitedShapes; TopTools_MapOfShape visitedShapes;
map< SMESH_subMesh*, set< int > > visitedEdgeSM2Faces; map< SMESH_subMesh*, set< int > > visitedEdgeSM2Faces;
@ -712,9 +737,8 @@ bool NETGENPlugin_Mesher::fillNgMesh(const netgen::OCCGeometry& occgeom,
seg.edgenr = ngMesh.GetNSeg() + 1; // segment id seg.edgenr = ngMesh.GetNSeg() + 1; // segment id
ngMesh.AddSegment (seg); ngMesh.AddSegment (seg);
netgen::Point3d ngP1(p1.node->X(), p1.node->Y(), p1.node->Z()); SMESH_TNodeXYZ np1( p1.node ), np2( p2.node );
netgen::Point3d ngP2(p2.node->X(), p2.node->Y(), p2.node->Z()); RestrictLocalSize( ngMesh, 0.5*(np1+np2), (np1-np2).Modulus() );
ngMesh.RestrictLocalH( netgen::Center( ngP1,ngP2), Dist(ngP1,ngP2));
#ifdef DUMP_SEGMENTS #ifdef DUMP_SEGMENTS
cout << "Segment: " << seg.edgenr << " on SMESH face " << helper.GetMeshDS()->ShapeToIndex( face ) << endl cout << "Segment: " << seg.edgenr << " on SMESH face " << helper.GetMeshDS()->ShapeToIndex( face ) << endl
@ -860,9 +884,23 @@ bool NETGENPlugin_Mesher::fillNgMesh(const netgen::OCCGeometry& occgeom,
case TopAbs_VERTEX: { // VERTEX case TopAbs_VERTEX: { // VERTEX
// -------------------------- // --------------------------
SMDS_NodeIteratorPtr nodeIt = smDS->GetNodes(); // issue 0021405. Add node only if a VERTEX is shared by a not meshed EDGE,
if ( nodeIt->more() ) // else netgen removes a free node and nodeVector becomes invalid
ngNodeId( nodeIt->next(), ngMesh, nodeNgIdMap ); PShapeIteratorPtr ansIt = helper.GetAncestors( sm->GetSubShape(),
*sm->GetFather(),
TopAbs_EDGE );
bool toAdd = false;
while ( const TopoDS_Shape* e = ansIt->next() )
{
SMESH_subMesh* eSub = helper.GetMesh()->GetSubMesh( *e );
if (( toAdd = eSub->IsEmpty() )) break;
}
if ( toAdd )
{
SMDS_NodeIteratorPtr nodeIt = smDS->GetNodes();
if ( nodeIt->more() )
ngNodeId( nodeIt->next(), ngMesh, nodeNgIdMap );
}
break; break;
} }
default:; default:;
@ -1662,8 +1700,7 @@ namespace
TopoDS_Iterator vIt( edge ); TopoDS_Iterator vIt( edge );
if ( !vIt.More() ) return; if ( !vIt.More() ) return;
gp_Pnt p = BRep_Tool::Pnt( TopoDS::Vertex( vIt.Value() )); gp_Pnt p = BRep_Tool::Pnt( TopoDS::Vertex( vIt.Value() ));
netgen::Point3d pi(p.X(), p.Y(), p.Z()); NETGENPlugin_Mesher::RestrictLocalSize( mesh, p.XYZ(), size );
mesh.RestrictLocalH(pi, size);
} }
else else
{ {
@ -1672,12 +1709,12 @@ namespace
{ {
Standard_Real u = u1 + delta*i; Standard_Real u = u1 + delta*i;
gp_Pnt p = curve->Value(u); gp_Pnt p = curve->Value(u);
NETGENPlugin_Mesher::RestrictLocalSize( mesh, p.XYZ(), size );
netgen::Point3d pi(p.X(), p.Y(), p.Z()); netgen::Point3d pi(p.X(), p.Y(), p.Z());
mesh.RestrictLocalH(pi, size);
double resultSize = mesh.GetH(pi); double resultSize = mesh.GetH(pi);
if ( resultSize - size > 0.1*size ) if ( resultSize - size > 0.1*size )
// netgen does restriction iff oldH/newH > 1.2 (localh.cpp:136) // netgen does restriction iff oldH/newH > 1.2 (localh.cpp:136)
mesh.RestrictLocalH(pi, resultSize/1.201); NETGENPlugin_Mesher::RestrictLocalSize( mesh, p.XYZ(), resultSize/1.201 );
} }
} }
} }
@ -1808,7 +1845,7 @@ bool NETGENPlugin_Mesher::Compute()
if ( mparams.maxh == 0.0 ) if ( mparams.maxh == 0.0 )
mparams.maxh = occgeo.boundingbox.Diam(); mparams.maxh = occgeo.boundingbox.Diam();
if ( _simpleHyp || mparams.minh == 0.0 ) if ( _simpleHyp || ( mparams.minh == 0.0 && _fineness != NETGENPlugin_Hypothesis::UserDefined))
mparams.minh = GetDefaultMinSize( _shape, mparams.maxh ); mparams.minh = GetDefaultMinSize( _shape, mparams.maxh );
#ifdef NETGEN_NEW #ifdef NETGEN_NEW
occgeo.face_maxh = mparams.maxh; occgeo.face_maxh = mparams.maxh;
@ -1870,8 +1907,7 @@ bool NETGENPlugin_Mesher::Compute()
const TopoDS_Shape& shape = ShapesWithLocalSize.FindKey(key); const TopoDS_Shape& shape = ShapesWithLocalSize.FindKey(key);
const TopoDS_Vertex& v = TopoDS::Vertex(shape); const TopoDS_Vertex& v = TopoDS::Vertex(shape);
gp_Pnt p = BRep_Tool::Pnt(v); gp_Pnt p = BRep_Tool::Pnt(v);
netgen::Point3d pi(p.X(), p.Y(), p.Z()); NETGENPlugin_Mesher::RestrictLocalSize( *ngMesh, p.XYZ(), hi );
ngMesh->RestrictLocalH(pi, hi);
} }
} }

3
src/NETGENPlugin/NETGENPlugin_Mesher.hxx
View File

@ -97,6 +97,8 @@ class NETGENPLUGIN_EXPORT NETGENPlugin_Mesher
static double GetDefaultMinSize(const TopoDS_Shape& shape, static double GetDefaultMinSize(const TopoDS_Shape& shape,
const double maxSize); const double maxSize);
static void RestrictLocalSize(netgen::Mesh& ngMesh, const gp_XYZ& p, const double size);
static int FillSMesh(const netgen::OCCGeometry& occgeom, static int FillSMesh(const netgen::OCCGeometry& occgeom,
const netgen::Mesh& ngMesh, const netgen::Mesh& ngMesh,
const NETGENPlugin_ngMeshInfo& initState, const NETGENPlugin_ngMeshInfo& initState,
@ -134,6 +136,7 @@ class NETGENPLUGIN_EXPORT NETGENPlugin_Mesher
const TopoDS_Shape& _shape; const TopoDS_Shape& _shape;
bool _isVolume; bool _isVolume;
bool _optimize; bool _optimize;
int _fineness;
const NETGENPlugin_SimpleHypothesis_2D * _simpleHyp; const NETGENPlugin_SimpleHypothesis_2D * _simpleHyp;
std::map< int, std::pair<int,int> > _faceDescriptors; std::map< int, std::pair<int,int> > _faceDescriptors;

11
src/NETGENPlugin/NETGENPlugin_NETGEN_2D_ONLY.cxx
View File

@ -292,8 +292,8 @@ static TError addSegmentsToMesh(netgen::Mesh& ngMesh,
edgeID = geom.emap.FindIndex( edge ); edgeID = geom.emap.FindIndex( edge );
posID = posShapeID; posID = posShapeID;
isInternalWire = ( edge.Orientation() == TopAbs_INTERNAL ); isInternalWire = ( edge.Orientation() == TopAbs_INTERNAL );
if ( onVertex ) // param on curve is different on each of two edges // if ( onVertex ) // param on curve is different on each of two edges
seg.epgeominfo[ iEnd ].dist = helper.GetNodeU( edge, pnt.node ); // seg.epgeominfo[ iEnd ].dist = helper.GetNodeU( edge, pnt.node );
} }
seg.epgeominfo[ iEnd ].edgenr = edgeID; // = geom.emap.FindIndex(edge); seg.epgeominfo[ iEnd ].edgenr = edgeID; // = geom.emap.FindIndex(edge);
} }
@ -301,15 +301,14 @@ static TError addSegmentsToMesh(netgen::Mesh& ngMesh,
ngMesh.AddSegment (seg); ngMesh.AddSegment (seg);
{ {
// restrict size of elements near the segment // restrict size of elements near the segment
netgen::Point3d ngP1(n->X(), n->Y(), n->Z()); SMESH_TNodeXYZ np1( n ), np2( uvPtVec[ i+1 ].node );
n = uvPtVec[ i+1 ].node;
netgen::Point3d ngP2(n->X(), n->Y(), n->Z());
// get an average size of adjacent segments to avoid sharp change of // get an average size of adjacent segments to avoid sharp change of
// element size (regression on issue 0020452, note 0010898) // element size (regression on issue 0020452, note 0010898)
int iPrev = SMESH_MesherHelper::WrapIndex( i-1, nbSegments ); int iPrev = SMESH_MesherHelper::WrapIndex( i-1, nbSegments );
int iNext = SMESH_MesherHelper::WrapIndex( i+1, nbSegments ); int iNext = SMESH_MesherHelper::WrapIndex( i+1, nbSegments );
double avgH = ( segLen[ iPrev ] + segLen[ i ] + segLen[ iNext ]) / 3; double avgH = ( segLen[ iPrev ] + segLen[ i ] + segLen[ iNext ]) / 3;
ngMesh.RestrictLocalH( netgen::Center( ngP1,ngP2), avgH );
NETGENPlugin_Mesher::RestrictLocalSize( ngMesh, 0.5*(np1+np2), avgH );
} }
#ifdef DUMP_SEGMENTS #ifdef DUMP_SEGMENTS
cout << "Segment: " << seg.edgenr << endl cout << "Segment: " << seg.edgenr << endl