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

Comparing ray/src/common/testBSDF.c (file contents):
Revision 1.2 by greg, Fri Jun 5 18:58:37 2015 UTC vs.
Revision 1.15 by greg, Sat Dec 28 18:05:14 2019 UTC

# Line 8 | Line 8 | static const char RCSid[] = "$Id$";
8   */
9  
10   #define _USE_MATH_DEFINES
11 – #include <stdio.h>
11   #include <stdlib.h>
12   #include <math.h>
13   #include <ctype.h>
# Line 20 | Line 19 | Usage(const char *prog)
19   {
20          printf("Usage: %s [bsdf_directory]\n", prog);
21          printf("Input commands:\n");
22 <        printf("  l bsdf.xml\t\t\t Load (make active) given BSDF input file\n");
23 <        printf("  q theta_i phi_i theta_o phi_o\t Query BSDF for given path\n");
24 <        printf("  s N theta phi\t\t\t Generate N ray directions at given incidence\n");
25 <        printf("  h theta phi\t\t\t Report hemispherical total at given incidence\n");
26 <        printf("  r theta phi\t\t\t Report hemispherical reflection at given incidence\n");
27 <        printf("  t theta phi\t\t\t Report hemispherical transmission at given incidence\n");
22 >        printf("  L bsdf.xml\t\t\t Load (make active) given BSDF input file\n");
23 >        printf("  i\t\t\t\t Report general information (metadata)\n");
24 >        printf("  c\t\t\t\t Report diffuse and specular components\n");
25 >        printf("  q theta_i phi_i theta_o phi_o\t Query BSDF for given path (CIE-XYZ)\n");
26 >        printf("  s[r|t][s|d] N theta phi\t Generate N ray directions & colors at given incidence\n");
27 >        printf("  h[s|d] theta phi\t\t Report hemispherical scattering at given incidence\n");
28 >        printf("  r[s|d] theta phi\t\t Report hemispherical reflection at given incidence\n");
29 >        printf("  t[s|d] theta phi\t\t Report hemispherical transmission at given incidence\n");
30          printf("  a theta phi [t2 p2]\t\t Report resolution (in proj. steradians) for given direction(s)\n");
31          printf("  ^D\t\t\t\t Quit program\n");
32   }
# Line 42 | Line 43 | vec_from_deg(FVECT v, double theta, double phi)
43          v[2] = cos(theta);
44   }
45  
46 + static void
47 + printXYZ(const char *intro, const SDValue *vp)
48 + {
49 +        if (vp->cieY <= 1e-9) {
50 +                printf("%s0 0 0\n", intro);
51 +                return;
52 +        }
53 +        printf("%s%.3e %.3e %.3e\n", intro,
54 +                        vp->spec.cx/vp->spec.cy*vp->cieY,
55 +                        vp->cieY,
56 +                        (1.-vp->spec.cx-vp->spec.cy)/
57 +                        vp->spec.cy*vp->cieY);
58 + }
59 +
60   int
61   main(int argc, char *argv[])
62   {
63 <        const char      *directory = "";
63 >        const char      *directory = NULL;
64          char            inp[512], path[512];
65          const SDData    *bsdf = NULL;
66  
# Line 56 | Line 71 | main(int argc, char *argv[])
71          if (argc == 2)
72                  directory = argv[1];
73          
74 <        SDretainSet = SDretainBSDFs;            /* keep BSDFs loaded */
74 >        SDretainSet = SDretainBSDFs;            /* keep BSDFs in memory */
75  
76                                                  /* loop on command */
77 <        while (fgets(inp, sizeof(inp), stdin) != NULL) {
77 >        while (fgets(inp, sizeof(inp), stdin)) {
78                  int     sflags = SDsampAll;
79                  char    *cp = inp;
80                  char    *cp2;
81                  FVECT   vin, vout;
82                  double  proja[2];
83 <                int     n;
83 >                int     n, i;
84                  SDValue val;
85                  
86                  while (isspace(*cp)) cp++;
87  
88 <                switch (*cp) {
74 <                case 'l':
88 >                switch (toupper(*cp)) {
89                  case 'L':                       /* load/activate BSDF input */
90                          cp2 = cp = sskip2(cp, 1);
91                          if (!*cp)
92                                  break;
93                          while (*cp) cp++;
94                          while (isspace(*--cp)) *cp = '\0';
95 <                        if (directory[0])
95 >                        if (directory)
96                                  sprintf(path, "%s/%s", directory, cp2);
97                          else
98                                  strcpy(path, cp2);
99 <                        if (bsdf != NULL)
99 >                        if (bsdf)
100                                  SDfreeCache(bsdf);
101                          bsdf = SDcacheFile(path);
102                          continue;
103 <                case 'q':
103 >                case 'I':                       /* report general info. */
104 >                        if (!bsdf)
105 >                                goto noBSDFerr;
106 >                        printf("Material: '%s'\n", bsdf->matn);
107 >                        printf("Manufacturer: '%s'\n", bsdf->makr);
108 >                        printf("Width, Height, Thickness (m): %.4e, %.4e, %.4e\n",
109 >                                        bsdf->dim[0], bsdf->dim[1], bsdf->dim[2]);
110 >                        if (bsdf->mgf)
111 >                                printf("Has geometry: %lu bytes\n",
112 >                                                (unsigned long)strlen(bsdf->mgf));
113 >                        else
114 >                                printf("Has geometry: no\n");
115 >                        break;
116 >                case 'C':                       /* report constant values */
117 >                        if (!bsdf)
118 >                                goto noBSDFerr;
119 >                        if (bsdf->rf)
120 >                                printf("Peak front hemispherical reflectance: %.3e\n",
121 >                                                bsdf->rLambFront.cieY +
122 >                                                bsdf->rf->maxHemi);
123 >                        if (bsdf->rb)
124 >                                printf("Peak back hemispherical reflectance: %.3e\n",
125 >                                                bsdf->rLambBack.cieY +
126 >                                                bsdf->rb->maxHemi);
127 >                        if (bsdf->tf)
128 >                                printf("Peak front hemispherical transmittance: %.3e\n",
129 >                                                bsdf->tLamb.cieY + bsdf->tf->maxHemi);
130 >                        if (bsdf->tb)
131 >                                printf("Peak back hemispherical transmittance: %.3e\n",
132 >                                                bsdf->tLamb.cieY + bsdf->tb->maxHemi);
133 >                        printXYZ("Diffuse Front Reflectance: ", &bsdf->rLambFront);
134 >                        printXYZ("Diffuse Back Reflectance: ", &bsdf->rLambBack);
135 >                        printXYZ("Diffuse Transmittance: ", &bsdf->tLamb);
136 >                        break;
137                  case 'Q':                       /* query BSDF value */
138 <                        if (bsdf == NULL)
138 >                        if (!bsdf)
139                                  goto noBSDFerr;
140                          if (!*sskip2(cp,4))
141                                  break;
142                          vec_from_deg(vin, atof(sskip2(cp,1)), atof(sskip2(cp,2)));
143                          vec_from_deg(vout, atof(sskip2(cp,3)), atof(sskip2(cp,4)));
144                          if (!SDreportError(SDevalBSDF(&val, vout, vin, bsdf), stderr))
145 <                                printf("%.3e\n", val.cieY);
146 <                        continue;
100 <                case 's':
145 >                                printXYZ("", &val);
146 >                        break;
147                  case 'S':                       /* sample BSDF */
148 <                        if (bsdf == NULL)
148 >                        if (!bsdf)
149                                  goto noBSDFerr;
150                          if (!*sskip2(cp,3))
151                                  break;
152 <                        n = atoi(sskip2(cp,1));
152 >                        if (toupper(cp[1]) == 'R') {
153 >                                sflags &= ~SDsampT;
154 >                                ++cp;
155 >                        } else if (toupper(cp[1]) == 'T') {
156 >                                sflags &= ~SDsampR;
157 >                                ++cp;
158 >                        }
159 >                        if (toupper(cp[1]) == 'S')
160 >                                sflags &= ~SDsampDf;
161 >                        else if (toupper(cp[1]) == 'D')
162 >                                sflags &= ~SDsampSp;
163 >                        i = n = atoi(sskip2(cp,1));
164                          vec_from_deg(vin, atof(sskip2(cp,2)), atof(sskip2(cp,3)));
165 <                        while (n-- > 0) {
166 <                                if (SDreportError(SDsampBSDF(&val, vin,
167 <                                                rand()*(1./(RAND_MAX+.5)),
165 >                        while (i-- > 0) {
166 >                                VCOPY(vout, vin);
167 >                                if (SDreportError(SDsampBSDF(&val, vout,
168 >                                                (i+rand()*(1./(RAND_MAX+.5)))/(double)n,
169                                                  sflags, bsdf), stderr))
170                                          break;
171 <                                printf("%.8f %.8f %.8f\n", vin[0], vin[1], vin[2]);
171 >                                printf("%.8f %.8f %.8f ", vout[0], vout[1], vout[2]);
172 >                                printXYZ("", &val);
173                          }
174 <                        continue;
175 <                case 'h':
117 <                case 'H':                       /* hemispherical totals */
118 <                case 'r':
174 >                        break;
175 >                case 'H':                       /* hemispherical values */
176                  case 'R':
120 –                case 't':
177                  case 'T':
178 <                        if (bsdf == NULL)
178 >                        if (!bsdf)
179                                  goto noBSDFerr;
180                          if (!*sskip2(cp,2))
181                                  break;
182 <                        if (tolower(*cp) == 'r')
183 <                                sflags ^= SDsampT;
184 <                        else if (tolower(*cp) == 't')
185 <                                sflags ^= SDsampR;
182 >                        if (toupper(cp[0]) == 'R')
183 >                                sflags &= ~SDsampT;
184 >                        else if (toupper(cp[0]) == 'T')
185 >                                sflags &= ~SDsampR;
186 >                        if (toupper(cp[1]) == 'S')
187 >                                sflags &= ~SDsampDf;
188 >                        else if (toupper(cp[1]) == 'D')
189 >                                sflags &= ~SDsampSp;
190                          vec_from_deg(vin, atof(sskip2(cp,1)), atof(sskip2(cp,2)));
191                          printf("%.4e\n", SDdirectHemi(vin, sflags, bsdf));
192 <                        continue;
193 <                case 'a':
194 <                case 'A':                       /* resolution in degrees */
135 <                        if (bsdf == NULL)
192 >                        break;
193 >                case 'A':                       /* resolution in proj. steradians */
194 >                        if (!bsdf)
195                                  goto noBSDFerr;
196                          if (!*sskip2(cp,2))
197                                  break;
198                          vec_from_deg(vin, atof(sskip2(cp,1)), atof(sskip2(cp,2)));
199                          if (*sskip2(cp,4)) {
200                                  vec_from_deg(vout, atof(sskip2(cp,3)), atof(sskip2(cp,4)));
201 <                                if (SDreportError(SDsizeBSDF(proja, vin, vout,
201 >                                if (SDreportError(SDsizeBSDF(proja, vout, vin,
202                                                  SDqueryMin+SDqueryMax, bsdf), stderr))
203                                          continue;
204                          } else if (SDreportError(SDsizeBSDF(proja, vin, NULL,
205                                                  SDqueryMin+SDqueryMax, bsdf), stderr))
206                                          continue;
207                          printf("%.4e %.4e\n", proja[0], proja[1]);
208 <                        continue;
208 >                        break;
209 >                default:
210 >                        Usage(argv[0]);
211 >                        break;
212                  }
213 <                Usage(argv[0]);
213 >                fflush(stdout);                 /* in case we're on remote */
214                  continue;
215   noBSDFerr:
216 <                fprintf(stderr, "%s: need to use 'l' command to load BSDF\n", argv[0]);
216 >                fprintf(stderr, "%s: First, use 'L' command to load BSDF\n", argv[0]);
217          }
218          return 0;
219   }

Diff Legend

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