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.8 by greg, Thu Feb 2 00:34:55 2017 UTC vs.
Revision 1.12 by greg, Mon May 15 22:15:40 2017 UTC

# Line 24 | Line 24 | Usage(const char *prog)
24          printf("  i\t\t\t\t Report general information (metadata)\n");
25          printf("  c\t\t\t\t Report diffuse and specular components\n");
26          printf("  q theta_i phi_i theta_o phi_o\t Query BSDF for given path (CIE-XYZ)\n");
27 <        printf("  s N theta phi\t\t\t Generate N ray directions at given incidence\n");
28 <        printf("  h theta phi\t\t\t Report hemispherical total at given incidence\n");
29 <        printf("  r theta phi\t\t\t Report hemispherical reflection at given incidence\n");
30 <        printf("  t theta phi\t\t\t Report hemispherical transmission at given incidence\n");
27 >        printf("  s[r|t][s|d] N theta phi\t Generate N ray directions & colors at given incidence\n");
28 >        printf("  h[s|d] theta phi\t\t Report hemispherical scattering at given incidence\n");
29 >        printf("  r[s|d] theta phi\t\t Report hemispherical reflection at given incidence\n");
30 >        printf("  t[s|d] theta phi\t\t Report hemispherical transmission at given incidence\n");
31          printf("  a theta phi [t2 p2]\t\t Report resolution (in proj. steradians) for given direction(s)\n");
32          printf("  ^D\t\t\t\t Quit program\n");
33   }
# Line 47 | Line 47 | vec_from_deg(FVECT v, double theta, double phi)
47   static void
48   printXYZ(const char *intro, const SDValue *vp)
49   {
50 +        if (vp->cieY <= 1e-9) {
51 +                printf("%s0 0 0\n", intro);
52 +                return;
53 +        }
54          printf("%s%.3e %.3e %.3e\n", intro,
55                          vp->spec.cx/vp->spec.cy*vp->cieY,
56                          vp->cieY,
# Line 77 | Line 81 | main(int argc, char *argv[])
81                  char    *cp2;
82                  FVECT   vin, vout;
83                  double  proja[2];
84 <                int     n;
84 >                int     n, i;
85                  SDValue val;
86                  
87                  while (isspace(*cp)) cp++;
# Line 142 | Line 146 | main(int argc, char *argv[])
146                                  goto noBSDFerr;
147                          if (!*sskip2(cp,3))
148                                  break;
149 <                        n = atoi(sskip2(cp,1));
149 >                        if (toupper(cp[1]) == 'R') {
150 >                                sflags &= ~SDsampT;
151 >                                ++cp;
152 >                        } else if (toupper(cp[1]) == 'T') {
153 >                                sflags &= ~SDsampR;
154 >                                ++cp;
155 >                        }
156 >                        if (toupper(cp[1]) == 'S')
157 >                                sflags &= ~SDsampDf;
158 >                        else if (toupper(cp[1]) == 'D')
159 >                                sflags &= ~SDsampSp;
160 >                        i = n = atoi(sskip2(cp,1));
161                          vec_from_deg(vin, atof(sskip2(cp,2)), atof(sskip2(cp,3)));
162 <                        while (n-- > 0) {
163 <                                if (SDreportError(SDsampBSDF(&val, vin,
164 <                                                rand()*(1./(RAND_MAX+.5)),
162 >                        while (i-- > 0) {
163 >                                VCOPY(vout, vin);
164 >                                if (SDreportError(SDsampBSDF(&val, vout,
165 >                                                (i+rand()*(1./(RAND_MAX+.5)))/(double)n,
166                                                  sflags, bsdf), stderr))
167                                          break;
168 <                                printf("%.8f %.8f %.8f\n", vin[0], vin[1], vin[2]);
168 >                                printf("%.8f %.8f %.8f ", vout[0], vout[1], vout[2]);
169 >                                printXYZ("", &val);
170                          }
171                          continue;
172 <                case 'H':                       /* hemispherical totals */
172 >                case 'H':                       /* hemispherical values */
173                  case 'R':
174                  case 'T':
175                          if (!bsdf)
176                                  goto noBSDFerr;
177                          if (!*sskip2(cp,2))
178                                  break;
179 <                        if (tolower(*cp) == 'r')
179 >                        if (toupper(cp[0]) == 'R')
180                                  sflags &= ~SDsampT;
181 <                        else if (tolower(*cp) == 't')
181 >                        else if (toupper(cp[0]) == 'T')
182                                  sflags &= ~SDsampR;
183 +                        if (toupper(cp[1]) == 'S')
184 +                                sflags &= ~SDsampDf;
185 +                        else if (toupper(cp[1]) == 'D')
186 +                                sflags &= ~SDsampSp;
187                          vec_from_deg(vin, atof(sskip2(cp,1)), atof(sskip2(cp,2)));
188                          printf("%.4e\n", SDdirectHemi(vin, sflags, bsdf));
189                          continue;
# Line 174 | Line 195 | main(int argc, char *argv[])
195                          vec_from_deg(vin, atof(sskip2(cp,1)), atof(sskip2(cp,2)));
196                          if (*sskip2(cp,4)) {
197                                  vec_from_deg(vout, atof(sskip2(cp,3)), atof(sskip2(cp,4)));
198 <                                if (SDreportError(SDsizeBSDF(proja, vin, vout,
198 >                                if (SDreportError(SDsizeBSDF(proja, vout, vin,
199                                                  SDqueryMin+SDqueryMax, bsdf), stderr))
200                                          continue;
201                          } else if (SDreportError(SDsizeBSDF(proja, vin, NULL,

Diff Legend

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