diff options
| author | Even Rouault <even.rouault@spatialys.com> | 2017-02-28 01:02:25 +0100 |
|---|---|---|
| committer | Even Rouault <even.rouault@spatialys.com> | 2017-02-28 17:02:53 +0100 |
| commit | 21d16428760469717fd25ebf9d03c286db935c66 (patch) | |
| tree | 0747dd9d46c106e11a4e85d0caaacc002f2557e9 /src/PJ_mod_ster.c | |
| parent | 8bf1acf2692e7628c2a7b41ad53d876f9b9b09e9 (diff) | |
| download | PROJ-21d16428760469717fd25ebf9d03c286db935c66.tar.gz PROJ-21d16428760469717fd25ebf9d03c286db935c66.zip | |
Enable cppcheck and fix related mostly false-positive warnings
Diffstat (limited to 'src/PJ_mod_ster.c')
| -rw-r--r-- | src/PJ_mod_ster.c | 3 |
1 files changed, 2 insertions, 1 deletions
diff --git a/src/PJ_mod_ster.c b/src/PJ_mod_ster.c index bcf2bf9c..f8c90bc4 100644 --- a/src/PJ_mod_ster.c +++ b/src/PJ_mod_ster.c @@ -46,7 +46,7 @@ static LP e_inverse (XY xy, PJ *P) { /* Ellipsoidal, inverse */ struct pj_opaque *Q = P->opaque; int nn; COMPLEX p, fxy, fpxy, dp; - double den, rh = 0.0, z, sinz = 0.0, cosz = 0.0, chi, phi = 0.0, dphi, esphi; + double den, rh = 0.0, z, sinz = 0.0, cosz = 0.0, chi, phi = 0.0, esphi; p.r = xy.x; p.i = xy.y; @@ -79,6 +79,7 @@ static LP e_inverse (XY xy, PJ *P) { /* Ellipsoidal, inverse */ chi = aasin(P->ctx, cosz * Q->schio + p.i * sinz * Q->cchio / rh); phi = chi; for (nn = 20; nn ;--nn) { + double dphi; esphi = P->e * sin(phi); dphi = 2. * atan(tan((M_HALFPI + chi) * .5) * pow((1. + esphi) / (1. - esphi), P->e * .5)) - M_HALFPI - phi; |