ForkMaster/netgen
1
0
Fork 0
mirror of https://github.com/NGSolve/netgen.git synced 2024-09-23 02:02:49 +05:00
Code Issues Packages Projects Releases Wiki Activity

Merge branch 'edge_midpoint_without_geo_edge' into 'master'

Browse Source 
if PointBetweenEdge gets invalid edgenr, return midpoint

See merge request jschoeberl/netgen!419
This commit is contained in:
Joachim Schöberl 2021-09-17 10:00:31 +00:00
commit 46a6062bdf
1 changed files with 5 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

9
libsrc/occ/occgeom.cpp
View File

@ -1442,12 +1442,13 @@ namespace netgen
{ {
double s0, s1; double s0, s1;
Point<3> hnewp = p1+secpoint*(p2-p1); newp = p1+secpoint*(p2-p1);
gp_Pnt pnt(hnewp(0), hnewp(1), hnewp(2)); if(ap1.edgenr > emap.Size() || ap1.edgenr == 0)
return;
gp_Pnt pnt(newp(0), newp(1), newp(2));
GeomAPI_ProjectPointOnCurve proj(pnt, BRep_Tool::Curve(TopoDS::Edge(emap(ap1.edgenr)), s0, s1)); GeomAPI_ProjectPointOnCurve proj(pnt, BRep_Tool::Curve(TopoDS::Edge(emap(ap1.edgenr)), s0, s1));
pnt = proj.NearestPoint(); pnt = proj.NearestPoint();
hnewp = Point<3> (pnt.X(), pnt.Y(), pnt.Z()); newp = Point<3> (pnt.X(), pnt.Y(), pnt.Z());
newp = hnewp;
newgi = ap1; newgi = ap1;
}; };