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.25 by greg, Wed Apr 12 00:31:40 2017 UTC vs.
Revision 2.33 by greg, Mon Oct 2 22:11:32 2017 UTC

# Line 13 | Line 13 | static const char RCSid[] = "$Id$";
13   #include "resolu.h"
14   #include "bsdfrep.h"
15  
16 + #ifndef NINCIDENT
17   #define NINCIDENT       37              /* number of samples/hemisphere */
18 <
18 > #endif
19 > #ifndef GRIDSTEP
20   #define GRIDSTEP        2               /* our grid step size */
21 + #endif
22   #define SAMPRES         (GRIDRES/GRIDSTEP)
23  
24   int     front_comp = 0;                 /* front component flags (SDsamp*) */
# Line 24 | Line 27 | double overall_min = 1./PI;            /* overall minimum BSDF v
27   double  min_log10;                      /* smallest log10 value for plotting */
28   double  overall_max = .0;               /* overall maximum BSDF value */
29  
30 < char    ourTempDir[TEMPLEN] = "";       /* our temporary directory */
30 > char    ourTempDir[TEMPLEN+1] = "";     /* our temporary directory */
31  
32 < const char      frpref[] = "frefl";
33 < const char      ftpref[] = "ftrans";
34 < const char      brpref[] = "brefl";
35 < const char      btpref[] = "btrans";
32 > const char      frpref[] = "rf";
33 > const char      ftpref[] = "tf";
34 > const char      brpref[] = "rb";
35 > const char      btpref[] = "tb";
36   const char      dsuffix[] = ".txt";
37  
38   const char      sph_fmat[] = "fBSDFmat";
# Line 42 | Line 45 | const double   sph_xoffset = 15.;
45  
46   #define FEQ(a,b)        ((a)-(b) <= 1e-7 && (b)-(a) <= 1e-7)
47  
48 < #define set_minlog()    (min_log10 = log10(overall_min + 1e-5) - .1)
48 > #define set_minlog()    overall_min = (overall_min < 1e-5) ? 1e-5 : overall_min; \
49 >                                min_log10 = log10(overall_min) - .1
50  
51   char    *progname;
52  
# Line 76 | Line 80 | cvt_sposition(FVECT sp, const FVECT iv, int inc_side)
80   static char *
81   tfile_name(const char *prefix, const char *suffix, int i)
82   {
83 <        static char     buf[128];
83 >        static char     buf[256];
84  
85          if (!ourTempDir[0]) {           /* create temporary directory */
86                  mktemp(strcpy(ourTempDir,TEMPLATE));
# Line 274 | Line 278 | put_mirror_arrow(const FVECT origin, const FVECT nrm)
278   {
279          const double    arrow_len = 1.2*bsdf_rad;
280          const double    tip_len = 0.2*bsdf_rad;
281 +        static int      cnt = 1;
282          FVECT           refl;
283          int             i;
284  
# Line 281 | Line 286 | put_mirror_arrow(const FVECT origin, const FVECT nrm)
286          refl[1] = 2.*nrm[2]*nrm[1];
287          refl[2] = 2.*nrm[2]*nrm[2] - 1.;
288  
289 <        printf("\n# Mirror arrow\n");
290 <        printf("\narrow_mat cylinder inc_dir\n0\n0\n7");
289 >        printf("\n# Mirror arrow #%d\n", cnt);
290 >        printf("\nshaft_mat cylinder inc_dir%d\n0\n0\n7", cnt);
291          printf("\n\t%f %f %f\n\t%f %f %f\n\t%f\n",
292                          origin[0], origin[1], origin[2]+arrow_len,
293                          origin[0], origin[1], origin[2],
294                          arrow_rad);
295 <        printf("\narrow_mat cylinder mir_dir\n0\n0\n7");
295 >        printf("\nshaft_mat cylinder mir_dir%d\n0\n0\n7", cnt);
296          printf("\n\t%f %f %f\n\t%f %f %f\n\t%f\n",
297                          origin[0], origin[1], origin[2],
298                          origin[0] + arrow_len*refl[0],
299                          origin[1] + arrow_len*refl[1],
300                          origin[2] + arrow_len*refl[2],
301                          arrow_rad);
302 <        printf("\narrow_mat cone mir_tip\n0\n0\n8");
302 >        printf("\ntip_mat cone mir_tip%d\n0\n0\n8", cnt);
303          printf("\n\t%f %f %f\n\t%f %f %f\n\t%f 0\n",
304                          origin[0] + (arrow_len-.5*tip_len)*refl[0],
305                          origin[1] + (arrow_len-.5*tip_len)*refl[1],
# Line 303 | Line 308 | put_mirror_arrow(const FVECT origin, const FVECT nrm)
308                          origin[1] + (arrow_len+.5*tip_len)*refl[1],
309                          origin[2] + (arrow_len+.5*tip_len)*refl[2],
310                          2.*arrow_rad);
311 +        ++cnt;
312   }
313  
314   /* Put out transmitted direction arrow for the given incident vector */
# Line 311 | Line 317 | put_trans_arrow(const FVECT origin)
317   {
318          const double    arrow_len = 1.2*bsdf_rad;
319          const double    tip_len = 0.2*bsdf_rad;
320 +        static int      cnt = 1;
321          int             i;
322  
323 <        printf("\n# Transmission arrow\n");
324 <        printf("\narrow_mat cylinder trans_dir\n0\n0\n7");
323 >        printf("\n# Transmission arrow #%d\n", cnt);
324 >        printf("\nshaft_mat cylinder trans_dir%d\n0\n0\n7", cnt);
325          printf("\n\t%f %f %f\n\t%f %f %f\n\t%f\n",
326                          origin[0], origin[1], origin[2],
327                          origin[0], origin[1], origin[2]-arrow_len,
328                          arrow_rad);
329 <        printf("\narrow_mat cone trans_tip\n0\n0\n8");
329 >        printf("\ntip_mat cone trans_tip%d\n0\n0\n8", cnt);
330          printf("\n\t%f %f %f\n\t%f %f %f\n\t%f 0\n",
331                          origin[0], origin[1], origin[2]-arrow_len+.5*tip_len,
332                          origin[0], origin[1], origin[2]-arrow_len-.5*tip_len,
333 <                        2.*arrow_rad);  
333 >                        2.*arrow_rad);
334 >        ++cnt;
335   }
336  
337   /* Compute rotation (x,y,z) => (xp,yp,zp) */
# Line 378 | Line 386 | put_BSDFs(void)
386          printf("\n# Gensurf output corresponding to %d incident directions\n",
387                          NINCIDENT);
388  
389 <        printf("\nvoid glow arrow_glow\n0\n0\n4 1 0 1 0\n");
390 <        printf("\nvoid mixfunc arrow_mat\n4 arrow_glow void 0.25 .\n0\n0\n");
389 >        printf("\nvoid glow tip_mat\n0\n0\n4 1 0 1 0\n");
390 >        printf("\nvoid mixfunc shaft_mat\n4 tip_mat void 0.25 .\n0\n0\n");
391  
384 –        if (front_comp & SDsampR)                       /* front reflection */
392          for (i = 0; i < NINCIDENT; i++) {
393                  get_ivector(ivec, i);
394                  nrm[0] = -ivec[0]; nrm[1] = -ivec[1]; nrm[2] = ivec[2];
# Line 471 | Line 478 | put_matBSDF(const char *XMLfile)
478          }
479          switch (XMLfile[0]) {           /* avoid RAYPATH search */
480          case '.':
481 +        case '~':
482          CASEDIRSEP:
483                  curdir = "";
484                  break;
# Line 512 | Line 520 | put_hemispheres(void)
520          if (front_comp) {
521                  printf(
522   "\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",
523 <                                sph_fmat, nsegs, sph_rad, 0.495*PI, sph_xoffset);
523 >                                sph_fmat, nsegs, sph_rad, 0.5*PI, sph_xoffset);
524                  printf("\nvoid brighttext front_text\n3 helvet.fnt . FRONT\n0\n");
525                  printf("12\n\t%f %f 0\n\t%f 0 0\n\t0 %f 0\n\t.01 1 -.1\n",
526                                  -.22*sph_rad + sph_xoffset, -1.4*sph_rad,
# Line 528 | Line 536 | put_hemispheres(void)
536          if (back_comp) {
537                  printf(
538   "\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",
539 <                                sph_bmat, nsegs, sph_rad, 0.495*PI, -sph_xoffset);
539 >                                sph_bmat, nsegs, sph_rad, 0.5*PI, -sph_xoffset);
540                  printf("\nvoid brighttext back_text\n3 helvet.fnt . BACK\n0\n");
541                  printf("12\n\t%f %f 0\n\t%f 0 0\n\t0 %f 0\n\t.01 1 -.1\n",
542                                  -.22*sph_rad - sph_xoffset, -1.4*sph_rad,
# Line 750 | Line 758 | main(int argc, char *argv[])
758                  strcpy(bsdf_manuf, myBSDF.makr);
759                  put_matBSDF(argv[1]);
760          } else {
761 <                FILE    *fp;
761 >                FILE    *fp[4];
762 >                if (argc > 5) {
763 >                        fprintf(stderr, "%s: more than 4 hemispheres!\n", progname);
764 >                        return(1);
765 >                }
766                  for (n = 1; n < argc; n++) {
767 <                        fp = fopen(argv[n], "rb");
768 <                        if (fp == NULL) {
767 >                        fp[n-1] = fopen(argv[n], "rb");
768 >                        if (fp[n-1] == NULL) {
769                                  fprintf(stderr, "%s: cannot open BSDF interpolant '%s'\n",
770                                                  progname, argv[n]);
771                                  return(1);
772                          }
773 <                        if (getheader(fp, rbf_headline, NULL) < 0) {
773 >                        if (getheader(fp[n-1], rbf_headline, NULL) < 0) {
774                                  fprintf(stderr, "%s: bad BSDF interpolant '%s'\n",
775                                                  progname, argv[n]);
776                                  return(1);
777                          }
766 –                        fclose(fp);
778                  }
779                  set_minlog();
780                  for (n = 1; n < argc; n++) {
781 <                        fp = fopen(argv[n], "rb");
782 <                        if (!load_bsdf_rep(fp))
781 >                        if (fseek(fp[n-1], 0L, SEEK_SET) < 0) {
782 >                                fprintf(stderr, "%s: cannot seek on '%s'\n",
783 >                                                progname, argv[n]);
784                                  return(1);
785 <                        fclose(fp);
785 >                        }
786 >                        if (!load_bsdf_rep(fp[n-1]))
787 >                                return(1);
788 >                        fclose(fp[n-1]);
789                          if (!build_wRBF())
790                                  return(1);
791                  }
# Line 781 | Line 796 | main(int argc, char *argv[])
796          put_scale();
797          if (inpXML && myBSDF.mgf)
798                  convert_mgf(myBSDF.mgf);
799 <        if (!put_BSDFs())
799 >        if (!put_BSDFs())               /* most of the output happens here */
800                  return(1);
801          cleanup_tmp();
802          return(0);

Diff Legend

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