ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/cal/rcalc.c
(Generate patch)

Comparing ray/src/cal/rcalc.c (file contents):
Revision 1.13 by greg, Wed Jun 9 16:00:13 2004 UTC vs.
Revision 1.32 by greg, Fri Mar 11 22:50:13 2022 UTC

# Line 2 | Line 2
2   static const char RCSid[] = "$Id$";
3   #endif
4   /*
5 < *  rcalc.c - record calculator program.
5 > * rcalc.c - record calculator program.
6   *
7 < *     9/11/87
7 > * 9/11/87
8   */
9  
10 < #include  <stdlib.h>
11 < #include  <fcntl.h>
12 < #include  <stdio.h>
13 < #include  <string.h>
14 < #include  <math.h>
15 < #include  <ctype.h>
10 > #include <stdlib.h>
11 > #include <math.h>
12 > #include <ctype.h>
13  
14 < #include  "platform.h"
15 < #include  "rterror.h"
16 < #include  "rtmisc.h"
17 < #include  "rtio.h"
18 < #include  "calcomp.h"
14 > #include "platform.h"
15 > #include "rterror.h"
16 > #include "rtmisc.h"
17 > #include "rtio.h"
18 > #include "calcomp.h"
19  
20 < #define  isnum(c)       (isdigit(c) || (c)=='-' || (c)=='.' \
21 <                                || (c)=='+' || (c)=='e' || (c)=='E')
20 > #define isnum(c) (isdigit(c) || ((c)=='-') | ((c)=='.') \
21 >                                | ((c)=='+') | ((c)=='e') | ((c)=='E'))
22  
23 < #define  isblnk(c)      (igneol ? isspace(c) : (c)==' '||(c)=='\t')
23 > #define isblnk(c) (igneol ? isspace(c) : ((c)==' ')|((c)=='\t'))
24  
25 < #define  INBSIZ         4096    /* longest record */
26 < #define  MAXCOL         32      /* number of columns recorded */
25 > #define INBSIZ   16384  /* longest record */
26 > #define MAXCOL   32 /* number of columns recorded */
27  
28                                  /* field type specifications */
29 < #define  F_NUL          0               /* empty */
30 < #define  F_TYP          0x7000          /* mask for type */
31 < #define  F_WID          0x0fff          /* mask for width */
32 < #define  T_LIT          0x1000          /* string literal */
33 < #define  T_STR          0x2000          /* string variable */
34 < #define  T_NUM          0x3000          /* numeric value */
29 > #define F_NUL    0       /* empty */
30 > #define F_TYP    0x7000  /* mask for type */
31 > #define F_WID    0x0fff  /* mask for width */
32 > #define T_LIT    0x1000  /* string literal */
33 > #define T_STR    0x2000  /* string variable */
34 > #define T_NUM    0x3000  /* numeric value */
35  
36 < struct strvar {                 /* string variable */
37 <        char  *name;
38 <        char  *val;
39 <        char  *preset;
40 <        struct strvar  *next;
36 > struct strvar {          /* string variable */
37 >        char *name;
38 >        char *val;
39 >        char *preset;
40 >        struct strvar *next;
41   };
42  
43 < struct field {                  /* record format structure */
44 <        int  type;                      /* type of field (& width) */
43 > struct field {           /* record format structure */
44 >        int type;                /* type of field (& width) */
45          union {
46 <                char  *sl;                      /* string literal */
47 <                struct strvar  *sv;             /* string variable */
48 <                char  *nv;                      /* numeric variable */
49 <                EPNODE  *ne;                    /* numeric expression */
50 <        } f;                            /* field contents */
51 <        struct field  *next;            /* next field in record */
46 >                char *sl;                /* string literal */
47 >                struct strvar *sv;       /* string variable */
48 >                char *nv;                /* numeric variable */
49 >                EPNODE *ne;              /* numeric expression */
50 >        } f;                     /* field contents */
51 >        struct field *next;      /* next field in record */
52   };
53  
54 < #define  savqstr(s)     strcpy(emalloc(strlen(s)+1),s)
55 < #define  freqstr(s)     efree(s)
54 > #define savqstr(s) strcpy(emalloc(strlen(s)+1),s)
55 > #define freqstr(s) efree(s)
56  
57   static int getinputrec(FILE *fp);
58 < static void scaninp(void), advinp(void), resetinp(void);
58 > static void scaninp(void), advinp(void), passinp(void), skipinp(void);
59   static void putrec(void), putout(void), nbsynch(void);
60   static int getrec(void);
61   static void execute(char *file);
# Line 72 | Line 69 | static void bchanset(int n, double v);
69   static struct strvar* getsvar(char *svname);
70   static double l_in(char *);
71  
72 < struct field  *inpfmt = NULL;   /* input record format */
73 < struct field  *outfmt = NULL;   /* output record structure */
74 < struct strvar  *svhead = NULL;  /* string variables */
72 > struct field *inpfmt = NULL;    /* input record format */
73 > struct field *outfmt = NULL;    /* output record structure */
74 > struct strvar *svhead = NULL;   /* string variables */
75  
76 < int  blnkeq = 1;                /* blanks compare equal? */
77 < int  igneol = 0;                /* ignore end of line? */
81 < int  passive = 0;               /* passive mode (transmit unmatched input) */
82 < char  sepchar = '\t';           /* input/output separator */
83 < int  noinput = 0;               /* no input records? */
84 < int  nbicols = 0;               /* number of binary input columns */
85 < int  bocols = 0;                /* produce binary output columns */
86 < char  inpbuf[INBSIZ];           /* input buffer */
87 < double  colval[MAXCOL];         /* input column values */
88 < unsigned long  colflg = 0;      /* column retrieved flags */
89 < int  colpos;                    /* output column position */
76 > long incnt = 0;                 /* limit number of input records? */
77 > long outcnt = 0;                /* limit number of output records? */
78  
79 < int  nowarn = 0;                /* non-fatal diagnostic output */
80 < int  unbuff = 0;                /* unbuffered output (flush each record) */
79 > int blnkeq = 1;                 /* blanks compare equal? */
80 > int igneol = 0;                 /* ignore end of line? */
81 > int passive = 0;                /* passive mode (transmit unmatched input) */
82 > char sepchar = '\t';            /* input/output separator */
83 > int noinput = 0;                /* no input records? */
84 > int itype = 'a';                /* input type (a/f/F/d/D) */
85 > int nbicols = 0;                /* number of binary input columns */
86 > int otype = 'a';                /* output format (a/f/F/d/D) */
87 > char inpbuf[INBSIZ];            /* input buffer */
88 > double colval[MAXCOL];          /* input column values */
89 > unsigned long colflg = 0;       /* column retrieved flags */
90 > int colpos;                     /* output column position */
91  
92 + int nowarn = 0;                 /* non-fatal diagnostic output */
93 + int unbuff = 0;                 /* unbuffered output (flush each record) */
94 +
95   struct {
96 <        FILE  *fin;                     /* input file */
97 <        int  chr;                       /* next character */
98 <        char  *beg;                     /* home position */
99 <        char  *pos;                     /* scan position */
100 <        char  *end;                     /* read position */
101 < } ipb;                          /* circular lookahead buffer */
96 >        FILE *fin;               /* input file */
97 >        int chr;                 /* next character */
98 >        char *beg;               /* home position */
99 >        char *pos;               /* scan position */
100 >        char *end;               /* read position */
101 > } ipb;                   /* circular lookahead buffer */
102  
103  
104   int
105   main(
106 < int  argc,
107 < char  *argv[]
106 > int argc,
107 > char *argv[]
108   )
109   {
110 <        int  i;
110 >        char *fpath;
111 >        int i;
112  
113          esupport |= E_VARIABLE|E_FUNCTION|E_INCHAN|E_OUTCHAN|E_RCONST;
114          esupport &= ~(E_REDEFW);
115  
116 < #ifdef  BIGGERLIB
116 > #ifdef BIGGERLIB
117          biggerlib();
118   #endif
119          varset("PI", ':', 3.14159265358979323846);
# Line 126 | Line 128 | char  *argv[]
128                          igneol = !igneol;
129                          break;
130                  case 'p':
131 <                        passive = !passive;
131 >                        passive = 1;
132                          break;
133 +                case 'P':
134 +                        passive = -1;
135 +                        break;
136                  case 't':
137                          sepchar = argv[i][2];
138                          break;
# Line 135 | Line 140 | char  *argv[]
140                          svpreset(argv[++i]);
141                          break;
142                  case 'f':
143 <                        fcompile(argv[++i]);
143 >                        fpath = getpath(argv[++i], getrlibpath(), 0);
144 >                        if (fpath == NULL) {
145 >                                eputs(argv[0]);
146 >                                eputs(": cannot find file '");
147 >                                eputs(argv[i]);
148 >                                eputs("'\n");
149 >                                quit(1);
150 >                        }
151 >                        fcompile(fpath);
152                          break;
153                  case 'e':
154                          scompile(argv[++i], NULL, 0);
# Line 146 | Line 159 | char  *argv[]
159                  case 'i':
160                          switch (argv[i][2]) {
161                          case '\0':
162 +                                itype = 'a';
163                                  nbicols = 0;
164                                  readfmt(argv[++i], 0);
165                                  break;
166 +                        case 'n':
167 +                                incnt = atol(argv[++i]);
168 +                                break;
169                          case 'a':
170 +                                itype = 'a';
171                                  nbicols = 0;
172                                  break;
173                          case 'd':
174 +                        case 'D':
175 +                                itype = argv[i][2];
176                                  if (isdigit(argv[i][3]))
177                                          nbicols = atoi(argv[i]+3);
178                                  else
# Line 164 | Line 184 | char  *argv[]
184                                  }
185                                  break;
186                          case 'f':
187 +                        case 'F':
188 +                                itype = argv[i][2];
189                                  if (isdigit(argv[i][3]))
190 <                                        nbicols = -atoi(argv[i]+3);
190 >                                        nbicols = atoi(argv[i]+3);
191                                  else
192 <                                        nbicols = -1;
193 <                                if (-nbicols*sizeof(float) > INBSIZ) {
192 >                                        nbicols = 1;
193 >                                if (nbicols*sizeof(float) > INBSIZ) {
194                                          eputs(argv[0]);
195                                          eputs(": too many input columns\n");
196                                          quit(1);
# Line 181 | Line 203 | char  *argv[]
203                  case 'o':
204                          switch (argv[i][2]) {
205                          case '\0':
206 <                                bocols = 0;
206 >                                otype = 'a';
207                                  readfmt(argv[++i], 1);
208                                  break;
209 +                        case 'n':
210 +                                outcnt = atol(argv[++i]);
211 +                                break;
212                          case 'a':
213 <                                bocols = 0;
213 >                                otype = 'a';
214                                  break;
215                          case 'd':
216 <                                bocols = 1;
192 <                                break;
216 >                        case 'D':
217                          case 'f':
218 <                                bocols = -1;
218 >                        case 'F':
219 >                                otype = argv[i][2];
220                                  break;
221 +                        default:
222 +                                goto userr;
223                          }
224                          break;
225                  case 'w':
# Line 205 | Line 232 | char  *argv[]
232                  userr:
233                          eputs("Usage: ");
234                          eputs(argv[0]);
235 < eputs(" [-b][-l][-n][-p][-w][-u][-tS][-s svar=sval][-e expr][-f source][-i infmt][-o outfmt] [file]\n");
235 > eputs(" [-b][-l][-n][-p|-P][-w][-u][-tS][-s svar=sval][-e expr][-f source][-i infmt][-o outfmt] [file]\n");
236                          quit(1);
237                  }
238 <
239 <        if (noinput) {          /* produce a single output record */
238 >        if (otype != 'a')
239 >                SET_FILE_BINARY(stdout);
240 > #ifdef getc_unlocked            /* avoid lock/unlock overhead */
241 >        flockfile(stdout);
242 > #endif
243 >        if (noinput) {   /* produce a single output record */
244 >                if (i < argc) {
245 >                        eputs(argv[0]);
246 >                        eputs(": file argument(s) incompatible with -n\n");
247 >                        quit(1);
248 >                }
249                  eclock++;
250                  putout();
251                  quit(0);
252          }
253 <
254 <        if (blnkeq)             /* for efficiency */
253 >        if (passive && (inpfmt == NULL) | (outfmt == NULL)) {
254 >                eputs(argv[0]);
255 >                eputs(": options -p and -P require -i and -o formats\n");
256 >                quit(1);
257 >        }
258 >        if (blnkeq)      /* for efficiency */
259                  nbsynch();
260  
261 <        if (i == argc)          /* from stdin */
261 >        if (i == argc)   /* from stdin */
262                  execute(NULL);
263 <        else                    /* from one or more files */
263 >        else             /* from one or more files */
264                  for ( ; i < argc; i++)
265                          execute(argv[i]);
266          
# Line 230 | Line 270 | eputs(" [-b][-l][-n][-p][-w][-u][-tS][-s svar=sval][-e
270  
271  
272   static void
273 < nbsynch(void)               /* non-blank starting synch character */
273 > nbsynch(void)    /* non-blank starting synch character */
274   {
275          if (inpfmt == NULL || (inpfmt->type & F_TYP) != T_LIT)
276                  return;
# Line 243 | Line 283 | nbsynch(void)               /* non-blank starting sync
283  
284   static int
285   getinputrec(            /* get next input record */
286 < FILE  *fp
286 > FILE *fp
287   )
288   {
289 <        if (inpfmt != NULL)
290 <                return(getrec());
291 <        if (nbicols > 0)
292 <                return(fread(inpbuf, sizeof(double),
293 <                                        nbicols, fp) == nbicols);
294 <        if (nbicols < 0)
295 <                return(fread(inpbuf, sizeof(float),
296 <                                        -nbicols, fp) == -nbicols);
289 >        if ((itype == 'd') | (itype == 'D')) {
290 >                if (getbinary(inpbuf, sizeof(double), nbicols, fp) != nbicols)
291 >                        return(0);
292 >                if (itype == 'D')
293 >                        swap64(inpbuf, nbicols);
294 >                return(1);
295 >        }
296 >        if ((itype == 'f') | (itype == 'F')) {
297 >                if (getbinary(inpbuf, sizeof(float), nbicols, fp) != nbicols)
298 >                        return(0);
299 >                if (itype == 'F')
300 >                        swap32(inpbuf, nbicols);
301 >                return(1);
302 >        }
303          return(fgets(inpbuf, INBSIZ, fp) != NULL);
304   }
305  
306  
307   static void
308 < execute(           /* process a file */
309 < char  *file
308 > execute(         /* process a file */
309 > char *file
310   )
311   {
312 <        int  conditional = vardefined("cond");
313 <        long  nrecs = 0;
314 <        long  nout = 0;
315 <        FILE  *fp;
312 >        const int conditional = vardefined("cond");
313 >        const int set_recno = (varlookup("recno") != NULL);
314 >        const int set_outno = (varlookup("outno") != NULL);
315 >        long nrecs = 0;
316 >        long nout = 0;
317 >        FILE *fp;
318          
319          if (file == NULL)
320                  fp = stdin;
# Line 275 | Line 323 | char  *file
323                  eputs(": cannot open\n");
324                  quit(1);
325          }
326 +        if (itype != 'a')
327 +                SET_FILE_BINARY(fp);
328 + #ifdef getc_unlocked            /* avoid lock/unlock overhead */
329 +        flockfile(fp);
330 + #endif
331          if (inpfmt != NULL)
332                  initinp(fp);
333          
334 <        while (getinputrec(fp)) {
335 <                varset("recno", '=', (double)++nrecs);
334 >        while (inpfmt != NULL ? getrec() : getinputrec(fp)) {
335 >                ++nrecs;
336 >                if (set_recno)
337 >                        varset("recno", '=', (double)nrecs);
338 >                if (set_outno)
339 >                        varset("outno", '=', (double)(nout+1));
340                  colflg = 0;
341                  eclock++;
342                  if (!conditional || varvalue("cond") > 0.0) {
286 –                        varset("outno", '=', (double)++nout);
343                          putout();
344 +                        ++nout;
345 +                        advinp();
346 +                } else if (inpfmt != NULL) {
347 +                        if (passive < 0)
348 +                                passinp();
349 +                        else
350 +                                advinp();
351                  }
352 +                if (incnt && nrecs >= incnt)
353 +                        break;
354 +                if (outcnt && nout >= outcnt)
355 +                        break;
356          }
357          fclose(fp);
358   }
359  
360  
361   static void
362 < putout(void)                /* produce an output record */
362 > putout(void)            /* produce an output record */
363   {
364  
365          colpos = 0;
366          if (outfmt != NULL)
367                  putrec();
368 <        else if (bocols)
302 <                chanout(bchanset);
303 <        else
368 >        else if (otype == 'a')
369                  chanout(chanset);
370 <        if (colpos && !bocols)
370 >        else
371 >                chanout(bchanset);
372 >        if (colpos && otype == 'a')
373                  putchar('\n');
374          if (unbuff)
375                  fflush(stdout);
# Line 312 | Line 379 | putout(void)                /* produce an output recor
379   static double
380   l_in(char *funame)      /* function call for $channel */
381   {
382 <        int  n;
383 <        register char  *cp;
382 >        int n;
383 >        char *cp;
384                          /* get argument as integer */
385          n = (int)(argument(1) + .5);
386          if (n != 0)     /* return channel value */
# Line 321 | Line 388 | l_in(char *funame)     /* function call for $channel */
388                          /* determine number of channels */
389          if (noinput || inpfmt != NULL)
390                  return(0);
391 <        if (nbicols > 0)
391 >        if (nbicols)
392                  return(nbicols);
326 –        if (nbicols < 0)
327 –                return(-nbicols);
393          cp = inpbuf;    /* need to count */
394          for (n = 0; *cp; )
395                  if (blnkeq && isspace(sepchar)) {
# Line 342 | Line 407 | l_in(char *funame)     /* function call for $channel */
407   }
408  
409   double
410 < chanvalue(            /* return value for column n */
411 < int  n
410 > chanvalue(       /* return value for column n */
411 > int n
412   )
413   {
414 <        int  i;
415 <        register char  *cp;
414 >        int i;
415 >        char *cp;
416  
417          if (noinput || inpfmt != NULL) {
418                  eputs("no column input\n");
# Line 357 | Line 422 | int  n
422                  eputs("illegal channel number\n");
423                  quit(1);
424          }
425 <        if (nbicols > 0) {
425 >        if (nbicols) {
426                  if (n > nbicols)
427                          return(0.0);
428 <                cp = inpbuf + (n-1)*sizeof(double);
429 <                return(*(double *)cp);
430 <        }
431 <        if (nbicols < 0) {
367 <                if (n > -nbicols)
368 <                        return(0.0);
428 >                if ((itype == 'd') | (itype == 'D')) {
429 >                        cp = inpbuf + (n-1)*sizeof(double);
430 >                        return(*(double *)cp);
431 >                }
432                  cp = inpbuf + (n-1)*sizeof(float);
433                  return(*(float *)cp);
434          }
# Line 383 | Line 446 | int  n
446                          while (*cp && *cp++ != sepchar)
447                                  ;
448  
449 <        while (isspace(*cp))            /* some atof()'s don't like tabs */
449 >        while (isspace(*cp))     /* some atof()'s don't like tabs */
450                  cp++;
451  
452          if (n <= MAXCOL) {
# Line 395 | Line 458 | int  n
458  
459  
460   void
461 < chanset(                   /* output column n */
462 < int  n,
463 < double  v
461 > chanset(                 /* output column n */
462 > int n,
463 > double v
464   )
465   {
466 <        if (colpos == 0)                /* no leading separator */
466 >        if (colpos == 0)                /* no leading separator */
467                  colpos = 1;
468          while (colpos < n) {
469                  putchar(sepchar);
# Line 411 | Line 474 | double  v
474  
475  
476   void
477 < bchanset(                   /* output binary channel n */
478 < int  n,
479 < double  v
477 > bchanset(                /* output binary channel n */
478 > int n,
479 > double v
480   )
481   {
482          static char     zerobuf[sizeof(double)];
483 +        const int       otlen = ((otype == 'd') | (otype == 'D')) ?
484 +                                        sizeof(double) : sizeof(float);
485 +        float   fval = v;
486  
487          while (++colpos < n)
488 <                fwrite(zerobuf,
489 <                        bocols>0 ? sizeof(double) : sizeof(float),
490 <                        1, stdout);
491 <        if (bocols > 0)
492 <                fwrite(&v, sizeof(double), 1, stdout);
493 <        else {
494 <                float   fval = v;
495 <                fwrite(&fval, sizeof(float), 1, stdout);
488 >                putbinary(zerobuf, otlen, 1, stdout);
489 >        switch (otype) {
490 >        case 'D':
491 >                swap64((char *)&v, 1);
492 >                /* fall through */
493 >        case 'd':
494 >                putbinary(&v, sizeof(double), 1, stdout);
495 >                break;
496 >        case 'F':
497 >                swap32((char *)&fval, 1);
498 >                /* fall through */
499 >        case 'f':
500 >                putbinary(&fval, sizeof(float), 1, stdout);
501 >                break;
502          }
503   }
504  
505  
506   static void
507 < readfmt(                   /* read record format */
508 < char  *spec,
509 < int  output
507 > readfmt(                 /* read record format */
508 > char *spec,
509 > int output
510   )
511   {
512 <        int  fd;
513 <        char  *inptr;
514 <        struct field  fmt;
515 <        int  res;
516 <        register struct field  *f;
512 >        int fd;
513 >        char *inptr;
514 >        struct field fmt;
515 >        int res;
516 >        struct field *f;
517                                                  /* check for inline format */
518          for (inptr = spec; *inptr; inptr++)
519                  if (*inptr == '$')
520                          break;
521 <        if (*inptr)                             /* inline */
521 >        if (*inptr)                      /* inline */
522                  inptr = spec;
523 <        else {                                  /* from file */
523 >        else {                           /* from file */
524                  if ((fd = open(spec, 0)) == -1) {
525                          eputs(spec);
526                          eputs(": cannot open\n");
527                          quit(1);
528                  }
529 <                res = read(fd, inpbuf+1, INBSIZ-1);
529 >                res = read(fd, inpbuf+2, INBSIZ-2);
530                  if (res <= 0 || res >= INBSIZ-1) {
531                          eputs(spec);
532                          if (res < 0)
# Line 466 | Line 538 | int  output
538                          quit(1);
539                  }
540                  close(fd);
541 <                (inptr=inpbuf+1)[res] = '\0';
541 >                (inptr=inpbuf+2)[res] = '\0';
542          }
543 <        f = &fmt;                               /* get fields */
543 >        f = &fmt;                        /* get fields */
544          while ((res = readfield(&inptr)) != F_NUL) {
545                  f->next = (struct field *)emalloc(sizeof(struct field));
546                  f = f->next;
# Line 500 | Line 572 | int  output
572  
573  
574   static int
575 < readfield(                   /* get next field in format */
576 < register char  **pp
575 > readfield(               /* get next field in format */
576 > char **pp
577   )
578   {
579 <        int  type = F_NUL;
580 <        int  width = 0;
581 <        register char  *cp;
579 >        int type = F_NUL;
580 >        int width = 0;
581 >        char *cp;
582          
583          cp = inpbuf;
584          while (cp < &inpbuf[INBSIZ-1] && **pp != '\0') {
# Line 567 | Line 639 | register char  **pp
639  
640  
641   struct strvar *
642 < getsvar(                         /* get string variable */
643 < char  *svname
642 > getsvar(                         /* get string variable */
643 > char *svname
644   )
645   {
646 <        register struct strvar  *sv;
646 >        struct strvar *sv;
647          
648          for (sv = svhead; sv != NULL; sv = sv->next)
649                  if (!strcmp(sv->name, svname))
# Line 586 | Line 658 | char  *svname
658  
659  
660   static void
661 < svpreset(                    /* preset a string variable */
662 < char  *eqn
661 > svpreset(                /* preset a string variable */
662 > char *eqn
663   )
664   {
665 <        register struct strvar  *sv;
666 <        register char  *val;
665 >        struct strvar *sv;
666 >        char *val;
667  
668          for (val = eqn; *val != '='; val++)
669                  if (!*val)
# Line 610 | Line 682 | char  *eqn
682   static void
683   clearrec(void)                  /* clear input record variables */
684   {
685 <        register struct field  *f;
685 >        struct field *f;
686  
687          for (f = inpfmt; f != NULL; f = f->next)
688                  switch (f->type & F_TYP) {
# Line 628 | Line 700 | clearrec(void)                 /* clear input record variables */
700  
701  
702   static int
703 < getrec(void)                            /* get next record from file */
703 > getrec(void)                    /* get next record from file */
704   {
705 <        int  eatline;
706 <        register struct field  *f;
707 <        
705 >        int eatline;
706 >        struct field *f;
707 >
708          while (ipb.chr != EOF) {
709 <                eatline = !igneol && ipb.chr != '\n';
638 <                if (blnkeq)             /* beware of nbsynch() */
709 >                if (blnkeq) {           /* beware of nbsynch() */
710                          while (isblnk(ipb.chr))
711 <                                resetinp();
711 >                                skipinp();
712 >                        if (ipb.chr == EOF)
713 >                                return(0);
714 >                }
715 >                eatline = !igneol & (ipb.chr != '\n');
716                  clearrec();             /* start with fresh record */
717                  for (f = inpfmt; f != NULL; f = f->next)
718 <                        if (getfield(f) == -1)
718 >                        if (!getfield(f))
719                                  break;
720 <                if (f == NULL) {
646 <                        advinp();       /* got one! */
720 >                if (f == NULL)          /* got one? */
721                          return(1);
722 <                }
723 <                resetinp();             /* eat false start */
650 <                if (eatline) {          /* eat rest of line */
722 >                skipinp();              /* else eat false start */
723 >                if (eatline) {          /* eat rest of line */
724                          while (ipb.chr != '\n') {
725                                  if (ipb.chr == EOF)
726                                          return(0);
727 <                                resetinp();
727 >                                skipinp();
728                          }
729 <                        resetinp();
729 >                        skipinp();
730                  }
731          }
732          return(0);
# Line 661 | Line 734 | getrec(void)                           /* get next record from file */
734  
735  
736   static int
737 < getfield(                             /* get next field */
738 < register struct field  *f
737 > getfield(                        /* get next field */
738 > struct field *f
739   )
740   {
741 <        static char  buf[RMAXWORD+1];            /* no recursion! */
742 <        int  delim, inword;
743 <        double  d;
744 <        char  *np;
745 <        register char  *cp;
741 >        static char buf[RMAXWORD+1];     /* no recursion! */
742 >        int delim, inword;
743 >        double d;
744 >        char *np;
745 >        char *cp;
746  
747          switch (f->type & F_TYP) {
748          case T_LIT:
# Line 677 | Line 750 | register struct field  *f
750                  do {
751                          if (blnkeq && isblnk(*cp)) {
752                                  if (!isblnk(ipb.chr))
753 <                                        return(-1);
753 >                                        return(0);
754                                  do
755                                          cp++;
756                                  while (isblnk(*cp));
# Line 688 | Line 761 | register struct field  *f
761                                  cp++;
762                                  scaninp();
763                          } else
764 <                                return(-1);
764 >                                return(0);
765                  } while (*cp);
766 <                return(0);
766 >                break;
767          case T_STR:
768                  if (f->next == NULL || (f->next->type & F_TYP) != T_LIT)
769                          delim = EOF;
# Line 698 | Line 771 | register struct field  *f
771                          delim = f->next->f.sl[0];
772                  cp = buf;
773                  do {
774 <                        if (ipb.chr == EOF || ipb.chr == '\n')
774 >                        if ((ipb.chr == EOF) | (ipb.chr == '\n'))
775                                  inword = 0;
776                          else if (blnkeq && delim != EOF)
777                                  inword = isblnk(delim) ?
# Line 715 | Line 788 | register struct field  *f
788                  if (f->f.sv->val == NULL)
789                          f->f.sv->val = savqstr(buf);    /* first setting */
790                  else if (strcmp(f->f.sv->val, buf))
791 <                        return(-1);                     /* doesn't match! */
792 <                return(0);
791 >                        return(0);                      /* doesn't match! */
792 >                break;
793          case T_NUM:
794                  if (f->next == NULL || (f->next->type & F_TYP) != T_LIT)
795                          delim = EOF;
# Line 746 | Line 819 | register struct field  *f
819                          varset(f->f.nv, '=', d);        /* first setting */
820                  else if ((d = (varvalue(f->f.nv)-d)/(d==0.?1.:d)) > .001
821                                  || d < -.001)
822 <                        return(-1);                     /* doesn't match! */
823 <                return(0);
822 >                        return(0);                      /* doesn't match! */
823 >                break;
824          }
825 <        return -1; /* pro forma return */
825 >        return(1);      /* success! */
826   }
827  
828  
829   static void
830 < putrec(void)                                /* output a record */
830 > putrec(void)                            /* output a record */
831   {
832 <        char  fmt[32];
833 <        register int  n;
834 <        register struct field  *f;
835 <        int  adlast, adnext;
832 >        char fmt[32], typ[16];
833 >        int n;
834 >        struct field *f;
835 >        int adlast, adnext;
836 >        double dv, av;
837          
838          adlast = 0;
839          for (f = outfmt; f != NULL; f = f->next) {
840 <                adnext =        blnkeq &&
840 >                adnext =        blnkeq &&
841                                  f->next != NULL &&
842                                  !( (f->next->type&F_TYP) == T_LIT &&
843                                          f->next->f.sl[0] == ' ' );
# Line 790 | Line 864 | putrec(void)                                /* output
864                          break;
865                  case T_NUM:
866                          n = f->type & F_WID;
867 +                        dv = evalue(f->f.ne);
868 +                        av = fabs(dv);
869 +                        if (n <= 9)
870 +                                strcpy(typ, "g");
871 +                        else
872 +                                sprintf(typ, ".%de", n-5);
873 +                        if (av < 1L<<31) {
874 +                                long    iv = (int)(av + .5);
875 +                                if (iv && fabs(av-iv) <= av*1e-14)
876 +                                        strcpy(typ, ".0f");
877 +                        }
878                          if (adlast && adnext)
879 <                                strcpy(fmt, "%g");
879 >                                sprintf(fmt, "%%%s", typ);
880                          else if (adlast)
881 <                                sprintf(fmt, "%%-%dg", n);
881 >                                sprintf(fmt, "%%-%d%s", n, typ);
882                          else
883 <                                sprintf(fmt, "%%%dg", n);
884 <                        printf(fmt, evalue(f->f.ne));
883 >                                sprintf(fmt, "%%%d%s", n, typ);
884 >                        printf(fmt, dv);
885                          adlast = 1;
886                          break;
887                  }
# Line 805 | Line 890 | putrec(void)                                /* output
890  
891  
892   static void
893 < initinp(FILE  *fp)                     /* prepare lookahead buffer */
893 > initinp(FILE *fp)                /* prepare lookahead buffer */
894  
895   {
896          ipb.fin = fp;
897          ipb.beg = ipb.end = inpbuf;
898 <        ipb.pos = inpbuf-1;             /* position before beginning */
898 >        ipb.pos = inpbuf-1;             /* position before beginning */
899          ipb.chr = '\0';
900          scaninp();
901   }
902  
903  
904   static void
905 < scaninp(void)                       /* scan next character */
905 > scaninp(void)                   /* scan next character */
906   {
907          if (ipb.chr == EOF)
908                  return;
909          if (++ipb.pos >= &inpbuf[INBSIZ])
910                  ipb.pos = inpbuf;
911 <        if (ipb.pos == ipb.end) {               /* new character */
911 >        if (ipb.pos == ipb.end) {        /* new character */
912                  if ((ipb.chr = getc(ipb.fin)) != EOF) {
913                          *ipb.end = ipb.chr;
914                          if (++ipb.end >= &inpbuf[INBSIZ])
# Line 837 | Line 922 | scaninp(void)                       /* scan next chara
922  
923  
924   static void
925 < advinp(void)                        /* move home to current position */
925 > advinp(void)                    /* move home to current position */
926   {
927          ipb.beg = ipb.pos;
928   }
929  
930  
931   static void
932 < resetinp(void)                      /* rewind position and advance 1 */
932 > passinp(void)                   /* pass beginning to current position */
933   {
934 <        if (ipb.beg == NULL)            /* full */
934 >        if (!passive | (ipb.beg == NULL)) {
935 >                advinp();
936 >                return;
937 >        }
938 >        while (ipb.beg != ipb.pos) {    /* transfer buffer unaltered */
939 >                putchar(*ipb.beg);
940 >                if (++ipb.beg >= &inpbuf[INBSIZ])
941 >                        ipb.beg = inpbuf;
942 >        }
943 > }
944 >
945 >
946 > static void
947 > skipinp(void)                   /* rewind position and advance 1 */
948 > {
949 >        if (ipb.beg == NULL)            /* full */
950                  ipb.beg = ipb.end;
951          ipb.pos = ipb.beg;
952          ipb.chr = *ipb.pos;
953          if (passive)                    /* transmit unmatched character? */
954 <                fputc(ipb.chr, stdout);
954 >                putchar(ipb.chr);
955          if (++ipb.beg >= &inpbuf[INBSIZ])
956                  ipb.beg = inpbuf;
957          scaninp();
# Line 859 | Line 959 | resetinp(void)                      /* rewind position
959  
960  
961   void
962 < eputs(char  *msg)
962 > eputs(char *msg)
963   {
964          fputs(msg, stderr);
965   }
966  
967  
968   void
969 < wputs(char  *msg)
969 > wputs(char *msg)
970   {
971          if (!nowarn)
972                  eputs(msg);
# Line 874 | Line 974 | wputs(char  *msg)
974  
975  
976   void
977 < quit(int  code)
977 > quit(int code)
978   {
979          exit(code);
980   }

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines