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

Comparing ray/src/cv/bsdfquery.c (file contents):
Revision 2.1 by greg, Thu Nov 21 23:57:28 2013 UTC vs.
Revision 2.7 by greg, Thu Aug 18 00:52:48 2016 UTC

# Line 14 | Line 14 | static const char RCSid[] = "$Id$";
14   #include <string.h>
15   #include <stdlib.h>
16   #include "rtmath.h"
17 + #include "rtio.h"
18   #include "bsdfrep.h"
19  
20   char    *progname;
# Line 27 | Line 28 | readIOdir(FVECT idir, FVECT odir, FILE *fp, int fmt)
28  
29          switch (fmt) {
30          case 'a':
31 <                if (fscanf(fp, "%lf %lf %lf %lf %lf %lf",
32 <                        dvec, dvec+1, dvec+2, dvec+3, dvec+4, dvec+5) != 6)
31 >                if (fscanf(fp, FVFORMAT, &idir[0], &idir[1], &idir[2]) != 3 ||
32 >                                fscanf(fp, FVFORMAT, &odir[0], &odir[1], &odir[2]) != 3)
33                          return(0);
33 –                VCOPY(idir, dvec);
34 –                VCOPY(odir, dvec+3);
34                  break;
35          case 'd':
36 <                if (fread(dvec, sizeof(double), 6, fp) != 6)
36 >                if (getbinary(dvec, sizeof(double), 6, fp) != 6)
37                          return(0);
38                  VCOPY(idir, dvec);
39                  VCOPY(odir, dvec+3);
40                  break;
41          case 'f':
42 <                if (fread(fvec, sizeof(float), 6, fp) != 6)
42 >                if (getbinary(fvec, sizeof(float), 6, fp) != 6)
43                          return(0);
44                  VCOPY(idir, fvec);
45                  VCOPY(odir, fvec+3);
# Line 57 | Line 56 | readIOdir(FVECT idir, FVECT odir, FILE *fp, int fmt)
56   int
57   main(int argc, char *argv[])
58   {
59 +        int     repXYZ = 0;
60          int     inpXML = -1;
61          int     inpfmt = 'a';
62          int     outfmt = 'a';
# Line 67 | Line 67 | main(int argc, char *argv[])
67          int     n;
68                                                  /* check arguments */
69          progname = argv[0];
70 <        if (argc > 2 && argv[1][0] == '-' && argv[1][1] == 'f' &&
71 <                        argv[1][2] && strchr("afd", argv[1][2]) != NULL) {
72 <                inpfmt = outfmt = argv[1][2];
73 <                if (argv[1][3] && strchr("afd", argv[1][3]) != NULL)
74 <                        outfmt = argv[1][3];
70 >        while (argc > 2 && argv[1][0] == '-') {
71 >                switch (argv[1][1]) {
72 >                case 'c':                       /* color output */
73 >                        repXYZ = 1;
74 >                        break;
75 >                case 'f':                       /* i/o format */
76 >                        if (!argv[1][2] || strchr("afd", argv[1][2]) == NULL)
77 >                                goto userr;
78 >                        inpfmt = outfmt = argv[1][2];
79 >                        if (argv[1][3] && strchr("afd", argv[1][3]) != NULL)
80 >                                outfmt = argv[1][3];
81 >                        break;
82 >                default:
83 >                        goto userr;
84 >                }
85                  ++argv; --argc;
86          }
87 <        if ((n = strlen(argv[1])-4) > 0) {
87 >        if (argc > 1 && (n = strlen(argv[1])-4) > 0) {
88                  if (!strcasecmp(argv[1]+n, ".xml"))
89                          inpXML = 1;
90                  else if (!strcasecmp(argv[1]+n, ".sir"))
91                          inpXML = 0;
92          }
93 <        if (inpXML < 0) {
94 <                fprintf(stderr, "Usage: %s [-fio] bsdf.{sir|xml}\n", progname);
85 <                return(1);
86 <        }
93 >        if ((argc != 2) | (inpXML < 0))
94 >                goto userr;
95                                                  /* load BSDF representation */
96          if (inpXML) {
97                  SDclearBSDF(&myBSDF, argv[1]);
# Line 102 | Line 110 | main(int argc, char *argv[])
110          }
111                                                  /* query BSDF values */
112          while (readIOdir(idir, odir, stdin, inpfmt)) {
113 <                double  bsdf;
106 <                float   fval;
113 >                SDValue sval;
114                  if (inpXML) {
108 –                        SDValue sval;
115                          if (SDreportError(SDevalBSDF(&sval, odir,
116                                                  idir, &myBSDF), stderr))
117                                  return(1);
112 –                        bsdf = sval.cieY;
118                  } else {
119                          int32   inpDir = encodedir(idir);
120                          if (inpDir != prevInpDir) {
# Line 127 | Line 132 | main(int argc, char *argv[])
132                                                  progname);
133                                  return(1);
134                          }
135 <                        bsdf = eval_rbfrep(rbf, odir)/(output_orient*odir[2]);
135 >                        if (SDreportError(eval_rbfcol(&sval, rbf, odir), stderr))
136 >                                return(1);
137                  }
138 +                if (repXYZ)                     /* ensure we have CIE (x,y) */
139 +                        c_ccvt(&sval.spec, C_CSXY);
140 +
141                  switch (outfmt) {               /* write to stdout */
142                  case 'a':
143 <                        printf("%.6e\n", bsdf);
143 >                        if (repXYZ) {
144 >                                double  cieX = sval.spec.cx/sval.spec.cy*sval.cieY;
145 >                                double  cieZ = (1. - sval.spec.cx - sval.spec.cy) /
146 >                                                sval.spec.cy * sval.cieY;
147 >                                printf("%.6e %.6e %.6e\n", cieX, sval.cieY, cieZ);
148 >                        } else
149 >                                printf("%.6e\n", sval.cieY);
150                          break;
151                  case 'd':
152 <                        fwrite(&bsdf, sizeof(double), 1, stdout);
152 >                        if (repXYZ) {
153 >                                double  cieXYZ[3];
154 >                                cieXYZ[0] = sval.spec.cx/sval.spec.cy*sval.cieY;
155 >                                cieXYZ[1] = sval.cieY;
156 >                                cieXYZ[2] = (1. - sval.spec.cx - sval.spec.cy) /
157 >                                                sval.spec.cy * sval.cieY;
158 >                                putbinary(cieXYZ, sizeof(double), 3, stdout);
159 >                        } else
160 >                                putbinary(&sval.cieY, sizeof(double), 1, stdout);
161                          break;
162                  case 'f':
163 <                        fval = bsdf;
164 <                        fwrite(&fval, sizeof(float), 1, stdout);
163 >                        if (repXYZ) {
164 >                                float   cieXYZ[3];
165 >                                cieXYZ[0] = sval.spec.cx/sval.spec.cy*sval.cieY;
166 >                                cieXYZ[1] = sval.cieY;
167 >                                cieXYZ[2] = (1. - sval.spec.cx - sval.spec.cy) /
168 >                                                sval.spec.cy * sval.cieY;
169 >                                putbinary(cieXYZ, sizeof(float), 3, stdout);
170 >                        } else {
171 >                                float   cieY = sval.cieY;
172 >                                putbinary(&cieY, sizeof(float), 1, stdout);
173 >                        }
174                          break;
175                  }
176          }
177 +        /* if (rbf != NULL) free(rbf); */
178          return(0);
179 + userr:
180 +        fprintf(stderr, "Usage: %s [-c][-fio] bsdf.{sir|xml}\n", progname);
181 +        return(1);
182   }

Diff Legend

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