diff options
Diffstat (limited to 'src/3rdparty/proj/PJ_wink2.c')
-rw-r--r-- | src/3rdparty/proj/PJ_wink2.c | 34 |
1 files changed, 0 insertions, 34 deletions
diff --git a/src/3rdparty/proj/PJ_wink2.c b/src/3rdparty/proj/PJ_wink2.c deleted file mode 100644 index 8eace009..00000000 --- a/src/3rdparty/proj/PJ_wink2.c +++ /dev/null @@ -1,34 +0,0 @@ -#define PROJ_PARMS__ \ - double cosphi1; -#define PJ_LIB__ -# include <projects.h> -PROJ_HEAD(wink2, "Winkel II") "\n\tPCyl., Sph., no inv.\n\tlat_1="; -#define MAX_ITER 10 -#define LOOP_TOL 1e-7 -#define TWO_D_PI 0.636619772367581343 -FORWARD(s_forward); /* spheroid */ - double k, V; - int i; - - xy.y = lp.phi * TWO_D_PI; - k = PI * sin(lp.phi); - lp.phi *= 1.8; - for (i = MAX_ITER; i ; --i) { - lp.phi -= V = (lp.phi + sin(lp.phi) - k) / - (1. + cos(lp.phi)); - if (fabs(V) < LOOP_TOL) - break; - } - if (!i) - lp.phi = (lp.phi < 0.) ? -HALFPI : HALFPI; - else - lp.phi *= 0.5; - xy.x = 0.5 * lp.lam * (cos(lp.phi) + P->cosphi1); - xy.y = FORTPI * (sin(lp.phi) + xy.y); - return (xy); -} -FREEUP; if (P) pj_dalloc(P); } -ENTRY0(wink2) - P->cosphi1 = cos(pj_param(P->params, "rlat_1").f); - P->es = 0.; P->inv = 0; P->fwd = s_forward; -ENDENTRY(P) |