From 565a4bd035b9d4a83955808efef20f1d8dfa24cf Mon Sep 17 00:00:00 2001
From: Frank Warmerdam <warmerdam@pobox.com>
Date: Thu, 18 Mar 1999 16:34:52 +0000
Subject: New

git-svn-id: http://svn.osgeo.org/metacrs/proj/trunk@776 4e78687f-474d-0410-85f9-8d5e500ac6b2
---
 src/PJ_hammer.c | 32 ++++++++++++++++++++++++++++++++
 1 file changed, 32 insertions(+)
 create mode 100644 src/PJ_hammer.c

(limited to 'src/PJ_hammer.c')

diff --git a/src/PJ_hammer.c b/src/PJ_hammer.c
new file mode 100644
index 00000000..813d6e5f
--- /dev/null
+++ b/src/PJ_hammer.c
@@ -0,0 +1,32 @@
+#ifndef lint
+static const char SCCSID[]="@(#)PJ_hammer.c	4.1	94/02/15	GIE	REL";
+#endif
+#define PROJ_PARMS__ \
+	double w; \
+	double m, rm;
+#define PJ_LIB__
+# include	<projects.h>
+PROJ_HEAD(hammer, "Hammer & Eckert-Greifendorff")
+	"\n\tMisc Sph, no inv.\n\tW= M=";
+FORWARD(s_forward); /* spheroid */
+	double cosphi, d;
+
+	d = sqrt(2./(1. + (cosphi = cos(lp.phi)) * cos(lp.lam *= P->w)));
+	xy.x = P->m * d * cosphi * sin(lp.lam);
+	xy.y = P->rm * d * sin(lp.phi);
+	return (xy);
+}
+FREEUP; if (P) pj_dalloc(P); }
+ENTRY0(hammer)
+	if (pj_param(P->params, "tW").i) {
+		if ((P->w = fabs(pj_param(P->params, "dW").f)) <= 0.) E_ERROR(-27);
+	} else
+		P->w = .5;
+	if (pj_param(P->params, "tM").i) {
+		if ((P->m = fabs(pj_param(P->params, "dM").f)) <= 0.) E_ERROR(-27);
+	} else
+		P->m = 1.;
+	P->rm = 1. / P->m;
+	P->m /= P->w;
+	P->es = 0.; P->fwd = s_forward;
+ENDENTRY(P)
-- 
cgit v1.2.3