diff options
| author | Kristian Evers <kristianevers@gmail.com> | 2016-10-14 14:06:06 +0200 |
|---|---|---|
| committer | Kristian Evers <kristianevers@gmail.com> | 2016-10-14 14:19:25 +0200 |
| commit | 217fb477d0359ae08149044f0181f822804fb293 (patch) | |
| tree | 081ce35b8accea54a7663495c6f7986072a97711 | |
| parent | 6f0e836a92a3c128ac18aaf597a12fcef6fa6470 (diff) | |
| download | PROJ-217fb477d0359ae08149044f0181f822804fb293.tar.gz PROJ-217fb477d0359ae08149044f0181f822804fb293.zip | |
Replace calls to fabls with fabs (not using long doubles anyway). Fixes #426.
| -rw-r--r-- | src/PJ_healpix.c | 12 |
1 files changed, 6 insertions, 6 deletions
diff --git a/src/PJ_healpix.c b/src/PJ_healpix.c index c5374df7..318f42bc 100644 --- a/src/PJ_healpix.c +++ b/src/PJ_healpix.c @@ -231,7 +231,7 @@ double auth_lat(PJ *P, double alpha, int inverse) { double qp = Q->qp; double ratio = q/qp; - if (fabsl(ratio) > 1) { + if (fabs(ratio) > 1) { /* Rounding error. */ ratio = pj_sign(ratio); } @@ -254,12 +254,12 @@ XY healpix_sphere(LP lp) { XY xy; /* equatorial region */ - if ( fabsl(phi) <= phi0) { + if ( fabs(phi) <= phi0) { xy.x = lam; xy.y = 3*M_PI/8*sin(phi); } else { double lamc; - double sigma = sqrt(3*(1 - fabsl(sin(phi)))); + double sigma = sqrt(3*(1 - fabs(sin(phi)))); double cn = floor(2*lam / M_PI + 2); if (cn >= 4) { cn = 3; @@ -282,17 +282,17 @@ LP healpix_sphere_inverse(XY xy) { double y0 = M_FORTPI; /* Equatorial region. */ - if (fabsl(y) <= y0) { + if (fabs(y) <= y0) { lp.lam = x; lp.phi = asin(8*y/(3*M_PI)); - } else if (fabsl(y) < M_HALFPI) { + } else if (fabs(y) < M_HALFPI) { double cn = floor(2*x/M_PI + 2); double xc, tau; if (cn >= 4) { cn = 3; } xc = -3*M_FORTPI + M_HALFPI*cn; - tau = 2.0 - 4*fabsl(y)/M_PI; + tau = 2.0 - 4*fabs(y)/M_PI; lp.lam = xc + (x - xc)/tau; lp.phi = pj_sign(y)*asin(1.0 - pow(tau, 2)/3.0); } else { |