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

Comparing src/cv/bsdf2rad.c (file contents):
Revision 2.18 by greg, Sun Apr 9 22:51:19 2017 UTC vs.
Revision 2.37 by greg, Fri Feb 12 00:53:56 2021 UTC

# Line 5 | Line 5 | static const char RCSid[] = "$Id$";
5   *  Plot 3-D BSDF output based on scattering interpolant or XML representation
6   */
7  
8 – #include <stdio.h>
9 – #include <string.h>
8   #include <stdlib.h>
9 + #include "rtio.h"
10   #include "paths.h"
11   #include "rtmath.h"
13 – #include "resolu.h"
12   #include "bsdfrep.h"
13  
14 + #ifndef NINCIDENT
15   #define NINCIDENT       37              /* number of samples/hemisphere */
16 <
16 > #endif
17 > #ifndef GRIDSTEP
18   #define GRIDSTEP        2               /* our grid step size */
19 + #endif
20   #define SAMPRES         (GRIDRES/GRIDSTEP)
21  
22   int     front_comp = 0;                 /* front component flags (SDsamp*) */
# Line 24 | Line 25 | double overall_min = 1./PI;            /* overall minimum BSDF v
25   double  min_log10;                      /* smallest log10 value for plotting */
26   double  overall_max = .0;               /* overall maximum BSDF value */
27  
28 < char    ourTempDir[TEMPLEN] = "";       /* our temporary directory */
28 > char    ourTempDir[TEMPLEN+1] = "";     /* our temporary directory */
29  
30 < const FVECT     Xaxis = {1., 0., 0.};
31 < const FVECT     Yaxis = {0., 1., 0.};
32 < const FVECT     Zaxis = {0., 0., 1.};
33 <
33 < const char      frpref[] = "frefl";
34 < const char      ftpref[] = "ftrans";
35 < const char      brpref[] = "brefl";
36 < const char      btpref[] = "btrans";
30 > const char      frpref[] = "rf";
31 > const char      ftpref[] = "tf";
32 > const char      brpref[] = "rb";
33 > const char      btpref[] = "tb";
34   const char      dsuffix[] = ".txt";
35  
36 < const char      sph_mat[] = "BSDFmat";
36 > const char      sph_fmat[] = "fBSDFmat";
37 > const char      sph_bmat[] = "bBSDFmat";
38   const double    sph_rad = 10.;
39   const double    sph_xoffset = 15.;
40  
41   #define bsdf_rad        (sph_rad*.25)
42   #define arrow_rad       (bsdf_rad*.015)
43  
44 < #define FEQ(a,b)        ((a)-(b) <= 1e-7 && (b)-(a) <= 1e-7)
44 > #define set_minlog()    overall_min = (overall_min < 1e-5) ? 1e-5 : overall_min; \
45 >                                min_log10 = log10(overall_min) - .1
46  
48 – #define set_minlog()    (min_log10 = log10(overall_min + 1e-5) - .1)
49 –
47   char    *progname;
48  
49   /* Get Fibonacci sphere vector (0 to NINCIDENT-1) */
50 < static void
50 > static RREAL *
51   get_ivector(FVECT iv, int i)
52   {
53          const double    phistep = PI*(3. - 2.236067978);
# Line 60 | Line 57 | get_ivector(FVECT iv, int i)
57          r = sqrt(1. - iv[2]*iv[2]);
58          iv[0] = r * cos((i+1.)*phistep);
59          iv[1] = r * sin((i+1.)*phistep);
60 +
61 +        return(iv);
62   }
63  
64 + /* Convert incident vector into sphere position */
65 + static RREAL *
66 + cvt_sposition(FVECT sp, const FVECT iv, int inc_side)
67 + {
68 +        sp[0] = -iv[0]*sph_rad + inc_side*sph_xoffset;
69 +        sp[1] = -iv[1]*sph_rad;
70 +        sp[2] = iv[2]*sph_rad;
71 +
72 +        return(sp);
73 + }
74 +
75   /* Get temporary file name */
76   static char *
77   tfile_name(const char *prefix, const char *suffix, int i)
78   {
79 <        static char     buf[128];
79 >        static char     buf[256];
80  
81          if (!ourTempDir[0]) {           /* create temporary directory */
82                  mktemp(strcpy(ourTempDir,TEMPLATE));
# Line 248 | Line 258 | build_wRBF(void)
258                  RBFNODE *rbf;
259                  get_ivector(ivec, i);
260                  if (input_orient < 0) {
261 <                        ivec[0] = -ivec[0]; ivec[1] = -ivec[1]; ivec[2] = -ivec[2];
261 >                        ivec[0] = -ivec[0]; ivec[2] = -ivec[2];
262                  }
263                  rbf = advect_rbf(ivec, 15000);
264                  if (!plotRBF(tfile_name(pref, dsuffix, i), rbf))
# Line 260 | 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 >        static int      cnt = 1;
278 >        FVECT           refl;
279          int             i;
280  
281 <        for (i = 3; i--; ) origin[i] = ivec[i]*sph_rad;
282 <        origin[0] -= inc_side*sph_xoffset;
281 >        refl[0] = 2.*nrm[2]*nrm[0];
282 >        refl[1] = 2.*nrm[2]*nrm[1];
283 >        refl[2] = 2.*nrm[2]*nrm[2] - 1.;
284  
285 <        refl[0] = 2.*ivec[2]*ivec[0];
286 <        refl[1] = 2.*ivec[2]*ivec[1];
275 <        refl[2] = 2.*ivec[2]*ivec[2] - 1.;
276 <
277 <        printf("\n# Mirror arrow\n");
278 <        printf("\narrow_mat cylinder inc_dir\n0\n0\n7");
285 >        printf("\n# Mirror arrow #%d\n", cnt);
286 >        printf("\nshaft_mat cylinder inc_dir%d\n0\n0\n7", cnt);
287          printf("\n\t%f %f %f\n\t%f %f %f\n\t%f\n",
288                          origin[0], origin[1], origin[2]+arrow_len,
289                          origin[0], origin[1], origin[2],
290                          arrow_rad);
291 <        printf("\narrow_mat cylinder mir_dir\n0\n0\n7");
291 >        printf("\nshaft_mat cylinder mir_dir%d\n0\n0\n7", cnt);
292          printf("\n\t%f %f %f\n\t%f %f %f\n\t%f\n",
293                          origin[0], origin[1], origin[2],
294                          origin[0] + arrow_len*refl[0],
295                          origin[1] + arrow_len*refl[1],
296                          origin[2] + arrow_len*refl[2],
297                          arrow_rad);
298 <        printf("\narrow_mat cone mir_tip\n0\n0\n8");
298 >        printf("\ntip_mat cone mir_tip%d\n0\n0\n8", cnt);
299          printf("\n\t%f %f %f\n\t%f %f %f\n\t%f 0\n",
300                          origin[0] + (arrow_len-.5*tip_len)*refl[0],
301                          origin[1] + (arrow_len-.5*tip_len)*refl[1],
# Line 296 | Line 304 | put_mirror_arrow(const FVECT ivec, int inc_side)
304                          origin[1] + (arrow_len+.5*tip_len)*refl[1],
305                          origin[2] + (arrow_len+.5*tip_len)*refl[2],
306                          2.*arrow_rad);
307 +        ++cnt;
308   }
309  
310   /* Put out transmitted direction arrow for the given incident vector */
311   static void
312 < put_trans_arrow(const FVECT ivec, int inc_side)
312 > put_trans_arrow(const FVECT origin)
313   {
314          const double    arrow_len = 1.2*bsdf_rad;
315          const double    tip_len = 0.2*bsdf_rad;
316 <        FVECT           origin;
316 >        static int      cnt = 1;
317          int             i;
318  
319 <        for (i = 3; i--; ) origin[i] = ivec[i]*sph_rad;
320 <        origin[0] -= inc_side*sph_xoffset;
312 <
313 <        printf("\n# Transmission arrow\n");
314 <        printf("\narrow_mat cylinder trans_dir\n0\n0\n7");
319 >        printf("\n# Transmission arrow #%d\n", cnt);
320 >        printf("\nshaft_mat cylinder trans_dir%d\n0\n0\n7", cnt);
321          printf("\n\t%f %f %f\n\t%f %f %f\n\t%f\n",
322                          origin[0], origin[1], origin[2],
323                          origin[0], origin[1], origin[2]-arrow_len,
324                          arrow_rad);
325 <        printf("\narrow_mat cone trans_tip\n0\n0\n8");
325 >        printf("\ntip_mat cone trans_tip%d\n0\n0\n8", cnt);
326          printf("\n\t%f %f %f\n\t%f %f %f\n\t%f 0\n",
327                          origin[0], origin[1], origin[2]-arrow_len+.5*tip_len,
328                          origin[0], origin[1], origin[2]-arrow_len-.5*tip_len,
329 <                        2.*arrow_rad);  
329 >                        2.*arrow_rad);
330 >        ++cnt;
331   }
332  
333   /* Compute rotation (x,y,z) => (xp,yp,zp) */
# Line 339 | Line 346 | addrot(char *xf, const FVECT xp, const FVECT yp, const
346                  return(4);
347          }
348          theta = atan2(yp[2], zp[2]);
349 <        if (!FEQ(theta,0.0)) {
349 >        if (!FABSEQ(theta,0.0)) {
350                  sprintf(xf, " -rx %f", theta*(180./PI));
351                  while (*xf) ++xf;
352                  n += 2;
353          }
354          theta = Asin(-xp[2]);
355 <        if (!FEQ(theta,0.0)) {
355 >        if (!FABSEQ(theta,0.0)) {
356                  sprintf(xf, " -ry %f", theta*(180./PI));
357                  while (*xf) ++xf;
358                  n += 2;
359          }
360          theta = atan2(xp[1], xp[0]);
361 <        if (!FEQ(theta,0.0)) {
361 >        if (!FABSEQ(theta,0.0)) {
362                  sprintf(xf, " -rz %f", theta*(180./PI));
363                  /* while (*xf) ++xf; */
364                  n += 2;
# Line 364 | Line 371 | static int
371   put_BSDFs(void)
372   {
373          const double    scalef = bsdf_rad/(log10(overall_max) - min_log10);
374 <        FVECT           ivec;
374 >        FVECT           ivec, sorg, nrm, upv;
375          RREAL           vMtx[3][3];
376          char            *fname;
377          char            cmdbuf[256];
378 <        char            xfargs[128];
379 <        int             nxfa;
378 >        char            rotargs[64];
379 >        int             nrota;
380          int             i;
381  
382          printf("\n# Gensurf output corresponding to %d incident directions\n",
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");
385 >        printf("\nvoid glow tip_mat\n0\n0\n4 1 0 1 0\n");
386 >        printf("\nvoid mixfunc shaft_mat\n4 tip_mat void 0.25 .\n0\n0\n");
387  
388 <        if (front_comp & SDsampR)
389 <                for (i = 0; i < NINCIDENT; i++) {
390 <                        get_ivector(ivec, i);
391 <                        put_mirror_arrow(ivec, 1);
392 <                        sprintf(xfargs, "-s %f -t %f %f %f", bsdf_rad,
393 <                                        ivec[0]*sph_rad - sph_xoffset,
394 <                                        ivec[1]*sph_rad, ivec[2]*sph_rad);
395 <                        nxfa = 6;
388 >        for (i = 0; i < NINCIDENT; i++) {
389 >                get_ivector(ivec, i);
390 >                nrm[0] = -ivec[0]; nrm[1] = -ivec[1]; nrm[2] = ivec[2];
391 >                upv[0] = nrm[0]*nrm[1]*(nrm[2] - 1.);
392 >                upv[1] = nrm[0]*nrm[0] + nrm[1]*nrm[1]*nrm[2];
393 >                upv[2] = -nrm[1]*(nrm[0]*nrm[0] + nrm[1]*nrm[1]);
394 >                if (SDcompXform(vMtx, nrm, upv) != SDEnone)
395 >                        continue;
396 >                nrota = addrot(rotargs, vMtx[0], vMtx[1], vMtx[2]);
397 >                if (front_comp) {
398 >                        cvt_sposition(sorg, ivec, 1);
399                          printf("\nvoid colorfunc scale_pat\n");
400 <                        printf("%d bsdf_red bsdf_grn bsdf_blu bsdf2rad.cal\n\t%s\n0\n0\n",
401 <                                        4+nxfa, xfargs);
400 >                        printf("10 bsdf_red bsdf_grn bsdf_blu bsdf2rad.cal\n");
401 >                        printf("\t-s %f -t %f %f %f\n0\n0\n",
402 >                                        bsdf_rad, sorg[0], sorg[1], sorg[2]);
403                          printf("\nscale_pat glow scale_mat\n0\n0\n4 1 1 1 0\n");
404 <                        SDcompXform(vMtx, ivec, Yaxis);
405 <                        nxfa = addrot(xfargs, vMtx[0], vMtx[1], vMtx[2]);
406 <                        sprintf(xfargs+strlen(xfargs), " -s %f -t %f %f %f",
396 <                                        scalef, ivec[0]*sph_rad - sph_xoffset,
397 <                                        ivec[1]*sph_rad, ivec[2]*sph_rad);
398 <                        nxfa += 6;
404 >                }
405 >                if (front_comp & SDsampR) {
406 >                        put_mirror_arrow(sorg, nrm);
407                          fname = tfile_name(frpref, dsuffix, i);
408 <                        sprintf(cmdbuf, "gensurf scale_mat %s%d %s %s %s %d %d | xform -mx -my %s",
409 <                                        frpref, i+1, fname, fname, fname, SAMPRES-1, SAMPRES-1,
410 <                                        xfargs);
408 >                        sprintf(cmdbuf,
409 >                "gensurf scale_mat %s%d %s %s %s %d %d | xform %s -s %f -t %f %f %f",
410 >                                        frpref, i, fname, fname, fname, SAMPRES-1, SAMPRES-1,
411 >                                        rotargs, scalef, sorg[0], sorg[1], sorg[2]);
412                          if (!run_cmd(cmdbuf))
413                                  return(0);
414                  }
415 <        if (front_comp & SDsampT)
416 <                for (i = 0; i < NINCIDENT; i++) {
408 <                        get_ivector(ivec, i);
409 <                        put_trans_arrow(ivec, 1);
410 <                        sprintf(xfargs, "-s %f -t %f %f %f", bsdf_rad,
411 <                                        ivec[0]*sph_rad - sph_xoffset,
412 <                                        ivec[1]*sph_rad, ivec[2]*sph_rad);
413 <                        nxfa = 6;
414 <                        printf("\nvoid colorfunc scale_pat\n");
415 <                        printf("%d bsdf_red bsdf_grn bsdf_blu bsdf2rad.cal\n\t%s\n0\n0\n",
416 <                                        4+nxfa, xfargs);
417 <                        printf("\nscale_pat glow scale_mat\n0\n0\n4 1 1 1 0\n");
418 <                        SDcompXform(vMtx, ivec, Yaxis);
419 <                        nxfa = addrot(xfargs, vMtx[0], vMtx[1], vMtx[2]);
420 <                        sprintf(xfargs+strlen(xfargs), " -s %f -t %f %f %f",
421 <                                        scalef, ivec[0]*sph_rad - sph_xoffset,
422 <                                        ivec[1]*sph_rad, ivec[2]*sph_rad);
423 <                        nxfa += 6;
415 >                if (front_comp & SDsampT) {
416 >                        put_trans_arrow(sorg);
417                          fname = tfile_name(ftpref, dsuffix, i);
418 <                        sprintf(cmdbuf, "gensurf scale_mat %s%d %s %s %s %d %d | xform -I -mx -my %s",
419 <                                        ftpref, i+1, fname, fname, fname, SAMPRES-1, SAMPRES-1,
420 <                                        xfargs);
418 >                        sprintf(cmdbuf,
419 >                "gensurf scale_mat %s%d %s %s %s %d %d | xform -I %s -s %f -t %f %f %f",
420 >                                        ftpref, i, fname, fname, fname, SAMPRES-1, SAMPRES-1,
421 >                                        rotargs, scalef, sorg[0], sorg[1], sorg[2]);
422                          if (!run_cmd(cmdbuf))
423                                  return(0);
424                  }
425 <        if (back_comp & SDsampR)
426 <                for (i = 0; i < NINCIDENT; i++) {
433 <                        get_ivector(ivec, i);
434 <                        put_mirror_arrow(ivec, -1);
435 <                        fname = tfile_name(brpref, dsuffix, i);
436 <                        sprintf(xfargs, "-s %f -t %f %f %f", bsdf_rad,
437 <                                        ivec[0]*sph_rad + sph_xoffset,
438 <                                        ivec[1]*sph_rad, ivec[2]*sph_rad);
439 <                        nxfa = 6;
425 >                if (back_comp) {
426 >                        cvt_sposition(sorg, ivec, -1);
427                          printf("\nvoid colorfunc scale_pat\n");
428 <                        printf("%d bsdf_red bsdf_grn bsdf_blu bsdf2rad.cal\n\t%s\n0\n0\n",
429 <                                        4+nxfa, xfargs);
428 >                        printf("10 bsdf_red bsdf_grn bsdf_blu bsdf2rad.cal\n");
429 >                        printf("\t-s %f -t %f %f %f\n0\n0\n",
430 >                                        bsdf_rad, sorg[0], sorg[1], sorg[2]);
431                          printf("\nscale_pat glow scale_mat\n0\n0\n4 1 1 1 0\n");
432 <                        SDcompXform(vMtx, ivec, Yaxis);
433 <                        nxfa = addrot(xfargs, vMtx[0], vMtx[1], vMtx[2]);
434 <                        sprintf(xfargs+strlen(xfargs), " -s %f -t %f %f %f",
447 <                                        scalef, ivec[0]*sph_rad + sph_xoffset,
448 <                                        ivec[1]*sph_rad, ivec[2]*sph_rad);
449 <                        nxfa += 6;
432 >                }
433 >                if (back_comp & SDsampR) {
434 >                        put_mirror_arrow(sorg, nrm);
435                          fname = tfile_name(brpref, dsuffix, i);
436 <                        sprintf(cmdbuf, "gensurf scale_mat %s%d %s %s %s %d %d | xform -I -ry 180 -mx -my %s",
437 <                                        brpref, i+1, fname, fname, fname, SAMPRES-1, SAMPRES-1,
438 <                                        xfargs);
436 >                        sprintf(cmdbuf,
437 >                "gensurf scale_mat %s%d %s %s %s %d %d | xform -I -ry 180 %s -s %f -t %f %f %f",
438 >                                        brpref, i, fname, fname, fname, SAMPRES-1, SAMPRES-1,
439 >                                        rotargs, scalef, sorg[0], sorg[1], sorg[2]);
440                          if (!run_cmd(cmdbuf))
441                                  return(0);
442                  }
443 <        if (back_comp & SDsampT)
444 <                for (i = 0; i < NINCIDENT; i++) {
459 <                        get_ivector(ivec, i);
460 <                        put_trans_arrow(ivec, -1);
443 >                if (back_comp & SDsampT) {
444 >                        put_trans_arrow(sorg);
445                          fname = tfile_name(btpref, dsuffix, i);
446 <                        sprintf(xfargs, "-s %f -t %f %f %f", bsdf_rad,
447 <                                        ivec[0]*sph_rad + sph_xoffset,
448 <                                        ivec[1]*sph_rad, ivec[2]*sph_rad);
449 <                        nxfa = 6;
466 <                        printf("\nvoid colorfunc scale_pat\n");
467 <                        printf("%d bsdf_red bsdf_grn bsdf_blu bsdf2rad.cal\n\t%s\n0\n0\n",
468 <                                        4+nxfa, xfargs);
469 <                        printf("\nscale_pat glow scale_mat\n0\n0\n4 1 1 1 0\n");
470 <                        SDcompXform(vMtx, ivec, Yaxis);
471 <                        nxfa = addrot(xfargs, vMtx[0], vMtx[1], vMtx[2]);
472 <                        sprintf(xfargs+strlen(xfargs), " -s %f -t %f %f %f",
473 <                                        scalef, ivec[0]*sph_rad + sph_xoffset,
474 <                                        ivec[1]*sph_rad, ivec[2]*sph_rad);
475 <                        nxfa += 6;
476 <                        fname = tfile_name(btpref, dsuffix, i);
477 <                        sprintf(cmdbuf, "gensurf scale_mat %s%d %s %s %s %d %d | xform -ry 180 -mx -my %s",
478 <                                        btpref, i+1, fname, fname, fname, SAMPRES-1, SAMPRES-1,
479 <                                        xfargs);
446 >                        sprintf(cmdbuf,
447 >                "gensurf scale_mat %s%d %s %s %s %d %d | xform -ry 180 %s -s %f -t %f %f %f",
448 >                                        btpref, i, fname, fname, fname, SAMPRES-1, SAMPRES-1,
449 >                                        rotargs, scalef, sorg[0], sorg[1], sorg[2]);
450                          if (!run_cmd(cmdbuf))
451                                  return(0);
452                  }
453 +        }
454          return(1);
455   }
456  
# Line 493 | Line 464 | put_matBSDF(const char *XMLfile)
464                  printf("\n# Simplified material because we have no XML input\n");
465                  printf("\nvoid brightfunc latlong\n2 latlong bsdf2rad.cal\n0\n0\n");
466                  if ((front_comp|back_comp) & SDsampT)
467 <                        printf("\nlatlong trans %s\n0\n0\n7 .75 .75 .75 0 0 .5 .8\n",
468 <                                        sph_mat);
467 >                        printf("\nlatlong trans %s\n0\n0\n7 .75 .75 .75 0 .04 .5 .8\n",
468 >                                        sph_fmat);
469                  else
470                          printf("\nlatlong plastic %s\n0\n0\n5 .5 .5 .5 0 0\n",
471 <                                        sph_mat);
471 >                                        sph_fmat);
472 >                printf("\ninherit alias %s %s\n", sph_bmat, sph_fmat);
473                  return;
474          }
475          switch (XMLfile[0]) {           /* avoid RAYPATH search */
476          case '.':
477 +        case '~':
478          CASEDIRSEP:
479                  curdir = "";
480                  break;
# Line 510 | Line 483 | put_matBSDF(const char *XMLfile)
483                  exit(1);
484                  break;
485          }
486 <        printf("\n# Actual BSDF material for rendering the hemispheres\n");
487 <        printf("\nvoid BSDF BSDFmat\n6 0 \"%s%s\" 0 1 0 .\n0\n0\n",
486 >        printf("\n# Actual BSDF materials for rendering the hemispheres\n");
487 >        printf("\nvoid BSDF BSDF_f\n6 0 \"%s%s\" upx upy upz bsdf2rad.cal\n0\n0\n",
488                          curdir, XMLfile);
489          printf("\nvoid plastic black\n0\n0\n5 0 0 0 0 0\n");
490 <        printf("\nvoid mixfunc %s\n4 BSDFmat black latlong bsdf2rad.cal\n0\n0\n",
491 <                        sph_mat);
490 >        printf("\nvoid mixfunc %s\n4 BSDF_f black latlong bsdf2rad.cal\n0\n0\n",
491 >                        sph_fmat);
492 >        printf("\nvoid BSDF BSDF_b\n8 0 \"%s%s\" upx upy upz bsdf2rad.cal -ry 180\n0\n0\n",
493 >                        curdir, XMLfile);
494 >        printf("\nvoid mixfunc %s\n4 BSDF_b black latlong bsdf2rad.cal\n0\n0\n",
495 >                        sph_bmat);
496   }
497  
498   /* Put out overhead parallel light source */
# Line 523 | Line 500 | static void
500   put_source(void)
501   {
502          printf("\n# Overhead parallel light source\n");
503 <        printf("\nvoid light bright\n0\n0\n3 2000 2000 2000\n");
503 >        printf("\nvoid light bright\n0\n0\n3 2500 2500 2500\n");
504          printf("\nbright source light\n0\n0\n4 0 0 1 2\n");
505          printf("\n# Material used for labels\n");
506          printf("\nvoid trans vellum\n0\n0\n7 1 1 1 0 0 .5 0\n");
# Line 533 | Line 510 | put_source(void)
510   static void
511   put_hemispheres(void)
512   {
513 +        const int       nsegs = 131;
514 +
515          printf("\n# Hemisphere(s) for showing BSDF appearance (if XML file)\n");
537 –        printf("\nvoid antimatter anti_sph\n2 void %s\n0\n0\n", sph_mat);
516          if (front_comp) {
517 <                printf("\n%s sphere Front\n0\n0\n4 %f 0 0 %f\n",
518 <                                sph_mat, -sph_xoffset, sph_rad);
519 <                printf("\n!genbox anti_sph sph_eraser %f %f %f | xform -t %f %f %f\n",
542 <                                2.02*sph_rad, 2.02*sph_rad, 1.02*sph_rad,
543 <                                -1.01*sph_rad - sph_xoffset, -1.01*sph_rad, -1.01*sph_rad);
517 >                printf(
518 > "\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",
519 >                                sph_fmat, nsegs, sph_rad, 0.5*PI, sph_xoffset);
520                  printf("\nvoid brighttext front_text\n3 helvet.fnt . FRONT\n0\n");
521                  printf("12\n\t%f %f 0\n\t%f 0 0\n\t0 %f 0\n\t.01 1 -.1\n",
522 <                                -.22*sph_rad - sph_xoffset, -1.4*sph_rad,
522 >                                -.22*sph_rad + sph_xoffset, -1.4*sph_rad,
523                                  .35/5.*sph_rad, -1.6*.35/5.*sph_rad);
524                  printf("\nfront_text alias front_label_mat vellum\n");
525                  printf("\nfront_label_mat polygon front_label\n0\n0\n12");
526                  printf("\n\t%f %f 0\n\t%f %f 0\n\t%f %f 0\n\t%f %f 0\n",
527 <                                -.25*sph_rad - sph_xoffset, -1.3*sph_rad,
528 <                                -.25*sph_rad - sph_xoffset, (-1.4-1.6*.35/5.-.1)*sph_rad,
529 <                                .25*sph_rad - sph_xoffset, (-1.4-1.6*.35/5.-.1)*sph_rad,
530 <                                .25*sph_rad - sph_xoffset, -1.3*sph_rad );
527 >                                -.25*sph_rad + sph_xoffset, -1.3*sph_rad,
528 >                                -.25*sph_rad + sph_xoffset, (-1.4-1.6*.35/5.-.1)*sph_rad,
529 >                                .25*sph_rad + sph_xoffset, (-1.4-1.6*.35/5.-.1)*sph_rad,
530 >                                .25*sph_rad + sph_xoffset, -1.3*sph_rad );
531          }
532          if (back_comp) {
533 <                printf("\n%s bubble Back\n0\n0\n4 %f 0 0 %f\n",
534 <                                sph_mat, sph_xoffset, sph_rad);
535 <                printf("\n!genbox anti_sph sph_eraser %f %f %f | xform -t %f %f %f\n",
560 <                                2.02*sph_rad, 2.02*sph_rad, 1.02*sph_rad,
561 <                                -1.01*sph_rad + sph_xoffset, -1.01*sph_rad, -1.01*sph_rad);
533 >                printf(
534 > "\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",
535 >                                sph_bmat, nsegs, sph_rad, 0.5*PI, -sph_xoffset);
536                  printf("\nvoid brighttext back_text\n3 helvet.fnt . BACK\n0\n");
537                  printf("12\n\t%f %f 0\n\t%f 0 0\n\t0 %f 0\n\t.01 1 -.1\n",
538 <                                -.22*sph_rad + sph_xoffset, -1.4*sph_rad,
538 >                                -.22*sph_rad - sph_xoffset, -1.4*sph_rad,
539                                  .35/4.*sph_rad, -1.6*.35/4.*sph_rad);
540                  printf("\nback_text alias back_label_mat vellum\n");
541                  printf("\nback_label_mat polygon back_label\n0\n0\n12");
542                  printf("\n\t%f %f 0\n\t%f %f 0\n\t%f %f 0\n\t%f %f 0\n",
543 <                                -.25*sph_rad + sph_xoffset, -1.3*sph_rad,
544 <                                -.25*sph_rad + sph_xoffset, (-1.4-1.6*.35/4.-.1)*sph_rad,
545 <                                .25*sph_rad + sph_xoffset, (-1.4-1.6*.35/4.-.1)*sph_rad,
546 <                                .25*sph_rad + sph_xoffset, -1.3*sph_rad );
543 >                                -.25*sph_rad - sph_xoffset, -1.3*sph_rad,
544 >                                -.25*sph_rad - sph_xoffset, (-1.4-1.6*.35/4.-.1)*sph_rad,
545 >                                .25*sph_rad - sph_xoffset, (-1.4-1.6*.35/4.-.1)*sph_rad,
546 >                                .25*sph_rad - sph_xoffset, -1.3*sph_rad );
547          }
548   }
549  
# Line 578 | Line 552 | static void
552   put_scale(void)
553   {
554          const double    max_log10 = log10(overall_max);
555 <        const double    leg_width = 2.*.75*(sph_xoffset - sph_rad);
555 >        const double    leg_width = 2.*.75*(fabs(sph_xoffset) - sph_rad);
556          const double    leg_height = 2.*sph_rad;
557          const int       text_lines = 6;
558          const int       text_digits = 8;
# Line 677 | Line 651 | convert_mgf(const char *mgfdata)
651                  sprintf(cmdbuf, "xform -t %f %f %f -s %f -t %f %f 0 %s",
652                                  -.5*(xmin+xmax), -.5*(ymin+ymax), -zmax,
653                                  1.5*sph_rad/max_dim,
654 <                                -sph_xoffset, -2.5*sph_rad,
654 >                                sph_xoffset, -2.5*sph_rad,
655                                  radfn);
656                  if (!run_cmd(cmdbuf))
657                          return;
# Line 687 | Line 661 | convert_mgf(const char *mgfdata)
661                  sprintf(cmdbuf, "xform -t %f %f %f -s %f -ry 180 -t %f %f 0 %s",
662                                  -.5*(xmin+xmax), -.5*(ymin+ymax), -zmin,
663                                  1.5*sph_rad/max_dim,
664 <                                sph_xoffset, -2.5*sph_rad,
664 >                                -sph_xoffset, -2.5*sph_rad,
665                                  radfn);
666                  if (!run_cmd(cmdbuf))
667                          return;
# Line 698 | Line 672 | convert_mgf(const char *mgfdata)
672   static int
673   rbf_headline(char *s, void *p)
674   {
675 <        char    fmt[64];
675 >        char    fmt[MAXFMTLEN];
676  
677          if (formatval(fmt, s)) {
678                  if (strcmp(fmt, BSDFREP_FMT))
# Line 732 | Line 706 | int
706   main(int argc, char *argv[])
707   {
708          int     inpXML = -1;
709 +        double  myLim[2];
710          SDData  myBSDF;
711 <        int     n;
711 >        int     a, n;
712                                                  /* check arguments */
713          progname = argv[0];
714 <        if (argc > 1 && (n = strlen(argv[1])-4) > 0) {
715 <                if (!strcasecmp(argv[1]+n, ".xml"))
714 >        a = 1;
715 >        myLim[0] = -1; myLim[1] = -2;           /* specified BSDF range? */
716 >        if (argc > a+3 && argv[a][0] == '-' && argv[a][1] == 'r') {
717 >                myLim[0] = atof(argv[++a]);
718 >                myLim[1] = atof(argv[++a]);
719 >                ++a;
720 >        }
721 >        if (argc > a && (n = strlen(argv[a])-4) > 0) {
722 >                if (!strcasecmp(argv[a]+n, ".xml"))
723                          inpXML = 1;
724 <                else if (!strcasecmp(argv[1]+n, ".sir"))
724 >                else if (!strcasecmp(argv[a]+n, ".sir"))
725                          inpXML = 0;
726          }
727 <        if (inpXML < 0 || inpXML & (argc > 2)) {
728 <                fprintf(stderr, "Usage: %s bsdf.xml > output.rad\n", progname);
729 <                fprintf(stderr, "   Or: %s hemi1.sir hemi2.sir .. > output.rad\n", progname);
727 >        if (inpXML < 0 || inpXML & (argc > a+1)) {
728 >                fprintf(stderr, "Usage: %s [-r min max] bsdf.xml > output.rad\n", progname);
729 >                fprintf(stderr, "   Or: %s [-r min max] hemi1.sir hemi2.sir .. > output.rad\n", progname);
730                  return(1);
731          }
732          fputs("# ", stdout);                    /* copy our command */
733          printargs(argc, argv, stdout);
734                                                  /* evaluate BSDF */
735          if (inpXML) {
736 <                SDclearBSDF(&myBSDF, argv[1]);
737 <                if (SDreportError(SDloadFile(&myBSDF, argv[1]), stderr))
736 >                SDclearBSDF(&myBSDF, argv[a]);
737 >                if (SDreportError(SDloadFile(&myBSDF, argv[a]), stderr))
738                          return(1);
739                  if (myBSDF.rf != NULL) front_comp |= SDsampR;
740                  if (myBSDF.tf != NULL) front_comp |= SDsampT;
# Line 760 | Line 742 | main(int argc, char *argv[])
742                  if (myBSDF.tb != NULL) back_comp |= SDsampT;
743                  if (!front_comp & !back_comp) {
744                          fprintf(stderr, "%s: nothing to plot in '%s'\n",
745 <                                        progname, argv[1]);
745 >                                        progname, argv[a]);
746                          return(1);
747                  }
748 <                if (front_comp & SDsampR && myBSDF.rLambFront.cieY < overall_min*PI)
749 <                        overall_min = myBSDF.rLambFront.cieY/PI;
750 <                if (back_comp & SDsampR && myBSDF.rLambBack.cieY < overall_min*PI)
751 <                        overall_min = myBSDF.rLambBack.cieY/PI;
752 <                if ((front_comp|back_comp) & SDsampT &&
753 <                                myBSDF.tLamb.cieY < overall_min*PI)
754 <                        overall_min = myBSDF.tLamb.cieY/PI;
748 >                if (myLim[0] >= 0)
749 >                        overall_min = myLim[0];
750 >                else {
751 >                        if (front_comp & SDsampR && myBSDF.rLambFront.cieY < overall_min*PI)
752 >                                overall_min = myBSDF.rLambFront.cieY/PI;
753 >                        if (back_comp & SDsampR && myBSDF.rLambBack.cieY < overall_min*PI)
754 >                                overall_min = myBSDF.rLambBack.cieY/PI;
755 >                        if ((front_comp|back_comp) & SDsampT &&
756 >                                        myBSDF.tLamb.cieY < overall_min*PI)
757 >                                overall_min = myBSDF.tLamb.cieY/PI;
758 >                }
759                  set_minlog();
760                  if (!build_wBSDF(&myBSDF))
761                          return(1);
# Line 778 | Line 764 | main(int argc, char *argv[])
764                  else
765                          strcpy(bsdf_name, myBSDF.name);
766                  strcpy(bsdf_manuf, myBSDF.makr);
767 <                put_matBSDF(argv[1]);
767 >                put_matBSDF(argv[a]);
768          } else {
769 <                FILE    *fp;
770 <                for (n = 1; n < argc; n++) {
771 <                        fp = fopen(argv[n], "rb");
772 <                        if (fp == NULL) {
769 >                FILE    *fp[4];
770 >                if (argc > a+4) {
771 >                        fprintf(stderr, "%s: more than 4 hemispheres!\n", progname);
772 >                        return(1);
773 >                }
774 >                for (n = a; n < argc; n++) {
775 >                        fp[n-a] = fopen(argv[n], "rb");
776 >                        if (fp[n-a] == NULL) {
777                                  fprintf(stderr, "%s: cannot open BSDF interpolant '%s'\n",
778                                                  progname, argv[n]);
779                                  return(1);
780                          }
781 <                        if (getheader(fp, rbf_headline, NULL) < 0) {
781 >                        if (getheader(fp[n-a], rbf_headline, NULL) < 0) {
782                                  fprintf(stderr, "%s: bad BSDF interpolant '%s'\n",
783                                                  progname, argv[n]);
784                                  return(1);
785                          }
796 –                        fclose(fp);
786                  }
787 +                if (myLim[0] >= 0)
788 +                        overall_min = myLim[0];
789                  set_minlog();
790 <                for (n = 1; n < argc; n++) {
791 <                        fp = fopen(argv[n], "rb");
792 <                        if (!load_bsdf_rep(fp))
790 >                for (n = a; n < argc; n++) {
791 >                        if (fseek(fp[n-a], 0L, SEEK_SET) < 0) {
792 >                                fprintf(stderr, "%s: cannot seek on '%s'\n",
793 >                                                progname, argv[n]);
794                                  return(1);
795 <                        fclose(fp);
795 >                        }
796 >                        if (!load_bsdf_rep(fp[n-a]))
797 >                                return(1);
798 >                        fclose(fp[n-a]);
799                          if (!build_wRBF())
800                                  return(1);
801                  }
802                  put_matBSDF(NULL);
803          }
804 +        if (myLim[1] > myLim[0])        /* override maximum BSDF? */
805 +                overall_max = myLim[1];
806          put_source();                   /* before hemispheres & labels */
807          put_hemispheres();
808          put_scale();
809          if (inpXML && myBSDF.mgf)
810                  convert_mgf(myBSDF.mgf);
811 <        if (!put_BSDFs())
811 >        if (!put_BSDFs())               /* most of the output happens here */
812                  return(1);
813          cleanup_tmp();
814          return(0);

Diff Legend

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