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.1 by greg, Mon Oct 15 17:52:52 2018 UTC vs.
Revision 2.9 by greg, Tue Oct 16 16:52:05 2018 UTC

# Line 33 | Line 33 | double max_lim = 0.25;         /* difference limit if non-neg
33  
34   int     lin1cnt=0, lin2cnt=0;   /* file line position */
35  
36 + const char      nsuffix[10][3] = {              /* 1st, 2nd, 3rd, etc. */
37 +                        "th","st","nd","rd","th","th","th","th","th","th"
38 +                };
39 + #define num_sfx(n)      nsuffix[(n)%10]
40 +
41                                  /* file types */
42   const char      *file_type[] = {
43                          "Unrecognized",
# Line 40 | Line 45 | const char     *file_type[] = {
45                          "ascii",
46                          COLRFMT,
47                          CIEFMT,
48 <                        "BINARY_float",
49 <                        "BINARY_double",
48 >                        "float",
49 >                        "double",
50                          "BSDF_RBFmesh",
51                          "Radiance_octree",
52                          "Radiance_tmesh",
53                          "BINARY_unknown",
54                          NULL    /* terminator */
55                  };
56 <
56 >                                /* keep consistent with above */
57   enum {TYP_UNKNOWN, TYP_TEXT, TYP_ASCII, TYP_RGBE, TYP_XYZE, TYP_FLOAT,
58                  TYP_DOUBLE, TYP_RBFMESH, TYP_OCTREE, TYP_TMESH, TYP_BINARY};
59                  
# Line 74 | Line 79 | char           *progname = NULL;
79   const char      stdin_name[] = "<stdin>";
80   const char      *f1name=NULL, *f2name=NULL;
81   FILE            *f1in=NULL, *f2in=NULL;
82 +
83                                  /* running real differences */
84   double  diff2sum = 0;
85 < int     nsum = 0;
85 > long    nsum = 0;
86  
87   /* Report usage and exit */
88   static void
# Line 84 | Line 90 | usage()
90   {
91          fputs("Usage: ", stderr);
92          fputs(progname, stderr);
93 <        fputs(" [-h][-s|-w|-v][-rel min_test][-rms epsilon][-max epsilon] file1 file2\n",
93 >        fputs(" [-h][-s|-w|-v][-rel min_test][-rms epsilon][-max epsilon] reference test\n",
94                          stderr);
95 <        exit(1);
95 >        exit(2);
96   }
97  
98   /* Get type ID from name (or 0 if not found) */
# Line 116 | Line 122 | real_check(double r1, double r2)
122          }
123          if (max_lim >= 0 && diff2 > max_lim*max_lim) {
124                  if (report != REP_QUIET)
125 <                        fprintf(stderr,
126 <                        "%s: %sdifference between %.8g and %.8g exceeds epsilon\n",
125 >                        printf(
126 >                        "%s: %sdifference between %.8g and %.8g exceeds epsilon of %.8g\n",
127                                          progname,
128                                          (rel_min > 0) ? "relative " : "",
129 <                                        r1, r2);
129 >                                        r1, r2, max_lim);
130                  return(0);
131          }
132          diff2sum += diff2;
# Line 128 | Line 134 | real_check(double r1, double r2)
134          return(1);
135   }
136  
137 + /* Compare two color values for equivalence */
138 + static int
139 + color_check(COLOR c1, COLOR c2)
140 + {
141 +        int     p;
142 +
143 +        if (!real_check(colval(c1,RED)+colval(c1,GRN)+colval(c1,BLU)*(1./3.),
144 +                        colval(c2,RED)+colval(c2,GRN)+colval(c2,BLU))*(1./3.))
145 +                return(0);
146 +
147 +        p = (colval(c1,GRN) > colval(c1,RED)) ? GRN : RED;
148 +        if (colval(c1,BLU) > colval(c1,p)) p = BLU;
149 +        
150 +        return(real_check(colval(c1,p), colval(c2,p)));
151 + }
152 +
153   /* Compare two strings for equivalence */
154   static int
155   equiv_string(char *s1, char *s2)
# Line 138 | Line 160 | equiv_string(char *s1, char *s2)
160                                  /* skip whitespace at beginning */
161          while (isspace(*s1)) s1++;
162          while (isspace(*s2)) s2++;
141 –        if (!strcmp(s1, s2))    /* quick check */
142 –                return(1);
163          while (*s1) {           /* check each word */
164 <                int     inquote = *s1;
164 >                int     inquote;
165 >                if (!*s2)       /* unexpected EOL in s2? */
166 >                        return(0);
167 >                inquote = *s1;
168                  if ((inquote != '\'') & (inquote != '"'))
169                          inquote = 0;
170                  if (inquote) {  /* quoted text must match exactly */
# Line 205 | Line 228 | setheadvar(char *val, void *p)
228          int     kln, vln;
229          int     n;
230  
231 +        adv_linecnt(htp);       /* side-effect is to count lines */
232          if (!isalpha(*val))     /* key must start line */
233                  return(0);
234          key = val++;
# Line 220 | Line 244 | setheadvar(char *val, void *p)
244                  val++;
245          if (!*val)              /* nothing set? */
246                  return(0);
223 –        vln = strlen(val);      /* eliminate space and newline at end */
224 –        while (--vln > 0 && isspace(val[vln]))
225 –                ;
226 –        val[++vln] = '\0';
247                                  /* check if key to ignore */
248          for (n = 0; hdr_ignkey[n]; n++)
249                  if (!strcmp(key, hdr_ignkey[n]))
250                          return(0);
251 +        vln = strlen(val);      /* eliminate space and newline at end */
252 +        while (isspace(val[--vln]))
253 +                ;
254 +        val[++vln] = '\0';
255          if (!(tep = lu_find(htp, key)))
256 <                return(-1);
256 >                return(-1);     /* memory allocation error */
257          if (!tep->key)
258                  tep->key = strcpy(malloc(kln+1), key);
259          if (tep->data)
260                  free(tep->data);
261          tep->data = strcpy(malloc(vln+1), val);
238 –        adv_linecnt(htp);
262          return(1);
263   }
264  
# Line 246 | Line 269 | match_val(const LUENT *ep1, void *p2)
269          const LUENT     *ep2 = lu_find((LUTAB *)p2, ep1->key);
270          if (!ep2 || !ep2->data) {
271                  if (report != REP_QUIET)
272 <                        fprintf(stderr, "%s: Variable '%s' missing in '%s'\n",
272 >                        printf("%s: variable '%s' missing in '%s'\n",
273                                          progname, ep1->key, f2name);
274                  return(-1);
275          }
276          if (!equiv_string((char *)ep1->data, (char *)ep2->data)) {
277 <                if (report != REP_QUIET)
278 <                        fprintf(stderr, "%s: Header variables '%s' have different values\n",
277 >                if (report != REP_QUIET) {
278 >                        printf("%s: header variable '%s' has different values\n",
279                                          progname, ep1->key);
280 +                        if (report >= REP_VERBOSE) {
281 +                                printf("%s: %s=%s\n", f1name,
282 +                                                ep1->key, (char *)ep1->data);
283 +                                printf("%s: %s=%s\n", f2name,
284 +                                                ep2->key, (char *)ep2->data);
285 +                        }
286 +                }
287                  return(-1);
288          }
289          return(1);              /* good match */
# Line 266 | Line 296 | headers_match(LUTAB *hp1, LUTAB *hp2)
296          int     ne = lu_doall(hp1, match_val, hp2);
297          if (ne < 0)
298                  return(0);      /* something didn't match! */
299 <        ne = lu_doall(hp2, NULL, NULL) - ne;
300 <        if (ne) {
301 <                if (report != REP_QUIET)
272 <                        fprintf(stderr, "%s: '%s' has %d extra header variable(s)\n",
299 >                                /* non-fatal if second header has extra */
300 >        if (report >= REP_WARN && (ne = lu_doall(hp2, NULL, NULL) - ne))
301 >                printf("%s: warning - '%s' has %d extra header setting(s)\n",
302                                          progname, f2name, ne);
274 –                return(0);
275 –        }
303          return(1);              /* good match */
304   }
305  
# Line 284 | Line 311 | input_is_binary(FILE *fin)
311          int     c = 0;
312  
313          while ((c = getc(fin)) != EOF) {
314 +                ++n;
315                  if (!c | (c > 127))
316                          break;                  /* non-ascii character */
317 <                if (++n >= 10240)
317 >                if (n >= 10240)
318                          break;                  /* enough to be confident */
319          }
320          if (!n)
# Line 316 | Line 344 | identify_type(const char *name, FILE *fin, LUTAB *htp)
344                  char    sbuf[32];
345                  if (!fgets(sbuf, sizeof(sbuf), fin))
346                          goto badeof;
347 +                adv_linecnt(htp);               /* for #?ID string */
348                  if (report >= REP_WARN && strncmp(sbuf, "RADIANCE", 8)) {
349 <                        fputs(name, stderr);
350 <                        fputs(": warning - unexpected header ID: ", stderr);
351 <                        fputs(sbuf, stderr);
349 >                        fputs(name, stdout);
350 >                        fputs(": warning - unexpected header ID: ", stdout);
351 >                        fputs(sbuf, stdout);
352                  }
324 –                adv_linecnt(htp);               /* for #?ID string */
353                  if (getheader(fin, setheadvar, htp) < 0) {
354                          fputs(name, stderr);
355 <                        fputs(": Unknown error reading header\n", stderr);
355 >                        fputs(": unknown error reading header\n", stderr);
356                          return(-1);
357                  }
358                  adv_linecnt(htp);               /* for trailing emtpy line */
# Line 342 | Line 370 | identify_type(const char *name, FILE *fin, LUTAB *htp)
370          return(TYP_TEXT);
371   badeof:
372          if (report != REP_QUIET) {
373 <                fputs(name, stderr);
374 <                fputs(": Unexpected end-of-file\n", stderr);
373 >                fputs(name, stdout);
374 >                fputs(": unexpected end-of-file\n", stdout);
375          }
376          return(-1);
377   }
# Line 356 | Line 384 | good_RMS()
384                  return(1);
385          if (diff2sum/(double)nsum > rms_lim*rms_lim) {
386                  if (report != REP_QUIET)
387 <                        fprintf(stderr,
388 <        "%s: %sRMS difference between '%s' and '%s' of %.5g exceeds limit\n",
387 >                        printf(
388 >        "%s: %sRMS difference between '%s' and '%s' of %.5g exceeds limit of %.5g\n",
389                                          progname,
390                                          (rel_min > 0) ? "relative " : "",
391                                          f1name, f2name,
392 <                                        sqrt(diff2sum/(double)nsum));
392 >                                        sqrt(diff2sum/(double)nsum), rms_lim);
393                  return(0);
394          }
395          if (report >= REP_VERBOSE)
396 <                fprintf(stderr,
369 <                        "%s: %sRMS difference of reals in '%s' and '%s' is %.5g\n",
396 >                printf("%s: %sRMS difference of reals in '%s' and '%s' is %.5g\n",
397                                  progname, (rel_min > 0) ? "relative " : "",
398                                  f1name, f2name, sqrt(diff2sum/(double)nsum));
399          return(1);
# Line 376 | Line 403 | good_RMS()
403   static int
404   compare_binary()
405   {
406 +        int     c1=0, c2=0;
407 +
408          if (report >= REP_VERBOSE) {
409 <                fputs(progname, stderr);
410 <                fputs(": comparing inputs as binary\n", stderr);
409 >                fputs(progname, stdout);
410 >                fputs(": comparing inputs as binary\n", stdout);
411          }
412          for ( ; ; ) {                           /* exact byte matching */
413 <                int     c1 = getc(f1in);
414 <                int     c2 = getc(f2in);
413 >                c1 = getc(f1in);
414 >                c2 = getc(f2in);
415                  if (c1 == EOF) {
416                          if (c2 == EOF)
417                                  return(1);      /* success! */
418                          if (report != REP_QUIET) {
419 <                                fputs(f1name, stderr);
420 <                                fputs(": Unexpected end-of-file\n", stderr);
419 >                                fputs(f1name, stdout);
420 >                                fputs(": unexpected end-of-file\n", stdout);
421                          }
422                          return(0);
423                  }
424                  if (c2 == EOF) {
425                          if (report != REP_QUIET) {
426 <                                fputs(f2name, stderr);
427 <                                fputs(": Unexpected end-of-file\n", stderr);
426 >                                fputs(f2name, stdout);
427 >                                fputs(": unexpected end-of-file\n", stdout);
428                          }
429                          return(0);
430                  }
# Line 404 | Line 433 | compare_binary()
433          }
434          if (report == REP_QUIET)
435                  return(0);
436 <        fprintf(stderr, "%s: binary files '%s' and '%s' differ at offset %ld|%ld\n",
436 >        printf("%s: binary files '%s' and '%s' differ at byte offset %ld|%ld\n",
437                          progname, f1name, f2name, ftell(f1in), ftell(f2in));
438 +        if (report >= REP_VERBOSE)
439 +                printf("%s: byte in '%s' is 0x%X, byte in '%s' is 0x%X\n",
440 +                                progname, f1name, c1, f2name, c2);
441          return(0);
442   }
443  
# Line 416 | Line 448 | compare_text()
448          char    l1buf[4096], l2buf[4096];
449  
450          if (report >= REP_VERBOSE) {
451 <                fputs(progname, stderr);
452 <                fputs(": comparing inputs as ASCII text\n", stderr);
451 >                fputs(progname, stdout);
452 >                fputs(": comparing inputs as ASCII text\n", stdout);
453          }
454                                                  /* compare a line at a time */
455          while (fgets(l1buf, sizeof(l1buf), f1in)) {
# Line 431 | Line 463 | compare_text()
463                  }
464                  if (feof(f2in)) {
465                          if (report != REP_QUIET) {
466 <                                fputs(f2name, stderr);
467 <                                fputs(": Unexpected end-of-file\n", stderr);
466 >                                fputs(f2name, stdout);
467 >                                fputs(": unexpected end-of-file\n", stdout);
468                          }
469                          return(0);
470                  }
471                                                  /* compare non-empty lines */
472                  if (!equiv_string(l1buf, l2buf)) {
473                          if (report != REP_QUIET) {
474 <                                fprintf(stderr, "%s: inputs '%s' and '%s' differ at line %d|%d\n",
474 >                                printf("%s: inputs '%s' and '%s' differ at line %d|%d\n",
475                                                  progname, f1name, f2name,
476                                                  lin1cnt, lin2cnt);
477                                  if (report >= REP_VERBOSE) {
478 <                                        fputs("------------- Mismatch -------------\n", stderr);
479 <                                        fprintf(stderr, "%s(%d):\t%s", f1name,
478 >                                        fputs("------------- Mismatch -------------\n", stdout);
479 >                                        printf("%s@%d:\t%s", f1name,
480                                                          lin1cnt, l1buf);
481 <                                        fprintf(stderr, "%s(%d):\t%s", f2name,
481 >                                        printf("%s@%d:\t%s", f2name,
482                                                          lin2cnt, l2buf);
483                                  }
484                          }
# Line 458 | Line 490 | compare_text()
490                  if (!*sskip2(l2buf,0))
491                          continue;
492                  if (report != REP_QUIET) {
493 <                        fputs(f1name, stderr);
494 <                        fputs(": Unexpected end-of-file\n", stderr);
493 >                        fputs(f1name, stdout);
494 >                        fputs(": unexpected end-of-file\n", stdout);
495                  }
496                  return(0);
497          }
# Line 474 | Line 506 | compare_hdr()
506          COLOR   *scan1, *scan2;
507          int     x, y;
508  
509 +        if (report >= REP_VERBOSE) {
510 +                fputs(progname, stdout);
511 +                fputs(": comparing inputs as HDR images\n", stdout);
512 +        }
513          fgetsresolu(&rs1, f1in);
514          fgetsresolu(&rs2, f2in);
515          if (rs1.rt != rs2.rt) {
516                  if (report != REP_QUIET)
517 <                        fprintf(stderr,
517 >                        printf(
518                          "%s: Images '%s' and '%s' have different pixel ordering\n",
519                                          progname, f1name, f2name);
520                  return(0);
521          }
522          if ((rs1.xr != rs2.xr) | (rs1.yr != rs2.yr)) {
523                  if (report != REP_QUIET)
524 <                        fprintf(stderr,
524 >                        printf(
525                          "%s: Images '%s' and '%s' are different sizes\n",
526                                          progname, f1name, f2name);
527                  return(0);
# Line 500 | Line 536 | compare_hdr()
536                  if ((freadscan(scan1, scanlen(&rs1), f1in) < 0) |
537                                  (freadscan(scan2, scanlen(&rs2), f2in) < 0)) {
538                          if (report != REP_QUIET)
539 <                                fprintf(stderr, "%s: Unexpected end-of-file\n",
539 >                                printf("%s: unexpected end-of-file\n",
540                                                  progname);
541                          free(scan1);
542                          free(scan2);
543                          return(0);
544                  }
545 <                for (x = scanlen(&rs1); x--; ) {
546 <                        if (real_check(colval(scan1[x],RED),
511 <                                                colval(scan2[x],RED)) &
512 <                                        real_check(colval(scan1[x],GRN),
513 <                                                colval(scan2[x],GRN)) &
514 <                                        real_check(colval(scan1[x],BLU),
515 <                                                colval(scan2[x],BLU)))
545 >                for (x = 0; x < scanlen(&rs1); x++) {
546 >                        if (color_check(scan1[x], scan2[x]))
547                                  continue;
548                          if (report != REP_QUIET) {
549 <                                fprintf(stderr,
549 >                                printf(
550                                  "%s: pixels at scanline %d offset %d differ\n",
551                                          progname, y, x);
552                                  if (report >= REP_VERBOSE) {
553 <                                        fprintf(stderr, "%s: (R,G,B)=(%g,%g,%g)\n",
553 >                                        printf("%s: (R,G,B)=(%g,%g,%g)\n",
554                                                  f1name, colval(scan1[x],RED),
555                                                  colval(scan1[x],GRN),
556                                                  colval(scan1[x],BLU));
557 <                                        fprintf(stderr, "%s: (R,G,B)=(%g,%g,%g)\n",
557 >                                        printf("%s: (R,G,B)=(%g,%g,%g)\n",
558                                                  f2name, colval(scan2[x],RED),
559                                                  colval(scan2[x],GRN),
560                                                  colval(scan2[x],BLU));
# Line 539 | Line 570 | compare_hdr()
570          return(good_RMS());                     /* final check of RMS */
571   }
572  
542 – const char      nsuffix[10][3] = {              /* 1st, 2nd, 3rd, etc. */
543 –                        "th","st","nd","rd","th","th","th","th","th","th"
544 –                };
545 –
573   /* Compare two inputs that are known to be 32-bit floating-point data */
574   static int
575   compare_float()
# Line 550 | Line 577 | compare_float()
577          long    nread = 0;
578          float   f1, f2;
579  
580 +        if (report >= REP_VERBOSE) {
581 +                fputs(progname, stdout);
582 +                fputs(": comparing inputs as 32-bit IEEE floats\n", stdout);
583 +        }
584          while (getbinary(&f1, sizeof(f1), 1, f1in)) {
585                  if (!getbinary(&f2, sizeof(f2), 1, f2in))
586                          goto badeof;
# Line 557 | Line 588 | compare_float()
588                  if (real_check(f1, f2))
589                          continue;
590                  if (report != REP_QUIET)
591 <                        fprintf(stderr, "%s: %ld%s float values differ\n",
592 <                                        progname, nread, nsuffix[nread%10]);
591 >                        printf("%s: %ld%s float values differ\n",
592 >                                        progname, nread, num_sfx(nread));
593                  return(0);
594          }
595          if (!getbinary(&f2, sizeof(f2), 1, f2in))
596                  return(good_RMS());             /* final check of RMS */
597   badeof:
598          if (report != REP_QUIET)
599 <                fprintf(stderr, "%s: Unexpected end-of-file\n", progname);
599 >                printf("%s: unexpected end-of-file\n", progname);
600          return(0);
601   }
602  
# Line 576 | Line 607 | compare_double()
607          long    nread = 0;
608          double  f1, f2;
609  
610 +        if (report >= REP_VERBOSE) {
611 +                fputs(progname, stdout);
612 +                fputs(": comparing inputs as 64-bit IEEE doubles\n", stdout);
613 +        }
614          while (getbinary(&f1, sizeof(f1), 1, f1in)) {
615                  if (!getbinary(&f2, sizeof(f2), 1, f2in))
616                          goto badeof;
# Line 583 | Line 618 | compare_double()
618                  if (real_check(f1, f2))
619                          continue;
620                  if (report != REP_QUIET)
621 <                        fprintf(stderr, "%s: %ld%s float values differ\n",
622 <                                        progname, nread, nsuffix[nread%10]);
621 >                        printf("%s: %ld%s double values differ\n",
622 >                                        progname, nread, num_sfx(nread));
623                  return(0);
624          }
625          if (!getbinary(&f2, sizeof(f2), 1, f2in))
626                  return(good_RMS());             /* final check of RMS */
627   badeof:
628          if (report != REP_QUIET)
629 <                fprintf(stderr, "%s: Unexpected end-of-file\n", progname);
629 >                printf("%s: unexpected end-of-file\n", progname);
630          return(0);
631   }
632  
# Line 638 | Line 673 | main(int argc, char *argv[])
673                  }
674                  break;
675          }
676 <        if (a != argc-2)
676 >        if (a != argc-2)                        /* make sure of two inputs */
677                  usage();
678 <        if (!f1name) f1name = argv[a];
644 <        if (!f2name) f2name = argv[a+1];
645 <
646 <        if (!strcmp(f1name, f2name)) {          /* inputs are same? */
678 >        if (!strcmp(argv[a], argv[a+1])) {      /* inputs are same? */
679                  if (report >= REP_WARN)
680 <                        fprintf(stderr, "%s: warning - identical inputs given\n",
680 >                        printf("%s: warning - identical inputs given\n",
681                                          progname);
682                  return(0);
683          }
684 +        if (!f1name) f1name = argv[a];
685 +        if (!f2name) f2name = argv[a+1];
686                                                  /* open inputs */
687          SET_FILE_BINARY(stdin);                 /* in case we're using it */
688          if (!f1in && !(f1in = fopen(f1name, "rb"))) {
689                  fprintf(stderr, "%s: cannot open for reading\n", f1name);
690 <                return(1);
690 >                return(2);
691          }
692          if (!strcmp(f2name, "-")) {
693                  f2in = stdin;
694                  f2name = stdin_name;
695          } else if (!(f2in = fopen(f2name, "rb"))) {
696                  fprintf(stderr, "%s: cannot open for reading\n", f2name);
697 <                return(1);
697 >                return(2);
698          }
699                                                  /* load headers */
700          if ((typ1 = identify_type(f1name, f1in, &hdr1)) < 0)
701 <                return(1);
701 >                return(2);
702          if ((typ2 = identify_type(f2name, f2in, &hdr2)) < 0)
703 <                return(1);
703 >                return(2);
704          if (typ1 != typ2) {
705                  if (report != REP_QUIET)
706 <                        fprintf(stderr, "%s: '%s' is %s and '%s' is %s\n",
706 >                        printf("%s: '%s' is %s and '%s' is %s\n",
707                                          progname, f1name, file_type[typ1],
708                                          f2name, file_type[typ2]);
709                  return(1);
# Line 680 | Line 714 | main(int argc, char *argv[])
714          lu_done(&hdr1); lu_done(&hdr2);
715          if (!ign_header & (report >= REP_WARN)) {
716                  if (typ1 == TYP_UNKNOWN)
717 <                        fprintf(stderr, "%s: warning - unrecognized format, comparing as binary\n",
717 >                        printf("%s: warning - unrecognized format, comparing as binary\n",
718                                          progname);
719                  if (lin1cnt != lin2cnt)
720 <                        fprintf(stderr, "%s: warning - headers are different lengths\n",
720 >                        printf("%s: warning - headers are different lengths\n",
721                                          progname);
722          }
723          if (report >= REP_VERBOSE)
724 <                fprintf(stderr, "%s: input file type is %s\n",
724 >                printf("%s: input file type is %s\n",
725                                  progname, file_type[typ1]);
726  
727          switch (typ1) {                         /* compare based on type */

Diff Legend

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