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
|
/* print row coefficients of Tseries structure */
#include "proj.h"
#include "proj_internal.h"
#include <stdio.h>
#include <string.h>
#define NF 20 /* length of final format string */
#define CUT 60 /* check length of line */
/* FIXME: put the declaration in a header. Also used in gen_cheb.c */
void p_series(Tseries *T, FILE *file, char *fmt);
void p_series(Tseries *T, FILE *file, char *fmt) {
int i, j, n, L;
char format[NF+1];
*format = ' ';
strncpy(format + 1, fmt, NF - 3);
strcat(format, "%n");
fprintf(file, "u: %d\n", T->mu+1);
for (i = 0; i <= T->mu; ++i)
if (T->cu[i].m) {
fprintf(file, "%d %d%n", i, T->cu[i].m, &L);
n = 0;
for (j = 0; j < T->cu[i].m; ++j) {
if ((L += n) > CUT)
fprintf(file, "\n %n", &L);
fprintf(file, format, T->cu[i].c[j], &n);
}
fputc('\n', file);
}
fprintf(file, "v: %d\n", T->mv+1);
for (i = 0; i <= T->mv; ++i)
if (T->cv[i].m) {
fprintf(file, "%d %d%n", i, T->cv[i].m, &L);
n = 0;
for (j = 0; j < T->cv[i].m; ++j) {
if ((L += n) > 60)
fprintf(file, "\n %n", &L);
fprintf(file, format, T->cv[i].c[j], &n);
}
fputc('\n', file);
}
}
|
