/Proj4/p_series.c
C | 42 lines | 40 code | 1 blank | 1 comment | 8 complexity | 860ad58f1d23bbb53c73ed71572c4dcf MD5 | raw file
- /* print row coefficients of Tseries structure */
- #ifndef lint
- static const char SCCSID[]="@(#)p_series.c 4.6 95/08/19 GIE REL";
- #endif
- #include "projects.h"
- #include <stdio.h>
- #include <string.h>
- #define NF 20 /* length of final format string */
- #define CUT 60 /* check length of line */
- 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) {
- L = fprintf(file, "%d %d", i, T->cu[i].m);
- n = 0;
- for (j = 0; j < T->cu[i].m; ++j) {
- if ((L += n) > CUT)
- L = fprintf(file, "\n ");
- 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) {
- L = fprintf(file, "%d %d", i, T->cv[i].m);
- n = 0;
- for (j = 0; j < T->cv[i].m; ++j) {
- if ((L += n) > 60)
- L = fprintf(file, "\n ");
- fprintf(file, format, T->cv[i].c[j], &n);
- }
- fputc('\n', file);
- }
- }