From 35c0972b0c958d45ad1e3a375be2d3a8b951586a Mon Sep 17 00:00:00 2001 From: Glynn Clements Date: Fri, 14 May 2010 18:55:26 +0100 Subject: Fix inverse stereo projection on an ellipsoid. The problem is that the setup function doesn't initialise all of the values that it needs for the ellipsoid, inverse, equatorial case. In the forward and spherical inverse cases, the equatorial case (|phi0| < 1e-10) is a simplification of the oblique case, and requires fewer parameters. However, for the ellipsoid inverse case (e_inverse), the oblique calculation is used, but the setup function doesn't initialise all of the parameters. --- src/PJ_stere.c | 2 -- 1 file changed, 2 deletions(-) (limited to 'src') diff --git a/src/PJ_stere.c b/src/PJ_stere.c index 52ec1df1..a0524627 100644 --- a/src/PJ_stere.c +++ b/src/PJ_stere.c @@ -192,8 +192,6 @@ setup(PJ *P) { /* general initialization */ } break; case EQUIT: - P->akm1 = 2. * P->k0; - break; case OBLIQ: t = sin(P->phi0); X = 2. * atan(ssfn_(P->phi0, t, P->e)) - HALFPI; -- cgit v1.2.3