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.7 by schorsch, Sun Jul 27 22:12:01 2003 UTC vs.
Revision 1.31 by greg, Mon Feb 21 23:00:55 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  <stdio.h>
12 < #include  <string.h>
13 < #include  <math.h>
14 < #include  <ctype.h>
10 > #include <stdlib.h>
11 > #include <math.h>
12 > #include <ctype.h>
13  
14 < #include  "platform.h"
15 < #include  "calcomp.h"
16 < #include  "rterror.h"
14 > #include "platform.h"
15 > #include "rterror.h"
16 > #include "rtmisc.h"
17 > #include "rtio.h"
18 > #include "calcomp.h"
19  
20 < #ifdef  CPM
21 < #define  getc           agetc   /* text files only, right? */
22 < #endif
20 > #define isnum(c) (isdigit(c) || ((c)=='-') | ((c)=='.') \
21 >                                | ((c)=='+') | ((c)=='e') | ((c)=='E'))
22  
23 < #define  isnum(c)       (isdigit(c) || (c)=='-' || (c)=='.' \
25 <                                || (c)=='+' || (c)=='e' || (c)=='E')
23 > #define isblnk(c) (igneol ? isspace(c) : ((c)==' ')|((c)=='\t'))
24  
25 < #define  isblnk(c)      (igneol ? isspace(c) : (c)==' '||(c)=='\t')
25 > #define INBSIZ   16384  /* longest record */
26 > #define MAXCOL   32 /* number of columns recorded */
27  
29 – #define  INBSIZ         4096    /* longest record */
30 – #define  MAXCOL         32      /* number of columns recorded */
31 –
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 void scaninp(void), advinp(void), resetinp(void);
57 > static int getinputrec(FILE *fp);
58 > static void scaninp(void), advinp(void), skipinp(void);
59   static void putrec(void), putout(void), nbsynch(void);
60   static int getrec(void);
61   static void execute(char *file);
# Line 70 | Line 67 | static int getfield(struct field *f);
67   static void chanset(int n, double v);
68   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? */
80 < char  sepchar = '\t';           /* input/output separator */
81 < int  noinput = 0;               /* no input records? */
82 < int  nbicols = 0;               /* number of binary input columns */
83 < int  bocols = 0;                /* produce binary output columns */
84 < char  inpbuf[INBSIZ];           /* input buffer */
85 < double  colval[MAXCOL];         /* input column values */
86 < unsigned long  colflg = 0;      /* column retrieved flags */
87 < 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);
120 +        funset("in", 1, '=', &l_in);
121  
122          for (i = 1; i < argc && argv[i][0] == '-'; i++)
123                  switch (argv[i][1]) {
# Line 122 | Line 127 | char  *argv[]
127                  case 'l':
128                          igneol = !igneol;
129                          break;
130 +                case 'p':
131 +                        passive = !passive;
132 +                        break;
133                  case 't':
134                          sepchar = argv[i][2];
135                          break;
# Line 129 | Line 137 | char  *argv[]
137                          svpreset(argv[++i]);
138                          break;
139                  case 'f':
140 <                        fcompile(argv[++i]);
140 >                        fpath = getpath(argv[++i], getrlibpath(), 0);
141 >                        if (fpath == NULL) {
142 >                                eputs(argv[0]);
143 >                                eputs(": cannot find file '");
144 >                                eputs(argv[i]);
145 >                                eputs("'\n");
146 >                                quit(1);
147 >                        }
148 >                        fcompile(fpath);
149                          break;
150                  case 'e':
151                          scompile(argv[++i], NULL, 0);
# Line 140 | Line 156 | char  *argv[]
156                  case 'i':
157                          switch (argv[i][2]) {
158                          case '\0':
159 +                                itype = 'a';
160                                  nbicols = 0;
161                                  readfmt(argv[++i], 0);
162                                  break;
163 +                        case 'n':
164 +                                incnt = atol(argv[++i]);
165 +                                break;
166                          case 'a':
167 +                                itype = 'a';
168                                  nbicols = 0;
169                                  break;
170                          case 'd':
171 +                        case 'D':
172 +                                itype = argv[i][2];
173                                  if (isdigit(argv[i][3]))
174                                          nbicols = atoi(argv[i]+3);
175                                  else
176                                          nbicols = 1;
177 +                                if (nbicols*sizeof(double) > INBSIZ) {
178 +                                        eputs(argv[0]);
179 +                                        eputs(": too many input columns\n");
180 +                                        quit(1);
181 +                                }
182                                  break;
183                          case 'f':
184 +                        case 'F':
185 +                                itype = argv[i][2];
186                                  if (isdigit(argv[i][3]))
187 <                                        nbicols = -atoi(argv[i]+3);
187 >                                        nbicols = atoi(argv[i]+3);
188                                  else
189 <                                        nbicols = -1;
189 >                                        nbicols = 1;
190 >                                if (nbicols*sizeof(float) > INBSIZ) {
191 >                                        eputs(argv[0]);
192 >                                        eputs(": too many input columns\n");
193 >                                        quit(1);
194 >                                }
195                                  break;
196                          default:
197                                  goto userr;
# Line 165 | Line 200 | char  *argv[]
200                  case 'o':
201                          switch (argv[i][2]) {
202                          case '\0':
203 <                                bocols = 0;
203 >                                otype = 'a';
204                                  readfmt(argv[++i], 1);
205                                  break;
206 +                        case 'n':
207 +                                outcnt = atol(argv[++i]);
208 +                                break;
209                          case 'a':
210 <                                bocols = 0;
210 >                                otype = 'a';
211                                  break;
212                          case 'd':
213 <                                bocols = 1;
176 <                                break;
213 >                        case 'D':
214                          case 'f':
215 <                                bocols = -1;
215 >                        case 'F':
216 >                                otype = argv[i][2];
217                                  break;
218 +                        default:
219 +                                goto userr;
220                          }
221                          break;
222                  case 'w':
# Line 189 | Line 229 | char  *argv[]
229                  userr:
230                          eputs("Usage: ");
231                          eputs(argv[0]);
232 < eputs(" [-b][-l][-n][-w][-u][-tS][-s svar=sval][-e expr][-f source][-i infmt][-o outfmt] [file]\n");
232 > eputs(" [-b][-l][-n][-p][-w][-u][-tS][-s svar=sval][-e expr][-f source][-i infmt][-o outfmt] [file]\n");
233                          quit(1);
234                  }
235 <
236 <        if (noinput) {          /* produce a single output record */
235 >        if (otype != 'a')
236 >                SET_FILE_BINARY(stdout);
237 > #ifdef getc_unlocked            /* avoid lock/unlock overhead */
238 >        flockfile(stdout);
239 > #endif
240 >        if (noinput) {   /* produce a single output record */
241 >                if (i < argc) {
242 >                        eputs(argv[0]);
243 >                        eputs(": file argument(s) incompatible with -n\n");
244 >                        quit(1);
245 >                }
246                  eclock++;
247                  putout();
248                  quit(0);
249          }
250 <
202 <        if (blnkeq)             /* for efficiency */
250 >        if (blnkeq)      /* for efficiency */
251                  nbsynch();
252  
253 <        if (i == argc)          /* from stdin */
253 >        if (i == argc)   /* from stdin */
254                  execute(NULL);
255 <        else                    /* from one or more files */
255 >        else             /* from one or more files */
256                  for ( ; i < argc; i++)
257                          execute(argv[i]);
258          
259          quit(0);
260 +        return 0; /* pro forma return */
261   }
262  
263  
264   static void
265 < nbsynch(void)               /* non-blank starting synch character */
265 > nbsynch(void)    /* non-blank starting synch character */
266   {
267          if (inpfmt == NULL || (inpfmt->type & F_TYP) != T_LIT)
268                  return;
# Line 224 | Line 273 | nbsynch(void)               /* non-blank starting sync
273   }
274  
275  
276 < int
276 > static int
277   getinputrec(            /* get next input record */
278 < FILE  *fp
278 > FILE *fp
279   )
280   {
281          if (inpfmt != NULL)
282                  return(getrec());
283 <        if (nbicols > 0)
284 <                return(fread(inpbuf, sizeof(double),
285 <                                        nbicols, fp) == nbicols);
286 <        if (nbicols < 0)
287 <                return(fread(inpbuf, sizeof(float),
288 <                                        -nbicols, fp) == -nbicols);
283 >        if ((itype == 'd') | (itype == 'D')) {
284 >                if (getbinary(inpbuf, sizeof(double), nbicols, fp) != nbicols)
285 >                        return(0);
286 >                if (itype == 'D')
287 >                        swap64(inpbuf, nbicols);
288 >                return(1);
289 >        }
290 >        if ((itype == 'f') | (itype == 'F')) {
291 >                if (getbinary(inpbuf, sizeof(float), nbicols, fp) != nbicols)
292 >                        return(0);
293 >                if (itype == 'F')
294 >                        swap32(inpbuf, nbicols);
295 >                return(1);
296 >        }
297          return(fgets(inpbuf, INBSIZ, fp) != NULL);
298   }
299  
300  
301   static void
302 < execute(           /* process a file */
303 < char  *file
302 > execute(         /* process a file */
303 > char *file
304   )
305   {
306 <        int  conditional = vardefined("cond");
307 <        long  nrecs = 0;
308 <        long  nout = 0;
309 <        FILE  *fp;
306 >        const int conditional = vardefined("cond");
307 >        const int set_recno = (varlookup("recno") != NULL);
308 >        const int set_outno = (varlookup("outno") != NULL);
309 >        long nrecs = 0;
310 >        long nout = 0;
311 >        FILE *fp;
312          
313          if (file == NULL)
314                  fp = stdin;
# Line 258 | Line 317 | char  *file
317                  eputs(": cannot open\n");
318                  quit(1);
319          }
320 +        if (itype != 'a')
321 +                SET_FILE_BINARY(fp);
322 + #ifdef getc_unlocked            /* avoid lock/unlock overhead */
323 +        flockfile(fp);
324 + #endif
325          if (inpfmt != NULL)
326                  initinp(fp);
327          
328          while (getinputrec(fp)) {
329 <                varset("recno", '=', (double)++nrecs);
329 >                ++nrecs;
330 >                if (set_recno)
331 >                        varset("recno", '=', (double)nrecs);
332 >                if (set_outno)
333 >                        varset("outno", '=', (double)(nout+1));
334                  colflg = 0;
335                  eclock++;
336                  if (!conditional || varvalue("cond") > 0.0) {
269 –                        varset("outno", '=', (double)++nout);
337                          putout();
338 +                        ++nout;
339                  }
340 +                if (incnt && nrecs >= incnt)
341 +                        break;
342 +                if (outcnt && nout >= outcnt)
343 +                        break;
344          }
345          fclose(fp);
346   }
347  
348  
349   static void
350 < putout(void)                /* produce an output record */
350 > putout(void)            /* produce an output record */
351   {
352  
353          colpos = 0;
354          if (outfmt != NULL)
355                  putrec();
356 <        else if (bocols)
285 <                chanout(bchanset);
286 <        else
356 >        else if (otype == 'a')
357                  chanout(chanset);
358 <        if (colpos && !bocols)
358 >        else
359 >                chanout(bchanset);
360 >        if (colpos && otype == 'a')
361                  putchar('\n');
362          if (unbuff)
363                  fflush(stdout);
364   }
365  
366  
367 + static double
368 + l_in(char *funame)      /* function call for $channel */
369 + {
370 +        int n;
371 +        char *cp;
372 +                        /* get argument as integer */
373 +        n = (int)(argument(1) + .5);
374 +        if (n != 0)     /* return channel value */
375 +                return(chanvalue(n));
376 +                        /* determine number of channels */
377 +        if (noinput || inpfmt != NULL)
378 +                return(0);
379 +        if (nbicols)
380 +                return(nbicols);
381 +        cp = inpbuf;    /* need to count */
382 +        for (n = 0; *cp; )
383 +                if (blnkeq && isspace(sepchar)) {
384 +                        while (isspace(*cp))
385 +                                cp++;
386 +                        n += *cp != '\0';
387 +                        while (*cp && !isspace(*cp))
388 +                                cp++;
389 +                } else {
390 +                        n += *cp != '\n';
391 +                        while (*cp && *cp++ != sepchar)
392 +                                ;
393 +                }
394 +        return(n);
395 + }
396 +
397   double
398 < chanvalue(            /* return value for column n */
399 < int  n
398 > chanvalue(       /* return value for column n */
399 > int n
400   )
401   {
402 <        int  i;
403 <        register char  *cp;
402 >        int i;
403 >        char *cp;
404  
405          if (noinput || inpfmt != NULL) {
406                  eputs("no column input\n");
# Line 308 | Line 410 | int  n
410                  eputs("illegal channel number\n");
411                  quit(1);
412          }
413 <        if (nbicols > 0) {
413 >        if (nbicols) {
414                  if (n > nbicols)
415                          return(0.0);
416 <                cp = inpbuf + (n-1)*sizeof(double);
417 <                return(*(double *)cp);
418 <        }
419 <        if (nbicols < 0) {
318 <                if (n > -nbicols)
319 <                        return(0.0);
416 >                if ((itype == 'd') | (itype == 'D')) {
417 >                        cp = inpbuf + (n-1)*sizeof(double);
418 >                        return(*(double *)cp);
419 >                }
420                  cp = inpbuf + (n-1)*sizeof(float);
421                  return(*(float *)cp);
422          }
# Line 334 | Line 434 | int  n
434                          while (*cp && *cp++ != sepchar)
435                                  ;
436  
437 <        while (isspace(*cp))            /* some atof()'s don't like tabs */
437 >        while (isspace(*cp))     /* some atof()'s don't like tabs */
438                  cp++;
439  
440          if (n <= MAXCOL) {
# Line 346 | Line 446 | int  n
446  
447  
448   void
449 < chanset(                   /* output column n */
450 < int  n,
451 < double  v
449 > chanset(                 /* output column n */
450 > int n,
451 > double v
452   )
453   {
454 <        if (colpos == 0)                /* no leading separator */
454 >        if (colpos == 0)                /* no leading separator */
455                  colpos = 1;
456          while (colpos < n) {
457                  putchar(sepchar);
# Line 362 | Line 462 | double  v
462  
463  
464   void
465 < bchanset(                   /* output binary channel n */
466 < int  n,
467 < double  v
465 > bchanset(                /* output binary channel n */
466 > int n,
467 > double v
468   )
469   {
470          static char     zerobuf[sizeof(double)];
471 +        const int       otlen = ((otype == 'd') | (otype == 'D')) ?
472 +                                        sizeof(double) : sizeof(float);
473 +        float   fval = v;
474  
475          while (++colpos < n)
476 <                fwrite(zerobuf,
477 <                        bocols>0 ? sizeof(double) : sizeof(float),
478 <                        1, stdout);
479 <        if (bocols > 0)
480 <                fwrite(&v, sizeof(double), 1, stdout);
481 <        else {
482 <                float   fval = v;
483 <                fwrite(&fval, sizeof(float), 1, stdout);
476 >                putbinary(zerobuf, otlen, 1, stdout);
477 >        switch (otype) {
478 >        case 'D':
479 >                swap64((char *)&v, 1);
480 >                /* fall through */
481 >        case 'd':
482 >                putbinary(&v, sizeof(double), 1, stdout);
483 >                break;
484 >        case 'F':
485 >                swap32((char *)&fval, 1);
486 >                /* fall through */
487 >        case 'f':
488 >                putbinary(&fval, sizeof(float), 1, stdout);
489 >                break;
490          }
491   }
492  
493  
494   static void
495 < readfmt(                   /* read record format */
496 < char  *spec,
497 < int  output
495 > readfmt(                 /* read record format */
496 > char *spec,
497 > int output
498   )
499   {
500 <        int  fd;
501 <        char  *inptr;
502 <        struct field  fmt;
503 <        int  res;
504 <        register struct field  *f;
500 >        int fd;
501 >        char *inptr;
502 >        struct field fmt;
503 >        int res;
504 >        struct field *f;
505                                                  /* check for inline format */
506          for (inptr = spec; *inptr; inptr++)
507                  if (*inptr == '$')
508                          break;
509 <        if (*inptr)                             /* inline */
509 >        if (*inptr)                      /* inline */
510                  inptr = spec;
511 <        else {                                  /* from file */
511 >        else {                           /* from file */
512                  if ((fd = open(spec, 0)) == -1) {
513                          eputs(spec);
514                          eputs(": cannot open\n");
515                          quit(1);
516                  }
517 <                res = read(fd, inpbuf+1, INBSIZ-1);
517 >                res = read(fd, inpbuf+2, INBSIZ-2);
518                  if (res <= 0 || res >= INBSIZ-1) {
519                          eputs(spec);
520                          if (res < 0)
# Line 417 | Line 526 | int  output
526                          quit(1);
527                  }
528                  close(fd);
529 <                (inptr=inpbuf+1)[res] = '\0';
529 >                (inptr=inpbuf+2)[res] = '\0';
530          }
531 <        f = &fmt;                               /* get fields */
531 >        f = &fmt;                        /* get fields */
532          while ((res = readfield(&inptr)) != F_NUL) {
533                  f->next = (struct field *)emalloc(sizeof(struct field));
534                  f = f->next;
# Line 451 | Line 560 | int  output
560  
561  
562   static int
563 < readfield(                   /* get next field in format */
564 < register char  **pp
563 > readfield(               /* get next field in format */
564 > char **pp
565   )
566   {
567 <        int  type = F_NUL;
568 <        int  width = 0;
569 <        register char  *cp;
567 >        int type = F_NUL;
568 >        int width = 0;
569 >        char *cp;
570          
571          cp = inpbuf;
572          while (cp < &inpbuf[INBSIZ-1] && **pp != '\0') {
# Line 518 | Line 627 | register char  **pp
627  
628  
629   struct strvar *
630 < getsvar(                         /* get string variable */
631 < char  *svname
630 > getsvar(                         /* get string variable */
631 > char *svname
632   )
633   {
634 <        register struct strvar  *sv;
634 >        struct strvar *sv;
635          
636          for (sv = svhead; sv != NULL; sv = sv->next)
637                  if (!strcmp(sv->name, svname))
# Line 537 | Line 646 | char  *svname
646  
647  
648   static void
649 < svpreset(                    /* preset a string variable */
650 < char  *eqn
649 > svpreset(                /* preset a string variable */
650 > char *eqn
651   )
652   {
653 <        register struct strvar  *sv;
654 <        register char  *val;
653 >        struct strvar *sv;
654 >        char *val;
655  
656          for (val = eqn; *val != '='; val++)
657                  if (!*val)
# Line 561 | Line 670 | char  *eqn
670   static void
671   clearrec(void)                  /* clear input record variables */
672   {
673 <        register struct field  *f;
673 >        struct field *f;
674  
675          for (f = inpfmt; f != NULL; f = f->next)
676                  switch (f->type & F_TYP) {
# Line 579 | Line 688 | clearrec(void)                 /* clear input record variables */
688  
689  
690   static int
691 < getrec(void)                                /* get next record from file */
691 > getrec(void)                            /* get next record from file */
692   {
693 <        int  eatline;
694 <        register struct field  *f;
695 <        
693 >        int eatline;
694 >        struct field *f;
695 >
696          while (ipb.chr != EOF) {
697 <                eatline = !igneol && ipb.chr != '\n';
589 <                if (blnkeq)             /* beware of nbsynch() */
697 >                if (blnkeq) {           /* beware of nbsynch() */
698                          while (isblnk(ipb.chr))
699 <                                scaninp();
699 >                                skipinp();
700 >                        if (ipb.chr == EOF)
701 >                                return(0);
702 >                }
703 >                eatline = !igneol & (ipb.chr != '\n');
704                  clearrec();             /* start with fresh record */
705                  for (f = inpfmt; f != NULL; f = f->next)
706 <                        if (getfield(f) == -1)
706 >                        if (!getfield(f))
707                                  break;
708                  if (f == NULL) {
709 <                        advinp();
709 >                        advinp(); /* got one! */
710                          return(1);
711                  }
712 <                resetinp();
713 <                if (eatline) {          /* eat rest of line */
712 >                skipinp();              /* eat false start */
713 >                if (eatline) {   /* eat rest of line */
714                          while (ipb.chr != '\n') {
715                                  if (ipb.chr == EOF)
716                                          return(0);
717 <                                scaninp();
717 >                                skipinp();
718                          }
719 <                        scaninp();
608 <                        advinp();
719 >                        skipinp();
720                  }
721          }
722          return(0);
# Line 613 | Line 724 | getrec(void)                                /* get nex
724  
725  
726   static int
727 < getfield(                             /* get next field */
728 < register struct field  *f
727 > getfield(                        /* get next field */
728 > struct field *f
729   )
730   {
731 <        static char  buf[RMAXWORD+1];            /* no recursion! */
732 <        int  delim, inword;
733 <        double  d;
734 <        char  *np;
735 <        register char  *cp;
731 >        static char buf[RMAXWORD+1];     /* no recursion! */
732 >        int delim, inword;
733 >        double d;
734 >        char *np;
735 >        char *cp;
736  
737          switch (f->type & F_TYP) {
738          case T_LIT:
# Line 629 | Line 740 | register struct field  *f
740                  do {
741                          if (blnkeq && isblnk(*cp)) {
742                                  if (!isblnk(ipb.chr))
743 <                                        return(-1);
743 >                                        return(0);
744                                  do
745                                          cp++;
746                                  while (isblnk(*cp));
# Line 640 | Line 751 | register struct field  *f
751                                  cp++;
752                                  scaninp();
753                          } else
754 <                                return(-1);
754 >                                return(0);
755                  } while (*cp);
756 <                return(0);
756 >                break;
757          case T_STR:
758                  if (f->next == NULL || (f->next->type & F_TYP) != T_LIT)
759                          delim = EOF;
# Line 650 | Line 761 | register struct field  *f
761                          delim = f->next->f.sl[0];
762                  cp = buf;
763                  do {
764 <                        if (ipb.chr == EOF)
764 >                        if ((ipb.chr == EOF) | (ipb.chr == '\n'))
765                                  inword = 0;
766                          else if (blnkeq && delim != EOF)
767                                  inword = isblnk(delim) ?
# Line 667 | Line 778 | register struct field  *f
778                  if (f->f.sv->val == NULL)
779                          f->f.sv->val = savqstr(buf);    /* first setting */
780                  else if (strcmp(f->f.sv->val, buf))
781 <                        return(-1);                     /* doesn't match! */
782 <                return(0);
781 >                        return(0);                      /* doesn't match! */
782 >                break;
783          case T_NUM:
784                  if (f->next == NULL || (f->next->type & F_TYP) != T_LIT)
785                          delim = EOF;
# Line 698 | Line 809 | register struct field  *f
809                          varset(f->f.nv, '=', d);        /* first setting */
810                  else if ((d = (varvalue(f->f.nv)-d)/(d==0.?1.:d)) > .001
811                                  || d < -.001)
812 <                        return(-1);                     /* doesn't match! */
813 <                return(0);
812 >                        return(0);                      /* doesn't match! */
813 >                break;
814          }
815 <        return -1; /* pro forma return */
815 >        return(1);      /* success! */
816   }
817  
818  
819   static void
820 < putrec(void)                                /* output a record */
820 > putrec(void)                            /* output a record */
821   {
822 <        char  fmt[32];
823 <        register int  n;
824 <        register struct field  *f;
825 <        int  adlast, adnext;
822 >        char fmt[32], typ[16];
823 >        int n;
824 >        struct field *f;
825 >        int adlast, adnext;
826 >        double dv, av;
827          
828          adlast = 0;
829          for (f = outfmt; f != NULL; f = f->next) {
830 <                adnext =        blnkeq &&
830 >                adnext =        blnkeq &&
831                                  f->next != NULL &&
832                                  !( (f->next->type&F_TYP) == T_LIT &&
833                                          f->next->f.sl[0] == ' ' );
# Line 742 | Line 854 | putrec(void)                                /* output
854                          break;
855                  case T_NUM:
856                          n = f->type & F_WID;
857 +                        dv = evalue(f->f.ne);
858 +                        av = fabs(dv);
859 +                        if (n <= 9)
860 +                                strcpy(typ, "g");
861 +                        else
862 +                                sprintf(typ, ".%de", n-5);
863 +                        if (av < 1L<<31) {
864 +                                long    iv = (int)(av + .5);
865 +                                if (iv && fabs(av-iv) <= av*1e-14)
866 +                                        strcpy(typ, ".0f");
867 +                        }
868                          if (adlast && adnext)
869 <                                strcpy(fmt, "%g");
869 >                                sprintf(fmt, "%%%s", typ);
870                          else if (adlast)
871 <                                sprintf(fmt, "%%-%dg", n);
871 >                                sprintf(fmt, "%%-%d%s", n, typ);
872                          else
873 <                                sprintf(fmt, "%%%dg", n);
874 <                        printf(fmt, evalue(f->f.ne));
873 >                                sprintf(fmt, "%%%d%s", n, typ);
874 >                        printf(fmt, dv);
875                          adlast = 1;
876                          break;
877                  }
# Line 757 | Line 880 | putrec(void)                                /* output
880  
881  
882   static void
883 < initinp(FILE  *fp)                     /* prepare lookahead buffer */
883 > initinp(FILE *fp)                /* prepare lookahead buffer */
884  
885   {
886          ipb.fin = fp;
887          ipb.beg = ipb.end = inpbuf;
888 <        ipb.pos = inpbuf-1;             /* position before beginning */
888 >        ipb.pos = inpbuf-1;             /* position before beginning */
889          ipb.chr = '\0';
890          scaninp();
891   }
892  
893  
894   static void
895 < scaninp(void)                       /* scan next character */
895 > scaninp(void)                   /* scan next character */
896   {
897          if (ipb.chr == EOF)
898                  return;
899          if (++ipb.pos >= &inpbuf[INBSIZ])
900                  ipb.pos = inpbuf;
901 <        if (ipb.pos == ipb.end) {               /* new character */
901 >        if (ipb.pos == ipb.end) {        /* new character */
902                  if ((ipb.chr = getc(ipb.fin)) != EOF) {
903                          *ipb.end = ipb.chr;
904                          if (++ipb.end >= &inpbuf[INBSIZ])
# Line 789 | Line 912 | scaninp(void)                       /* scan next chara
912  
913  
914   static void
915 < advinp(void)                        /* move home to current position */
915 > advinp(void)                    /* move home to current position */
916   {
917          ipb.beg = ipb.pos;
918   }
919  
920  
921   static void
922 < resetinp(void)                      /* rewind position and advance 1 */
922 > skipinp(void)            /* rewind position and advance 1 */
923   {
924 <        if (ipb.beg == NULL)            /* full */
924 >        if (ipb.beg == NULL)    /* full */
925                  ipb.beg = ipb.end;
926          ipb.pos = ipb.beg;
927          ipb.chr = *ipb.pos;
928 +        if (passive)                    /* transmit unmatched character? */
929 +                putchar(ipb.chr);
930          if (++ipb.beg >= &inpbuf[INBSIZ])
931                  ipb.beg = inpbuf;
932          scaninp();
# Line 809 | Line 934 | resetinp(void)                      /* rewind position
934  
935  
936   void
937 < eputs(char  *msg)
937 > eputs(char *msg)
938   {
939          fputs(msg, stderr);
940   }
941  
942  
943   void
944 < wputs(char  *msg)
944 > wputs(char *msg)
945   {
946          if (!nowarn)
947                  eputs(msg);
# Line 824 | Line 949 | wputs(char  *msg)
949  
950  
951   void
952 < quit(int  code)
952 > quit(int code)
953   {
954          exit(code);
955   }

Diff Legend

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