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

Comparing ray/src/cv/bsdf2rad.c (file contents):
Revision 2.22 by greg, Tue Apr 11 02:21:37 2017 UTC vs.
Revision 2.26 by greg, Wed Apr 12 03:35:01 2017 UTC

# Line 26 | Line 26 | double overall_max = .0;               /* overall maximum BSDF valu
26  
27   char    ourTempDir[TEMPLEN] = "";       /* our temporary directory */
28  
29 < const char      frpref[] = "frefl";
30 < const char      ftpref[] = "ftrans";
31 < const char      brpref[] = "brefl";
32 < const char      btpref[] = "btrans";
29 > const char      frpref[] = "rf";
30 > const char      ftpref[] = "tf";
31 > const char      brpref[] = "rb";
32 > const char      btpref[] = "tb";
33   const char      dsuffix[] = ".txt";
34  
35 < const char      sph_mat[] = "BSDFmat";
35 > const char      sph_fmat[] = "fBSDFmat";
36 > const char      sph_bmat[] = "bBSDFmat";
37   const double    sph_rad = 10.;
38   const double    sph_xoffset = 15.;
39  
# Line 269 | Line 270 | build_wRBF(void)
270  
271   /* Put out mirror arrow for the given incident vector */
272   static void
273 < put_mirror_arrow(const FVECT ivec, int inc_side)
273 > put_mirror_arrow(const FVECT origin, const FVECT nrm)
274   {
275          const double    arrow_len = 1.2*bsdf_rad;
276          const double    tip_len = 0.2*bsdf_rad;
277 <        FVECT           origin, refl;
277 >        FVECT           refl;
278          int             i;
279  
280 <        cvt_sposition(origin, ivec, inc_side);
280 >        refl[0] = 2.*nrm[2]*nrm[0];
281 >        refl[1] = 2.*nrm[2]*nrm[1];
282 >        refl[2] = 2.*nrm[2]*nrm[2] - 1.;
283  
281 –        refl[0] = -2.*ivec[2]*ivec[0];
282 –        refl[1] = -2.*ivec[2]*ivec[1];
283 –        refl[2] = 2.*ivec[2]*ivec[2] - 1.;
284 –
284          printf("\n# Mirror arrow\n");
285          printf("\narrow_mat cylinder inc_dir\n0\n0\n7");
286          printf("\n\t%f %f %f\n\t%f %f %f\n\t%f\n",
# Line 308 | Line 307 | put_mirror_arrow(const FVECT ivec, int inc_side)
307  
308   /* Put out transmitted direction arrow for the given incident vector */
309   static void
310 < put_trans_arrow(const FVECT ivec, int inc_side)
310 > put_trans_arrow(const FVECT origin)
311   {
312          const double    arrow_len = 1.2*bsdf_rad;
313          const double    tip_len = 0.2*bsdf_rad;
315 –        FVECT           origin;
314          int             i;
315  
318 –        cvt_sposition(origin, ivec, inc_side);
319 –
316          printf("\n# Transmission arrow\n");
317          printf("\narrow_mat cylinder trans_dir\n0\n0\n7");
318          printf("\n\t%f %f %f\n\t%f %f %f\n\t%f\n",
# Line 371 | Line 367 | static int
367   put_BSDFs(void)
368   {
369          const double    scalef = bsdf_rad/(log10(overall_max) - min_log10);
370 <        FVECT           ivec, sorg, upv;
370 >        FVECT           ivec, sorg, nrm, upv;
371          RREAL           vMtx[3][3];
372          char            *fname;
373          char            cmdbuf[256];
374 <        char            xfargs[128];
375 <        int             nxfa;
374 >        char            rotargs[64];
375 >        int             nrota;
376          int             i;
377  
378          printf("\n# Gensurf output corresponding to %d incident directions\n",
379                          NINCIDENT);
380  
381          printf("\nvoid glow arrow_glow\n0\n0\n4 1 0 1 0\n");
382 <        printf("\nvoid mixfunc arrow_mat\n4 arrow_glow void 0.5 .\n0\n0\n");
382 >        printf("\nvoid mixfunc arrow_mat\n4 arrow_glow void 0.25 .\n0\n0\n");
383  
384 <        if (front_comp & SDsampR)                       /* front reflection */
385 <                for (i = 0; i < NINCIDENT; i++) {
386 <                        get_ivector(ivec, i);
387 <                        put_mirror_arrow(ivec, 1);
384 >        for (i = 0; i < NINCIDENT; i++) {
385 >                get_ivector(ivec, i);
386 >                nrm[0] = -ivec[0]; nrm[1] = -ivec[1]; nrm[2] = ivec[2];
387 >                upv[0] = nrm[0]*nrm[1]*(nrm[2] - 1.);
388 >                upv[1] = nrm[0]*nrm[0] + nrm[1]*nrm[1]*nrm[2];
389 >                upv[2] = -nrm[1]*(nrm[0]*nrm[0] + nrm[1]*nrm[1]);
390 >                if (SDcompXform(vMtx, nrm, upv) != SDEnone)
391 >                        continue;
392 >                nrota = addrot(rotargs, vMtx[0], vMtx[1], vMtx[2]);
393 >                if (front_comp) {
394                          cvt_sposition(sorg, ivec, 1);
393 –                        ivec[0] = -ivec[0]; ivec[1] = -ivec[1]; /* normal */
394 –                        upv[0] = ivec[0]*ivec[1]*(ivec[2] - 1.);
395 –                        upv[1] = ivec[0]*ivec[0] + ivec[1]*ivec[1]*ivec[2];
396 –                        upv[2] = -ivec[1]*(ivec[0]*ivec[0] + ivec[1]*ivec[1]);
397 –                        sprintf(xfargs, "-s %f -t %f %f %f", bsdf_rad,
398 –                                        sorg[0], sorg[1], sorg[2]);
399 –                        nxfa = 6;
395                          printf("\nvoid colorfunc scale_pat\n");
396 <                        printf("%d bsdf_red bsdf_grn bsdf_blu bsdf2rad.cal\n\t%s\n0\n0\n",
397 <                                        4+nxfa, xfargs);
396 >                        printf("10 bsdf_red bsdf_grn bsdf_blu bsdf2rad.cal\n");
397 >                        printf("\t-s %f -t %f %f %f\n0\n0\n",
398 >                                        bsdf_rad, sorg[0], sorg[1], sorg[2]);
399                          printf("\nscale_pat glow scale_mat\n0\n0\n4 1 1 1 0\n");
400 <                        SDcompXform(vMtx, ivec, upv);
401 <                        nxfa = addrot(xfargs, vMtx[0], vMtx[1], vMtx[2]);
402 <                        sprintf(xfargs+strlen(xfargs), " -s %f -t %f %f %f",
407 <                                        scalef, sorg[0], sorg[1], sorg[2]);
408 <                        nxfa += 6;
400 >                }
401 >                if (front_comp & SDsampR) {
402 >                        put_mirror_arrow(sorg, nrm);
403                          fname = tfile_name(frpref, dsuffix, i);
404 <                        sprintf(cmdbuf, "gensurf scale_mat %s%d %s %s %s %d %d | xform %s",
404 >                        sprintf(cmdbuf,
405 >                "gensurf scale_mat %s%d %s %s %s %d %d | xform %s -s %f -t %f %f %f",
406                                          frpref, i, fname, fname, fname, SAMPRES-1, SAMPRES-1,
407 <                                        xfargs);
407 >                                        rotargs, scalef, sorg[0], sorg[1], sorg[2]);
408                          if (!run_cmd(cmdbuf))
409                                  return(0);
410                  }
411 <        if (front_comp & SDsampT)                       /* front transmission */
412 <                for (i = 0; i < NINCIDENT; i++) {
418 <                        get_ivector(ivec, i);
419 <                        put_trans_arrow(ivec, 1);
420 <                        cvt_sposition(sorg, ivec, 1);
421 <                        ivec[0] = -ivec[0]; ivec[1] = -ivec[1]; /* normal */
422 <                        upv[0] = ivec[0]*ivec[1]*(ivec[2] - 1.);
423 <                        upv[1] = ivec[0]*ivec[0] + ivec[1]*ivec[1]*ivec[2];
424 <                        upv[2] = -ivec[1]*(ivec[0]*ivec[0] + ivec[1]*ivec[1]);
425 <                        sprintf(xfargs, "-s %f -t %f %f %f", bsdf_rad,
426 <                                        sorg[0], sorg[1], sorg[2]);
427 <                        nxfa = 6;
428 <                        printf("\nvoid colorfunc scale_pat\n");
429 <                        printf("%d bsdf_red bsdf_grn bsdf_blu bsdf2rad.cal\n\t%s\n0\n0\n",
430 <                                        4+nxfa, xfargs);
431 <                        printf("\nscale_pat glow scale_mat\n0\n0\n4 1 1 1 0\n");
432 <                        SDcompXform(vMtx, ivec, upv);
433 <                        nxfa = addrot(xfargs, vMtx[0], vMtx[1], vMtx[2]);
434 <                        sprintf(xfargs+strlen(xfargs), " -s %f -t %f %f %f",
435 <                                        scalef, sorg[0], sorg[1], sorg[2]);
436 <                        nxfa += 6;
411 >                if (front_comp & SDsampT) {
412 >                        put_trans_arrow(sorg);
413                          fname = tfile_name(ftpref, dsuffix, i);
414 <                        sprintf(cmdbuf, "gensurf scale_mat %s%d %s %s %s %d %d | xform -I %s",
414 >                        sprintf(cmdbuf,
415 >                "gensurf scale_mat %s%d %s %s %s %d %d | xform -I %s -s %f -t %f %f %f",
416                                          ftpref, i, fname, fname, fname, SAMPRES-1, SAMPRES-1,
417 <                                        xfargs);
417 >                                        rotargs, scalef, sorg[0], sorg[1], sorg[2]);
418                          if (!run_cmd(cmdbuf))
419                                  return(0);
420                  }
421 <        if (back_comp & SDsampR)                        /* rear reflection */
445 <                for (i = 0; i < NINCIDENT; i++) {
446 <                        get_ivector(ivec, i);
447 <                        put_mirror_arrow(ivec, -1);
421 >                if (back_comp) {
422                          cvt_sposition(sorg, ivec, -1);
449 –                        ivec[0] = -ivec[0]; ivec[1] = -ivec[1]; /* normal */
450 –                        upv[0] = ivec[0]*ivec[1]*(ivec[2] - 1.);
451 –                        upv[1] = ivec[0]*ivec[0] + ivec[1]*ivec[1]*ivec[2];
452 –                        upv[2] = -ivec[1]*(ivec[0]*ivec[0] + ivec[1]*ivec[1]);
453 –                        sprintf(xfargs, "-s %f -t %f %f %f", bsdf_rad,
454 –                                        sorg[0], sorg[1], sorg[2]);
455 –                        nxfa = 6;
423                          printf("\nvoid colorfunc scale_pat\n");
424 <                        printf("%d bsdf_red bsdf_grn bsdf_blu bsdf2rad.cal\n\t%s\n0\n0\n",
425 <                                        4+nxfa, xfargs);
424 >                        printf("10 bsdf_red bsdf_grn bsdf_blu bsdf2rad.cal\n");
425 >                        printf("\t-s %f -t %f %f %f\n0\n0\n",
426 >                                        bsdf_rad, sorg[0], sorg[1], sorg[2]);
427                          printf("\nscale_pat glow scale_mat\n0\n0\n4 1 1 1 0\n");
428 <                        SDcompXform(vMtx, ivec, upv);
429 <                        nxfa = addrot(xfargs, vMtx[0], vMtx[1], vMtx[2]);
430 <                        sprintf(xfargs+strlen(xfargs), " -s %f -t %f %f %f",
463 <                                        scalef, sorg[0], sorg[1], sorg[2]);
464 <                        nxfa += 6;
428 >                }
429 >                if (back_comp & SDsampR) {
430 >                        put_mirror_arrow(sorg, nrm);
431                          fname = tfile_name(brpref, dsuffix, i);
432 <                        sprintf(cmdbuf, "gensurf scale_mat %s%d %s %s %s %d %d | xform -I -ry 180 %s",
432 >                        sprintf(cmdbuf,
433 >                "gensurf scale_mat %s%d %s %s %s %d %d | xform -I -ry 180 %s -s %f -t %f %f %f",
434                                          brpref, i, fname, fname, fname, SAMPRES-1, SAMPRES-1,
435 <                                        xfargs);
435 >                                        rotargs, scalef, sorg[0], sorg[1], sorg[2]);
436                          if (!run_cmd(cmdbuf))
437                                  return(0);
438                  }
439 <        if (back_comp & SDsampT)                        /* rear transmission */
440 <                for (i = 0; i < NINCIDENT; i++) {
474 <                        get_ivector(ivec, i);
475 <                        put_trans_arrow(ivec, -1);
476 <                        cvt_sposition(sorg, ivec, -1);
477 <                        ivec[0] = -ivec[0]; ivec[1] = -ivec[1]; /* normal */
478 <                        upv[0] = ivec[0]*ivec[1]*(ivec[2] - 1.);
479 <                        upv[1] = ivec[0]*ivec[0] + ivec[1]*ivec[1]*ivec[2];
480 <                        upv[2] = -ivec[1]*(ivec[0]*ivec[0] + ivec[1]*ivec[1]);
481 <                        sprintf(xfargs, "-s %f -t %f %f %f", bsdf_rad,
482 <                                        sorg[0], sorg[1], sorg[2]);
483 <                        nxfa = 6;
484 <                        printf("\nvoid colorfunc scale_pat\n");
485 <                        printf("%d bsdf_red bsdf_grn bsdf_blu bsdf2rad.cal\n\t%s\n0\n0\n",
486 <                                        4+nxfa, xfargs);
487 <                        printf("\nscale_pat glow scale_mat\n0\n0\n4 1 1 1 0\n");
488 <                        SDcompXform(vMtx, ivec, upv);
489 <                        nxfa = addrot(xfargs, vMtx[0], vMtx[1], vMtx[2]);
490 <                        sprintf(xfargs+strlen(xfargs), " -s %f -t %f %f %f",
491 <                                        scalef, sorg[0], sorg[1], sorg[2]);
492 <                        nxfa += 6;
439 >                if (back_comp & SDsampT) {
440 >                        put_trans_arrow(sorg);
441                          fname = tfile_name(btpref, dsuffix, i);
442 <                        sprintf(cmdbuf, "gensurf scale_mat %s%d %s %s %s %d %d | xform -ry 180 %s",
442 >                        sprintf(cmdbuf,
443 >                "gensurf scale_mat %s%d %s %s %s %d %d | xform -ry 180 %s -s %f -t %f %f %f",
444                                          btpref, i, fname, fname, fname, SAMPRES-1, SAMPRES-1,
445 <                                        xfargs);
445 >                                        rotargs, scalef, sorg[0], sorg[1], sorg[2]);
446                          if (!run_cmd(cmdbuf))
447                                  return(0);
448                  }
449 +        }
450          return(1);
451   }
452  
# Line 511 | Line 461 | put_matBSDF(const char *XMLfile)
461                  printf("\nvoid brightfunc latlong\n2 latlong bsdf2rad.cal\n0\n0\n");
462                  if ((front_comp|back_comp) & SDsampT)
463                          printf("\nlatlong trans %s\n0\n0\n7 .75 .75 .75 0 .04 .5 .8\n",
464 <                                        sph_mat);
464 >                                        sph_fmat);
465                  else
466                          printf("\nlatlong plastic %s\n0\n0\n5 .5 .5 .5 0 0\n",
467 <                                        sph_mat);
467 >                                        sph_fmat);
468 >                printf("\ninherit alias %s %s\n", sph_bmat, sph_fmat);
469                  return;
470          }
471          switch (XMLfile[0]) {           /* avoid RAYPATH search */
# Line 527 | Line 478 | put_matBSDF(const char *XMLfile)
478                  exit(1);
479                  break;
480          }
481 <        printf("\n# Actual BSDF material for rendering the hemispheres\n");
482 <        printf("\nvoid BSDF BSDFmat\n6 0 \"%s%s\" upx upy upz bsdf2rad.cal\n0\n0\n",
481 >        printf("\n# Actual BSDF materials for rendering the hemispheres\n");
482 >        printf("\nvoid BSDF BSDF_f\n6 0 \"%s%s\" upx upy upz bsdf2rad.cal\n0\n0\n",
483                          curdir, XMLfile);
484          printf("\nvoid plastic black\n0\n0\n5 0 0 0 0 0\n");
485 <        printf("\nvoid mixfunc %s\n4 BSDFmat black latlong bsdf2rad.cal\n0\n0\n",
486 <                        sph_mat);
485 >        printf("\nvoid mixfunc %s\n4 BSDF_f black latlong bsdf2rad.cal\n0\n0\n",
486 >                        sph_fmat);
487 >        printf("\nvoid BSDF BSDF_b\n8 0 \"%s%s\" upx upy upz bsdf2rad.cal -ry 180\n0\n0\n",
488 >                        curdir, XMLfile);
489 >        printf("\nvoid mixfunc %s\n4 BSDF_b black latlong bsdf2rad.cal\n0\n0\n",
490 >                        sph_bmat);
491   }
492  
493   /* Put out overhead parallel light source */
# Line 556 | Line 511 | put_hemispheres(void)
511          if (front_comp) {
512                  printf(
513   "\n!genrev %s Front \"R*sin(A*t)\" \"R*cos(A*t)\" %d -e \"R:%g;A:%f\" -s | xform -t %g 0 0\n",
514 <                                sph_mat, nsegs, sph_rad, 0.495*PI, sph_xoffset);
514 >                                sph_fmat, nsegs, sph_rad, 0.495*PI, sph_xoffset);
515                  printf("\nvoid brighttext front_text\n3 helvet.fnt . FRONT\n0\n");
516                  printf("12\n\t%f %f 0\n\t%f 0 0\n\t0 %f 0\n\t.01 1 -.1\n",
517                                  -.22*sph_rad + sph_xoffset, -1.4*sph_rad,
# Line 572 | Line 527 | put_hemispheres(void)
527          if (back_comp) {
528                  printf(
529   "\n!genrev %s Back \"R*cos(A*t)\" \"R*sin(A*t)\" %d -e \"R:%g;A:%f\" -s | xform -t %g 0 0\n",
530 <                                sph_mat, nsegs, sph_rad, 0.495*PI, -sph_xoffset);
530 >                                sph_bmat, nsegs, sph_rad, 0.495*PI, -sph_xoffset);
531                  printf("\nvoid brighttext back_text\n3 helvet.fnt . BACK\n0\n");
532                  printf("12\n\t%f %f 0\n\t%f 0 0\n\t0 %f 0\n\t.01 1 -.1\n",
533                                  -.22*sph_rad - sph_xoffset, -1.4*sph_rad,

Diff Legend

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