Skip to content
This repository
Fetching contributors…

Cannot retrieve contributors at this time

file 71 lines (69 sloc) 2.291 kb
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71
/* generates 'T' option output */
#define PJ_LIB__
#include "projects.h"
#include <stdio.h>
#include <string.h>
#include <errno.h>
#include "emess.h"
#ifndef COEF_LINE_MAX
#define COEF_LINE_MAX 60
#endif
void
gen_cheb(int inverse, projUV (*proj)(projUV), char *s, PJ *P, int iargc, char **iargv) {
int NU = 15, NV = 15, i, res = -1, errin = 0, pwr;
char *arg, fmt[15];
projUV low, upp, resid;
Tseries *F;
extern void p_series(Tseries *, FILE *, char *);
double (*input)(const char *, char **);

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");
}
Something went wrong with that request. Please try again.