diff options
| author | Kristian Evers <kristianevers@gmail.com> | 2016-04-18 12:32:21 +0200 |
|---|---|---|
| committer | Kristian Evers <kristianevers@gmail.com> | 2016-04-18 12:32:21 +0200 |
| commit | 8531aac07f11691fafa2704d357373459b6e9075 (patch) | |
| tree | 1e630d7f17b94e598ac66df1e0729e0a5afdcf1d /src | |
| parent | dbba67bd749d6cf8080658fecc51fde90502f4d1 (diff) | |
| download | PROJ-8531aac07f11691fafa2704d357373459b6e9075.tar.gz PROJ-8531aac07f11691fafa2704d357373459b6e9075.zip | |
Converted gins8. Expanded tabs.
Diffstat (limited to 'src')
| -rw-r--r-- | src/PJ_aea.c | 1 | ||||
| -rw-r--r-- | src/PJ_gins8.c | 87 |
2 files changed, 72 insertions, 16 deletions
diff --git a/src/PJ_aea.c b/src/PJ_aea.c index c5f808a5..d99217bb 100644 --- a/src/PJ_aea.c +++ b/src/PJ_aea.c @@ -360,7 +360,6 @@ int pj_eqc_selftest (void) {return 10000;} int pj_eqdc_selftest (void) {return 10000;} int pj_etmerc_selftest (void) {return 10000;} int pj_geocent_selftest (void) {return 10000;} -int pj_gins8_selftest (void) {return 10000;} int pj_gn_sinu_selftest (void) {return 10000;} int pj_gnom_selftest (void) {return 10000;} int pj_goode_selftest (void) {return 10000;} diff --git a/src/PJ_gins8.c b/src/PJ_gins8.c index b7f38ad9..3eae7efa 100644 --- a/src/PJ_gins8.c +++ b/src/PJ_gins8.c @@ -1,18 +1,75 @@ #define PJ_LIB__ -# include <projects.h> +#include <projects.h> + PROJ_HEAD(gins8, "Ginsburg VIII (TsNIIGAiK)") "\n\tPCyl, Sph., no inv."; -#define Cl 0.000952426 -#define Cp 0.162388 -#define C12 0.08333333333333333 -FORWARD(s_forward); /* spheroid */ - double t = lp.phi * lp.phi; - (void) P; - - xy.y = lp.phi * (1. + t * C12); - xy.x = lp.lam * (1. - Cp * t); - t = lp.lam * lp.lam; - xy.x *= (0.87 - Cl * t * t); - return (xy); + +#define Cl 0.000952426 +#define Cp 0.162388 +#define C12 0.08333333333333333 + + +static XY s_forward (LP lp, PJ *P) { /* Spheroidal, forward */ + XY xy = {0.0,0.0}; + double t = lp.phi * lp.phi; + (void) P; + + xy.y = lp.phi * (1. + t * C12); + xy.x = lp.lam * (1. - Cp * t); + t = lp.lam * lp.lam; + xy.x *= (0.87 - Cl * t * t); + + return xy; +} + + +static void *freeup_new (PJ *P) { /* Destructor */ + if (0==P) + return 0; + + return pj_dealloc(P); +} + +static void freeup (PJ *P) { + freeup_new (P); + return; +} + + +PJ *PROJECTION(gins8) { + P->es = 0.0; + P->inv = 0; + P->fwd = s_forward; + + return P; } -FREEUP; if (P) pj_dalloc(P); } -ENTRY0(gins8) P->es = 0.; P->inv = 0; P->fwd = s_forward; ENDENTRY(P) + + +#ifdef PJ_OMIT_SELFTEST +int pj_gins8_selftest (void) {return 0;} +#else + +int pj_gins8_selftest (void) { + double tolerance_lp = 1e-10; + double tolerance_xy = 1e-7; + + char s_args[] = {"+proj=gins8 +a=6400000 +lat_1=0.5 +lat_2=2"}; + + LP fwd_in[] = { + { 2, 1}, + { 2,-1}, + {-2, 1}, + {-2,-1} + }; + + XY s_fwd_expect[] = { + { 194350.25093959007, 111703.90763533533}, + { 194350.25093959007, -111703.90763533533}, + {-194350.25093959007, 111703.90763533533}, + {-194350.25093959007, -111703.90763533533}, + }; + + return pj_generic_selftest (0, s_args, tolerance_xy, tolerance_lp, 4, 4, fwd_in, 0, s_fwd_expect, 0, 0, 0); +} + + +#endif |
