ForkMaster/smesh
1
0
Fork 0
mirror of https://git.salome-platform.org/gitpub/modules/smesh.git synced 2025-02-07 21:14:17 +05:00
Code Issues Packages Projects Releases Wiki Activity

Fix size of the returned vector in planar case

Browse Source
This commit is contained in:
eap 2017-10-03 20:52:49 +03:00
parent ba27d23aec
commit 4b2ed66f4e
1 changed files with 6 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
src/SMESH/SMESH_MeshEditor.cxx
View File

@ -13209,6 +13209,7 @@ void SMESH_MeshEditor::MakePolyLine( TListOfPolySegments& theSegments,
// get cutting planes // get cutting planes
std::vector< bool > isVectorOK( theSegments.size(), true ); std::vector< bool > isVectorOK( theSegments.size(), true );
const double planarCoef = 0.333; // plane height in planar case
for ( size_t iSeg = 0; iSeg < theSegments.size(); ++iSeg ) for ( size_t iSeg = 0; iSeg < theSegments.size(); ++iSeg )
{ {
@ -13231,7 +13232,7 @@ void SMESH_MeshEditor::MakePolyLine( TListOfPolySegments& theSegments,
if ( polySeg.myMidProjPoint.Distance( pMid ) < Precision::Confusion() ) if ( polySeg.myMidProjPoint.Distance( pMid ) < Precision::Confusion() )
{ {
SMESH_MeshAlgos::FaceNormal( face, const_cast< gp_XYZ& >( polySeg.myVector.XYZ() )); SMESH_MeshAlgos::FaceNormal( face, const_cast< gp_XYZ& >( polySeg.myVector.XYZ() ));
polySeg.myMidProjPoint = pMid + polySeg.myVector.XYZ() * ( p1 - p2 ).Modulus() * 0.333; polySeg.myMidProjPoint = pMid + polySeg.myVector.XYZ() * ( p1 - p2 ).Modulus() * planarCoef;
} }
else else
{ {
@ -13288,6 +13289,7 @@ void SMESH_MeshEditor::MakePolyLine( TListOfPolySegments& theSegments,
// return a vector // return a vector
gp_XYZ pMid = 0.5 * ( path.myPoints[0] + path.myPoints.back() );
if ( isVectorOK[ iSeg ]) if ( isVectorOK[ iSeg ])
{ {
// find the most distance point of a path // find the most distance point of a path
@ -13301,8 +13303,10 @@ void SMESH_MeshEditor::MakePolyLine( TListOfPolySegments& theSegments,
theSegments[iSeg].myMidProjPoint = path.myPoints[iP]; theSegments[iSeg].myMidProjPoint = path.myPoints[iP];
} }
} }
if ( maxDist < Precision::Confusion() ) // planar case
theSegments[iSeg].myMidProjPoint =
pMid + theSegments[iSeg].myVector.XYZ().Normalized() * path.myLength * planarCoef;
} }
gp_XYZ pMid = 0.5 * ( path.myPoints[0] + path.myPoints.back() );
theSegments[iSeg].myVector = gp_Vec( pMid, theSegments[iSeg].myMidProjPoint ); theSegments[iSeg].myVector = gp_Vec( pMid, theSegments[iSeg].myMidProjPoint );
} }