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

modified spline projection

Browse Source
This commit is contained in:
Joachim Schöberl 2016-09-07 09:02:57 +02:00
parent 6c44a273ce
commit 0eaedadf10
1 changed files with 15 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

16
libsrc/gprim/spline.cpp
View File

@ -197,11 +197,25 @@ namespace netgen
{
double t_old = -1;
/*
if(proj_latest_t > 0. && proj_latest_t < 1.)
t = proj_latest_t;
else
t = 0.5;
*/
double tmin = 1;
double dist_min2 = Dist2 (GetPoint(tmin), point);
for (double ti = 0; ti < 0.99; ti += 0.25)
{
double di = Dist2(GetPoint(ti), point);
if (di < dist_min2)
{
tmin = ti;
dist_min2 = di;
}
}
t = tmin;
Point<D> phi;
Vec<D> phip,phipp,phimp;