aboutsummaryrefslogtreecommitdiff
path: root/src/PJ_larr.c
diff options
context:
space:
mode:
authorThomas Knudsen <knudsen.thomas@gmail.com>2016-05-20 00:44:56 +0200
committerThomas Knudsen <knudsen.thomas@gmail.com>2016-05-20 00:44:56 +0200
commit8b4c079485a68593b810ebfc14848ef6bc7864fa (patch)
tree5f8bf20027529b6365caba608d5372e19fce4fb3 /src/PJ_larr.c
parent1384b978373cead004211d8548518373d1f4950c (diff)
parent605dc6f3b47bf72c61712b275284212c66e541d0 (diff)
downloadPROJ-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.c77
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