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.16 by greg, Sun Apr 9 21:39:26 2017 UTC vs.
Revision 2.32 by greg, Thu Aug 3 19:50:12 2017 UTC

# Line 13 | Line 13 | static const char RCSid[] = "$Id$";
13   #include "resolu.h"
14   #include "bsdfrep.h"
15  
16 < #define NINCIDENT       25              /* number of samples/hemisphere */
16 > #define NINCIDENT       37              /* number of samples/hemisphere */
17  
18   #define GRIDSTEP        2               /* our grid step size */
19   #define SAMPRES         (GRIDRES/GRIDSTEP)
# 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 <
33 < const char      frpref[] = "frefl";
34 < const char      ftpref[] = "ftrans";
35 < const char      brpref[] = "brefl";
36 < 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 45 | Line 42 | const double   sph_xoffset = 15.;
42  
43   #define FEQ(a,b)        ((a)-(b) <= 1e-7 && (b)-(a) <= 1e-7)
44  
45 < #define set_minlog()    (min_log10 = log10(overall_min + 1e-5) - .1)
45 > #define set_minlog()    overall_min = (overall_min < 1e-5) ? 1e-5 : overall_min; \
46 >                                min_log10 = log10(overall_min) - .1
47  
48   char    *progname;
49  
50   /* Get Fibonacci sphere vector (0 to NINCIDENT-1) */
51 < static void
51 > static RREAL *
52   get_ivector(FVECT iv, int i)
53   {
56 –        const double    zstep = 1./(2*NINCIDENT-1);
54          const double    phistep = PI*(3. - 2.236067978);
55          double          r;
56  
57 <        iv[2] = 1. - (i+.5)*zstep;
57 >        iv[2] = 1. - (i+.5)*(1./NINCIDENT);
58          r = sqrt(1. - iv[2]*iv[2]);
59          iv[0] = r * cos((i+1.)*phistep);
60          iv[1] = r * sin((i+1.)*phistep);
61 +
62 +        return(iv);
63   }
64  
65 + /* Convert incident vector into sphere position */
66 + static RREAL *
67 + cvt_sposition(FVECT sp, const FVECT iv, int inc_side)
68 + {
69 +        sp[0] = -iv[0]*sph_rad + inc_side*sph_xoffset;
70 +        sp[1] = -iv[1]*sph_rad;
71 +        sp[2] = iv[2]*sph_rad;
72 +
73 +        return(sp);
74 + }
75 +
76   /* Get temporary file name */
77   static char *
78   tfile_name(const char *prefix, const char *suffix, int i)
# Line 249 | Line 259 | build_wRBF(void)
259                  RBFNODE *rbf;
260                  get_ivector(ivec, i);
261                  if (input_orient < 0) {
262 <                        ivec[0] = -ivec[0]; ivec[1] = -ivec[1]; ivec[2] = -ivec[2];
262 >                        ivec[0] = -ivec[0]; ivec[2] = -ivec[2];
263                  }
264                  rbf = advect_rbf(ivec, 15000);
265                  if (!plotRBF(tfile_name(pref, dsuffix, i), rbf))
# Line 261 | Line 271 | build_wRBF(void)
271  
272   /* Put out mirror arrow for the given incident vector */
273   static void
274 < put_mirror_arrow(const FVECT ivec, int inc_side)
274 > put_mirror_arrow(const FVECT origin, const FVECT nrm)
275   {
276          const double    arrow_len = 1.2*bsdf_rad;
277          const double    tip_len = 0.2*bsdf_rad;
278 <        FVECT           origin, refl;
278 >        static int      cnt = 1;
279 >        FVECT           refl;
280          int             i;
281  
282 <        for (i = 3; i--; ) origin[i] = ivec[i]*sph_rad;
283 <        origin[0] -= inc_side*sph_xoffset;
282 >        refl[0] = 2.*nrm[2]*nrm[0];
283 >        refl[1] = 2.*nrm[2]*nrm[1];
284 >        refl[2] = 2.*nrm[2]*nrm[2] - 1.;
285  
286 <        refl[0] = 2.*ivec[2]*ivec[0];
287 <        refl[1] = 2.*ivec[2]*ivec[1];
276 <        refl[2] = 2.*ivec[2]*ivec[2] - 1.;
277 <
278 <        printf("\n# Mirror arrow\n");
279 <        printf("\narrow_mat cylinder inc_dir\n0\n0\n7");
286 >        printf("\n# Mirror arrow #%d\n", cnt);
287 >        printf("\nshaft_mat cylinder inc_dir%d\n0\n0\n7", cnt);
288          printf("\n\t%f %f %f\n\t%f %f %f\n\t%f\n",
289                          origin[0], origin[1], origin[2]+arrow_len,
290                          origin[0], origin[1], origin[2],
291                          arrow_rad);
292 <        printf("\narrow_mat cylinder mir_dir\n0\n0\n7");
292 >        printf("\nshaft_mat cylinder mir_dir%d\n0\n0\n7", cnt);
293          printf("\n\t%f %f %f\n\t%f %f %f\n\t%f\n",
294                          origin[0], origin[1], origin[2],
295                          origin[0] + arrow_len*refl[0],
296                          origin[1] + arrow_len*refl[1],
297                          origin[2] + arrow_len*refl[2],
298                          arrow_rad);
299 <        printf("\narrow_mat cone mir_tip\n0\n0\n8");
299 >        printf("\ntip_mat cone mir_tip%d\n0\n0\n8", cnt);
300          printf("\n\t%f %f %f\n\t%f %f %f\n\t%f 0\n",
301                          origin[0] + (arrow_len-.5*tip_len)*refl[0],
302                          origin[1] + (arrow_len-.5*tip_len)*refl[1],
# Line 297 | Line 305 | put_mirror_arrow(const FVECT ivec, int inc_side)
305                          origin[1] + (arrow_len+.5*tip_len)*refl[1],
306                          origin[2] + (arrow_len+.5*tip_len)*refl[2],
307                          2.*arrow_rad);
308 +        ++cnt;
309   }
310  
311   /* Put out transmitted direction arrow for the given incident vector */
312   static void
313 < put_trans_arrow(const FVECT ivec, int inc_side)
313 > put_trans_arrow(const FVECT origin)
314   {
315          const double    arrow_len = 1.2*bsdf_rad;
316          const double    tip_len = 0.2*bsdf_rad;
317 <        FVECT           origin;
317 >        static int      cnt = 1;
318          int             i;
319  
320 <        for (i = 3; i--; ) origin[i] = ivec[i]*sph_rad;
321 <        origin[0] -= inc_side*sph_xoffset;
313 <
314 <        printf("\n# Transmission arrow\n");
315 <        printf("\narrow_mat cylinder trans_dir\n0\n0\n7");
320 >        printf("\n# Transmission arrow #%d\n", cnt);
321 >        printf("\nshaft_mat cylinder trans_dir%d\n0\n0\n7", cnt);
322          printf("\n\t%f %f %f\n\t%f %f %f\n\t%f\n",
323                          origin[0], origin[1], origin[2],
324                          origin[0], origin[1], origin[2]-arrow_len,
325                          arrow_rad);
326 <        printf("\narrow_mat cone trans_tip\n0\n0\n8");
326 >        printf("\ntip_mat cone trans_tip%d\n0\n0\n8", cnt);
327          printf("\n\t%f %f %f\n\t%f %f %f\n\t%f 0\n",
328                          origin[0], origin[1], origin[2]-arrow_len+.5*tip_len,
329                          origin[0], origin[1], origin[2]-arrow_len-.5*tip_len,
330 <                        2.*arrow_rad);  
330 >                        2.*arrow_rad);
331 >        ++cnt;
332   }
333  
334   /* Compute rotation (x,y,z) => (xp,yp,zp) */
# Line 365 | Line 372 | static int
372   put_BSDFs(void)
373   {
374          const double    scalef = bsdf_rad/(log10(overall_max) - min_log10);
375 <        FVECT           ivec;
375 >        FVECT           ivec, sorg, nrm, upv;
376          RREAL           vMtx[3][3];
377          char            *fname;
378          char            cmdbuf[256];
379 <        char            xfargs[128];
380 <        int             nxfa;
379 >        char            rotargs[64];
380 >        int             nrota;
381          int             i;
382  
383          printf("\n# Gensurf output corresponding to %d incident directions\n",
384                          NINCIDENT);
385  
386 <        printf("\nvoid glow arrow_glow\n0\n0\n4 1 0 1 0\n");
387 <        printf("\nvoid mixfunc arrow_mat\n4 arrow_glow void .5 .\n0\n0\n");
386 >        printf("\nvoid glow tip_mat\n0\n0\n4 1 0 1 0\n");
387 >        printf("\nvoid mixfunc shaft_mat\n4 tip_mat void 0.25 .\n0\n0\n");
388  
389 <        if (front_comp & SDsampR)
390 <                for (i = 0; i < NINCIDENT; i++) {
391 <                        get_ivector(ivec, i);
392 <                        put_mirror_arrow(ivec, 1);
393 <                        sprintf(xfargs, "-s %f -t %f %f %f", bsdf_rad,
394 <                                        ivec[0]*sph_rad - sph_xoffset,
395 <                                        ivec[1]*sph_rad, ivec[2]*sph_rad);
396 <                        nxfa = 6;
389 >        for (i = 0; i < NINCIDENT; i++) {
390 >                get_ivector(ivec, i);
391 >                nrm[0] = -ivec[0]; nrm[1] = -ivec[1]; nrm[2] = ivec[2];
392 >                upv[0] = nrm[0]*nrm[1]*(nrm[2] - 1.);
393 >                upv[1] = nrm[0]*nrm[0] + nrm[1]*nrm[1]*nrm[2];
394 >                upv[2] = -nrm[1]*(nrm[0]*nrm[0] + nrm[1]*nrm[1]);
395 >                if (SDcompXform(vMtx, nrm, upv) != SDEnone)
396 >                        continue;
397 >                nrota = addrot(rotargs, vMtx[0], vMtx[1], vMtx[2]);
398 >                if (front_comp) {
399 >                        cvt_sposition(sorg, ivec, 1);
400                          printf("\nvoid colorfunc scale_pat\n");
401 <                        printf("%d bsdf_red bsdf_grn bsdf_blu bsdf2rad.cal\n\t%s\n0\n0\n",
402 <                                        4+nxfa, xfargs);
401 >                        printf("10 bsdf_red bsdf_grn bsdf_blu bsdf2rad.cal\n");
402 >                        printf("\t-s %f -t %f %f %f\n0\n0\n",
403 >                                        bsdf_rad, sorg[0], sorg[1], sorg[2]);
404                          printf("\nscale_pat glow scale_mat\n0\n0\n4 1 1 1 0\n");
405 <                        SDcompXform(vMtx, ivec, Yaxis);
406 <                        nxfa = addrot(xfargs, vMtx[0], vMtx[1], vMtx[2]);
407 <                        sprintf(xfargs+strlen(xfargs), " -s %f -t %f %f %f",
397 <                                        scalef, ivec[0]*sph_rad - sph_xoffset,
398 <                                        ivec[1]*sph_rad, ivec[2]*sph_rad);
399 <                        nxfa += 6;
405 >                }
406 >                if (front_comp & SDsampR) {
407 >                        put_mirror_arrow(sorg, nrm);
408                          fname = tfile_name(frpref, dsuffix, i);
409 <                        sprintf(cmdbuf, "gensurf scale_mat %s%d %s %s %s %d %d | xform -mx -my %s",
410 <                                        frpref, i+1, fname, fname, fname, SAMPRES-1, SAMPRES-1,
411 <                                        xfargs);
409 >                        sprintf(cmdbuf,
410 >                "gensurf scale_mat %s%d %s %s %s %d %d | xform %s -s %f -t %f %f %f",
411 >                                        frpref, i, fname, fname, fname, SAMPRES-1, SAMPRES-1,
412 >                                        rotargs, scalef, sorg[0], sorg[1], sorg[2]);
413                          if (!run_cmd(cmdbuf))
414                                  return(0);
415                  }
416 <        if (front_comp & SDsampT)
417 <                for (i = 0; i < NINCIDENT; i++) {
409 <                        get_ivector(ivec, i);
410 <                        put_trans_arrow(ivec, 1);
411 <                        sprintf(xfargs, "-s %f -t %f %f %f", bsdf_rad,
412 <                                        ivec[0]*sph_rad - sph_xoffset,
413 <                                        ivec[1]*sph_rad, ivec[2]*sph_rad);
414 <                        nxfa = 6;
415 <                        printf("\nvoid colorfunc scale_pat\n");
416 <                        printf("%d bsdf_red bsdf_grn bsdf_blu bsdf2rad.cal\n\t%s\n0\n0\n",
417 <                                        4+nxfa, xfargs);
418 <                        printf("\nscale_pat glow scale_mat\n0\n0\n4 1 1 1 0\n");
419 <                        SDcompXform(vMtx, ivec, Yaxis);
420 <                        nxfa = addrot(xfargs, vMtx[0], vMtx[1], vMtx[2]);
421 <                        sprintf(xfargs+strlen(xfargs), " -s %f -t %f %f %f",
422 <                                        scalef, ivec[0]*sph_rad - sph_xoffset,
423 <                                        ivec[1]*sph_rad, ivec[2]*sph_rad);
424 <                        nxfa += 6;
416 >                if (front_comp & SDsampT) {
417 >                        put_trans_arrow(sorg);
418                          fname = tfile_name(ftpref, dsuffix, i);
419 <                        sprintf(cmdbuf, "gensurf scale_mat %s%d %s %s %s %d %d | xform -I -mx -my %s",
420 <                                        ftpref, i+1, fname, fname, fname, SAMPRES-1, SAMPRES-1,
421 <                                        xfargs);
419 >                        sprintf(cmdbuf,
420 >                "gensurf scale_mat %s%d %s %s %s %d %d | xform -I %s -s %f -t %f %f %f",
421 >                                        ftpref, i, fname, fname, fname, SAMPRES-1, SAMPRES-1,
422 >                                        rotargs, scalef, sorg[0], sorg[1], sorg[2]);
423                          if (!run_cmd(cmdbuf))
424                                  return(0);
425                  }
426 <        if (back_comp & SDsampR)
427 <                for (i = 0; i < NINCIDENT; i++) {
434 <                        get_ivector(ivec, i);
435 <                        put_mirror_arrow(ivec, -1);
436 <                        fname = tfile_name(brpref, dsuffix, i);
437 <                        sprintf(xfargs, "-s %f -t %f %f %f", bsdf_rad,
438 <                                        ivec[0]*sph_rad + sph_xoffset,
439 <                                        ivec[1]*sph_rad, ivec[2]*sph_rad);
440 <                        nxfa = 6;
426 >                if (back_comp) {
427 >                        cvt_sposition(sorg, ivec, -1);
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);
429 >                        printf("10 bsdf_red bsdf_grn bsdf_blu bsdf2rad.cal\n");
430 >                        printf("\t-s %f -t %f %f %f\n0\n0\n",
431 >                                        bsdf_rad, sorg[0], sorg[1], sorg[2]);
432                          printf("\nscale_pat glow scale_mat\n0\n0\n4 1 1 1 0\n");
433 <                        SDcompXform(vMtx, ivec, Yaxis);
434 <                        nxfa = addrot(xfargs, vMtx[0], vMtx[1], vMtx[2]);
435 <                        sprintf(xfargs+strlen(xfargs), " -s %f -t %f %f %f",
448 <                                        scalef, ivec[0]*sph_rad + sph_xoffset,
449 <                                        ivec[1]*sph_rad, ivec[2]*sph_rad);
450 <                        nxfa += 6;
433 >                }
434 >                if (back_comp & SDsampR) {
435 >                        put_mirror_arrow(sorg, nrm);
436                          fname = tfile_name(brpref, dsuffix, i);
437 <                        sprintf(cmdbuf, "gensurf scale_mat %s%d %s %s %s %d %d | xform -I -ry 180 -mx -my %s",
438 <                                        brpref, i+1, fname, fname, fname, SAMPRES-1, SAMPRES-1,
439 <                                        xfargs);
437 >                        sprintf(cmdbuf,
438 >                "gensurf scale_mat %s%d %s %s %s %d %d | xform -I -ry 180 %s -s %f -t %f %f %f",
439 >                                        brpref, i, fname, fname, fname, SAMPRES-1, SAMPRES-1,
440 >                                        rotargs, scalef, sorg[0], sorg[1], sorg[2]);
441                          if (!run_cmd(cmdbuf))
442                                  return(0);
443                  }
444 <        if (back_comp & SDsampT)
445 <                for (i = 0; i < NINCIDENT; i++) {
460 <                        get_ivector(ivec, i);
461 <                        put_trans_arrow(ivec, -1);
444 >                if (back_comp & SDsampT) {
445 >                        put_trans_arrow(sorg);
446                          fname = tfile_name(btpref, dsuffix, i);
447 <                        sprintf(xfargs, "-s %f -t %f %f %f", bsdf_rad,
448 <                                        ivec[0]*sph_rad + sph_xoffset,
449 <                                        ivec[1]*sph_rad, ivec[2]*sph_rad);
450 <                        nxfa = 6;
467 <                        printf("\nvoid colorfunc scale_pat\n");
468 <                        printf("%d bsdf_red bsdf_grn bsdf_blu bsdf2rad.cal\n\t%s\n0\n0\n",
469 <                                        4+nxfa, xfargs);
470 <                        printf("\nscale_pat glow scale_mat\n0\n0\n4 1 1 1 0\n");
471 <                        SDcompXform(vMtx, ivec, Yaxis);
472 <                        nxfa = addrot(xfargs, vMtx[0], vMtx[1], vMtx[2]);
473 <                        sprintf(xfargs+strlen(xfargs), " -s %f -t %f %f %f",
474 <                                        scalef, ivec[0]*sph_rad + sph_xoffset,
475 <                                        ivec[1]*sph_rad, ivec[2]*sph_rad);
476 <                        nxfa += 6;
477 <                        fname = tfile_name(btpref, dsuffix, i);
478 <                        sprintf(cmdbuf, "gensurf scale_mat %s%d %s %s %s %d %d | xform -ry 180 -mx -my %s",
479 <                                        btpref, i+1, fname, fname, fname, SAMPRES-1, SAMPRES-1,
480 <                                        xfargs);
447 >                        sprintf(cmdbuf,
448 >                "gensurf scale_mat %s%d %s %s %s %d %d | xform -ry 180 %s -s %f -t %f %f %f",
449 >                                        btpref, i, fname, fname, fname, SAMPRES-1, SAMPRES-1,
450 >                                        rotargs, scalef, sorg[0], sorg[1], sorg[2]);
451                          if (!run_cmd(cmdbuf))
452                                  return(0);
453                  }
454 +        }
455          return(1);
456   }
457  
# Line 494 | Line 465 | put_matBSDF(const char *XMLfile)
465                  printf("\n# Simplified material because we have no XML input\n");
466                  printf("\nvoid brightfunc latlong\n2 latlong bsdf2rad.cal\n0\n0\n");
467                  if ((front_comp|back_comp) & SDsampT)
468 <                        printf("\nlatlong trans %s\n0\n0\n7 .75 .75 .75 0 0 .5 .8\n",
469 <                                        sph_mat);
468 >                        printf("\nlatlong trans %s\n0\n0\n7 .75 .75 .75 0 .04 .5 .8\n",
469 >                                        sph_fmat);
470                  else
471                          printf("\nlatlong plastic %s\n0\n0\n5 .5 .5 .5 0 0\n",
472 <                                        sph_mat);
472 >                                        sph_fmat);
473 >                printf("\ninherit alias %s %s\n", sph_bmat, sph_fmat);
474                  return;
475          }
476          switch (XMLfile[0]) {           /* avoid RAYPATH search */
477          case '.':
478 +        case '~':
479          CASEDIRSEP:
480                  curdir = "";
481                  break;
# Line 511 | Line 484 | put_matBSDF(const char *XMLfile)
484                  exit(1);
485                  break;
486          }
487 <        printf("\n# Actual BSDF material for rendering the hemispheres\n");
488 <        printf("\nvoid BSDF BSDFmat\n6 0 \"%s%s\" 0 1 0 .\n0\n0\n",
487 >        printf("\n# Actual BSDF materials for rendering the hemispheres\n");
488 >        printf("\nvoid BSDF BSDF_f\n6 0 \"%s%s\" upx upy upz bsdf2rad.cal\n0\n0\n",
489                          curdir, XMLfile);
490          printf("\nvoid plastic black\n0\n0\n5 0 0 0 0 0\n");
491 <        printf("\nvoid mixfunc %s\n4 BSDFmat black latlong bsdf2rad.cal\n0\n0\n",
492 <                        sph_mat);
491 >        printf("\nvoid mixfunc %s\n4 BSDF_f black latlong bsdf2rad.cal\n0\n0\n",
492 >                        sph_fmat);
493 >        printf("\nvoid BSDF BSDF_b\n8 0 \"%s%s\" upx upy upz bsdf2rad.cal -ry 180\n0\n0\n",
494 >                        curdir, XMLfile);
495 >        printf("\nvoid mixfunc %s\n4 BSDF_b black latlong bsdf2rad.cal\n0\n0\n",
496 >                        sph_bmat);
497   }
498  
499   /* Put out overhead parallel light source */
# Line 524 | Line 501 | static void
501   put_source(void)
502   {
503          printf("\n# Overhead parallel light source\n");
504 <        printf("\nvoid light bright\n0\n0\n3 1000 1000 1000\n");
504 >        printf("\nvoid light bright\n0\n0\n3 2500 2500 2500\n");
505          printf("\nbright source light\n0\n0\n4 0 0 1 2\n");
506          printf("\n# Material used for labels\n");
507          printf("\nvoid trans vellum\n0\n0\n7 1 1 1 0 0 .5 0\n");
# Line 534 | Line 511 | put_source(void)
511   static void
512   put_hemispheres(void)
513   {
514 +        const int       nsegs = 131;
515 +
516          printf("\n# Hemisphere(s) for showing BSDF appearance (if XML file)\n");
538 –        printf("\nvoid antimatter anti_sph\n2 void %s\n0\n0\n", sph_mat);
517          if (front_comp) {
518 <                printf("\n%s sphere Front\n0\n0\n4 %f 0 0 %f\n",
519 <                                sph_mat, -sph_xoffset, sph_rad);
520 <                printf("\n!genbox anti_sph sph_eraser %f %f %f | xform -t %f %f %f\n",
543 <                                2.02*sph_rad, 2.02*sph_rad, 1.02*sph_rad,
544 <                                -1.01*sph_rad - sph_xoffset, -1.01*sph_rad, -1.01*sph_rad);
518 >                printf(
519 > "\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",
520 >                                sph_fmat, nsegs, sph_rad, 0.5*PI, sph_xoffset);
521                  printf("\nvoid brighttext front_text\n3 helvet.fnt . FRONT\n0\n");
522                  printf("12\n\t%f %f 0\n\t%f 0 0\n\t0 %f 0\n\t.01 1 -.1\n",
523 <                                -.22*sph_rad - sph_xoffset, -1.2*sph_rad,
523 >                                -.22*sph_rad + sph_xoffset, -1.4*sph_rad,
524                                  .35/5.*sph_rad, -1.6*.35/5.*sph_rad);
525                  printf("\nfront_text alias front_label_mat vellum\n");
526                  printf("\nfront_label_mat polygon front_label\n0\n0\n12");
527                  printf("\n\t%f %f 0\n\t%f %f 0\n\t%f %f 0\n\t%f %f 0\n",
528 <                                -.25*sph_rad - sph_xoffset, -1.1*sph_rad,
529 <                                -.25*sph_rad - sph_xoffset, (-1.2-1.6*.35/5.-.1)*sph_rad,
530 <                                .25*sph_rad - sph_xoffset, (-1.2-1.6*.35/5.-.1)*sph_rad,
531 <                                .25*sph_rad - sph_xoffset, -1.1*sph_rad );
528 >                                -.25*sph_rad + sph_xoffset, -1.3*sph_rad,
529 >                                -.25*sph_rad + sph_xoffset, (-1.4-1.6*.35/5.-.1)*sph_rad,
530 >                                .25*sph_rad + sph_xoffset, (-1.4-1.6*.35/5.-.1)*sph_rad,
531 >                                .25*sph_rad + sph_xoffset, -1.3*sph_rad );
532          }
533          if (back_comp) {
534 <                printf("\n%s bubble Back\n0\n0\n4 %f 0 0 %f\n",
535 <                                sph_mat, sph_xoffset, sph_rad);
536 <                printf("\n!genbox anti_sph sph_eraser %f %f %f | xform -t %f %f %f\n",
561 <                                2.02*sph_rad, 2.02*sph_rad, 1.02*sph_rad,
562 <                                -1.01*sph_rad + sph_xoffset, -1.01*sph_rad, -1.01*sph_rad);
534 >                printf(
535 > "\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",
536 >                                sph_bmat, nsegs, sph_rad, 0.5*PI, -sph_xoffset);
537                  printf("\nvoid brighttext back_text\n3 helvet.fnt . BACK\n0\n");
538                  printf("12\n\t%f %f 0\n\t%f 0 0\n\t0 %f 0\n\t.01 1 -.1\n",
539 <                                -.22*sph_rad + sph_xoffset, -1.2*sph_rad,
539 >                                -.22*sph_rad - sph_xoffset, -1.4*sph_rad,
540                                  .35/4.*sph_rad, -1.6*.35/4.*sph_rad);
541                  printf("\nback_text alias back_label_mat vellum\n");
542                  printf("\nback_label_mat polygon back_label\n0\n0\n12");
543                  printf("\n\t%f %f 0\n\t%f %f 0\n\t%f %f 0\n\t%f %f 0\n",
544 <                                -.25*sph_rad + sph_xoffset, -1.1*sph_rad,
545 <                                -.25*sph_rad + sph_xoffset, (-1.2-1.6*.35/4.-.1)*sph_rad,
546 <                                .25*sph_rad + sph_xoffset, (-1.2-1.6*.35/4.-.1)*sph_rad,
547 <                                .25*sph_rad + sph_xoffset, -1.1*sph_rad );
544 >                                -.25*sph_rad - sph_xoffset, -1.3*sph_rad,
545 >                                -.25*sph_rad - sph_xoffset, (-1.4-1.6*.35/4.-.1)*sph_rad,
546 >                                .25*sph_rad - sph_xoffset, (-1.4-1.6*.35/4.-.1)*sph_rad,
547 >                                .25*sph_rad - sph_xoffset, -1.3*sph_rad );
548          }
549   }
550  
# Line 579 | Line 553 | static void
553   put_scale(void)
554   {
555          const double    max_log10 = log10(overall_max);
556 <        const double    leg_width = 2.*.75*(sph_xoffset - sph_rad);
556 >        const double    leg_width = 2.*.75*(fabs(sph_xoffset) - sph_rad);
557          const double    leg_height = 2.*sph_rad;
558          const int       text_lines = 6;
559 <        const int       text_digits = 7;
559 >        const int       text_digits = 8;
560          char            fmt[16];
561          int             i;
562  
563          printf("\n# BSDF legend with falsecolor scale\n");
564          printf("\nvoid colorfunc lscale\n10 sca_red(Py) sca_grn(Py) sca_blu(Py)");
565          printf("\n\tbsdf2rad.cal -s %f -t 0 %f 0\n0\n0\n", leg_height, -.5*leg_height);
566 <        sprintf(fmt, "%%%dg", text_digits-2);
566 >        sprintf(fmt, "%%.%df", text_digits-3);
567          for (i = 0; i < text_lines; i++) {
568                  char    vbuf[16];
569                  sprintf(vbuf, fmt, pow(10., (i+.5)/text_lines*(max_log10-min_log10)+min_log10));
570                  printf("\nlscale brighttext lscale\n");
571                  printf("3 helvet.fnt . %s\n0\n12\n", vbuf);
572                  printf("\t%f %f 0\n", -.45*leg_width, ((i+.9)/text_lines-.5)*leg_height);
573 <                printf("\t%f 0 0\n", .9*leg_width/strlen(vbuf));
573 >                printf("\t%f 0 0\n", .8*leg_width/strlen(vbuf));
574                  printf("\t0 %f 0\n", -.9/text_lines*leg_height);
575                  printf("\t.01 1 -.1\n");
576          }
# Line 623 | Line 597 | put_scale(void)
597                  return;
598          printf("\nvoid brighttext BSDFname\n3 helvet.fnt . \"%s\"\n0\n12\n", bsdf_name);
599          printf("\t%f %f 0\n", -.95*leg_width, -.6*leg_height);
600 <        printf("\t%f 0 0\n", .4/strlen(bsdf_name)*leg_width);
600 >        printf("\t%f 0 0\n", 1.8/strlen(bsdf_name)*leg_width);
601          printf("\t0 %f 0\n", -.1*leg_height);
602          printf("\t.01 1 -.1\n");
603          printf("\nBSDFname alias name_mat vellum\n");
# Line 678 | Line 652 | convert_mgf(const char *mgfdata)
652                  sprintf(cmdbuf, "xform -t %f %f %f -s %f -t %f %f 0 %s",
653                                  -.5*(xmin+xmax), -.5*(ymin+ymax), -zmax,
654                                  1.5*sph_rad/max_dim,
655 <                                -sph_xoffset, -2.5*sph_rad,
655 >                                sph_xoffset, -2.5*sph_rad,
656                                  radfn);
657                  if (!run_cmd(cmdbuf))
658                          return;
# Line 688 | Line 662 | convert_mgf(const char *mgfdata)
662                  sprintf(cmdbuf, "xform -t %f %f %f -s %f -ry 180 -t %f %f 0 %s",
663                                  -.5*(xmin+xmax), -.5*(ymin+ymax), -zmin,
664                                  1.5*sph_rad/max_dim,
665 <                                sph_xoffset, -2.5*sph_rad,
665 >                                -sph_xoffset, -2.5*sph_rad,
666                                  radfn);
667                  if (!run_cmd(cmdbuf))
668                          return;
# Line 781 | Line 755 | main(int argc, char *argv[])
755                  strcpy(bsdf_manuf, myBSDF.makr);
756                  put_matBSDF(argv[1]);
757          } else {
758 <                FILE    *fp;
758 >                FILE    *fp[4];
759 >                if (argc > 5) {
760 >                        fprintf(stderr, "%s: more than 4 hemispheres!\n", progname);
761 >                        return(1);
762 >                }
763                  for (n = 1; n < argc; n++) {
764 <                        fp = fopen(argv[n], "rb");
765 <                        if (fp == NULL) {
764 >                        fp[n-1] = fopen(argv[n], "rb");
765 >                        if (fp[n-1] == NULL) {
766                                  fprintf(stderr, "%s: cannot open BSDF interpolant '%s'\n",
767                                                  progname, argv[n]);
768                                  return(1);
769                          }
770 <                        if (getheader(fp, rbf_headline, NULL) < 0) {
770 >                        if (getheader(fp[n-1], rbf_headline, NULL) < 0) {
771                                  fprintf(stderr, "%s: bad BSDF interpolant '%s'\n",
772                                                  progname, argv[n]);
773                                  return(1);
774                          }
797 –                        fclose(fp);
775                  }
776                  set_minlog();
777                  for (n = 1; n < argc; n++) {
778 <                        fp = fopen(argv[n], "rb");
779 <                        if (!load_bsdf_rep(fp))
778 >                        if (fseek(fp[n-1], 0L, SEEK_SET) < 0) {
779 >                                fprintf(stderr, "%s: cannot seek on '%s'\n",
780 >                                                progname, argv[n]);
781                                  return(1);
782 <                        fclose(fp);
782 >                        }
783 >                        if (!load_bsdf_rep(fp[n-1]))
784 >                                return(1);
785 >                        fclose(fp[n-1]);
786                          if (!build_wRBF())
787                                  return(1);
788                  }
# Line 812 | Line 793 | main(int argc, char *argv[])
793          put_scale();
794          if (inpXML && myBSDF.mgf)
795                  convert_mgf(myBSDF.mgf);
796 <        if (!put_BSDFs())
796 >        if (!put_BSDFs())               /* most of the output happens here */
797                  return(1);
798          cleanup_tmp();
799          return(0);

Diff Legend

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