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

Comparing ray/src/cv/bsdf2klems.c (file contents):
Revision 2.4 by greg, Thu Apr 25 03:12:08 2013 UTC vs.
Revision 2.9 by greg, Sun Aug 11 14:32:34 2013 UTC

# Line 174 | Line 174 | eval_bsdf(const char *fname)
174                      sum = 0;                    /* average over patches */
175                      for (n = npsamps; n-- > 0; ) {
176                          fo_getvec(vout, j+(n+frandom())/npsamps, abp);
177 <                        fi_getvec(vin, i+(n+frandom())/npsamps, abp);
177 >                        fi_getvec(vin, i+urand(n), abp);
178                          ec = SDevalBSDF(&sv, vout, vin, &bsd);
179                          if (ec != SDEnone)
180                                  goto err;
# Line 195 | Line 195 | eval_bsdf(const char *fname)
195                      sum = 0;                    /* average over patches */
196                      for (n = npsamps; n-- > 0; ) {
197                          bo_getvec(vout, j+(n+frandom())/npsamps, abp);
198 <                        bi_getvec(vin, i+(n+frandom())/npsamps, abp);
198 >                        bi_getvec(vin, i+urand(n), abp);
199                          ec = SDevalBSDF(&sv, vout, vin, &bsd);
200                          if (ec != SDEnone)
201                                  goto err;
# Line 216 | Line 216 | eval_bsdf(const char *fname)
216                      sum = 0;                    /* average over patches */
217                      for (n = npsamps; n-- > 0; ) {
218                          bo_getvec(vout, j+(n+frandom())/npsamps, abp);
219 <                        fi_getvec(vin, i+(n+frandom())/npsamps, abp);
219 >                        fi_getvec(vin, i+urand(n), abp);
220                          ec = SDevalBSDF(&sv, vout, vin, &bsd);
221                          if (ec != SDEnone)
222                                  goto err;
# Line 229 | Line 229 | eval_bsdf(const char *fname)
229              data_epilogue();
230          }
231                                                  /* back transmission */
232 <        if (bsd.tb != NULL) {
232 >        if ((bsd.tb != NULL) | (bsd.tf != NULL)) {
233              input_orient = -1; output_orient = 1;
234              data_prologue();
235              for (j = 0; j < abp->nangles; j++) {
236                  for (i = 0; i < abp->nangles; i++) {
237 <                    sum = 0;                    /* average over patches */
237 >                    sum = 0;            /* average over patches */
238                      for (n = npsamps; n-- > 0; ) {
239                          fo_getvec(vout, j+(n+frandom())/npsamps, abp);
240 <                        bi_getvec(vin, i+(n+frandom())/npsamps, abp);
240 >                        bi_getvec(vin, i+urand(n), abp);
241                          ec = SDevalBSDF(&sv, vout, vin, &bsd);
242                          if (ec != SDEnone)
243                                  goto err;
# Line 261 | Line 261 | static void
261   eval_function(char *funame)
262   {
263          ANGLE_BASIS     *abp = get_basis(kbasis);
264 +        int             assignD = (fundefined(funame) < 6);
265          double          iovec[6];
266          double          sum;
267          int             i, j, n;
# Line 281 | Line 282 | eval_function(char *funame)
282                      else
283                          bi_getvec(iovec, i+urand(n), abp);
284  
285 +                    if (assignD) {
286 +                        varset("Dx", '=', -iovec[3]);
287 +                        varset("Dy", '=', -iovec[4]);
288 +                        varset("Dz", '=', -iovec[5]);
289 +                        ++eclock;
290 +                    }
291                      sum += funvalue(funame, 6, iovec);
292                  }
293                  printf("\t%.3e\n", sum/npsamps);
# Line 326 | Line 333 | eval_rbf(void)
333                  }
334                  bsdfarr[j*abp->nangles + i] = sum*output_orient/npsamps;
335              }
336 +            if (rbf != NULL)
337 +                free(rbf);
338          }
339          n = 0;                                  /* write out our matrix */
340          for (j = 0; j < abp->nangles; j++) {
# Line 384 | Line 393 | main(int argc, char *argv[])
393                          fprintf(stderr,
394          "%s: need single function with 6 arguments: bsdf(ix,iy,iz,ox,oy,oz)\n",
395                                          progname);
396 +                        fprintf(stderr, "\tor 3 arguments using Dx,Dy,Dz: bsdf(ix,iy,iz)\n",
397 +                                        progname);
398                          goto userr;
399                  }
400 +                ++eclock;
401                  xml_header(argc, argv);                 /* start XML output */
402                  xml_prologue(NULL);
403                  if (dofwd) {

Diff Legend

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