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.19 by greg, Mon Apr 10 01:31:37 2017 UTC vs.
Revision 2.22 by greg, Tue Apr 11 02:21:37 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 FVECT     Xaxis = {1., 0., 0.};
30 – const FVECT     Yaxis = {0., 1., 0.};
31 – const FVECT     Zaxis = {0., 0., 1.};
32 –
29   const char      frpref[] = "frefl";
30   const char      ftpref[] = "ftrans";
31   const char      brpref[] = "brefl";
# Line 68 | Line 64 | get_ivector(FVECT iv, int i)
64   static RREAL *
65   cvt_sposition(FVECT sp, const FVECT iv, int inc_side)
66   {
67 <        sp[0] = -iv[0]*sph_rad - inc_side*sph_xoffset;
67 >        sp[0] = -iv[0]*sph_rad + inc_side*sph_xoffset;
68          sp[1] = -iv[1]*sph_rad;
69          sp[2] = iv[2]*sph_rad;
70  
# Line 375 | Line 371 | static int
371   put_BSDFs(void)
372   {
373          const double    scalef = bsdf_rad/(log10(overall_max) - min_log10);
374 <        FVECT           ivec, sorg;
374 >        FVECT           ivec, sorg, upv;
375          RREAL           vMtx[3][3];
376          char            *fname;
377          char            cmdbuf[256];
# Line 387 | Line 383 | put_BSDFs(void)
383                          NINCIDENT);
384  
385          printf("\nvoid glow arrow_glow\n0\n0\n4 1 0 1 0\n");
386 <        printf("\nvoid mixfunc arrow_mat\n4 arrow_glow void .5 .\n0\n0\n");
386 >        printf("\nvoid mixfunc arrow_mat\n4 arrow_glow void 0.5 .\n0\n0\n");
387  
388          if (front_comp & SDsampR)                       /* front reflection */
389                  for (i = 0; i < NINCIDENT; i++) {
# Line 395 | Line 391 | put_BSDFs(void)
391                          put_mirror_arrow(ivec, 1);
392                          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;
# Line 402 | Line 401 | put_BSDFs(void)
401                          printf("%d bsdf_red bsdf_grn bsdf_blu bsdf2rad.cal\n\t%s\n0\n0\n",
402                                          4+nxfa, xfargs);
403                          printf("\nscale_pat glow scale_mat\n0\n0\n4 1 1 1 0\n");
404 <                        SDcompXform(vMtx, ivec, Yaxis);
404 >                        SDcompXform(vMtx, ivec, upv);
405                          nxfa = addrot(xfargs, vMtx[0], vMtx[1], vMtx[2]);
406                          sprintf(xfargs+strlen(xfargs), " -s %f -t %f %f %f",
407                                          scalef, sorg[0], sorg[1], sorg[2]);
408                          nxfa += 6;
409                          fname = tfile_name(frpref, dsuffix, i);
410                          sprintf(cmdbuf, "gensurf scale_mat %s%d %s %s %s %d %d | xform %s",
411 <                                        frpref, i+1, fname, fname, fname, SAMPRES-1, SAMPRES-1,
411 >                                        frpref, i, fname, fname, fname, SAMPRES-1, SAMPRES-1,
412                                          xfargs);
413                          if (!run_cmd(cmdbuf))
414                                  return(0);
# Line 420 | Line 419 | put_BSDFs(void)
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;
# Line 427 | Line 429 | put_BSDFs(void)
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, Yaxis);
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;
437                          fname = tfile_name(ftpref, dsuffix, i);
438                          sprintf(cmdbuf, "gensurf scale_mat %s%d %s %s %s %d %d | xform -I %s",
439 <                                        ftpref, i+1, fname, fname, fname, SAMPRES-1, SAMPRES-1,
439 >                                        ftpref, i, fname, fname, fname, SAMPRES-1, SAMPRES-1,
440                                          xfargs);
441                          if (!run_cmd(cmdbuf))
442                                  return(0);
# Line 445 | Line 447 | put_BSDFs(void)
447                          put_mirror_arrow(ivec, -1);
448                          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;
# Line 452 | Line 457 | put_BSDFs(void)
457                          printf("%d bsdf_red bsdf_grn bsdf_blu bsdf2rad.cal\n\t%s\n0\n0\n",
458                                          4+nxfa, xfargs);
459                          printf("\nscale_pat glow scale_mat\n0\n0\n4 1 1 1 0\n");
460 <                        SDcompXform(vMtx, ivec, Yaxis);
460 >                        SDcompXform(vMtx, ivec, upv);
461                          nxfa = addrot(xfargs, vMtx[0], vMtx[1], vMtx[2]);
462                          sprintf(xfargs+strlen(xfargs), " -s %f -t %f %f %f",
463                                          scalef, sorg[0], sorg[1], sorg[2]);
464                          nxfa += 6;
465                          fname = tfile_name(brpref, dsuffix, i);
466                          sprintf(cmdbuf, "gensurf scale_mat %s%d %s %s %s %d %d | xform -I -ry 180 %s",
467 <                                        brpref, i+1, fname, fname, fname, SAMPRES-1, SAMPRES-1,
467 >                                        brpref, i, fname, fname, fname, SAMPRES-1, SAMPRES-1,
468                                          xfargs);
469                          if (!run_cmd(cmdbuf))
470                                  return(0);
# Line 470 | Line 475 | put_BSDFs(void)
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;
# Line 477 | Line 485 | put_BSDFs(void)
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, Yaxis);
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;
493                          fname = tfile_name(btpref, dsuffix, i);
494                          sprintf(cmdbuf, "gensurf scale_mat %s%d %s %s %s %d %d | xform -ry 180 %s",
495 <                                        btpref, i+1, fname, fname, fname, SAMPRES-1, SAMPRES-1,
495 >                                        btpref, i, fname, fname, fname, SAMPRES-1, SAMPRES-1,
496                                          xfargs);
497                          if (!run_cmd(cmdbuf))
498                                  return(0);
# Line 502 | Line 510 | put_matBSDF(const char *XMLfile)
510                  printf("\n# Simplified material because we have no XML input\n");
511                  printf("\nvoid brightfunc latlong\n2 latlong bsdf2rad.cal\n0\n0\n");
512                  if ((front_comp|back_comp) & SDsampT)
513 <                        printf("\nlatlong trans %s\n0\n0\n7 .75 .75 .75 0 0 .5 .8\n",
513 >                        printf("\nlatlong trans %s\n0\n0\n7 .75 .75 .75 0 .04 .5 .8\n",
514                                          sph_mat);
515                  else
516                          printf("\nlatlong plastic %s\n0\n0\n5 .5 .5 .5 0 0\n",
# Line 520 | Line 528 | put_matBSDF(const char *XMLfile)
528                  break;
529          }
530          printf("\n# Actual BSDF material for rendering the hemispheres\n");
531 <        printf("\nvoid BSDF BSDFmat\n6 0 \"%s%s\" 0 1 0 .\n0\n0\n",
531 >        printf("\nvoid BSDF BSDFmat\n6 0 \"%s%s\" upx upy upz bsdf2rad.cal\n0\n0\n",
532                          curdir, XMLfile);
533          printf("\nvoid plastic black\n0\n0\n5 0 0 0 0 0\n");
534          printf("\nvoid mixfunc %s\n4 BSDFmat black latlong bsdf2rad.cal\n0\n0\n",
# Line 532 | Line 540 | static void
540   put_source(void)
541   {
542          printf("\n# Overhead parallel light source\n");
543 <        printf("\nvoid light bright\n0\n0\n3 2000 2000 2000\n");
543 >        printf("\nvoid light bright\n0\n0\n3 2500 2500 2500\n");
544          printf("\nbright source light\n0\n0\n4 0 0 1 2\n");
545          printf("\n# Material used for labels\n");
546          printf("\nvoid trans vellum\n0\n0\n7 1 1 1 0 0 .5 0\n");
# Line 542 | Line 550 | put_source(void)
550   static void
551   put_hemispheres(void)
552   {
553 +        const int       nsegs = 131;
554 +
555          printf("\n# Hemisphere(s) for showing BSDF appearance (if XML file)\n");
546 –        printf("\nvoid antimatter anti_sph\n2 void %s\n0\n0\n", sph_mat);
556          if (front_comp) {
557 <                printf("\n%s sphere Front\n0\n0\n4 %f 0 0 %f\n",
558 <                                sph_mat, -sph_xoffset, sph_rad);
559 <                printf("\n!genbox anti_sph sph_eraser %f %f %f | xform -t %f %f %f\n",
551 <                                2.02*sph_rad, 2.02*sph_rad, 1.02*sph_rad,
552 <                                -1.01*sph_rad - sph_xoffset, -1.01*sph_rad, -1.01*sph_rad);
557 >                printf(
558 > "\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",
559 >                                sph_mat, nsegs, sph_rad, 0.495*PI, sph_xoffset);
560                  printf("\nvoid brighttext front_text\n3 helvet.fnt . FRONT\n0\n");
561                  printf("12\n\t%f %f 0\n\t%f 0 0\n\t0 %f 0\n\t.01 1 -.1\n",
562 <                                -.22*sph_rad - sph_xoffset, -1.4*sph_rad,
562 >                                -.22*sph_rad + sph_xoffset, -1.4*sph_rad,
563                                  .35/5.*sph_rad, -1.6*.35/5.*sph_rad);
564                  printf("\nfront_text alias front_label_mat vellum\n");
565                  printf("\nfront_label_mat polygon front_label\n0\n0\n12");
566                  printf("\n\t%f %f 0\n\t%f %f 0\n\t%f %f 0\n\t%f %f 0\n",
567 <                                -.25*sph_rad - sph_xoffset, -1.3*sph_rad,
568 <                                -.25*sph_rad - sph_xoffset, (-1.4-1.6*.35/5.-.1)*sph_rad,
569 <                                .25*sph_rad - sph_xoffset, (-1.4-1.6*.35/5.-.1)*sph_rad,
570 <                                .25*sph_rad - sph_xoffset, -1.3*sph_rad );
567 >                                -.25*sph_rad + sph_xoffset, -1.3*sph_rad,
568 >                                -.25*sph_rad + sph_xoffset, (-1.4-1.6*.35/5.-.1)*sph_rad,
569 >                                .25*sph_rad + sph_xoffset, (-1.4-1.6*.35/5.-.1)*sph_rad,
570 >                                .25*sph_rad + sph_xoffset, -1.3*sph_rad );
571          }
572          if (back_comp) {
573 <                printf("\n%s bubble Back\n0\n0\n4 %f 0 0 %f\n",
574 <                                sph_mat, sph_xoffset, sph_rad);
575 <                printf("\n!genbox anti_sph sph_eraser %f %f %f | xform -t %f %f %f\n",
569 <                                2.02*sph_rad, 2.02*sph_rad, 1.02*sph_rad,
570 <                                -1.01*sph_rad + sph_xoffset, -1.01*sph_rad, -1.01*sph_rad);
573 >                printf(
574 > "\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",
575 >                                sph_mat, nsegs, sph_rad, 0.495*PI, -sph_xoffset);
576                  printf("\nvoid brighttext back_text\n3 helvet.fnt . BACK\n0\n");
577                  printf("12\n\t%f %f 0\n\t%f 0 0\n\t0 %f 0\n\t.01 1 -.1\n",
578 <                                -.22*sph_rad + sph_xoffset, -1.4*sph_rad,
578 >                                -.22*sph_rad - sph_xoffset, -1.4*sph_rad,
579                                  .35/4.*sph_rad, -1.6*.35/4.*sph_rad);
580                  printf("\nback_text alias back_label_mat vellum\n");
581                  printf("\nback_label_mat polygon back_label\n0\n0\n12");
582                  printf("\n\t%f %f 0\n\t%f %f 0\n\t%f %f 0\n\t%f %f 0\n",
583 <                                -.25*sph_rad + sph_xoffset, -1.3*sph_rad,
584 <                                -.25*sph_rad + sph_xoffset, (-1.4-1.6*.35/4.-.1)*sph_rad,
585 <                                .25*sph_rad + sph_xoffset, (-1.4-1.6*.35/4.-.1)*sph_rad,
586 <                                .25*sph_rad + sph_xoffset, -1.3*sph_rad );
583 >                                -.25*sph_rad - sph_xoffset, -1.3*sph_rad,
584 >                                -.25*sph_rad - sph_xoffset, (-1.4-1.6*.35/4.-.1)*sph_rad,
585 >                                .25*sph_rad - sph_xoffset, (-1.4-1.6*.35/4.-.1)*sph_rad,
586 >                                .25*sph_rad - sph_xoffset, -1.3*sph_rad );
587          }
588   }
589  
# Line 587 | Line 592 | static void
592   put_scale(void)
593   {
594          const double    max_log10 = log10(overall_max);
595 <        const double    leg_width = 2.*.75*(sph_xoffset - sph_rad);
595 >        const double    leg_width = 2.*.75*(fabs(sph_xoffset) - sph_rad);
596          const double    leg_height = 2.*sph_rad;
597          const int       text_lines = 6;
598          const int       text_digits = 8;
# Line 686 | Line 691 | convert_mgf(const char *mgfdata)
691                  sprintf(cmdbuf, "xform -t %f %f %f -s %f -t %f %f 0 %s",
692                                  -.5*(xmin+xmax), -.5*(ymin+ymax), -zmax,
693                                  1.5*sph_rad/max_dim,
694 <                                -sph_xoffset, -2.5*sph_rad,
694 >                                sph_xoffset, -2.5*sph_rad,
695                                  radfn);
696                  if (!run_cmd(cmdbuf))
697                          return;
# Line 696 | Line 701 | convert_mgf(const char *mgfdata)
701                  sprintf(cmdbuf, "xform -t %f %f %f -s %f -ry 180 -t %f %f 0 %s",
702                                  -.5*(xmin+xmax), -.5*(ymin+ymax), -zmin,
703                                  1.5*sph_rad/max_dim,
704 <                                sph_xoffset, -2.5*sph_rad,
704 >                                -sph_xoffset, -2.5*sph_rad,
705                                  radfn);
706                  if (!run_cmd(cmdbuf))
707                          return;

Diff Legend

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