diff options
| author | Thomas Knudsen <knudsen.thomas@gmail.com> | 2016-05-20 00:44:56 +0200 |
|---|---|---|
| committer | Thomas Knudsen <knudsen.thomas@gmail.com> | 2016-05-20 00:44:56 +0200 |
| commit | 8b4c079485a68593b810ebfc14848ef6bc7864fa (patch) | |
| tree | 5f8bf20027529b6365caba608d5372e19fce4fb3 /src/PJ_larr.c | |
| parent | 1384b978373cead004211d8548518373d1f4950c (diff) | |
| parent | 605dc6f3b47bf72c61712b275284212c66e541d0 (diff) | |
| download | PROJ-8b4c079485a68593b810ebfc14848ef6bc7864fa.tar.gz PROJ-8b4c079485a68593b810ebfc14848ef6bc7864fa.zip | |
Merge pull request #373 from busstoptaktik/sdfe-refactor-macros--and-repair-generic-constructor-bug
First steps toward simplified macros/internals
Diffstat (limited to 'src/PJ_larr.c')
| -rw-r--r-- | src/PJ_larr.c | 77 |
1 files changed, 67 insertions, 10 deletions
diff --git a/src/PJ_larr.c b/src/PJ_larr.c index 118bef00..f55289b1 100644 --- a/src/PJ_larr.c +++ b/src/PJ_larr.c @@ -1,13 +1,70 @@ -#define PROJ_PARMS__ #define PJ_LIB__ -#include <projects.h> +#include <projects.h> + PROJ_HEAD(larr, "Larrivee") "\n\tMisc Sph, no inv."; + #define SIXTH .16666666666666666 -FORWARD(s_forward); /* sphere */ - (void) P; - xy.x = 0.5 * lp.lam * (1. + sqrt(cos(lp.phi))); - xy.y = lp.phi / (cos(0.5 * lp.phi) * cos(SIXTH * lp.lam)); - return (xy); -} -FREEUP; if (P) pj_dalloc(P); } -ENTRY0(larr) P->fwd = s_forward; P->inv = 0; P->es = 0.; ENDENTRY(P) + + +static XY s_forward (LP lp, PJ *P) { /* Spheroidal, forward */ + XY xy = {0.0,0.0}; + (void) P; + + xy.x = 0.5 * lp.lam * (1. + sqrt(cos(lp.phi))); + xy.y = lp.phi / (cos(0.5 * lp.phi) * cos(SIXTH * lp.lam)); + 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(larr) { + + P->es = 0; + P->fwd = s_forward; + + return P; +} + + +#ifdef PJ_OMIT_SELFTEST +int pj_larr_selftest (void) {return 0;} +#else + +int pj_larr_selftest (void) { + double tolerance_lp = 1e-10; + double tolerance_xy = 1e-7; + + char s_args[] = {"+proj=larr +a=6400000 +lat_1=0.5 +lat_2=2"}; + + LP fwd_in[] = { + { 2, 1}, + { 2,-1}, + {-2, 1}, + {-2,-1} + }; + + XY s_fwd_expect[] = { + {223393.637624200899, 111707.215961255497}, + {223393.637624200899, -111707.215961255497}, + {-223393.637624200899, 111707.215961255497}, + {-223393.637624200899, -111707.215961255497}, + }; + + return pj_generic_selftest (0, s_args, tolerance_xy, tolerance_lp, 4, 4, fwd_in, 0, s_fwd_expect, 0, 0, 0); +} + + +#endif |
