ForkMaster/netgen
1
0
Fork 0
mirror of https://github.com/NGSolve/netgen.git synced 2024-12-24 21:10:33 +05:00
Code Issues Packages Projects Releases Wiki Activity

store cumulated angle

Browse Source
This commit is contained in:
Christopher Lackner 2021-04-15 22:48:16 +02:00
parent 21ef833bbd
commit 087a830a67
1 changed files with 4 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
libsrc/csg/extrusion.cpp
View File

@ -42,6 +42,7 @@ namespace netgen
} }
} }
double cum_angle = 0.;
for(auto i : Range(path->GetSplines())) for(auto i : Range(path->GetSplines()))
{ {
const auto& sp = path->GetSpline(i); const auto& sp = path->GetSpline(i);
@ -49,7 +50,8 @@ namespace netgen
t1.Normalize(); t1.Normalize();
auto t2 = sp.GetTangent(1.); auto t2 = sp.GetTangent(1.);
t2.Normalize(); t2.Normalize();
angles.Append(acos(t1 * t2)); cum_angle += acos(t1 * t2);
angles.Append(cum_angle);
} }
profile->GetCoeff(profile_spline_coeff); profile->GetCoeff(profile_spline_coeff);
@ -678,11 +680,7 @@ namespace netgen
double h, int& zone) const double h, int& zone) const
{ {
Surface::ToPlane(p3d, p2d, h, zone); Surface::ToPlane(p3d, p2d, h, zone);
double angle = 0; double angle = angles[tangential_plane_seg] - angles[latest_seg];
for(int i = latest_seg; i < tangential_plane_seg; i++)
angle += angles[i];
for(int i = tangential_plane_seg; i < latest_seg; i++)
angle -= angles[i];
if(fabs(angle) > 3.14/2.) if(fabs(angle) > 3.14/2.)
zone = -1; zone = -1;
} }