New

git-svn-id: http://svn.osgeo.org/metacrs/proj/trunk@776 4e78687f-474d-0410-85f9-8d5e500ac6b2

1 files changed, 74 insertions, 0 deletions

diff --git a/src/gen_cheb.c b/src/gen_cheb.c
new file mode 100644
index 00000000..14ac3ce9
--- /dev/null
+++ b/src/gen_cheb.c
@@ -0,0 +1,74 @@
+/* generates 'T' option output */
+#ifndef lint
+static const char SCCSID[]="@(#)gen_cheb.c	4.9	95/09/23	GIE	REL";
+#endif
+#define PJ_LIB__
+#include <stdio.h>
+#include <string.h>
+#include <errno.h>
+#include <projects.h>
+#include "emess.h"
+#ifndef COEF_LINE_MAX
+#define COEF_LINE_MAX 60
+#endif
+	void
+gen_cheb(int inverse, UV (*proj)(), char *s, PJ *P, int iargc, char **iargv) {
+	int NU = 15, NV = 15, i, res = -1, errin = 0, pwr;
+	char *arg, fmt[15];
+	UV low, upp, resid;
+	Tseries *F;
+	extern void p_series(Tseries *, FILE *, char *);
+	double (*input)();
+
+	input = inverse ? strtod : dmstor;
+	if (*s) low.u = input(s, &s); else ++errin;
+	if (*s == ',') upp.u = input(s+1, &s); else ++errin;
+	if (*s == ',') low.v = input(s+1, &s); else ++errin;
+	if (*s == ',') upp.v = input(s+1, &s); else ++errin;
+	if (errin)
+		emess(16,"null or absent -T parameters");
+	if (*s == ',') if (*++s != ',') res = strtol(s, &s, 10);
+	if (*s == ',') if (*++s != ',') NU = strtol(s, &s, 10);
+	if (*s == ',') if (*++s != ',') NV = strtol(s, &s, 10);
+	pwr = s && *s && !strcmp(s, ",P");
+	(void)printf("#proj_%s\n#    run-line:\n",
+		pwr ? "Power" : "Chebyshev");
+	if (iargc > 0) { /* proj execution audit trail */
+		int n = 0, L;
+
+		for( i = 0 ; iargc ; --iargc) {
+			arg = *iargv++;
+			if (*arg != '+') {
+				if (!n) { putchar('#'); ++n; }
+				(void)printf(" %s%n",arg, &L);
+				if ((n += L) > 50) { putchar('\n'); n = 0; }
+			}
+		}
+		if (n) putchar('\n');
+	}
+	(void)printf("# projection parameters\n");
+	pj_pr_list(P);
+	if (low.u == upp.u || low.v >= upp.v)
+		emess(16,"approx. argument range error");
+	if (low.u > upp.u)
+		low.u -= TWOPI;
+	if (NU < 2 || NV < 2)
+		emess(16,"approx. work dimensions (%d %d) too small",NU,NV);
+	if (!(F = mk_cheby(low, upp, pow(10., (double)res)*.5, &resid, proj,
+		NU, NV, pwr)))
+		emess(16,"generation of approx failed\nreason: %s\n",
+			pj_strerrno(errno));
+	(void)printf("%c,%.12g,%.12g,%.12g,%.12g,%.12g\n",inverse?'I':'F',
+		P->lam0*RAD_TO_DEG,
+		low.u*(inverse?1.:RAD_TO_DEG),upp.u*(inverse?1.:RAD_TO_DEG),
+		low.v*(inverse?1.:RAD_TO_DEG),upp.v*(inverse?1.:RAD_TO_DEG));
+	if (pwr)
+		strcpy(fmt, "%.15g");
+	else if (res <= 0)
+		(void)sprintf(fmt,"%%.%df",-res+1);
+	else
+		(void)strcpy(fmt,"%.0f");
+	p_series(F, stdout, fmt);
+	(void)printf("# |u,v| sums %g %g\n#end_proj_%s\n",
+		resid.u, resid.v, pwr ? "Power" : "Chebyshev");
+}