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.17 by greg, Sat Nov 17 20:22:17 2018 UTC vs.
Revision 2.29 by greg, Thu Mar 3 00:01:48 2022 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 +                        DEPTH16FMT,
55 +                        NORMAL32FMT,
56                          "float",
57                          "double",
58                          "BSDF_RBFmesh",
59                          "Radiance_octree",
60                          "Radiance_tmesh",
61 +                        "8-bit_indexed_name",
62 +                        "16-bit_indexed_name",
63 +                        "24-bit_indexed_name",
64                          "BINARY_unknown",
65                          NULL    /* terminator */
66                  };
67                                  /* keep consistent with above */
68 < enum {TYP_UNKNOWN, TYP_TEXT, TYP_ASCII, TYP_RGBE, TYP_XYZE, TYP_FLOAT,
69 <                TYP_DOUBLE, TYP_RBFMESH, TYP_OCTREE, TYP_TMESH, TYP_BINARY};
68 > enum {TYP_UNKNOWN, TYP_TEXT, TYP_ASCII, TYP_RGBE, TYP_XYZE,
69 >                TYP_DEPTH, TYP_NORM, TYP_FLOAT, TYP_DOUBLE,
70 >                TYP_RBFMESH, TYP_OCTREE, TYP_TMESH,
71 >                TYP_ID8, TYP_ID16, TYP_ID24, TYP_BINARY};
72                  
73   #define has_header(t)   (!( 1L<<(t) & (1L<<TYP_TEXT | 1L<<TYP_BINARY) ))
74  
# Line 64 | Line 77 | const char     *hdr_ignkey[] = {
77                          "SOFTWARE",
78                          "CAPDATE",
79                          "GMT",
80 +                        "FRAME",
81                          NULL    /* terminator */
82                  };
83                                  /* header variable settings */
# Line 89 | Line 103 | char           *progname = NULL;
103   const char      stdin_name[] = "<stdin>";
104   const char      *f1name=NULL, *f2name=NULL;
105   FILE            *f1in=NULL, *f2in=NULL;
106 + int             f1swap=0, f2swap=0;
107  
108                                  /* running real differences */
109   double  diff2sum = 0;
# Line 100 | Line 115 | usage()
115   {
116          fputs("Usage: ", stderr);
117          fputs(progname, stderr);
118 <        fputs(" [-h][-s|-w|-v][-rel min_test][-rms epsilon][-max epsilon] reference test\n",
118 >        fputs(" [-h][-c#][-s|-w|-v][-rel min_test][-rms epsilon][-max epsilon] reference test\n",
119                          stderr);
120          exit(2);
121   }
# Line 119 | Line 134 | read_line(LINEBUF *bp, FILE *fp)
134          }
135          while (fgets(bp->str + bp->len, bp->siz - bp->len, fp)) {
136                  bp->len += strlen(bp->str + bp->len);
137 <                if (bp->str[bp->len-1] == '\n')
137 >                if (bp->str[bp->len-1] == '\n') {
138 >                        if (bp->len > 1 && bp->str[bp->len-2] == '\r') {
139 >                                bp->str[--bp->len] = '\0';
140 >                                bp->str[bp->len-1] = '\n';
141 >                        }
142 >                        if (escape_newlines && bp->len > 1 &&
143 >                                        bp->str[bp->len-2] == '\\') {
144 >                                bp->str[--bp->len] = '\0';
145 >                                bp->str[bp->len-1] = ' ';
146 >                                continue;
147 >                        }
148                          break;          /* found EOL */
149 +                }
150                  if (bp->len < bp->siz - 4)
151                          continue;       /* at EOF? */
152                  if (bp->siz >= MAXBUF) {
# Line 138 | Line 164 | read_line(LINEBUF *bp, FILE *fp)
164                  if (!bp->str)
165                          goto memerr;
166          }
167 +        if (comment_c) {                /* elide comment? */
168 +                char    *cp = sskip2(bp->str,0);
169 +                if (*cp == comment_c) {
170 +                        *cp++ = '\n';
171 +                        *cp = '\0';
172 +                        bp->len = cp - bp->str;
173 +                }
174 +        }
175          return(bp->len);
176   memerr:
177          fprintf(stderr,
# Line 176 | Line 210 | real_check(double r1, double r2)
210  
211          if (rel_min > 0) {      /* doing relative differences? */
212                  double  av2 = .25*(r1*r1 + 2.*fabs(r1*r2) + r2*r2);
213 <                if (av2 > rel_min*rel_min)
214 <                        diff2 /= av2;
213 >                if (av2 < rel_min*rel_min)
214 >                        av2 = rel_min*rel_min;
215 >                diff2 /= av2;
216          }
217          if (max_lim >= 0 && diff2 > max_lim*max_lim) {
218                  if (report != REP_QUIET)
# Line 209 | Line 244 | color_check(COLOR c1, COLOR c2)
244          return(real_check(colval(c1,p), colval(c2,p)));
245   }
246  
247 + /* Compare two normal directions for equivalence */
248 + static int
249 + norm_check(FVECT nv1, FVECT nv2)
250 + {
251 +        double  max2 = nv1[2]*nv2[2];
252 +        int     imax = 2;
253 +        int     i = 2;
254 +                                        /* identify largest component */
255 +        while (i--) {
256 +                double  tm2 = nv1[i]*nv2[i];
257 +                if (tm2 > max2) {
258 +                        imax = i;
259 +                        max2 = tm2;
260 +                }
261 +        }
262 +        i = 3;                          /* compare smaller components */
263 +        while (i--) {
264 +                if (i == imax)
265 +                        continue;
266 +                if (!real_check(nv1[i], nv2[i]))
267 +                        return(0);
268 +        }
269 +        return(1);
270 + }
271 +
272   /* Compare two strings for equivalence */
273   static int
274   equiv_string(char *s1, char *s2)
# Line 281 | Line 341 | equiv_string(char *s1, char *s2)
341   static int
342   setheadvar(char *val, void *p)
343   {
344 +        char    newval[128];
345          LUTAB   *htp = (LUTAB *)p;
346          LUENT   *tep;
347          char    *key;
# Line 290 | Line 351 | setheadvar(char *val, void *p)
351          adv_linecnt(htp);       /* side-effect is to count lines */
352          if (!isalpha(*val))     /* key must start line */
353                  return(0);
354 +                                /* check if we need to swap binary data */
355 +        if ((n = isbigendian(val)) >= 0) {
356 +                if (nativebigendian() == n)
357 +                        return(0);
358 +                f1swap += (htp == &hdr1);
359 +                f2swap += (htp == &hdr2);
360 +                return(0);
361 +        }
362          key = val++;
363          while (*val && !isspace(*val) & (*val != '='))
364                  val++;
# Line 315 | Line 384 | setheadvar(char *val, void *p)
384                  return(-1);     /* memory allocation error */
385          if (!tep->key)
386                  tep->key = strcpy(malloc(kln+1), key);
387 <        if (tep->data)
387 >        if (tep->data) {        /* check for special cases */
388 >                if (!strcmp(key, "EXPOSURE")) {
389 >                        sprintf(newval, "%.6e", atof(tep->data)*atof(val));
390 >                        vln = strlen(val = newval);
391 >                }
392                  free(tep->data);
393 +        }
394          tep->data = strcpy(malloc(vln+1), val);
395          return(1);
396   }
# Line 425 | Line 499 | identify_type(const char *name, FILE *fin, LUTAB *htp)
499          }
500          if (c)
501                  return(TYP_BINARY);
428        SET_FILE_TEXT(fin);                     /* originally set to binary */
502          return(TYP_TEXT);
503   badeof:
504          if (report != REP_QUIET) {
# Line 508 | Line 581 | compare_text()
581  
582          if (report >= REP_VERBOSE) {
583                  fputs(progname, stdout);
584 <                fputs(": comparing inputs as ASCII text\n", stdout);
584 >                fputs(": comparing inputs as ASCII text", stdout);
585 >                if (escape_newlines)
586 >                        fputs(", allowing escaped newlines", stdout);
587 >                if (comment_c) {
588 >                        fputs(", ignoring comments starting with '", stdout);
589 >                        fputc(comment_c, stdout);
590 >                        fputc('\'', stdout);
591 >                }
592 >                fputc('\n', stdout);
593          }
594 +        SET_FILE_TEXT(f1in);                    /* originally set to binary */
595 +        SET_FILE_TEXT(f2in);
596          init_line(&l1buf); init_line(&l2buf);   /* compare a line at a time */
597          while (read_line(&l1buf, f1in)) {
598                  lin1cnt++;
# Line 564 | Line 647 | compare_text()
647          return(good_RMS());                     /* final check for reals */
648   }
649  
650 + /* Check image/map resolutions */
651 + static int
652 + check_resolu(const char *class, RESOLU *r1p, RESOLU *r2p)
653 + {
654 +        if (r1p->rt != r2p->rt) {
655 +                if (report != REP_QUIET)
656 +                        printf(
657 +                        "%s: %ss '%s' and '%s' have different pixel ordering\n",
658 +                                        progname, class, f1name, f2name);
659 +                return(0);
660 +        }
661 +        if ((r1p->xr != r2p->xr) | (r1p->yr != r2p->yr)) {
662 +                if (report != REP_QUIET)
663 +                        printf(
664 +                        "%s: %ss '%s' and '%s' are different sizes\n",
665 +                                        progname, class, f1name, f2name);
666 +                return(0);
667 +        }
668 +        return(1);
669 + }
670 +
671   /* Compare two inputs that are known to be RGBE or XYZE images */
672   static int
673   compare_hdr()
# Line 578 | Line 682 | compare_hdr()
682          }
683          fgetsresolu(&rs1, f1in);
684          fgetsresolu(&rs2, f2in);
685 <        if (rs1.rt != rs2.rt) {
582 <                if (report != REP_QUIET)
583 <                        printf(
584 <                        "%s: Images '%s' and '%s' have different pixel ordering\n",
585 <                                        progname, f1name, f2name);
685 >        if (!check_resolu("HDR image", &rs1, &rs2))
686                  return(0);
587        }
588        if ((rs1.xr != rs2.xr) | (rs1.yr != rs2.yr)) {
589                if (report != REP_QUIET)
590                        printf(
591                        "%s: Images '%s' and '%s' are different sizes\n",
592                                        progname, f1name, f2name);
593                return(0);
594        }
687          scan1 = (COLOR *)malloc(scanlen(&rs1)*sizeof(COLOR));
688          scan2 = (COLOR *)malloc(scanlen(&rs2)*sizeof(COLOR));
689          if (!scan1 | !scan2) {
# Line 636 | Line 728 | compare_hdr()
728          return(good_RMS());                     /* final check of RMS */
729   }
730  
731 + /* Set reference depth based on header variable */
732 + static int
733 + set_refdepth(DEPTHCODEC *dcp, LUTAB *htp)
734 + {
735 +        static char     depthvar[] = DEPTHSTR;
736 +        const char      *drval;
737 +
738 +        depthvar[LDEPTHSTR-1] = '\0';
739 +        drval = (const char *)lu_find(htp, depthvar)->data;
740 +        if (!drval)
741 +                return(0);
742 +        dcp->refdepth = atof(drval);
743 +        if (dcp->refdepth <= 0) {
744 +                if (report != REP_QUIET) {
745 +                        fputs(dcp->inpname, stderr);
746 +                        fputs(": bad reference depth '", stderr);
747 +                        fputs(drval, stderr);
748 +                        fputs("'\n", stderr);
749 +                }
750 +                return(-1);
751 +        }
752 +        return(1);
753 + }
754 +
755 + /* Compare two encoded depth maps */
756 + static int
757 + compare_depth()
758 + {
759 +        long            nread = 0;
760 +        DEPTHCODEC      dc1, dc2;
761 +
762 +        if (report >= REP_VERBOSE) {
763 +                fputs(progname, stdout);
764 +                fputs(": comparing inputs as depth maps\n", stdout);
765 +        }
766 +        set_dc_defaults(&dc1);
767 +        dc1.hdrflags = HF_RESIN;
768 +        dc1.finp = f1in;
769 +        dc1.inpname = f1name;
770 +        set_dc_defaults(&dc2);
771 +        dc2.hdrflags = HF_RESIN;
772 +        dc2.finp = f2in;
773 +        dc2.inpname = f2name;
774 +        if (report != REP_QUIET) {
775 +                dc1.hdrflags |= HF_STDERR;
776 +                dc2.hdrflags |= HF_STDERR;
777 +        }
778 +        if (!process_dc_header(&dc1, 0, NULL))
779 +                return(0);
780 +        if (!process_dc_header(&dc2, 0, NULL))
781 +                return(0);
782 +        if (!check_resolu("Depth map", &dc1.res, &dc2.res))
783 +                return(0);
784 +        if (set_refdepth(&dc1, &hdr1) < 0)
785 +                return(0);
786 +        if (set_refdepth(&dc2, &hdr2) < 0)
787 +                return(0);
788 +        while (nread < dc1.res.xr*dc1.res.yr) {
789 +                double  d1 = decode_depth_next(&dc1);
790 +                double  d2 = decode_depth_next(&dc2);
791 +                if ((d1 < 0) | (d2 < 0)) {
792 +                        if (report != REP_QUIET)
793 +                                printf("%s: unexpected end-of-file\n",
794 +                                                progname);
795 +                        return(0);
796 +                }
797 +                ++nread;
798 +                if (real_check(d1, d2))
799 +                        continue;
800 +                if (report != REP_QUIET)
801 +                        printf("%s: %ld%s depth values differ\n",
802 +                                        progname, nread, num_sfx(nread));
803 +                return(0);
804 +        }
805 +        return(good_RMS());             /* final check of RMS */
806 + }
807 +
808 + /* Compare two encoded normal maps */
809 + static int
810 + compare_norm()
811 + {
812 +        long            nread = 0;
813 +        NORMCODEC       nc1, nc2;
814 +
815 +        if (report >= REP_VERBOSE) {
816 +                fputs(progname, stdout);
817 +                fputs(": comparing inputs as normal maps\n", stdout);
818 +        }
819 +        set_nc_defaults(&nc1);
820 +        nc1.hdrflags = HF_RESIN;
821 +        nc1.finp = f1in;
822 +        nc1.inpname = f1name;
823 +        set_nc_defaults(&nc2);
824 +        nc2.hdrflags = HF_RESIN;
825 +        nc2.finp = f2in;
826 +        nc2.inpname = f2name;
827 +        if (report != REP_QUIET) {
828 +                nc1.hdrflags |= HF_STDERR;
829 +                nc2.hdrflags |= HF_STDERR;
830 +        }
831 +        if (!process_nc_header(&nc1, 0, NULL))
832 +                return(0);
833 +        if (!process_nc_header(&nc2, 0, NULL))
834 +                return(0);
835 +        if (!check_resolu("Normal map", &nc1.res, &nc2.res))
836 +                return(0);
837 +        while (nread < nc1.res.xr*nc1.res.yr) {
838 +                FVECT   nv1, nv2;
839 +                int     rv1 = decode_normal_next(nv1, &nc1);
840 +                int     rv2 = decode_normal_next(nv2, &nc2);
841 +                if ((rv1 < 0) | (rv2 < 0)) {
842 +                        if (report != REP_QUIET)
843 +                                printf("%s: unexpected end-of-file\n",
844 +                                                progname);
845 +                        return(0);
846 +                }
847 +                ++nread;
848 +                if (rv1 == rv2 && (!rv1 || norm_check(nv1, nv2)))
849 +                        continue;
850 +                if (report != REP_QUIET)
851 +                        printf("%s: %ld%s normal vectors differ\n",
852 +                                        progname, nread, num_sfx(nread));
853 +                return(0);
854 +        }
855 +        return(good_RMS());             /* final check of RMS */
856 + }
857 +
858   /* Compare two inputs that are known to be 32-bit floating-point data */
859   static int
860   compare_float()
# Line 651 | Line 870 | compare_float()
870                  if (!getbinary(&f2, sizeof(f2), 1, f2in))
871                          goto badeof;
872                  ++nread;
873 +                if (f1swap) swap32((char *)&f1, 1);
874 +                if (f2swap) swap32((char *)&f2, 1);
875                  if (real_check(f1, f2))
876                          continue;
877                  if (report != REP_QUIET)
# Line 681 | Line 902 | compare_double()
902                  if (!getbinary(&f2, sizeof(f2), 1, f2in))
903                          goto badeof;
904                  ++nread;
905 +                if (f1swap) swap64((char *)&f1, 1);
906 +                if (f2swap) swap64((char *)&f2, 1);
907                  if (real_check(f1, f2))
908                          continue;
909                  if (report != REP_QUIET)
# Line 709 | Line 932 | main(int argc, char *argv[])
932                  case 'h':                       /* ignore header info. */
933                          ign_header = !ign_header;
934                          continue;
935 +                case 'n':                       /* allow newline escapes */
936 +                        escape_newlines = !escape_newlines;
937 +                        continue;
938 +                case 'c':                       /* ignore comments */
939 +                        comment_c = argv[a][2];
940 +                        continue;
941                  case 's':                       /* silent operation */
942                          report = REP_QUIET;
943                          continue;
# Line 769 | Line 998 | main(int argc, char *argv[])
998                  return(2);
999          if (typ1 != typ2) {
1000                  if (report != REP_QUIET)
1001 <                        printf("%s: '%s' is %s and '%s' is %s\n",
1001 >                        printf("%s: '%s' format is %s and '%s' is %s\n",
1002                                          progname, f1name, file_type[typ1],
1003                                          f2name, file_type[typ2]);
1004                  return(1);
# Line 777 | Line 1006 | main(int argc, char *argv[])
1006          ign_header |= !has_header(typ1);        /* check headers if indicated */
1007          if (!ign_header && !headers_match())
1008                  return(1);
780        lu_done(&hdr1); lu_done(&hdr2);         /* done with header info. */
1009          if (!ign_header & (report >= REP_WARN)) {
1010                  if (lin1cnt != lin2cnt)
1011                          printf("%s: warning - headers are different lengths\n",
# Line 786 | Line 1014 | main(int argc, char *argv[])
1014                          printf("%s: warning - unrecognized format\n",
1015                                          progname);
1016          }
1017 <        if (report >= REP_VERBOSE)
1018 <                printf("%s: input file type is %s\n",
1019 <                                progname, file_type[typ1]);
1020 <
1017 >        if (report >= REP_VERBOSE) {
1018 >                printf("%s: data format is %s\n", progname, file_type[typ1]);
1019 >                if ((typ1 == TYP_FLOAT) | (typ1 == TYP_DOUBLE)) {
1020 >                        if (f1swap)
1021 >                                printf("%s: input '%s' is byte-swapped\n",
1022 >                                                progname, f1name);
1023 >                        if (f2swap)
1024 >                                printf("%s: input '%s' is byte-swapped\n",
1025 >                                                progname, f2name);
1026 >                }
1027 >        }
1028          switch (typ1) {                         /* compare based on type */
1029          case TYP_BINARY:
1030          case TYP_TMESH:
1031          case TYP_OCTREE:
1032          case TYP_RBFMESH:
1033 +        case TYP_ID8:
1034 +        case TYP_ID16:
1035 +        case TYP_ID24:
1036          case TYP_UNKNOWN:
1037                  return( !compare_binary() );
1038          case TYP_TEXT:
# Line 803 | Line 1041 | main(int argc, char *argv[])
1041          case TYP_RGBE:
1042          case TYP_XYZE:
1043                  return( !compare_hdr() );
1044 +        case TYP_DEPTH:
1045 +                return( !compare_depth() );
1046 +        case TYP_NORM:
1047 +                return( !compare_norm() );
1048          case TYP_FLOAT:
1049                  return( !compare_float() );
1050          case TYP_DOUBLE:

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines