diff options
Diffstat (limited to 'src/PJ_august.c')
| -rw-r--r-- | src/PJ_august.c | 58 |
1 files changed, 55 insertions, 3 deletions
diff --git a/src/PJ_august.c b/src/PJ_august.c index 496fc3c3..83463e9f 100644 --- a/src/PJ_august.c +++ b/src/PJ_august.c @@ -1,8 +1,13 @@ #define PJ_LIB__ #include <projects.h> + + PROJ_HEAD(august, "August Epicycloidal") "\n\tMisc Sph, no inv."; #define M 1.333333333333333 -FORWARD(s_forward); /* spheroid */ + + +static XY s_forward (LP lp, PJ *P) { /* Spheroidal, forward */ + XY xy = {0.0,0.0}; double t, c1, c, x1, x12, y1, y12; (void) P; @@ -15,5 +20,52 @@ FORWARD(s_forward); /* spheroid */ xy.y = M * y1 * (3. + 3. * x12 - y12); return (xy); } -FREEUP; if (P) pj_dalloc(P); } -ENTRY0(august) P->inv = 0; P->fwd = s_forward; P->es = 0.; ENDENTRY(P) + + + +static void *freeup_new (PJ *P) { /* Destructor */ + return pj_dealloc(P); +} + +static void freeup (PJ *P) { + freeup_new (P); + return; +} + + +PJ *PROJECTION(august) { + P->inv = 0; + P->fwd = s_forward; + P->es = 0.; + return P; +} + +#ifdef PJ_OMIT_SELFTEST +int pj_august_selftest (void) {return 0;} +#else + +int pj_august_selftest (void) { + double tolerance_lp = 1e-10; + double tolerance_xy = 1e-7; + + char s_args[] = {"+proj=august +a=6400000 +lat_1=0 +lat_2=2"}; + + LP fwd_in[] = { + { 2, 1}, + { 2,-1}, + {-2, 1}, + {-2,-1} + }; + + XY s_fwd_expect[] = { + {223404.97818097242, 111722.34028976287}, + {223404.97818097242, -111722.34028976287}, + {-223404.97818097242, 111722.34028976287}, + {-223404.97818097242, -111722.34028976287}, + }; + + return pj_generic_selftest (0, s_args, tolerance_xy, tolerance_lp, 4, 4, fwd_in, 0, s_fwd_expect, 0, 0, 0); +} + + +#endif |