diff options
| author | Kristian Evers <kristianevers@gmail.com> | 2016-07-13 19:31:05 +0200 |
|---|---|---|
| committer | Kristian Evers <kristianevers@gmail.com> | 2016-07-13 19:31:05 +0200 |
| commit | 44e2818ac9dcb9e25d0057d106d3aa07d35c05bd (patch) | |
| tree | 2a8f8c01c28f2489aed4e79d8dd27d6f55e54dc7 /src | |
| parent | 196a1bb848dbd8f080d8ad6769f478c17c96d409 (diff) | |
| download | PROJ-44e2818ac9dcb9e25d0057d106d3aa07d35c05bd.tar.gz PROJ-44e2818ac9dcb9e25d0057d106d3aa07d35c05bd.zip | |
Fixed a few errors that occured during refactoring
Diffstat (limited to 'src')
| -rw-r--r-- | src/PJ_healpix.c | 10 |
1 files changed, 5 insertions, 5 deletions
diff --git a/src/PJ_healpix.c b/src/PJ_healpix.c index e1b48daa..1db484c1 100644 --- a/src/PJ_healpix.c +++ b/src/PJ_healpix.c @@ -162,8 +162,8 @@ int in_image(double x, double y, int proj, int north_square, int south_square) { double healpixVertsJit[][2] = { {-1.0*M_PI - EPS, M_FORTPI}, {-3.0*M_FORTPI, M_HALFPI + EPS}, - {M_HALFPI, M_FORTPI + EPS}, - {M_FORTPI, M_HALFPI + EPS}, + {-1.0*M_HALFPI, M_FORTPI + EPS}, + {-1.0*M_FORTPI, M_HALFPI + EPS}, {0.0, M_FORTPI + EPS}, {M_FORTPI, M_HALFPI + EPS}, {M_HALFPI, M_FORTPI + EPS}, @@ -192,7 +192,7 @@ int in_image(double x, double y, int proj, int north_square, int south_square) { rhealpixVertsJit[0][0] = -1.0*M_PI - EPS; rhealpixVertsJit[0][1] = M_FORTPI + EPS; rhealpixVertsJit[1][0] = -1.0*M_PI + north_square*M_HALFPI- EPS; - rhealpixVertsJit[1][2] = M_FORTPI + EPS; + rhealpixVertsJit[1][1] = M_FORTPI + EPS; rhealpixVertsJit[2][0] = -1.0*M_PI + north_square*M_HALFPI- EPS; rhealpixVertsJit[2][1] = 3*M_FORTPI + EPS; rhealpixVertsJit[3][0] = -1.0*M_PI + (north_square + 1.0)*M_HALFPI + EPS; @@ -266,7 +266,7 @@ XY healpix_sphere(LP lp) { if (cn >= 4) { cn = 3; } - lamc = -3*M_FORTPI + (M_HALFPI)*cn; + lamc = -3*M_FORTPI + M_HALFPI*cn; xy.x = lamc + (lam - lamc)*sigma; xy.y = pj_sign(phi)*M_FORTPI*(2 - sigma); } @@ -293,7 +293,7 @@ LP healpix_sphere_inverse(XY xy) { if (cn >= 4) { cn = 3; } - xc = -3.0*M_FORTPI + (M_HALFPI)*cn; + xc = -3.0*M_FORTPI + M_HALFPI*cn; tau = 2.0 - 4.0*fabsl(y)/M_PI; lp.lam = xc + (x - xc)/tau; lp.phi = pj_sign(y)*asin(1.0 - pow(tau , 2.0)/3.0); |