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

Comparing ray/src/util/radcompare.c (file contents):
Revision 2.14 by greg, Fri Oct 19 23:56:15 2018 UTC vs.
Revision 2.33 by greg, Tue Nov 21 19:22:38 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                  };
# Line 45 | Line 51 | const char     *file_type[] = {
51                          "ascii",
52                          COLRFMT,
53                          CIEFMT,
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                                  /* keep consistent with above */
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};
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 64 | Line 78 | const char     *hdr_ignkey[] = {
78                          "SOFTWARE",
79                          "CAPDATE",
80                          "GMT",
81 +                        "FRAME",
82                          NULL    /* terminator */
83                  };
84                                  /* header variable settings */
# Line 89 | Line 104 | 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;
# Line 100 | 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] reference test\n",
119 >        fputs(" [-h][-c#][-s|-w|-v][-rel min_test][-rms epsilon][-max epsilon] reference test\n",
120                          stderr);
121          exit(2);
122   }
# Line 109 | Line 125 | usage()
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);
# Line 117 | Line 135 | read_line(LINEBUF *bp, FILE *fp)
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')
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 <                        break;          /* don't go to extremes */
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,
# Line 167 | 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)
# Line 190 | Line 235 | 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.))
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;
# Line 200 | Line 245 | color_check(COLOR c1, COLOR c2)
245          return(real_check(colval(c1,p), colval(c2,p)));
246   }
247  
248 + /* Compare two color spectra for equivalence */
249 + static int
250 + spec_check(COLORV *sc1, COLORV *sc2)
251 + {
252 +        int     p, k;
253 +
254 +        if (!real_check(scolor_mean(sc1), scolor_mean(sc2)))
255 +                return(0);
256 +
257 +        p = 0;                          /* find max. component */
258 +        for (k = NCSAMP; --k; )
259 +                if (sc1[k] > sc1[p])
260 +                        p = k;
261 +
262 +        return(real_check(sc1[p], sc2[p]));
263 + }
264 +
265 + /* Compare two normal directions for equivalence */
266 + static int
267 + norm_check(FVECT nv1, FVECT nv2)
268 + {
269 +        double  max2 = nv1[2]*nv2[2];
270 +        int     imax = 2;
271 +        int     i = 2;
272 +                                        /* identify largest component */
273 +        while (i--) {
274 +                double  tm2 = nv1[i]*nv2[i];
275 +                if (tm2 > max2) {
276 +                        imax = i;
277 +                        max2 = tm2;
278 +                }
279 +        }
280 +        i = 3;                          /* compare smaller components */
281 +        while (i--) {
282 +                if (i == imax)
283 +                        continue;
284 +                if (!real_check(nv1[i], nv2[i]))
285 +                        return(0);
286 +        }
287 +        return(1);
288 + }
289 +
290   /* Compare two strings for equivalence */
291   static int
292   equiv_string(char *s1, char *s2)
# Line 272 | Line 359 | equiv_string(char *s1, char *s2)
359   static int
360   setheadvar(char *val, void *p)
361   {
362 +        char    newval[128];
363          LUTAB   *htp = (LUTAB *)p;
364          LUENT   *tep;
365          char    *key;
# Line 281 | Line 369 | setheadvar(char *val, void *p)
369          adv_linecnt(htp);       /* side-effect is to count lines */
370          if (!isalpha(*val))     /* key must start line */
371                  return(0);
372 +                                /* check if we need to swap binary data */
373 +        if ((n = isbigendian(val)) >= 0) {
374 +                if (nativebigendian() == n)
375 +                        return(0);
376 +                f1swap += (htp == &hdr1);
377 +                f2swap += (htp == &hdr2);
378 +                return(0);
379 +        }
380          key = val++;
381          while (*val && !isspace(*val) & (*val != '='))
382                  val++;
# Line 306 | Line 402 | setheadvar(char *val, void *p)
402                  return(-1);     /* memory allocation error */
403          if (!tep->key)
404                  tep->key = strcpy(malloc(kln+1), key);
405 <        if (tep->data)
405 >        if (tep->data) {        /* check for special cases */
406 >                if (!strcmp(key, "EXPOSURE")) {
407 >                        sprintf(newval, "%.6e", atof(tep->data)*atof(val));
408 >                        vln = strlen(val = newval);
409 >                }
410                  free(tep->data);
411 +        }
412          tep->data = strcpy(malloc(vln+1), val);
413          return(1);
414   }
# Line 416 | Line 517 | identify_type(const char *name, FILE *fin, LUTAB *htp)
517          }
518          if (c)
519                  return(TYP_BINARY);
419 –        SET_FILE_TEXT(fin);                     /* originally set to binary */
520          return(TYP_TEXT);
521   badeof:
522          if (report != REP_QUIET) {
# Line 499 | Line 599 | compare_text()
599  
600          if (report >= REP_VERBOSE) {
601                  fputs(progname, stdout);
602 <                fputs(": comparing inputs as ASCII text\n", stdout);
602 >                fputs(": comparing inputs as ASCII text", stdout);
603 >                if (escape_newlines)
604 >                        fputs(", allowing escaped newlines", stdout);
605 >                if (comment_c) {
606 >                        fputs(", ignoring comments starting with '", stdout);
607 >                        fputc(comment_c, stdout);
608 >                        fputc('\'', stdout);
609 >                }
610 >                fputc('\n', stdout);
611          }
612 +        SET_FILE_TEXT(f1in);                    /* originally set to binary */
613 +        SET_FILE_TEXT(f2in);
614          init_line(&l1buf); init_line(&l2buf);   /* compare a line at a time */
615          while (read_line(&l1buf, f1in)) {
616                  lin1cnt++;
# Line 523 | Line 633 | compare_text()
633                                                  /* compare non-empty lines */
634                  if (!equiv_string(l1buf.str, l2buf.str)) {
635                          if (report != REP_QUIET) {
636 <                                printf("%s: inputs '%s' and '%s' differ at line %d|%d\n",
637 <                                                progname, f1name, f2name,
638 <                                                lin1cnt, lin2cnt);
636 >                                printf("%s: inputs '%s' and '%s' differ at line %d",
637 >                                                progname, f1name, f2name, lin1cnt);
638 >                                if (lin1cnt != lin2cnt)
639 >                                        printf("|%d\n", lin2cnt);
640 >                                else
641 >                                        putchar('\n');
642                                  if ( report >= REP_VERBOSE &&
643                                                  (l1buf.len < 256) &
644                                                  (l2buf.len < 256) ) {
# Line 555 | Line 668 | compare_text()
668          return(good_RMS());                     /* final check for reals */
669   }
670  
671 + /* Set resolution based on NROWS, NCOLS in header */
672 + static int
673 + set_resolu(RESOLU *rs, LUTAB *htp)
674 + {
675 +        const char      *val;
676 +
677 +        rs->rt = PIXSTANDARD;
678 +        val = (const char *)lu_find(htp, "NROWS")->data;
679 +        if (!val) return(0);
680 +        rs->yr = atoi(val);
681 +        if (rs->yr <= 0) return(-1);
682 +        val = (const char *)lu_find(htp, "NCOLS")->data;
683 +        if (!val) return(0);
684 +        rs->xr = atoi(val);
685 +        if (rs->xr <= 0) return(-1);
686 +        return(1);
687 + }
688 +
689 + /* Check image/map resolutions */
690 + static int
691 + check_resolu(const char *class, RESOLU *r1p, RESOLU *r2p)
692 + {
693 +        if (r1p->rt != r2p->rt) {
694 +                if (report != REP_QUIET)
695 +                        printf(
696 +                        "%s: %ss '%s' and '%s' have different pixel ordering\n",
697 +                                        progname, class, f1name, f2name);
698 +                return(0);
699 +        }
700 +        if ((r1p->xr != r2p->xr) | (r1p->yr != r2p->yr)) {
701 +                if (report != REP_QUIET)
702 +                        printf(
703 +                        "%s: %ss '%s' and '%s' are different sizes\n",
704 +                                        progname, class, f1name, f2name);
705 +                return(0);
706 +        }
707 +        return(1);
708 + }
709 +
710   /* Compare two inputs that are known to be RGBE or XYZE images */
711   static int
712   compare_hdr()
# Line 569 | Line 721 | compare_hdr()
721          }
722          fgetsresolu(&rs1, f1in);
723          fgetsresolu(&rs2, f2in);
724 <        if (rs1.rt != rs2.rt) {
573 <                if (report != REP_QUIET)
574 <                        printf(
575 <                        "%s: Images '%s' and '%s' have different pixel ordering\n",
576 <                                        progname, f1name, f2name);
724 >        if (!check_resolu("HDR image", &rs1, &rs2))
725                  return(0);
578 –        }
579 –        if ((rs1.xr != rs2.xr) | (rs1.yr != rs2.yr)) {
580 –                if (report != REP_QUIET)
581 –                        printf(
582 –                        "%s: Images '%s' and '%s' are different sizes\n",
583 –                                        progname, f1name, f2name);
584 –                return(0);
585 –        }
726          scan1 = (COLOR *)malloc(scanlen(&rs1)*sizeof(COLOR));
727          scan2 = (COLOR *)malloc(scanlen(&rs2)*sizeof(COLOR));
728          if (!scan1 | !scan2) {
# Line 603 | Line 743 | compare_hdr()
743                          if (color_check(scan1[x], scan2[x]))
744                                  continue;
745                          if (report != REP_QUIET) {
746 <                                printf(
607 <                                "%s: pixels at scanline %d offset %d differ\n",
746 >                                printf("%s: pixels at scanline %d offset %d differ\n",
747                                          progname, y, x);
748                                  if (report >= REP_VERBOSE) {
749                                          printf("%s: (R,G,B)=(%g,%g,%g)\n",
# Line 627 | Line 766 | compare_hdr()
766          return(good_RMS());                     /* final check of RMS */
767   }
768  
769 + /* Compare two inputs that are known to be spectral images */
770 + static int
771 + compare_spec()
772 + {
773 +        static char     NCstr[] = NCOMPSTR;
774 +        RESOLU          rs1, rs2;
775 +        COLORV          *scan1, *scan2;
776 +        const char      *val;
777 +        int             x, y;
778 +
779 +        if (report >= REP_VERBOSE) {
780 +                fputs(progname, stdout);
781 +                fputs(": comparing inputs as spectral images\n", stdout);
782 +        }
783 +        if (!(set_resolu(&rs1, &hdr1) || fgetsresolu(&rs1, f1in)))
784 +                return(0);
785 +        if (!(set_resolu(&rs2, &hdr2) || fgetsresolu(&rs2, f2in)))
786 +                return(0);
787 +        if (!check_resolu("Spectral image", &rs1, &rs2))
788 +                return(0);
789 +        NCstr[LNCOMPSTR-1] = '\0';
790 +        val = (const char *)lu_find(&hdr1, NCstr)->data;
791 +        if (!val || (NCSAMP = atoi(val)) < 3) {
792 +                if (report != REP_QUIET) {
793 +                        if (val)
794 +                                printf("%s: illegal # components (%d) for spectral image\n",
795 +                                        progname, NCSAMP);
796 +                        else
797 +                                printf("%s: missing %s in header for spectral image\n",
798 +                                        progname, NCstr);
799 +                }
800 +                return(0);
801 +        }
802 +        scan1 = (COLORV *)malloc(sizeof(COLORV)*NCSAMP*scanlen(&rs1));
803 +        scan2 = (COLORV *)malloc(sizeof(COLORV)*NCSAMP*scanlen(&rs2));
804 +        if (!scan1 | !scan2) {
805 +                fprintf(stderr, "%s: out of memory in compare_hdr()\n", progname);
806 +                exit(2);
807 +        }
808 +        for (y = 0; y < numscans(&rs1); y++) {
809 +                if ((freadsscan(scan1, NCSAMP, scanlen(&rs1), f1in) < 0) |
810 +                                (freadsscan(scan2, NCSAMP, scanlen(&rs2), f2in) < 0)) {
811 +                        if (report != REP_QUIET)
812 +                                printf("%s: unexpected end-of-file\n", progname);
813 +                        free(scan1);
814 +                        free(scan2);
815 +                        return(0);
816 +                }
817 +                for (x = 0; x < scanlen(&rs1); x++) {
818 +                        if (spec_check(scan1+x*NCSAMP, scan2+x*NCSAMP))
819 +                                continue;
820 +                        if (report != REP_QUIET) {
821 +                                int     k;
822 +                                printf("%s: spectra at scanline %d offset %d differ\n",
823 +                                        progname, y, x);
824 +                                if (report >= REP_VERBOSE) {
825 +                                        printf("%s: spectrum =", f1name);
826 +                                        for (k = 0; k < NCSAMP; k++)
827 +                                                printf(" %g", scan1[x*NCSAMP+k]);
828 +                                        fputc('\n', stdout);
829 +                                        printf("%s: spectrum =", f2name);
830 +                                        for (k = 0; k < NCSAMP; k++)
831 +                                                printf(" %g", scan2[x*NCSAMP+k]);
832 +                                        fputc('\n', stdout);
833 +                                }
834 +                        }
835 +                        free(scan1);
836 +                        free(scan2);
837 +                        return(0);
838 +                }
839 +        }
840 +        free(scan1);
841 +        free(scan2);
842 +        return(good_RMS());                     /* final check of RMS */
843 + }
844 +
845 + /* Set reference depth based on header variable */
846 + static int
847 + set_refdepth(DEPTHCODEC *dcp, LUTAB *htp)
848 + {
849 +        static char     depthvar[] = DEPTHSTR;
850 +        const char      *drval;
851 +
852 +        depthvar[LDEPTHSTR-1] = '\0';
853 +        drval = (const char *)lu_find(htp, depthvar)->data;
854 +        if (!drval)
855 +                return(0);
856 +        dcp->refdepth = atof(drval);
857 +        if (dcp->refdepth <= 0) {
858 +                if (report != REP_QUIET) {
859 +                        fputs(dcp->inpname, stderr);
860 +                        fputs(": bad reference depth '", stderr);
861 +                        fputs(drval, stderr);
862 +                        fputs("'\n", stderr);
863 +                }
864 +                return(-1);
865 +        }
866 +        return(1);
867 + }
868 +
869 + /* Compare two encoded depth maps */
870 + static int
871 + compare_depth()
872 + {
873 +        long            nread = 0;
874 +        DEPTHCODEC      dc1, dc2;
875 +
876 +        if (report >= REP_VERBOSE) {
877 +                fputs(progname, stdout);
878 +                fputs(": comparing inputs as depth maps\n", stdout);
879 +        }
880 +        set_dc_defaults(&dc1);
881 +        dc1.hdrflags = HF_RESIN;
882 +        dc1.finp = f1in;
883 +        dc1.inpname = f1name;
884 +        set_dc_defaults(&dc2);
885 +        dc2.hdrflags = HF_RESIN;
886 +        dc2.finp = f2in;
887 +        dc2.inpname = f2name;
888 +        if (report != REP_QUIET) {
889 +                dc1.hdrflags |= HF_STDERR;
890 +                dc2.hdrflags |= HF_STDERR;
891 +        }
892 +        if (!process_dc_header(&dc1, 0, NULL))
893 +                return(0);
894 +        if (!process_dc_header(&dc2, 0, NULL))
895 +                return(0);
896 +        if (!check_resolu("Depth map", &dc1.res, &dc2.res))
897 +                return(0);
898 +        if (set_refdepth(&dc1, &hdr1) < 0)
899 +                return(0);
900 +        if (set_refdepth(&dc2, &hdr2) < 0)
901 +                return(0);
902 +        while (nread < dc1.res.xr*dc1.res.yr) {
903 +                double  d1 = decode_depth_next(&dc1);
904 +                double  d2 = decode_depth_next(&dc2);
905 +                if ((d1 < 0) | (d2 < 0)) {
906 +                        if (report != REP_QUIET)
907 +                                printf("%s: unexpected end-of-file\n",
908 +                                                progname);
909 +                        return(0);
910 +                }
911 +                ++nread;
912 +                if (real_check(d1, d2))
913 +                        continue;
914 +                if (report != REP_QUIET)
915 +                        printf("%s: %ld%s depth values differ\n",
916 +                                        progname, nread, num_sfx(nread));
917 +                return(0);
918 +        }
919 +        return(good_RMS());             /* final check of RMS */
920 + }
921 +
922 + /* Compare two encoded normal maps */
923 + static int
924 + compare_norm()
925 + {
926 +        long            nread = 0;
927 +        NORMCODEC       nc1, nc2;
928 +
929 +        if (report >= REP_VERBOSE) {
930 +                fputs(progname, stdout);
931 +                fputs(": comparing inputs as normal maps\n", stdout);
932 +        }
933 +        set_nc_defaults(&nc1);
934 +        nc1.hdrflags = HF_RESIN;
935 +        nc1.finp = f1in;
936 +        nc1.inpname = f1name;
937 +        set_nc_defaults(&nc2);
938 +        nc2.hdrflags = HF_RESIN;
939 +        nc2.finp = f2in;
940 +        nc2.inpname = f2name;
941 +        if (report != REP_QUIET) {
942 +                nc1.hdrflags |= HF_STDERR;
943 +                nc2.hdrflags |= HF_STDERR;
944 +        }
945 +        if (!process_nc_header(&nc1, 0, NULL))
946 +                return(0);
947 +        if (!process_nc_header(&nc2, 0, NULL))
948 +                return(0);
949 +        if (!check_resolu("Normal map", &nc1.res, &nc2.res))
950 +                return(0);
951 +        while (nread < nc1.res.xr*nc1.res.yr) {
952 +                FVECT   nv1, nv2;
953 +                int     rv1 = decode_normal_next(nv1, &nc1);
954 +                int     rv2 = decode_normal_next(nv2, &nc2);
955 +                if ((rv1 < 0) | (rv2 < 0)) {
956 +                        if (report != REP_QUIET)
957 +                                printf("%s: unexpected end-of-file\n",
958 +                                                progname);
959 +                        return(0);
960 +                }
961 +                ++nread;
962 +                if (rv1 == rv2 && (!rv1 || norm_check(nv1, nv2)))
963 +                        continue;
964 +                if (report != REP_QUIET)
965 +                        printf("%s: %ld%s normal vectors differ\n",
966 +                                        progname, nread, num_sfx(nread));
967 +                return(0);
968 +        }
969 +        return(good_RMS());             /* final check of RMS */
970 + }
971 +
972   /* Compare two inputs that are known to be 32-bit floating-point data */
973   static int
974   compare_float()
# Line 642 | Line 984 | compare_float()
984                  if (!getbinary(&f2, sizeof(f2), 1, f2in))
985                          goto badeof;
986                  ++nread;
987 +                if (f1swap) swap32((char *)&f1, 1);
988 +                if (f2swap) swap32((char *)&f2, 1);
989                  if (real_check(f1, f2))
990                          continue;
991                  if (report != REP_QUIET)
# Line 672 | Line 1016 | compare_double()
1016                  if (!getbinary(&f2, sizeof(f2), 1, f2in))
1017                          goto badeof;
1018                  ++nread;
1019 +                if (f1swap) swap64((char *)&f1, 1);
1020 +                if (f2swap) swap64((char *)&f2, 1);
1021                  if (real_check(f1, f2))
1022                          continue;
1023                  if (report != REP_QUIET)
# Line 700 | Line 1046 | main(int argc, char *argv[])
1046                  case 'h':                       /* ignore header info. */
1047                          ign_header = !ign_header;
1048                          continue;
1049 +                case 'n':                       /* allow newline escapes */
1050 +                        escape_newlines = !escape_newlines;
1051 +                        continue;
1052 +                case 'c':                       /* ignore comments */
1053 +                        comment_c = argv[a][2];
1054 +                        continue;
1055                  case 's':                       /* silent operation */
1056                          report = REP_QUIET;
1057                          continue;
# Line 760 | Line 1112 | main(int argc, char *argv[])
1112                  return(2);
1113          if (typ1 != typ2) {
1114                  if (report != REP_QUIET)
1115 <                        printf("%s: '%s' is %s and '%s' is %s\n",
1115 >                        printf("%s: '%s' format is %s and '%s' is %s\n",
1116                                          progname, f1name, file_type[typ1],
1117                                          f2name, file_type[typ2]);
1118                  return(1);
# Line 768 | Line 1120 | main(int argc, char *argv[])
1120          ign_header |= !has_header(typ1);        /* check headers if indicated */
1121          if (!ign_header && !headers_match())
1122                  return(1);
771 –        lu_done(&hdr1); lu_done(&hdr2);
1123          if (!ign_header & (report >= REP_WARN)) {
773 –                if (typ1 == TYP_UNKNOWN)
774 –                        printf("%s: warning - unrecognized format, comparing as binary\n",
775 –                                        progname);
1124                  if (lin1cnt != lin2cnt)
1125                          printf("%s: warning - headers are different lengths\n",
1126                                          progname);
1127 +                if (typ1 == TYP_UNKNOWN)
1128 +                        printf("%s: warning - unrecognized format\n",
1129 +                                        progname);
1130          }
1131 <        if (report >= REP_VERBOSE)
1132 <                printf("%s: input file type is %s\n",
1133 <                                progname, file_type[typ1]);
1134 <
1131 >        if (report >= REP_VERBOSE) {
1132 >                printf("%s: data format is %s\n", progname, file_type[typ1]);
1133 >                if ((typ1 == TYP_FLOAT) | (typ1 == TYP_DOUBLE)) {
1134 >                        if (f1swap)
1135 >                                printf("%s: input '%s' is byte-swapped\n",
1136 >                                                progname, f1name);
1137 >                        if (f2swap)
1138 >                                printf("%s: input '%s' is byte-swapped\n",
1139 >                                                progname, f2name);
1140 >                }
1141 >        }
1142          switch (typ1) {                         /* compare based on type */
1143          case TYP_BINARY:
1144          case TYP_TMESH:
1145          case TYP_OCTREE:
1146          case TYP_RBFMESH:
1147 +        case TYP_ID8:
1148 +        case TYP_ID16:
1149 +        case TYP_ID24:
1150          case TYP_UNKNOWN:
1151                  return( !compare_binary() );
1152          case TYP_TEXT:
# Line 794 | Line 1155 | main(int argc, char *argv[])
1155          case TYP_RGBE:
1156          case TYP_XYZE:
1157                  return( !compare_hdr() );
1158 +        case TYP_SPEC:
1159 +                return( !compare_spec() );
1160 +        case TYP_DEPTH:
1161 +                return( !compare_depth() );
1162 +        case TYP_NORM:
1163 +                return( !compare_norm() );
1164          case TYP_FLOAT:
1165                  return( !compare_float() );
1166          case TYP_DOUBLE:

Diff Legend

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