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.34 by greg, Tue Nov 21 20:01:39 2023 UTC

# Line 8 | Line 8 | static const char RCSid[] = "$Id$";
8   */
9  
10   #include <stdlib.h>
11 – #include <math.h>
11   #include <ctype.h>
12 + #include "rtmath.h"
13   #include "platform.h"
14   #include "rtio.h"
15   #include "resolu.h"
16   #include "color.h"
17 + #include "depthcodec.h"
18 + #include "normcodec.h"
19   #include "lookup.h"
20                          /* Reporting levels */
21   #define REP_QUIET       0       /* no reporting */
# Line 25 | Line 27 | int    report = REP_WARN;      /* reporting level */
27  
28   int     ign_header = 0;         /* ignore header differences? */
29  
30 + int     escape_newlines = 0;    /* allow backslash to skip newlines */
31 +
32   double  rel_min = 1e-5;         /* positive for relative comparisons */
33  
34   double  rms_lim = 0.01;         /* RMS difference limit */
# Line 33 | Line 37 | double max_lim = 0.25;         /* difference limit if non-neg
37  
38   int     lin1cnt=0, lin2cnt=0;   /* file line position */
39  
40 + int     comment_c = '\0';       /* comment delimiter for text files */
41 +
42 + const char      nsuffix[10][3] = {              /* 1st, 2nd, 3rd, etc. */
43 +                        "th","st","nd","rd","th","th","th","th","th","th"
44 +                };
45 + #define num_sfx(n)      nsuffix[(n)%10]
46 +
47                                  /* file types */
48   const char      *file_type[] = {
49                          "Unrecognized",
# Line 40 | Line 51 | const char     *file_type[] = {
51                          "ascii",
52                          COLRFMT,
53                          CIEFMT,
54 <                        "BINARY_float",
55 <                        "BINARY_double",
54 >                        SPECFMT,
55 >                        DEPTH16FMT,
56 >                        NORMAL32FMT,
57 >                        "float",
58 >                        "double",
59                          "BSDF_RBFmesh",
60                          "Radiance_octree",
61                          "Radiance_tmesh",
62 +                        "8-bit_indexed_name",
63 +                        "16-bit_indexed_name",
64 +                        "24-bit_indexed_name",
65                          "BINARY_unknown",
66                          NULL    /* terminator */
67                  };
68 <
69 < enum {TYP_UNKNOWN, TYP_TEXT, TYP_ASCII, TYP_RGBE, TYP_XYZE, TYP_FLOAT,
70 <                TYP_DOUBLE, TYP_RBFMESH, TYP_OCTREE, TYP_TMESH, TYP_BINARY};
68 >                                /* keep consistent with above */
69 > enum {TYP_UNKNOWN, TYP_TEXT, TYP_ASCII, TYP_RGBE, TYP_XYZE, TYP_SPEC,
70 >                TYP_DEPTH, TYP_NORM, TYP_FLOAT, TYP_DOUBLE,
71 >                TYP_RBFMESH, TYP_OCTREE, TYP_TMESH,
72 >                TYP_ID8, TYP_ID16, TYP_ID24, TYP_BINARY};
73                  
74   #define has_header(t)   (!( 1L<<(t) & (1L<<TYP_TEXT | 1L<<TYP_BINARY) ))
75  
# Line 59 | Line 78 | const char     *hdr_ignkey[] = {
78                          "SOFTWARE",
79                          "CAPDATE",
80                          "GMT",
81 +                        "FRAME",
82                          NULL    /* terminator */
83                  };
84                                  /* header variable settings */
# Line 69 | Line 89 | LUTAB  hdr2 = LU_SINIT(free,free);
89   #define adv_linecnt(htp)        (lin1cnt += (htp == &hdr1), \
90                                          lin2cnt += (htp == &hdr2))
91  
92 + typedef struct {                /* dynamic line buffer */
93 +        char    *str;
94 +        int     len;
95 +        int     siz;
96 + } LINEBUF;
97 +
98 + #define init_line(bp)   ((bp)->str = NULL, (bp)->siz = 0)
99 +                                /* 100 MByte limit on line buffer */
100 + #define MAXBUF          (100L<<20)
101 +
102                                  /* input files */
103   char            *progname = NULL;
104   const char      stdin_name[] = "<stdin>";
105   const char      *f1name=NULL, *f2name=NULL;
106   FILE            *f1in=NULL, *f2in=NULL;
107 + int             f1swap=0, f2swap=0;
108 +
109                                  /* running real differences */
110   double  diff2sum = 0;
111 < int     nsum = 0;
111 > long    nsum = 0;
112  
113   /* Report usage and exit */
114   static void
# Line 84 | Line 116 | usage()
116   {
117          fputs("Usage: ", stderr);
118          fputs(progname, stderr);
119 <        fputs(" [-h][-s|-w|-v][-rel min_test][-rms epsilon][-max epsilon] file1 file2\n",
119 >        fputs(" [-h][-c#][-s|-w|-v][-rel min_test][-rms epsilon][-max epsilon] reference test\n",
120                          stderr);
121 <        exit(1);
121 >        exit(2);
122   }
123  
124 + /* Read a text line, increasing buffer size as necessary */
125 + static int
126 + read_line(LINEBUF *bp, FILE *fp)
127 + {
128 +        static int      doneWarn = 0;
129 +
130 +        bp->len = 0;
131 +        if (!bp->str) {
132 +                bp->str = (char *)malloc(bp->siz = 512);
133 +                if (!bp->str)
134 +                        goto memerr;
135 +        }
136 +        while (fgets(bp->str + bp->len, bp->siz - bp->len, fp)) {
137 +                bp->len += strlen(bp->str + bp->len);
138 +                if (bp->str[bp->len-1] == '\n') {
139 +                        if (bp->len > 1 && bp->str[bp->len-2] == '\r') {
140 +                                bp->str[--bp->len] = '\0';
141 +                                bp->str[bp->len-1] = '\n';
142 +                        }
143 +                        if (escape_newlines && bp->len > 1 &&
144 +                                        bp->str[bp->len-2] == '\\') {
145 +                                bp->str[--bp->len] = '\0';
146 +                                bp->str[bp->len-1] = ' ';
147 +                                continue;
148 +                        }
149 +                        break;          /* found EOL */
150 +                }
151 +                if (bp->len < bp->siz - 4)
152 +                        continue;       /* at EOF? */
153 +                if (bp->siz >= MAXBUF) {
154 +                        if ((report >= REP_WARN) & !doneWarn) {
155 +                                fprintf(stderr,
156 +                        "%s: warning - input line(s) past %ld MByte limit\n",
157 +                                        progname, MAXBUF>>20);
158 +                                doneWarn++;
159 +                        }
160 +                        break;          /* return MAXBUF partial line */
161 +                }
162 +                if ((bp->siz += bp->siz/2) > MAXBUF)
163 +                        bp->siz = MAXBUF;
164 +                bp->str = (char *)realloc(bp->str, bp->siz);
165 +                if (!bp->str)
166 +                        goto memerr;
167 +        }
168 +        if (comment_c) {                /* elide comment? */
169 +                char    *cp = sskip2(bp->str,0);
170 +                if (*cp == comment_c) {
171 +                        *cp++ = '\n';
172 +                        *cp = '\0';
173 +                        bp->len = cp - bp->str;
174 +                }
175 +        }
176 +        return(bp->len);
177 + memerr:
178 +        fprintf(stderr,
179 +                "%s: out of memory in read_line() allocating %d byte buffer\n",
180 +                        progname, bp->siz);
181 +        exit(2);
182 + }
183 +
184 + /* Free line buffer */
185 + static void
186 + free_line(LINEBUF *bp)
187 + {
188 +        if (bp->str) free(bp->str);
189 +        init_line(bp);
190 + }
191 +
192   /* Get type ID from name (or 0 if not found) */
193   static int
194   xlate_type(const char *nm)
# Line 111 | Line 211 | real_check(double r1, double r2)
211  
212          if (rel_min > 0) {      /* doing relative differences? */
213                  double  av2 = .25*(r1*r1 + 2.*fabs(r1*r2) + r2*r2);
214 <                if (av2 > rel_min*rel_min)
215 <                        diff2 /= av2;
214 >                if (av2 < rel_min*rel_min)
215 >                        av2 = rel_min*rel_min;
216 >                diff2 /= av2;
217          }
218          if (max_lim >= 0 && diff2 > max_lim*max_lim) {
219                  if (report != REP_QUIET)
220 <                        fprintf(stderr,
221 <                        "%s: %sdifference between %.8g and %.8g exceeds epsilon\n",
220 >                        printf(
221 >                        "%s: %sdifference between %.8g and %.8g exceeds epsilon of %.8g\n",
222                                          progname,
223                                          (rel_min > 0) ? "relative " : "",
224 <                                        r1, r2);
224 >                                        r1, r2, max_lim);
225                  return(0);
226          }
227          diff2sum += diff2;
# Line 128 | Line 229 | real_check(double r1, double r2)
229          return(1);
230   }
231  
232 + /* Compare two color values for equivalence */
233 + static int
234 + color_check(COLOR c1, COLOR c2)
235 + {
236 +        int     p;
237 +
238 +        if (!real_check((colval(c1,RED)+colval(c1,GRN)+colval(c1,BLU))*(1./3.),
239 +                        (colval(c2,RED)+colval(c2,GRN)+colval(c2,BLU))*(1./3.)))
240 +                return(0);
241 +
242 +        p = (colval(c1,GRN) > colval(c1,RED)) ? GRN : RED;
243 +        if (colval(c1,BLU) > colval(c1,p)) p = BLU;
244 +        
245 +        return(real_check(colval(c1,p), colval(c2,p)));
246 + }
247 +
248 + #if 1
249 + /* Compare two color spectra for equivalence */
250 + static int
251 + spec_check(COLORV *sc1, COLORV *sc2)
252 + {
253 +        int     p, k;
254 +
255 +        if (!real_check(scolor_mean(sc1), scolor_mean(sc2)))
256 +                return(0);
257 +
258 +        p = 0;                          /* find max. component */
259 +        for (k = NCSAMP; --k; )
260 +                if (sc1[k] > sc1[p])
261 +                        p = k;
262 +
263 +        return(real_check(sc1[p], sc2[p]));
264 + }
265 + #else
266 + /* Compare two color spectra for equivalence */
267 + static int
268 + spec_check(COLORV *sc1, COLORV *sc2)
269 + {
270 +        COLOR   c1, c2;
271 +        int     p;
272 +
273 +        if (!real_check(scolor_mean(sc1), scolor_mean(sc2)))
274 +                return(0);
275 +                                        /* do comparisons in RGB space */
276 +        scolor_rgb(c1, sc1);
277 +        scolor_rgb(c2, sc2);
278 +
279 +        p = (colval(c1,GRN) > colval(c1,RED)) ? GRN : RED;
280 +        if (colval(c1,BLU) > colval(c1,p)) p = BLU;
281 +
282 +        return(real_check(colval(c1,p), colval(c2,p)));
283 + }
284 + #endif
285 +
286 + /* Compare two normal directions for equivalence */
287 + static int
288 + norm_check(FVECT nv1, FVECT nv2)
289 + {
290 +        double  max2 = nv1[2]*nv2[2];
291 +        int     imax = 2;
292 +        int     i = 2;
293 +                                        /* identify largest component */
294 +        while (i--) {
295 +                double  tm2 = nv1[i]*nv2[i];
296 +                if (tm2 > max2) {
297 +                        imax = i;
298 +                        max2 = tm2;
299 +                }
300 +        }
301 +        i = 3;                          /* compare smaller components */
302 +        while (i--) {
303 +                if (i == imax)
304 +                        continue;
305 +                if (!real_check(nv1[i], nv2[i]))
306 +                        return(0);
307 +        }
308 +        return(1);
309 + }
310 +
311   /* Compare two strings for equivalence */
312   static int
313   equiv_string(char *s1, char *s2)
# Line 138 | Line 318 | equiv_string(char *s1, char *s2)
318                                  /* skip whitespace at beginning */
319          while (isspace(*s1)) s1++;
320          while (isspace(*s2)) s2++;
141 –        if (!strcmp(s1, s2))    /* quick check */
142 –                return(1);
321          while (*s1) {           /* check each word */
322 <                int     inquote = *s1;
322 >                int     inquote;
323 >                if (!*s2)       /* unexpected EOL in s2? */
324 >                        return(0);
325 >                inquote = *s1;
326                  if ((inquote != '\'') & (inquote != '"'))
327                          inquote = 0;
328                  if (inquote) {  /* quoted text must match exactly */
# Line 199 | Line 380 | equiv_string(char *s1, char *s2)
380   static int
381   setheadvar(char *val, void *p)
382   {
383 +        char    newval[128];
384          LUTAB   *htp = (LUTAB *)p;
385          LUENT   *tep;
386          char    *key;
387          int     kln, vln;
388          int     n;
389  
390 +        adv_linecnt(htp);       /* side-effect is to count lines */
391          if (!isalpha(*val))     /* key must start line */
392                  return(0);
393 +                                /* check if we need to swap binary data */
394 +        if ((n = isbigendian(val)) >= 0) {
395 +                if (nativebigendian() == n)
396 +                        return(0);
397 +                f1swap += (htp == &hdr1);
398 +                f2swap += (htp == &hdr2);
399 +                return(0);
400 +        }
401          key = val++;
402          while (*val && !isspace(*val) & (*val != '='))
403                  val++;
# Line 220 | Line 411 | setheadvar(char *val, void *p)
411                  val++;
412          if (!*val)              /* nothing set? */
413                  return(0);
223 –        vln = strlen(val);      /* eliminate space and newline at end */
224 –        while (--vln > 0 && isspace(val[vln]))
225 –                ;
226 –        val[++vln] = '\0';
414                                  /* check if key to ignore */
415          for (n = 0; hdr_ignkey[n]; n++)
416                  if (!strcmp(key, hdr_ignkey[n]))
417                          return(0);
418 +        vln = strlen(val);      /* eliminate space and newline at end */
419 +        while (isspace(val[--vln]))
420 +                ;
421 +        val[++vln] = '\0';
422          if (!(tep = lu_find(htp, key)))
423 <                return(-1);
423 >                return(-1);     /* memory allocation error */
424          if (!tep->key)
425                  tep->key = strcpy(malloc(kln+1), key);
426 <        if (tep->data)
426 >        if (tep->data) {        /* check for special cases */
427 >                if (!strcmp(key, "EXPOSURE")) {
428 >                        sprintf(newval, "%.6e", atof(tep->data)*atof(val));
429 >                        vln = strlen(val = newval);
430 >                }
431                  free(tep->data);
432 +        }
433          tep->data = strcpy(malloc(vln+1), val);
238 –        adv_linecnt(htp);
434          return(1);
435   }
436  
# Line 246 | Line 441 | match_val(const LUENT *ep1, void *p2)
441          const LUENT     *ep2 = lu_find((LUTAB *)p2, ep1->key);
442          if (!ep2 || !ep2->data) {
443                  if (report != REP_QUIET)
444 <                        fprintf(stderr, "%s: Variable '%s' missing in '%s'\n",
444 >                        printf("%s: variable '%s' missing in '%s'\n",
445                                          progname, ep1->key, f2name);
446                  return(-1);
447          }
448          if (!equiv_string((char *)ep1->data, (char *)ep2->data)) {
449 <                if (report != REP_QUIET)
450 <                        fprintf(stderr, "%s: Header variables '%s' have different values\n",
449 >                if (report != REP_QUIET) {
450 >                        printf("%s: header variable '%s' has different values\n",
451                                          progname, ep1->key);
452 +                        if (report >= REP_VERBOSE) {
453 +                                printf("%s: %s=%s\n", f1name,
454 +                                                ep1->key, (char *)ep1->data);
455 +                                printf("%s: %s=%s\n", f2name,
456 +                                                ep2->key, (char *)ep2->data);
457 +                        }
458 +                }
459                  return(-1);
460          }
461          return(1);              /* good match */
# Line 261 | Line 463 | match_val(const LUENT *ep1, void *p2)
463  
464   /* Compare two sets of header variables */
465   static int
466 < headers_match(LUTAB *hp1, LUTAB *hp2)
466 > headers_match()
467   {
468 <        int     ne = lu_doall(hp1, match_val, hp2);
468 >        int     ne = lu_doall(&hdr1, match_val, &hdr2);
469          if (ne < 0)
470                  return(0);      /* something didn't match! */
471 <        ne = lu_doall(hp2, NULL, NULL) - ne;
472 <        if (ne) {
473 <                if (report != REP_QUIET)
272 <                        fprintf(stderr, "%s: '%s' has %d extra header variable(s)\n",
471 >                                /* non-fatal if second header has extra */
472 >        if (report >= REP_WARN && (ne = lu_doall(&hdr2, NULL, NULL) - ne))
473 >                printf("%s: warning - '%s' has %d extra header setting(s)\n",
474                                          progname, f2name, ne);
274 –                return(0);
275 –        }
475          return(1);              /* good match */
476   }
477  
# Line 284 | Line 483 | input_is_binary(FILE *fin)
483          int     c = 0;
484  
485          while ((c = getc(fin)) != EOF) {
486 +                ++n;
487                  if (!c | (c > 127))
488                          break;                  /* non-ascii character */
489 <                if (++n >= 10240)
489 >                if (n >= 10240)
490                          break;                  /* enough to be confident */
491          }
492          if (!n)
# Line 316 | Line 516 | identify_type(const char *name, FILE *fin, LUTAB *htp)
516                  char    sbuf[32];
517                  if (!fgets(sbuf, sizeof(sbuf), fin))
518                          goto badeof;
519 +                adv_linecnt(htp);               /* for #?ID string */
520                  if (report >= REP_WARN && strncmp(sbuf, "RADIANCE", 8)) {
521 <                        fputs(name, stderr);
522 <                        fputs(": warning - unexpected header ID: ", stderr);
523 <                        fputs(sbuf, stderr);
521 >                        fputs(name, stdout);
522 >                        fputs(": warning - unexpected header ID: ", stdout);
523 >                        fputs(sbuf, stdout);
524                  }
324 –                adv_linecnt(htp);               /* for #?ID string */
525                  if (getheader(fin, setheadvar, htp) < 0) {
526                          fputs(name, stderr);
527 <                        fputs(": Unknown error reading header\n", stderr);
527 >                        fputs(": unknown error reading header\n", stderr);
528                          return(-1);
529                  }
530                  adv_linecnt(htp);               /* for trailing emtpy line */
# Line 338 | Line 538 | identify_type(const char *name, FILE *fin, LUTAB *htp)
538          }
539          if (c)
540                  return(TYP_BINARY);
341 –        SET_FILE_TEXT(fin);                     /* originally set to binary */
541          return(TYP_TEXT);
542   badeof:
543          if (report != REP_QUIET) {
544 <                fputs(name, stderr);
545 <                fputs(": Unexpected end-of-file\n", stderr);
544 >                fputs(name, stdout);
545 >                fputs(": unexpected end-of-file\n", stdout);
546          }
547          return(-1);
548   }
# Line 356 | Line 555 | good_RMS()
555                  return(1);
556          if (diff2sum/(double)nsum > rms_lim*rms_lim) {
557                  if (report != REP_QUIET)
558 <                        fprintf(stderr,
559 <        "%s: %sRMS difference between '%s' and '%s' of %.5g exceeds limit\n",
558 >                        printf(
559 >        "%s: %sRMS difference between '%s' and '%s' of %.5g exceeds limit of %.5g\n",
560                                          progname,
561                                          (rel_min > 0) ? "relative " : "",
562                                          f1name, f2name,
563 <                                        sqrt(diff2sum/(double)nsum));
563 >                                        sqrt(diff2sum/(double)nsum), rms_lim);
564                  return(0);
565          }
566          if (report >= REP_VERBOSE)
567 <                fprintf(stderr,
369 <                        "%s: %sRMS difference of reals in '%s' and '%s' is %.5g\n",
567 >                printf("%s: %sRMS difference of reals in '%s' and '%s' is %.5g\n",
568                                  progname, (rel_min > 0) ? "relative " : "",
569                                  f1name, f2name, sqrt(diff2sum/(double)nsum));
570          return(1);
# Line 376 | Line 574 | good_RMS()
574   static int
575   compare_binary()
576   {
577 +        int     c1=0, c2=0;
578 +
579          if (report >= REP_VERBOSE) {
580 <                fputs(progname, stderr);
581 <                fputs(": comparing inputs as binary\n", stderr);
580 >                fputs(progname, stdout);
581 >                fputs(": comparing inputs as binary\n", stdout);
582          }
583          for ( ; ; ) {                           /* exact byte matching */
584 <                int     c1 = getc(f1in);
585 <                int     c2 = getc(f2in);
584 >                c1 = getc(f1in);
585 >                c2 = getc(f2in);
586                  if (c1 == EOF) {
587                          if (c2 == EOF)
588                                  return(1);      /* success! */
589                          if (report != REP_QUIET) {
590 <                                fputs(f1name, stderr);
591 <                                fputs(": Unexpected end-of-file\n", stderr);
590 >                                fputs(f1name, stdout);
591 >                                fputs(": unexpected end-of-file\n", stdout);
592                          }
593                          return(0);
594                  }
595                  if (c2 == EOF) {
596                          if (report != REP_QUIET) {
597 <                                fputs(f2name, stderr);
598 <                                fputs(": Unexpected end-of-file\n", stderr);
597 >                                fputs(f2name, stdout);
598 >                                fputs(": unexpected end-of-file\n", stdout);
599                          }
600                          return(0);
601                  }
# Line 404 | Line 604 | compare_binary()
604          }
605          if (report == REP_QUIET)
606                  return(0);
607 <        fprintf(stderr, "%s: binary files '%s' and '%s' differ at offset %ld|%ld\n",
607 >        printf("%s: binary files '%s' and '%s' differ at byte offset %ld|%ld\n",
608                          progname, f1name, f2name, ftell(f1in), ftell(f2in));
609 +        if (report >= REP_VERBOSE)
610 +                printf("%s: byte in '%s' is 0x%X, byte in '%s' is 0x%X\n",
611 +                                progname, f1name, c1, f2name, c2);
612          return(0);
613   }
614  
# Line 413 | Line 616 | compare_binary()
616   static int
617   compare_text()
618   {
619 <        char    l1buf[4096], l2buf[4096];
619 >        LINEBUF l1buf, l2buf;
620  
621          if (report >= REP_VERBOSE) {
622 <                fputs(progname, stderr);
623 <                fputs(": comparing inputs as ASCII text\n", stderr);
622 >                fputs(progname, stdout);
623 >                fputs(": comparing inputs as ASCII text", stdout);
624 >                if (escape_newlines)
625 >                        fputs(", allowing escaped newlines", stdout);
626 >                if (comment_c) {
627 >                        fputs(", ignoring comments starting with '", stdout);
628 >                        fputc(comment_c, stdout);
629 >                        fputc('\'', stdout);
630 >                }
631 >                fputc('\n', stdout);
632          }
633 <                                                /* compare a line at a time */
634 <        while (fgets(l1buf, sizeof(l1buf), f1in)) {
633 >        SET_FILE_TEXT(f1in);                    /* originally set to binary */
634 >        SET_FILE_TEXT(f2in);
635 >        init_line(&l1buf); init_line(&l2buf);   /* compare a line at a time */
636 >        while (read_line(&l1buf, f1in)) {
637                  lin1cnt++;
638 <                if (!*sskip2(l1buf,0))
638 >                if (!*sskip2(l1buf.str,0))
639                          continue;               /* ignore empty lines */
640 <                while (fgets(l2buf, sizeof(l2buf), f2in)) {
640 >
641 >                while (read_line(&l2buf, f2in)) {
642                          lin2cnt++;
643 <                        if (*sskip2(l2buf,0))
643 >                        if (*sskip2(l2buf.str,0))
644                                  break;          /* found other non-empty line */
645                  }
646 <                if (feof(f2in)) {
646 >                if (!l2buf.len) {               /* input 2 EOF? */
647                          if (report != REP_QUIET) {
648 <                                fputs(f2name, stderr);
649 <                                fputs(": Unexpected end-of-file\n", stderr);
648 >                                fputs(f2name, stdout);
649 >                                fputs(": unexpected end-of-file\n", stdout);
650                          }
651 +                        free_line(&l1buf); free_line(&l2buf);
652                          return(0);
653                  }
654                                                  /* compare non-empty lines */
655 <                if (!equiv_string(l1buf, l2buf)) {
655 >                if (!equiv_string(l1buf.str, l2buf.str)) {
656                          if (report != REP_QUIET) {
657 <                                fprintf(stderr, "%s: inputs '%s' and '%s' differ at line %d|%d\n",
658 <                                                progname, f1name, f2name,
659 <                                                lin1cnt, lin2cnt);
660 <                                if (report >= REP_VERBOSE) {
661 <                                        fputs("------------- Mismatch -------------\n", stderr);
662 <                                        fprintf(stderr, "%s(%d):\t%s", f1name,
663 <                                                        lin1cnt, l1buf);
664 <                                        fprintf(stderr, "%s(%d):\t%s", f2name,
665 <                                                        lin2cnt, l2buf);
657 >                                printf("%s: inputs '%s' and '%s' differ at line %d",
658 >                                                progname, f1name, f2name, lin1cnt);
659 >                                if (lin1cnt != lin2cnt)
660 >                                        printf("|%d\n", lin2cnt);
661 >                                else
662 >                                        putchar('\n');
663 >                                if ( report >= REP_VERBOSE &&
664 >                                                (l1buf.len < 256) &
665 >                                                (l2buf.len < 256) ) {
666 >                                        fputs("------------- Mismatch -------------\n", stdout);
667 >                                        printf("%s@%d:\t%s", f1name,
668 >                                                        lin1cnt, l1buf.str);
669 >                                        printf("%s@%d:\t%s", f2name,
670 >                                                        lin2cnt, l2buf.str);
671                                  }
672                          }
673 +                        free_line(&l1buf); free_line(&l2buf);
674                          return(0);
675                  }
676          }
677 <                                                /* check for EOF on input 2 */
678 <        while (fgets(l2buf, sizeof(l2buf), f2in)) {
679 <                if (!*sskip2(l2buf,0))
677 >        free_line(&l1buf);                      /* check for EOF on input 2 */
678 >        while (read_line(&l2buf, f2in)) {
679 >                if (!*sskip2(l2buf.str,0))
680                          continue;
681                  if (report != REP_QUIET) {
682 <                        fputs(f1name, stderr);
683 <                        fputs(": Unexpected end-of-file\n", stderr);
682 >                        fputs(f1name, stdout);
683 >                        fputs(": unexpected end-of-file\n", stdout);
684                  }
685 +                free_line(&l2buf);
686                  return(0);
687          }
688 +        free_line(&l2buf);
689          return(good_RMS());                     /* final check for reals */
690   }
691  
692 + /* Set resolution based on NROWS, NCOLS in header */
693 + static int
694 + set_resolu(RESOLU *rs, LUTAB *htp)
695 + {
696 +        const char      *val;
697 +
698 +        rs->rt = PIXSTANDARD;
699 +        val = (const char *)lu_find(htp, "NROWS")->data;
700 +        if (!val) return(0);
701 +        rs->yr = atoi(val);
702 +        if (rs->yr <= 0) return(-1);
703 +        val = (const char *)lu_find(htp, "NCOLS")->data;
704 +        if (!val) return(0);
705 +        rs->xr = atoi(val);
706 +        if (rs->xr <= 0) return(-1);
707 +        return(1);
708 + }
709 +
710 + /* Check image/map resolutions */
711 + static int
712 + check_resolu(const char *class, RESOLU *r1p, RESOLU *r2p)
713 + {
714 +        if (r1p->rt != r2p->rt) {
715 +                if (report != REP_QUIET)
716 +                        printf(
717 +                        "%s: %ss '%s' and '%s' have different pixel ordering\n",
718 +                                        progname, class, f1name, f2name);
719 +                return(0);
720 +        }
721 +        if ((r1p->xr != r2p->xr) | (r1p->yr != r2p->yr)) {
722 +                if (report != REP_QUIET)
723 +                        printf(
724 +                        "%s: %ss '%s' and '%s' are different sizes\n",
725 +                                        progname, class, f1name, f2name);
726 +                return(0);
727 +        }
728 +        return(1);
729 + }
730 +
731   /* Compare two inputs that are known to be RGBE or XYZE images */
732   static int
733   compare_hdr()
# Line 474 | Line 736 | compare_hdr()
736          COLOR   *scan1, *scan2;
737          int     x, y;
738  
739 +        if (report >= REP_VERBOSE) {
740 +                fputs(progname, stdout);
741 +                fputs(": comparing inputs as HDR images\n", stdout);
742 +        }
743          fgetsresolu(&rs1, f1in);
744          fgetsresolu(&rs2, f2in);
745 <        if (rs1.rt != rs2.rt) {
480 <                if (report != REP_QUIET)
481 <                        fprintf(stderr,
482 <                        "%s: Images '%s' and '%s' have different pixel ordering\n",
483 <                                        progname, f1name, f2name);
745 >        if (!check_resolu("HDR image", &rs1, &rs2))
746                  return(0);
485 –        }
486 –        if ((rs1.xr != rs2.xr) | (rs1.yr != rs2.yr)) {
487 –                if (report != REP_QUIET)
488 –                        fprintf(stderr,
489 –                        "%s: Images '%s' and '%s' are different sizes\n",
490 –                                        progname, f1name, f2name);
491 –                return(0);
492 –        }
747          scan1 = (COLOR *)malloc(scanlen(&rs1)*sizeof(COLOR));
748          scan2 = (COLOR *)malloc(scanlen(&rs2)*sizeof(COLOR));
749          if (!scan1 | !scan2) {
# Line 500 | Line 754 | compare_hdr()
754                  if ((freadscan(scan1, scanlen(&rs1), f1in) < 0) |
755                                  (freadscan(scan2, scanlen(&rs2), f2in) < 0)) {
756                          if (report != REP_QUIET)
757 <                                fprintf(stderr, "%s: Unexpected end-of-file\n",
757 >                                printf("%s: unexpected end-of-file\n",
758                                                  progname);
759                          free(scan1);
760                          free(scan2);
761                          return(0);
762                  }
763 <                for (x = scanlen(&rs1); x--; ) {
764 <                        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)))
763 >                for (x = 0; x < scanlen(&rs1); x++) {
764 >                        if (color_check(scan1[x], scan2[x]))
765                                  continue;
766                          if (report != REP_QUIET) {
767 <                                fprintf(stderr,
519 <                                "%s: pixels at scanline %d offset %d differ\n",
767 >                                printf("%s: pixels at scanline %d offset %d differ\n",
768                                          progname, y, x);
769                                  if (report >= REP_VERBOSE) {
770 <                                        fprintf(stderr, "%s: (R,G,B)=(%g,%g,%g)\n",
770 >                                        printf("%s: (R,G,B)=(%g,%g,%g)\n",
771                                                  f1name, colval(scan1[x],RED),
772                                                  colval(scan1[x],GRN),
773                                                  colval(scan1[x],BLU));
774 <                                        fprintf(stderr, "%s: (R,G,B)=(%g,%g,%g)\n",
774 >                                        printf("%s: (R,G,B)=(%g,%g,%g)\n",
775                                                  f2name, colval(scan2[x],RED),
776                                                  colval(scan2[x],GRN),
777                                                  colval(scan2[x],BLU));
# Line 539 | Line 787 | compare_hdr()
787          return(good_RMS());                     /* final check of RMS */
788   }
789  
790 < const char      nsuffix[10][3] = {              /* 1st, 2nd, 3rd, etc. */
791 <                        "th","st","nd","rd","th","th","th","th","th","th"
792 <                };
790 > /* Compare two inputs that are known to be spectral images */
791 > static int
792 > compare_spec()
793 > {
794 >        static char     NCstr[] = NCOMPSTR;
795 >        RESOLU          rs1, rs2;
796 >        COLORV          *scan1, *scan2;
797 >        const char      *val;
798 >        int             x, y;
799  
800 +        if (report >= REP_VERBOSE) {
801 +                fputs(progname, stdout);
802 +                fputs(": comparing inputs as spectral images\n", stdout);
803 +        }
804 +        if (!(set_resolu(&rs1, &hdr1) || fgetsresolu(&rs1, f1in)))
805 +                return(0);
806 +        if (!(set_resolu(&rs2, &hdr2) || fgetsresolu(&rs2, f2in)))
807 +                return(0);
808 +        if (!check_resolu("Spectral image", &rs1, &rs2))
809 +                return(0);
810 +        NCstr[LNCOMPSTR-1] = '\0';
811 +        val = (const char *)lu_find(&hdr1, NCstr)->data;
812 +        if (!val || (NCSAMP = atoi(val)) < 3) {
813 +                if (report != REP_QUIET) {
814 +                        if (val)
815 +                                printf("%s: illegal # components (%d) for spectral image\n",
816 +                                        progname, NCSAMP);
817 +                        else
818 +                                printf("%s: missing %s in header for spectral image\n",
819 +                                        progname, NCstr);
820 +                }
821 +                return(0);
822 +        }
823 +        scan1 = (COLORV *)malloc(sizeof(COLORV)*NCSAMP*scanlen(&rs1));
824 +        scan2 = (COLORV *)malloc(sizeof(COLORV)*NCSAMP*scanlen(&rs2));
825 +        if (!scan1 | !scan2) {
826 +                fprintf(stderr, "%s: out of memory in compare_hdr()\n", progname);
827 +                exit(2);
828 +        }
829 +        for (y = 0; y < numscans(&rs1); y++) {
830 +                if ((freadsscan(scan1, NCSAMP, scanlen(&rs1), f1in) < 0) |
831 +                                (freadsscan(scan2, NCSAMP, scanlen(&rs2), f2in) < 0)) {
832 +                        if (report != REP_QUIET)
833 +                                printf("%s: unexpected end-of-file\n", progname);
834 +                        free(scan1);
835 +                        free(scan2);
836 +                        return(0);
837 +                }
838 +                for (x = 0; x < scanlen(&rs1); x++) {
839 +                        if (spec_check(scan1+x*NCSAMP, scan2+x*NCSAMP))
840 +                                continue;
841 +                        if (report != REP_QUIET) {
842 +                                int     k;
843 +                                printf("%s: spectra at scanline %d offset %d differ\n",
844 +                                        progname, y, x);
845 +                                if (report >= REP_VERBOSE) {
846 +                                        printf("%s: spectrum =", f1name);
847 +                                        for (k = 0; k < NCSAMP; k++)
848 +                                                printf(" %g", scan1[x*NCSAMP+k]);
849 +                                        fputc('\n', stdout);
850 +                                        printf("%s: spectrum =", f2name);
851 +                                        for (k = 0; k < NCSAMP; k++)
852 +                                                printf(" %g", scan2[x*NCSAMP+k]);
853 +                                        fputc('\n', stdout);
854 +                                }
855 +                        }
856 +                        free(scan1);
857 +                        free(scan2);
858 +                        return(0);
859 +                }
860 +        }
861 +        free(scan1);
862 +        free(scan2);
863 +        return(good_RMS());                     /* final check of RMS */
864 + }
865 +
866 + /* Set reference depth based on header variable */
867 + static int
868 + set_refdepth(DEPTHCODEC *dcp, LUTAB *htp)
869 + {
870 +        static char     depthvar[] = DEPTHSTR;
871 +        const char      *drval;
872 +
873 +        depthvar[LDEPTHSTR-1] = '\0';
874 +        drval = (const char *)lu_find(htp, depthvar)->data;
875 +        if (!drval)
876 +                return(0);
877 +        dcp->refdepth = atof(drval);
878 +        if (dcp->refdepth <= 0) {
879 +                if (report != REP_QUIET) {
880 +                        fputs(dcp->inpname, stderr);
881 +                        fputs(": bad reference depth '", stderr);
882 +                        fputs(drval, stderr);
883 +                        fputs("'\n", stderr);
884 +                }
885 +                return(-1);
886 +        }
887 +        return(1);
888 + }
889 +
890 + /* Compare two encoded depth maps */
891 + static int
892 + compare_depth()
893 + {
894 +        long            nread = 0;
895 +        DEPTHCODEC      dc1, dc2;
896 +
897 +        if (report >= REP_VERBOSE) {
898 +                fputs(progname, stdout);
899 +                fputs(": comparing inputs as depth maps\n", stdout);
900 +        }
901 +        set_dc_defaults(&dc1);
902 +        dc1.hdrflags = HF_RESIN;
903 +        dc1.finp = f1in;
904 +        dc1.inpname = f1name;
905 +        set_dc_defaults(&dc2);
906 +        dc2.hdrflags = HF_RESIN;
907 +        dc2.finp = f2in;
908 +        dc2.inpname = f2name;
909 +        if (report != REP_QUIET) {
910 +                dc1.hdrflags |= HF_STDERR;
911 +                dc2.hdrflags |= HF_STDERR;
912 +        }
913 +        if (!process_dc_header(&dc1, 0, NULL))
914 +                return(0);
915 +        if (!process_dc_header(&dc2, 0, NULL))
916 +                return(0);
917 +        if (!check_resolu("Depth map", &dc1.res, &dc2.res))
918 +                return(0);
919 +        if (set_refdepth(&dc1, &hdr1) < 0)
920 +                return(0);
921 +        if (set_refdepth(&dc2, &hdr2) < 0)
922 +                return(0);
923 +        while (nread < dc1.res.xr*dc1.res.yr) {
924 +                double  d1 = decode_depth_next(&dc1);
925 +                double  d2 = decode_depth_next(&dc2);
926 +                if ((d1 < 0) | (d2 < 0)) {
927 +                        if (report != REP_QUIET)
928 +                                printf("%s: unexpected end-of-file\n",
929 +                                                progname);
930 +                        return(0);
931 +                }
932 +                ++nread;
933 +                if (real_check(d1, d2))
934 +                        continue;
935 +                if (report != REP_QUIET)
936 +                        printf("%s: %ld%s depth values differ\n",
937 +                                        progname, nread, num_sfx(nread));
938 +                return(0);
939 +        }
940 +        return(good_RMS());             /* final check of RMS */
941 + }
942 +
943 + /* Compare two encoded normal maps */
944 + static int
945 + compare_norm()
946 + {
947 +        long            nread = 0;
948 +        NORMCODEC       nc1, nc2;
949 +
950 +        if (report >= REP_VERBOSE) {
951 +                fputs(progname, stdout);
952 +                fputs(": comparing inputs as normal maps\n", stdout);
953 +        }
954 +        set_nc_defaults(&nc1);
955 +        nc1.hdrflags = HF_RESIN;
956 +        nc1.finp = f1in;
957 +        nc1.inpname = f1name;
958 +        set_nc_defaults(&nc2);
959 +        nc2.hdrflags = HF_RESIN;
960 +        nc2.finp = f2in;
961 +        nc2.inpname = f2name;
962 +        if (report != REP_QUIET) {
963 +                nc1.hdrflags |= HF_STDERR;
964 +                nc2.hdrflags |= HF_STDERR;
965 +        }
966 +        if (!process_nc_header(&nc1, 0, NULL))
967 +                return(0);
968 +        if (!process_nc_header(&nc2, 0, NULL))
969 +                return(0);
970 +        if (!check_resolu("Normal map", &nc1.res, &nc2.res))
971 +                return(0);
972 +        while (nread < nc1.res.xr*nc1.res.yr) {
973 +                FVECT   nv1, nv2;
974 +                int     rv1 = decode_normal_next(nv1, &nc1);
975 +                int     rv2 = decode_normal_next(nv2, &nc2);
976 +                if ((rv1 < 0) | (rv2 < 0)) {
977 +                        if (report != REP_QUIET)
978 +                                printf("%s: unexpected end-of-file\n",
979 +                                                progname);
980 +                        return(0);
981 +                }
982 +                ++nread;
983 +                if (rv1 == rv2 && (!rv1 || norm_check(nv1, nv2)))
984 +                        continue;
985 +                if (report != REP_QUIET)
986 +                        printf("%s: %ld%s normal vectors differ\n",
987 +                                        progname, nread, num_sfx(nread));
988 +                return(0);
989 +        }
990 +        return(good_RMS());             /* final check of RMS */
991 + }
992 +
993   /* Compare two inputs that are known to be 32-bit floating-point data */
994   static int
995   compare_float()
# Line 550 | Line 997 | compare_float()
997          long    nread = 0;
998          float   f1, f2;
999  
1000 +        if (report >= REP_VERBOSE) {
1001 +                fputs(progname, stdout);
1002 +                fputs(": comparing inputs as 32-bit IEEE floats\n", stdout);
1003 +        }
1004          while (getbinary(&f1, sizeof(f1), 1, f1in)) {
1005                  if (!getbinary(&f2, sizeof(f2), 1, f2in))
1006                          goto badeof;
1007                  ++nread;
1008 +                if (f1swap) swap32((char *)&f1, 1);
1009 +                if (f2swap) swap32((char *)&f2, 1);
1010                  if (real_check(f1, f2))
1011                          continue;
1012                  if (report != REP_QUIET)
1013 <                        fprintf(stderr, "%s: %ld%s float values differ\n",
1014 <                                        progname, nread, nsuffix[nread%10]);
1013 >                        printf("%s: %ld%s float values differ\n",
1014 >                                        progname, nread, num_sfx(nread));
1015                  return(0);
1016          }
1017          if (!getbinary(&f2, sizeof(f2), 1, f2in))
1018                  return(good_RMS());             /* final check of RMS */
1019   badeof:
1020          if (report != REP_QUIET)
1021 <                fprintf(stderr, "%s: Unexpected end-of-file\n", progname);
1021 >                printf("%s: unexpected end-of-file\n", progname);
1022          return(0);
1023   }
1024  
# Line 576 | Line 1029 | compare_double()
1029          long    nread = 0;
1030          double  f1, f2;
1031  
1032 +        if (report >= REP_VERBOSE) {
1033 +                fputs(progname, stdout);
1034 +                fputs(": comparing inputs as 64-bit IEEE doubles\n", stdout);
1035 +        }
1036          while (getbinary(&f1, sizeof(f1), 1, f1in)) {
1037                  if (!getbinary(&f2, sizeof(f2), 1, f2in))
1038                          goto badeof;
1039                  ++nread;
1040 +                if (f1swap) swap64((char *)&f1, 1);
1041 +                if (f2swap) swap64((char *)&f2, 1);
1042                  if (real_check(f1, f2))
1043                          continue;
1044                  if (report != REP_QUIET)
1045 <                        fprintf(stderr, "%s: %ld%s float values differ\n",
1046 <                                        progname, nread, nsuffix[nread%10]);
1045 >                        printf("%s: %ld%s double values differ\n",
1046 >                                        progname, nread, num_sfx(nread));
1047                  return(0);
1048          }
1049          if (!getbinary(&f2, sizeof(f2), 1, f2in))
1050                  return(good_RMS());             /* final check of RMS */
1051   badeof:
1052          if (report != REP_QUIET)
1053 <                fprintf(stderr, "%s: Unexpected end-of-file\n", progname);
1053 >                printf("%s: unexpected end-of-file\n", progname);
1054          return(0);
1055   }
1056  
# Line 608 | Line 1067 | main(int argc, char *argv[])
1067                  case 'h':                       /* ignore header info. */
1068                          ign_header = !ign_header;
1069                          continue;
1070 +                case 'n':                       /* allow newline escapes */
1071 +                        escape_newlines = !escape_newlines;
1072 +                        continue;
1073 +                case 'c':                       /* ignore comments */
1074 +                        comment_c = argv[a][2];
1075 +                        continue;
1076                  case 's':                       /* silent operation */
1077                          report = REP_QUIET;
1078                          continue;
# Line 638 | Line 1103 | main(int argc, char *argv[])
1103                  }
1104                  break;
1105          }
1106 <        if (a != argc-2)
1106 >        if (a != argc-2)                        /* make sure of two inputs */
1107                  usage();
1108 <        if (!f1name) f1name = argv[a];
644 <        if (!f2name) f2name = argv[a+1];
645 <
646 <        if (!strcmp(f1name, f2name)) {          /* inputs are same? */
1108 >        if (!strcmp(argv[a], argv[a+1])) {      /* inputs are same? */
1109                  if (report >= REP_WARN)
1110 <                        fprintf(stderr, "%s: warning - identical inputs given\n",
1110 >                        printf("%s: warning - identical inputs given\n",
1111                                          progname);
1112                  return(0);
1113          }
1114 +        if (!f1name) f1name = argv[a];
1115 +        if (!f2name) f2name = argv[a+1];
1116                                                  /* open inputs */
1117          SET_FILE_BINARY(stdin);                 /* in case we're using it */
1118          if (!f1in && !(f1in = fopen(f1name, "rb"))) {
1119                  fprintf(stderr, "%s: cannot open for reading\n", f1name);
1120 <                return(1);
1120 >                return(2);
1121          }
1122          if (!strcmp(f2name, "-")) {
1123                  f2in = stdin;
1124                  f2name = stdin_name;
1125          } else if (!(f2in = fopen(f2name, "rb"))) {
1126                  fprintf(stderr, "%s: cannot open for reading\n", f2name);
1127 <                return(1);
1127 >                return(2);
1128          }
1129                                                  /* load headers */
1130          if ((typ1 = identify_type(f1name, f1in, &hdr1)) < 0)
1131 <                return(1);
1131 >                return(2);
1132          if ((typ2 = identify_type(f2name, f2in, &hdr2)) < 0)
1133 <                return(1);
1133 >                return(2);
1134          if (typ1 != typ2) {
1135                  if (report != REP_QUIET)
1136 <                        fprintf(stderr, "%s: '%s' is %s and '%s' is %s\n",
1136 >                        printf("%s: '%s' format is %s and '%s' is %s\n",
1137                                          progname, f1name, file_type[typ1],
1138                                          f2name, file_type[typ2]);
1139                  return(1);
1140          }
1141          ign_header |= !has_header(typ1);        /* check headers if indicated */
1142 <        if (!ign_header && !headers_match(&hdr1, &hdr2))
1142 >        if (!ign_header && !headers_match())
1143                  return(1);
680 –        lu_done(&hdr1); lu_done(&hdr2);
1144          if (!ign_header & (report >= REP_WARN)) {
682 –                if (typ1 == TYP_UNKNOWN)
683 –                        fprintf(stderr, "%s: warning - unrecognized format, comparing as binary\n",
684 –                                        progname);
1145                  if (lin1cnt != lin2cnt)
1146 <                        fprintf(stderr, "%s: warning - headers are different lengths\n",
1146 >                        printf("%s: warning - headers are different lengths\n",
1147                                          progname);
1148 +                if (typ1 == TYP_UNKNOWN)
1149 +                        printf("%s: warning - unrecognized format\n",
1150 +                                        progname);
1151          }
1152 <        if (report >= REP_VERBOSE)
1153 <                fprintf(stderr, "%s: input file type is %s\n",
1154 <                                progname, file_type[typ1]);
1155 <
1152 >        if (report >= REP_VERBOSE) {
1153 >                printf("%s: data format is %s\n", progname, file_type[typ1]);
1154 >                if ((typ1 == TYP_FLOAT) | (typ1 == TYP_DOUBLE)) {
1155 >                        if (f1swap)
1156 >                                printf("%s: input '%s' is byte-swapped\n",
1157 >                                                progname, f1name);
1158 >                        if (f2swap)
1159 >                                printf("%s: input '%s' is byte-swapped\n",
1160 >                                                progname, f2name);
1161 >                }
1162 >        }
1163          switch (typ1) {                         /* compare based on type */
1164          case TYP_BINARY:
1165          case TYP_TMESH:
1166          case TYP_OCTREE:
1167          case TYP_RBFMESH:
1168 +        case TYP_ID8:
1169 +        case TYP_ID16:
1170 +        case TYP_ID24:
1171          case TYP_UNKNOWN:
1172                  return( !compare_binary() );
1173          case TYP_TEXT:
# Line 703 | Line 1176 | main(int argc, char *argv[])
1176          case TYP_RGBE:
1177          case TYP_XYZE:
1178                  return( !compare_hdr() );
1179 +        case TYP_SPEC:
1180 +                return( !compare_spec() );
1181 +        case TYP_DEPTH:
1182 +                return( !compare_depth() );
1183 +        case TYP_NORM:
1184 +                return( !compare_norm() );
1185          case TYP_FLOAT:
1186                  return( !compare_float() );
1187          case TYP_DOUBLE:

Diff Legend

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