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

Comparing ray/src/util/radcompare.c (file contents):
Revision 2.6 by greg, Mon Oct 15 22:38:31 2018 UTC vs.
Revision 2.13 by greg, Fri Oct 19 23:50:32 2018 UTC

# Line 74 | Line 74 | LUTAB  hdr2 = LU_SINIT(free,free);
74   #define adv_linecnt(htp)        (lin1cnt += (htp == &hdr1), \
75                                          lin2cnt += (htp == &hdr2))
76  
77 + typedef struct {                /* dynamic line buffer */
78 +        char    *ptr;
79 +        int     len;
80 +        int     siz;
81 + } LINEBUF;
82 +
83 + #define init_line(bp)   ((bp)->ptr = NULL, (bp)->siz = 0)
84 +                                /* 100 MByte limit on line buffer */
85 + #define MAXBUF          (100L<<20)
86 +
87                                  /* input files */
88   char            *progname = NULL;
89   const char      stdin_name[] = "<stdin>";
# Line 82 | Line 92 | FILE           *f1in=NULL, *f2in=NULL;
92  
93                                  /* running real differences */
94   double  diff2sum = 0;
95 < int     nsum = 0;
95 > long    nsum = 0;
96  
97   /* Report usage and exit */
98   static void
# Line 92 | Line 102 | usage()
102          fputs(progname, stderr);
103          fputs(" [-h][-s|-w|-v][-rel min_test][-rms epsilon][-max epsilon] reference test\n",
104                          stderr);
105 <        exit(1);
105 >        exit(2);
106   }
107  
108 + /* Read a text line, increasing buffer size as necessary */
109 + static int
110 + read_line(LINEBUF *bp, FILE *fp)
111 + {
112 +        bp->len = 0;
113 +        if (!bp->ptr) {
114 +                bp->ptr = (char *)malloc(bp->siz = 512);
115 +                if (!bp->ptr)
116 +                        goto memerr;
117 +        }
118 +        while (fgets(bp->ptr + bp->len, bp->siz - bp->len, fp)) {
119 +                bp->len += strlen(bp->ptr + bp->len);
120 +                if (bp->ptr[bp->len-1] == '\n')
121 +                        break;          /* found EOL */
122 +                if (bp->len < bp->siz - 4)
123 +                        continue;       /* at EOF? */
124 +                if (bp->siz >= MAXBUF)
125 +                        break;          /* don't go to extremes */
126 +                if ((bp->siz += bp->siz/2) > MAXBUF)
127 +                        bp->siz = MAXBUF;
128 +                bp->ptr = (char *)realloc(bp->ptr, bp->siz);
129 +                if (!bp->ptr)
130 +                        goto memerr;
131 +        }
132 +        return(bp->len);
133 + memerr:
134 +        fprintf(stderr,
135 +                "%s: out of memory in read_line() allocating %d byte buffer\n",
136 +                        progname, bp->siz);
137 +        exit(2);
138 + }
139 +
140 + /* Free line buffer */
141 + static void
142 + free_line(LINEBUF *bp)
143 + {
144 +        if (bp->ptr) free(bp->ptr);
145 +        init_line(bp);
146 + }
147 +
148   /* Get type ID from name (or 0 if not found) */
149   static int
150   xlate_type(const char *nm)
# Line 134 | Line 184 | real_check(double r1, double r2)
184          return(1);
185   }
186  
187 + /* Compare two color values for equivalence */
188 + static int
189 + color_check(COLOR c1, COLOR c2)
190 + {
191 +        int     p;
192 +
193 +        if (!real_check(colval(c1,RED)+colval(c1,GRN)+colval(c1,BLU)*(1./3.),
194 +                        colval(c2,RED)+colval(c2,GRN)+colval(c2,BLU))*(1./3.))
195 +                return(0);
196 +
197 +        p = (colval(c1,GRN) > colval(c1,RED)) ? GRN : RED;
198 +        if (colval(c1,BLU) > colval(c1,p)) p = BLU;
199 +        
200 +        return(real_check(colval(c1,p), colval(c2,p)));
201 + }
202 +
203   /* Compare two strings for equivalence */
204   static int
205   equiv_string(char *s1, char *s2)
# Line 275 | Line 341 | match_val(const LUENT *ep1, void *p2)
341  
342   /* Compare two sets of header variables */
343   static int
344 < headers_match(LUTAB *hp1, LUTAB *hp2)
344 > headers_match()
345   {
346 <        int     ne = lu_doall(hp1, match_val, hp2);
346 >        int     ne = lu_doall(&hdr1, match_val, &hdr2);
347          if (ne < 0)
348                  return(0);      /* something didn't match! */
349                                  /* non-fatal if second header has extra */
350 <        if (report >= REP_WARN && (ne = lu_doall(hp2, NULL, NULL) - ne))
350 >        if (report >= REP_WARN && (ne = lu_doall(&hdr2, NULL, NULL) - ne))
351                  printf("%s: warning - '%s' has %d extra header setting(s)\n",
352                                          progname, f2name, ne);
353          return(1);              /* good match */
# Line 429 | Line 495 | compare_binary()
495   static int
496   compare_text()
497   {
498 <        char    l1buf[4096], l2buf[4096];
498 >        LINEBUF l1buf, l2buf;
499  
500          if (report >= REP_VERBOSE) {
501                  fputs(progname, stdout);
502                  fputs(": comparing inputs as ASCII text\n", stdout);
503          }
504 <                                                /* compare a line at a time */
505 <        while (fgets(l1buf, sizeof(l1buf), f1in)) {
504 >        init_line(&l1buf); init_line(&l2buf);   /* compare a line at a time */
505 >        while (read_line(&l1buf, f1in)) {
506                  lin1cnt++;
507 <                if (!*sskip2(l1buf,0))
507 >                if (!*sskip2(l1buf.ptr,0))
508                          continue;               /* ignore empty lines */
509 <                while (fgets(l2buf, sizeof(l2buf), f2in)) {
509 >
510 >                while (read_line(&l2buf, f2in)) {
511                          lin2cnt++;
512 <                        if (*sskip2(l2buf,0))
512 >                        if (*sskip2(l2buf.ptr,0))
513                                  break;          /* found other non-empty line */
514                  }
515 <                if (feof(f2in)) {
515 >                if (!l2buf.len) {               /* input 2 EOF? */
516                          if (report != REP_QUIET) {
517                                  fputs(f2name, stdout);
518                                  fputs(": unexpected end-of-file\n", stdout);
519                          }
520 +                        free_line(&l1buf); free_line(&l2buf);
521                          return(0);
522                  }
523                                                  /* compare non-empty lines */
524 <                if (!equiv_string(l1buf, l2buf)) {
524 >                if (!equiv_string(l1buf.ptr, l2buf.ptr)) {
525                          if (report != REP_QUIET) {
526                                  printf("%s: inputs '%s' and '%s' differ at line %d|%d\n",
527                                                  progname, f1name, f2name,
528                                                  lin1cnt, lin2cnt);
529 <                                if (report >= REP_VERBOSE) {
529 >                                if ( report >= REP_VERBOSE &&
530 >                                                (l1buf.len < 256) &
531 >                                                (l2buf.len < 256) ) {
532                                          fputs("------------- Mismatch -------------\n", stdout);
533                                          printf("%s@%d:\t%s", f1name,
534 <                                                        lin1cnt, l1buf);
534 >                                                        lin1cnt, l1buf.ptr);
535                                          printf("%s@%d:\t%s", f2name,
536 <                                                        lin2cnt, l2buf);
536 >                                                        lin2cnt, l2buf.ptr);
537                                  }
538                          }
539 +                        free_line(&l1buf); free_line(&l2buf);
540                          return(0);
541                  }
542          }
543 <                                                /* check for EOF on input 2 */
544 <        while (fgets(l2buf, sizeof(l2buf), f2in)) {
545 <                if (!*sskip2(l2buf,0))
543 >        free_line(&l1buf);                      /* check for EOF on input 2 */
544 >        while (read_line(&l2buf, f2in)) {
545 >                if (!*sskip2(l2buf.ptr,0))
546                          continue;
547                  if (report != REP_QUIET) {
548                          fputs(f1name, stdout);
549                          fputs(": unexpected end-of-file\n", stdout);
550                  }
551 +                free_line(&l2buf);
552                  return(0);
553          }
554 +        free_line(&l2buf);
555          return(good_RMS());                     /* final check for reals */
556   }
557  
# Line 527 | Line 600 | compare_hdr()
600                          return(0);
601                  }
602                  for (x = 0; x < scanlen(&rs1); x++) {
603 <                        if (real_check(colval(scan1[x],RED),
531 <                                                colval(scan2[x],RED)) &
532 <                                        real_check(colval(scan1[x],GRN),
533 <                                                colval(scan2[x],GRN)) &
534 <                                        real_check(colval(scan1[x],BLU),
535 <                                                colval(scan2[x],BLU)))
603 >                        if (color_check(scan1[x], scan2[x]))
604                                  continue;
605                          if (report != REP_QUIET) {
606                                  printf(
# Line 568 | Line 636 | compare_float()
636  
637          if (report >= REP_VERBOSE) {
638                  fputs(progname, stdout);
639 <                fputs(": comparing inputs as 32-bit IEEE float\n", stdout);
639 >                fputs(": comparing inputs as 32-bit IEEE floats\n", stdout);
640          }
641          while (getbinary(&f1, sizeof(f1), 1, f1in)) {
642                  if (!getbinary(&f2, sizeof(f2), 1, f2in))
# Line 598 | Line 666 | compare_double()
666  
667          if (report >= REP_VERBOSE) {
668                  fputs(progname, stdout);
669 <                fputs(": comparing inputs as 64-bit IEEE double\n", stdout);
669 >                fputs(": comparing inputs as 64-bit IEEE doubles\n", stdout);
670          }
671          while (getbinary(&f1, sizeof(f1), 1, f1in)) {
672                  if (!getbinary(&f2, sizeof(f2), 1, f2in))
# Line 662 | Line 730 | main(int argc, char *argv[])
730                  }
731                  break;
732          }
733 <        if (a != argc-2)
733 >        if (a != argc-2)                        /* make sure of two inputs */
734                  usage();
735 <        if (!f1name) f1name = argv[a];
668 <        if (!f2name) f2name = argv[a+1];
669 <
670 <        if (!strcmp(f1name, f2name)) {          /* inputs are same? */
735 >        if (!strcmp(argv[a], argv[a+1])) {      /* inputs are same? */
736                  if (report >= REP_WARN)
737                          printf("%s: warning - identical inputs given\n",
738                                          progname);
739                  return(0);
740          }
741 +        if (!f1name) f1name = argv[a];
742 +        if (!f2name) f2name = argv[a+1];
743                                                  /* open inputs */
744          SET_FILE_BINARY(stdin);                 /* in case we're using it */
745          if (!f1in && !(f1in = fopen(f1name, "rb"))) {
746                  fprintf(stderr, "%s: cannot open for reading\n", f1name);
747 <                return(1);
747 >                return(2);
748          }
749          if (!strcmp(f2name, "-")) {
750                  f2in = stdin;
751                  f2name = stdin_name;
752          } else if (!(f2in = fopen(f2name, "rb"))) {
753                  fprintf(stderr, "%s: cannot open for reading\n", f2name);
754 <                return(1);
754 >                return(2);
755          }
756                                                  /* load headers */
757          if ((typ1 = identify_type(f1name, f1in, &hdr1)) < 0)
758 <                return(1);
758 >                return(2);
759          if ((typ2 = identify_type(f2name, f2in, &hdr2)) < 0)
760 <                return(1);
760 >                return(2);
761          if (typ1 != typ2) {
762                  if (report != REP_QUIET)
763                          printf("%s: '%s' is %s and '%s' is %s\n",
# Line 699 | Line 766 | main(int argc, char *argv[])
766                  return(1);
767          }
768          ign_header |= !has_header(typ1);        /* check headers if indicated */
769 <        if (!ign_header && !headers_match(&hdr1, &hdr2))
769 >        if (!ign_header && !headers_match())
770                  return(1);
771          lu_done(&hdr1); lu_done(&hdr2);
772          if (!ign_header & (report >= REP_WARN)) {

Diff Legend

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