From 46bfa242b7cfae5d0ac09d9c44086c4ef7d3ca87 Mon Sep 17 00:00:00 2001 From: Howard Butler Date: Thu, 1 Sep 2016 11:08:25 -0500 Subject: apply patch for #406 to improve Aitoff performance in some situations --- src/PJ_aitoff.c | 3 +-- 1 file changed, 1 insertion(+), 2 deletions(-) diff --git a/src/PJ_aitoff.c b/src/PJ_aitoff.c index 5cd717b1..ef6e7d71 100644 --- a/src/PJ_aitoff.c +++ b/src/PJ_aitoff.c @@ -121,8 +121,7 @@ static LP s_inverse (XY xy, PJ *P) { /* Spheroidal, inverse */ f1 -= xy.x; f2 -= xy.y; dl = (f2 * f1p - f1 * f2p) / (dp = f1p * f2l - f2p * f1l); dp = (f1 * f2l - f2 * f1l) / dp; - while (dl > M_PI) dl -= M_PI; /* set to interval [-M_PI, M_PI] */ - while (dl < -M_PI) dl += M_PI; /* set to interval [-M_PI, M_PI] */ + dl = fmod(dl, M_PI); /* set to interval [-M_PI, M_PI] */ lp.phi -= dp; lp.lam -= dl; } while ((fabs(dp) > EPSILON || fabs(dl) > EPSILON) && (iter++ < MAXITER)); if (lp.phi > M_PI_2) lp.phi -= 2.*(lp.phi-M_PI_2); /* correct if symmetrical solution for Aitoff */ -- cgit v1.2.3