ForkMaster/smesh
1
0
Fork 0
mirror of https://git.salome-platform.org/gitpub/modules/smesh.git synced 2024-11-13 17:18:35 +05:00
Code Issues Packages Projects Releases Wiki Activity

Fix failure of HEXOTIC/doc/salome/examples/hexoticdemo.py

Browse Source
This commit is contained in:
eap 2020-05-27 15:46:46 +03:00
parent 19720c1ef8
commit 2fcf0d1b8e
2 changed files with 57 additions and 64 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

93
src/SMESHUtils/SMESH_ControlPnt.cxx
View File

@ -58,9 +58,10 @@ namespace SMESHUtils
const double& theSize, const double& theSize,
std::vector<ControlPnt>& thePoints ); std::vector<ControlPnt>& thePoints );
std::vector<gp_Pnt> computePointsForSplitting( const gp_Pnt& p1, void computePointsForSplitting( const gp_Pnt& p1,
const gp_Pnt& p2, const gp_Pnt& p2,
const gp_Pnt& p3 ); const gp_Pnt& p3,
gp_Pnt midPoints[3]);
gp_Pnt tangencyPoint(const gp_Pnt& p1, gp_Pnt tangencyPoint(const gp_Pnt& p1,
const gp_Pnt& p2, const gp_Pnt& p2,
const gp_Pnt& Center); const gp_Pnt& Center);
@ -164,8 +165,7 @@ void SMESHUtils::createPointsSampleFromFace( const TopoDS_Face& theFace,
// Get triangles // Get triangles
int nbTriangles = aTri->NbTriangles(); int nbTriangles = aTri->NbTriangles();
Poly_Array1OfTriangle triangles(1,nbTriangles); const Poly_Array1OfTriangle& triangles = aTri->Triangles();
triangles=aTri->Triangles();
// GetNodes // GetNodes
int nbNodes = aTri->NbNodes(); int nbNodes = aTri->NbNodes();
@ -173,9 +173,10 @@ void SMESHUtils::createPointsSampleFromFace( const TopoDS_Face& theFace,
nodes = aTri->Nodes(); nodes = aTri->Nodes();
// Iterate on triangles and subdivide them // Iterate on triangles and subdivide them
for(int i=1; i<=nbTriangles; i++) thePoints.reserve( thePoints.size() + nbTriangles );
for ( int i = 1; i <= nbTriangles; i++ )
{ {
Poly_Triangle aTriangle = triangles.Value(i); const Poly_Triangle& aTriangle = triangles.Value(i);
gp_Pnt p1 = nodes.Value(aTriangle.Value(1)); gp_Pnt p1 = nodes.Value(aTriangle.Value(1));
gp_Pnt p2 = nodes.Value(aTriangle.Value(2)); gp_Pnt p2 = nodes.Value(aTriangle.Value(2));
gp_Pnt p3 = nodes.Value(aTriangle.Value(3)); gp_Pnt p3 = nodes.Value(aTriangle.Value(3));
@ -184,7 +185,7 @@ void SMESHUtils::createPointsSampleFromFace( const TopoDS_Face& theFace,
p2.Transform(aTrsf); p2.Transform(aTrsf);
p3.Transform(aTrsf); p3.Transform(aTrsf);
subdivideTriangle(p1, p2, p3, theSize, thePoints); subdivideTriangle( p1, p2, p3, theSize, thePoints );
} }
} }
@ -221,10 +222,8 @@ void SMESHUtils::createPointsSampleFromSolid( const TopoDS_Solid& theSolid,
// Step2 : for each face of theSolid: // Step2 : for each face of theSolid:
std::set<double> intersections; std::set<double> intersections;
std::set<double>::iterator it = intersections.begin();
TopExp_Explorer Ex; for ( TopExp_Explorer Ex( theSolid, TopAbs_FACE ); Ex.More(); Ex.Next() )
for (Ex.Init(theSolid,TopAbs_FACE); Ex.More(); Ex.Next())
{ {
// check if there is an intersection // check if there is an intersection
IntCurvesFace_Intersector anIntersector(TopoDS::Face(Ex.Current()), Precision::Confusion()); IntCurvesFace_Intersector anIntersector(TopoDS::Face(Ex.Current()), Precision::Confusion());
@ -232,13 +231,13 @@ void SMESHUtils::createPointsSampleFromSolid( const TopoDS_Solid& theSolid,
// get the intersection's parameter and store it // get the intersection's parameter and store it
int nbPoints = anIntersector.NbPnt(); int nbPoints = anIntersector.NbPnt();
for(int i = 0 ; i < nbPoints ; i++ ) for ( int i = 0 ; i < nbPoints; i++ )
{ {
it = intersections.insert( it, anIntersector.WParameter(i+1) ); intersections.insert( anIntersector.WParameter(i+1) );
} }
} }
// Step3 : go through the line chunk by chunk // Step3 : go through the line chunk by chunk
if ( intersections.begin() != intersections.end() ) if ( intersections.size() > 1 )
{ {
std::set<double>::iterator intersectionsIterator=intersections.begin(); std::set<double>::iterator intersectionsIterator=intersections.begin();
double first = *intersectionsIterator; double first = *intersectionsIterator;
@ -254,9 +253,9 @@ void SMESHUtils::createPointsSampleFromSolid( const TopoDS_Solid& theSolid,
double localStep = (second -first) / ceil( (second - first) / step ); double localStep = (second -first) / ceil( (second - first) / step );
for ( double z = Zmin + first; z < Zmin + second; z = z + localStep ) for ( double z = Zmin + first; z < Zmin + second; z = z + localStep )
{ {
thePoints.push_back(ControlPnt( x, y, z, theSize )); thePoints.emplace_back( x, y, z, theSize );
} }
thePoints.push_back(ControlPnt( x, y, Zmin + second, theSize )); thePoints.emplace_back( x, y, Zmin + second, theSize );
} }
first = second; first = second;
innerPoints = !innerPoints; innerPoints = !innerPoints;
@ -288,45 +287,45 @@ void SMESHUtils::subdivideTriangle( const gp_Pnt& p1,
// and the distance between two mass centers of two neighbouring triangles // and the distance between two mass centers of two neighbouring triangles
// sharing an edge is < 2 * 1/2 * S = S // sharing an edge is < 2 * 1/2 * S = S
// If the traingles share a Vertex and no Edge the distance of the mass centers // If the traingles share a Vertex and no Edge the distance of the mass centers
// to the Vertices is 2*D < S so the mass centers are distant of less than 2*S // to the Vertices is 2*D < S so the mass centers are distant of less than 2*S
double threshold = sqrt( 3. ) * theSize; double threshold = sqrt( 3. ) * theSize;
if ( (p1.Distance(p2) > threshold || if ( p1.Distance(p2) > threshold ||
p2.Distance(p3) > threshold || p2.Distance(p3) > threshold ||
p3.Distance(p1) > threshold)) p3.Distance(p1) > threshold )
{ try
std::vector<gp_Pnt> midPoints = computePointsForSplitting(p1, p2, p3); {
gp_Pnt midPoints[3];
computePointsForSplitting( p1, p2, p3, midPoints );
subdivideTriangle( midPoints[0], midPoints[1], midPoints[2], theSize, thePoints ); subdivideTriangle( midPoints[0], midPoints[1], midPoints[2], theSize, thePoints );
subdivideTriangle( midPoints[0], p2, midPoints[1], theSize, thePoints ); subdivideTriangle( midPoints[0], p2, midPoints[1], theSize, thePoints );
subdivideTriangle( midPoints[2], midPoints[1], p3, theSize, thePoints ); subdivideTriangle( midPoints[2], midPoints[1], p3, theSize, thePoints );
subdivideTriangle( p1, midPoints[0], midPoints[2], theSize, thePoints ); subdivideTriangle( p1, midPoints[0], midPoints[2], theSize, thePoints );
} return;
else }
{ catch (...)
double x = (p1.X() + p2.X() + p3.X()) / 3 ; {
double y = (p1.Y() + p2.Y() + p3.Y()) / 3 ; }
double z = (p1.Z() + p2.Z() + p3.Z()) / 3 ;
ControlPnt massCenter( x ,y ,z, theSize ); gp_Pnt massCenter = ( p1.XYZ() + p2.XYZ() + p3.XYZ() ) / 3.;
thePoints.push_back( massCenter ); thePoints.emplace_back( massCenter, theSize );
}
} }
//================================================================================ //================================================================================
/*! /*!
* \brief Returns the appropriate points for splitting a triangle * \brief Returns the appropriate points for splitting a triangle
* \brief the tangency points of the incircle are used in order to have mostly * the tangency points of the incircle are used in order to have mostly
* \brief well-shaped sub-triangles * well-shaped sub-triangles
*/ */
//================================================================================ //================================================================================
std::vector<gp_Pnt> SMESHUtils::computePointsForSplitting( const gp_Pnt& p1, void SMESHUtils::computePointsForSplitting( const gp_Pnt& p1,
const gp_Pnt& p2, const gp_Pnt& p2,
const gp_Pnt& p3 ) const gp_Pnt& p3,
gp_Pnt midPoints[3])
{ {
std::vector<gp_Pnt> midPoints;
//Change coordinates //Change coordinates
gp_Trsf Trsf_1; // Identity transformation gp_Trsf Trsf_1; // Identity transformation
gp_Ax3 reference_system(gp::Origin(), gp::DZ(), gp::DX()); // OXY gp_Ax3 reference_system(gp::Origin(), gp::DZ(), gp::DX()); // OXY
@ -359,15 +358,11 @@ std::vector<gp_Pnt> SMESHUtils::computePointsForSplitting( const gp_Pnt& p1,
gp_Pnt T2 = tangencyPoint( B, C, Center); gp_Pnt T2 = tangencyPoint( B, C, Center);
gp_Pnt T3 = tangencyPoint( C, A, Center); gp_Pnt T3 = tangencyPoint( C, A, Center);
gp_Pnt p1_2 = T1.Transformed(Trsf_1.Inverted()); midPoints[0] = T1.Transformed(Trsf_1.Inverted());
gp_Pnt p2_3 = T2.Transformed(Trsf_1.Inverted()); midPoints[1] = T2.Transformed(Trsf_1.Inverted());
gp_Pnt p3_1 = T3.Transformed(Trsf_1.Inverted()); midPoints[2] = T3.Transformed(Trsf_1.Inverted());
midPoints.push_back(p1_2); return;
midPoints.push_back(p2_3);
midPoints.push_back(p3_1);
return midPoints;
} }
//================================================================================ //================================================================================

28
src/SMESHUtils/SMESH_ControlPnt.hxx
View File

@ -42,11 +42,9 @@ namespace SMESHUtils
{ {
ControlPnt() ControlPnt()
: gp_Pnt(), size(0) {} : gp_Pnt(), size(0) {}
ControlPnt( const gp_Pnt& aPnt, double theSize) ControlPnt( const gp_Pnt& aPnt, double theSize=0)
: gp_Pnt( aPnt ), size( theSize ) {} : gp_Pnt( aPnt ), size( theSize ) {}
ControlPnt(double theX,double theY,double theZ) ControlPnt(double theX,double theY,double theZ, double theSize=0)
: gp_Pnt(theX, theY, theZ), size(0) {}
ControlPnt(double theX,double theY,double theZ, double theSize)
: gp_Pnt(theX, theY, theZ), size( theSize ) {} : gp_Pnt(theX, theY, theZ), size( theSize ) {}
double Size() const { return size; }; double Size() const { return size; };
@ -57,20 +55,20 @@ namespace SMESHUtils
// Functions to get sample point from shapes // Functions to get sample point from shapes
SMESHUtils_EXPORT void createControlPoints( const TopoDS_Shape& theShape, SMESHUtils_EXPORT void createControlPoints( const TopoDS_Shape& theShape,
const double& theSize, const double& theSize,
std::vector< ControlPnt >& thePoints ); std::vector< ControlPnt >& thePoints );
SMESHUtils_EXPORT void createPointsSampleFromEdge( const TopoDS_Edge& theEdge, SMESHUtils_EXPORT void createPointsSampleFromEdge( const TopoDS_Edge& theEdge,
const double& theSize, const double& theSize,
std::vector<ControlPnt>& thePoints ); std::vector<ControlPnt>& thePoints );
SMESHUtils_EXPORT void createPointsSampleFromFace( const TopoDS_Face& theFace, SMESHUtils_EXPORT void createPointsSampleFromFace( const TopoDS_Face& theFace,
const double& theSize, const double& theSize,
std::vector<ControlPnt>& thePoints ); std::vector<ControlPnt>& thePoints );
SMESHUtils_EXPORT void createPointsSampleFromSolid( const TopoDS_Solid& theSolid, SMESHUtils_EXPORT void createPointsSampleFromSolid( const TopoDS_Solid& theSolid,
const double& theSize, const double& theSize,
std::vector<ControlPnt>& thePoints ); std::vector<ControlPnt>& thePoints );
} }
#endif #endif