--- ray/src/cal/rcalc.c 2003/02/22 02:07:20 1.1 +++ ray/src/cal/rcalc.c 2003/08/01 14:14:23 1.8 @@ -1,5 +1,5 @@ #ifndef lint -static const char RCSid[] = "$Id: rcalc.c,v 1.1 2003/02/22 02:07:20 greg Exp $"; +static const char RCSid[] = "$Id: rcalc.c,v 1.8 2003/08/01 14:14:23 schorsch Exp $"; #endif /* * rcalc.c - record calculator program. @@ -7,20 +7,16 @@ static const char RCSid[] = "$Id: rcalc.c,v 1.1 2003/0 * 9/11/87 */ -#include - #include - +#include +#include #include - #include +#include "platform.h" #include "calcomp.h" +#include "rterror.h" -#ifdef CPM -#define getc agetc /* text files only, right? */ -#endif - #define isnum(c) (isdigit(c) || (c)=='-' || (c)=='.' \ || (c)=='+' || (c)=='e' || (c)=='E') @@ -58,8 +54,18 @@ struct field { /* record format struc #define savqstr(s) strcpy(emalloc(strlen(s)+1),s) #define freqstr(s) efree(s) -extern char *strcpy(), *emalloc(), *savestr(); -struct strvar *getsvar(); +static void scaninp(void), advinp(void), resetinp(void); +static void putrec(void), putout(void), nbsynch(void); +static int getrec(void); +static void execute(char *file); +static void initinp(FILE *fp); +static void svpreset(char *eqn); +static void readfmt(char *spec, int output); +static int readfield(char **pp); +static int getfield(struct field *f); +static void chanset(int n, double v); +static void bchanset(int n, double v); +static struct strvar* getsvar(char *svname); struct field *inpfmt = NULL; /* input record format */ struct field *outfmt = NULL; /* output record structure */ @@ -69,6 +75,8 @@ int blnkeq = 1; /* blanks compare equa int igneol = 0; /* ignore end of line? */ char sepchar = '\t'; /* input/output separator */ int noinput = 0; /* no input records? */ +int nbicols = 0; /* number of binary input columns */ +int bocols = 0; /* produce binary output columns */ char inpbuf[INBSIZ]; /* input buffer */ double colval[MAXCOL]; /* input column values */ unsigned long colflg = 0; /* column retrieved flags */ @@ -86,13 +94,16 @@ struct { } ipb; /* circular lookahead buffer */ -main(argc, argv) -int argc; -char *argv[]; +int +main( +int argc, +char *argv[] +) { int i; - esupport |= (E_VARIABLE|E_FUNCTION|E_INCHAN|E_OUTCHAN|E_RCONST); + esupport |= E_VARIABLE|E_FUNCTION|E_INCHAN|E_OUTCHAN|E_RCONST; + esupport &= ~(E_REDEFW); #ifdef BIGGERLIB biggerlib(); @@ -123,10 +134,46 @@ char *argv[]; noinput = 1; break; case 'i': - readfmt(argv[++i], 0); + switch (argv[i][2]) { + case '\0': + nbicols = 0; + readfmt(argv[++i], 0); + break; + case 'a': + nbicols = 0; + break; + case 'd': + if (isdigit(argv[i][3])) + nbicols = atoi(argv[i]+3); + else + nbicols = 1; + break; + case 'f': + if (isdigit(argv[i][3])) + nbicols = -atoi(argv[i]+3); + else + nbicols = -1; + break; + default: + goto userr; + } break; case 'o': - readfmt(argv[++i], 1); + switch (argv[i][2]) { + case '\0': + bocols = 0; + readfmt(argv[++i], 1); + break; + case 'a': + bocols = 0; + break; + case 'd': + bocols = 1; + break; + case 'f': + bocols = -1; + break; + } break; case 'w': nowarn = !nowarn; @@ -134,7 +181,8 @@ char *argv[]; case 'u': unbuff = !unbuff; break; - default: + default:; + userr: eputs("Usage: "); eputs(argv[0]); eputs(" [-b][-l][-n][-w][-u][-tS][-s svar=sval][-e expr][-f source][-i infmt][-o outfmt] [file]\n"); @@ -160,7 +208,8 @@ eputs(" [-b][-l][-n][-w][-u][-tS][-s svar=sval][-e exp } -nbsynch() /* non-blank starting synch character */ +static void +nbsynch(void) /* non-blank starting synch character */ { if (inpfmt == NULL || (inpfmt->type & F_TYP) != T_LIT) return; @@ -171,9 +220,28 @@ nbsynch() /* non-blank starting synch ch } -execute(file) /* process a file */ -char *file; +int +getinputrec( /* get next input record */ +FILE *fp +) { + if (inpfmt != NULL) + return(getrec()); + if (nbicols > 0) + return(fread(inpbuf, sizeof(double), + nbicols, fp) == nbicols); + if (nbicols < 0) + return(fread(inpbuf, sizeof(float), + -nbicols, fp) == -nbicols); + return(fgets(inpbuf, INBSIZ, fp) != NULL); +} + + +static void +execute( /* process a file */ +char *file +) +{ int conditional = vardefined("cond"); long nrecs = 0; long nout = 0; @@ -188,8 +256,8 @@ char *file; } if (inpfmt != NULL) initinp(fp); - - while (inpfmt != NULL ? getrec() : fgets(inpbuf, INBSIZ, fp) != NULL) { + + while (getinputrec(fp)) { varset("recno", '=', (double)++nrecs); colflg = 0; eclock++; @@ -202,16 +270,18 @@ char *file; } -putout() /* produce an output record */ +static void +putout(void) /* produce an output record */ { - extern int chanset(); colpos = 0; if (outfmt != NULL) putrec(); + else if (bocols) + chanout(bchanset); else chanout(chanset); - if (colpos) + if (colpos && !bocols) putchar('\n'); if (unbuff) fflush(stdout); @@ -219,8 +289,9 @@ putout() /* produce an output record */ double -chanvalue(n) /* return value for column n */ -int n; +chanvalue( /* return value for column n */ +int n +) { int i; register char *cp; @@ -233,6 +304,18 @@ int n; eputs("illegal channel number\n"); quit(1); } + if (nbicols > 0) { + if (n > nbicols) + return(0.0); + cp = inpbuf + (n-1)*sizeof(double); + return(*(double *)cp); + } + if (nbicols < 0) { + if (n > -nbicols) + return(0.0); + cp = inpbuf + (n-1)*sizeof(float); + return(*(float *)cp); + } if (n <= MAXCOL && colflg & 1L<<(n-1)) return(colval[n-1]); @@ -258,9 +341,11 @@ int n; } -chanset(n, v) /* output column n */ -int n; -double v; +void +chanset( /* output column n */ +int n, +double v +) { if (colpos == 0) /* no leading separator */ colpos = 1; @@ -272,10 +357,33 @@ double v; } -readfmt(spec, output) /* read record format */ -char *spec; -int output; +void +bchanset( /* output binary channel n */ +int n, +double v +) { + static char zerobuf[sizeof(double)]; + + while (++colpos < n) + fwrite(zerobuf, + bocols>0 ? sizeof(double) : sizeof(float), + 1, stdout); + if (bocols > 0) + fwrite(&v, sizeof(double), 1, stdout); + else { + float fval = v; + fwrite(&fval, sizeof(float), 1, stdout); + } +} + + +static void +readfmt( /* read record format */ +char *spec, +int output +) +{ int fd; char *inptr; struct field fmt; @@ -338,9 +446,10 @@ int output; } -int -readfield(pp) /* get next field in format */ -register char **pp; +static int +readfield( /* get next field in format */ +register char **pp +) { int type = F_NUL; int width = 0; @@ -405,8 +514,9 @@ register char **pp; struct strvar * -getsvar(svname) /* get string variable */ -char *svname; +getsvar( /* get string variable */ +char *svname +) { register struct strvar *sv; @@ -422,8 +532,10 @@ char *svname; } -svpreset(eqn) /* preset a string variable */ -char *eqn; +static void +svpreset( /* preset a string variable */ +char *eqn +) { register struct strvar *sv; register char *val; @@ -442,7 +554,8 @@ char *eqn; } -clearrec() /* clear input record variables */ +static void +clearrec(void) /* clear input record variables */ { register struct field *f; @@ -461,7 +574,8 @@ clearrec() /* clear input record variables */ } -getrec() /* get next record from file */ +static int +getrec(void) /* get next record from file */ { int eatline; register struct field *f; @@ -494,10 +608,12 @@ getrec() /* get next re } -getfield(f) /* get next field */ -register struct field *f; +static int +getfield( /* get next field */ +register struct field *f +) { - static char buf[MAXWORD+1]; /* no recursion! */ + static char buf[RMAXWORD+1]; /* no recursion! */ int delim, inword; double d; char *np; @@ -542,7 +658,7 @@ register struct field *f; *cp++ = ipb.chr; scaninp(); } - } while (inword && cp < &buf[MAXWORD]); + } while (inword && cp < &buf[RMAXWORD]); *cp = '\0'; if (f->f.sv->val == NULL) f->f.sv->val = savqstr(buf); /* first setting */ @@ -571,7 +687,7 @@ register struct field *f; *cp++ = ipb.chr; scaninp(); } - } while (inword && cp < &buf[MAXWORD]); + } while (inword && cp < &buf[RMAXWORD]); *cp = '\0'; d = np==NULL ? 0. : atof(np); if (!vardefined(f->f.nv)) @@ -581,10 +697,12 @@ register struct field *f; return(-1); /* doesn't match! */ return(0); } + return -1; /* pro forma return */ } -putrec() /* output a record */ +static void +putrec(void) /* output a record */ { char fmt[32]; register int n; @@ -634,8 +752,9 @@ putrec() /* output a re } -initinp(fp) /* prepare lookahead buffer */ -FILE *fp; +static void +initinp(FILE *fp) /* prepare lookahead buffer */ + { ipb.fin = fp; ipb.beg = ipb.end = inpbuf; @@ -645,7 +764,8 @@ FILE *fp; } -scaninp() /* scan next character */ +static void +scaninp(void) /* scan next character */ { if (ipb.chr == EOF) return; @@ -664,13 +784,15 @@ scaninp() /* scan next character } -advinp() /* move home to current position */ +static void +advinp(void) /* move home to current position */ { ipb.beg = ipb.pos; } -resetinp() /* rewind position and advance 1 */ +static void +resetinp(void) /* rewind position and advance 1 */ { if (ipb.beg == NULL) /* full */ ipb.beg = ipb.end; @@ -683,16 +805,14 @@ resetinp() /* rewind position and void -eputs(msg) -char *msg; +eputs(char *msg) { fputs(msg, stderr); } void -wputs(msg) -char *msg; +wputs(char *msg) { if (!nowarn) eputs(msg); @@ -700,8 +820,7 @@ char *msg; void -quit(code) -int code; +quit(int code) { exit(code); }