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) |
50 |
|
char *progname; |
51 |
|
|
52 |
|
/* Get Fibonacci sphere vector (0 to NINCIDENT-1) */ |
53 |
< |
static void |
53 |
> |
static RREAL * |
54 |
|
get_ivector(FVECT iv, int i) |
55 |
|
{ |
56 |
– |
const double zstep = 1./(2*NINCIDENT-1); |
56 |
|
const double phistep = PI*(3. - 2.236067978); |
57 |
|
double r; |
58 |
|
|
59 |
< |
iv[2] = 1. - (i+.5)*zstep; |
59 |
> |
iv[2] = 1. - (i+.5)*(1./NINCIDENT); |
60 |
|
r = sqrt(1. - iv[2]*iv[2]); |
61 |
|
iv[0] = r * cos((i+1.)*phistep); |
62 |
|
iv[1] = r * sin((i+1.)*phistep); |
63 |
+ |
|
64 |
+ |
return(iv); |
65 |
|
} |
66 |
|
|
67 |
+ |
/* Convert incident vector into sphere position */ |
68 |
+ |
static RREAL * |
69 |
+ |
cvt_sposition(FVECT sp, const FVECT iv, int inc_side) |
70 |
+ |
{ |
71 |
+ |
sp[0] = -iv[0]*sph_rad - inc_side*sph_xoffset; |
72 |
+ |
sp[1] = -iv[1]*sph_rad; |
73 |
+ |
sp[2] = iv[2]*sph_rad; |
74 |
+ |
|
75 |
+ |
return(sp); |
76 |
+ |
} |
77 |
+ |
|
78 |
|
/* Get temporary file name */ |
79 |
|
static char * |
80 |
|
tfile_name(const char *prefix, const char *suffix, int i) |
261 |
|
RBFNODE *rbf; |
262 |
|
get_ivector(ivec, i); |
263 |
|
if (input_orient < 0) { |
264 |
< |
ivec[0] = -ivec[0]; ivec[1] = -ivec[1]; ivec[2] = -ivec[2]; |
264 |
> |
ivec[0] = -ivec[0]; ivec[2] = -ivec[2]; |
265 |
|
} |
266 |
|
rbf = advect_rbf(ivec, 15000); |
267 |
|
if (!plotRBF(tfile_name(pref, dsuffix, i), rbf)) |
280 |
|
FVECT origin, refl; |
281 |
|
int i; |
282 |
|
|
283 |
< |
for (i = 3; i--; ) origin[i] = ivec[i]*sph_rad; |
272 |
< |
origin[0] -= inc_side*sph_xoffset; |
283 |
> |
cvt_sposition(origin, ivec, inc_side); |
284 |
|
|
285 |
< |
refl[0] = 2.*ivec[2]*ivec[0]; |
286 |
< |
refl[1] = 2.*ivec[2]*ivec[1]; |
285 |
> |
refl[0] = -2.*ivec[2]*ivec[0]; |
286 |
> |
refl[1] = -2.*ivec[2]*ivec[1]; |
287 |
|
refl[2] = 2.*ivec[2]*ivec[2] - 1.; |
288 |
|
|
289 |
|
printf("\n# Mirror arrow\n"); |
319 |
|
FVECT origin; |
320 |
|
int i; |
321 |
|
|
322 |
< |
for (i = 3; i--; ) origin[i] = ivec[i]*sph_rad; |
312 |
< |
origin[0] -= inc_side*sph_xoffset; |
322 |
> |
cvt_sposition(origin, ivec, inc_side); |
323 |
|
|
324 |
|
printf("\n# Transmission arrow\n"); |
325 |
|
printf("\narrow_mat cylinder trans_dir\n0\n0\n7"); |
375 |
|
put_BSDFs(void) |
376 |
|
{ |
377 |
|
const double scalef = bsdf_rad/(log10(overall_max) - min_log10); |
378 |
< |
FVECT ivec; |
378 |
> |
FVECT ivec, sorg; |
379 |
|
RREAL vMtx[3][3]; |
380 |
|
char *fname; |
381 |
|
char cmdbuf[256]; |
389 |
|
printf("\nvoid glow arrow_glow\n0\n0\n4 1 0 1 0\n"); |
390 |
|
printf("\nvoid mixfunc arrow_mat\n4 arrow_glow void .5 .\n0\n0\n"); |
391 |
|
|
392 |
< |
if (front_comp & SDsampR) |
392 |
> |
if (front_comp & SDsampR) /* front reflection */ |
393 |
|
for (i = 0; i < NINCIDENT; i++) { |
394 |
|
get_ivector(ivec, i); |
395 |
|
put_mirror_arrow(ivec, 1); |
396 |
+ |
cvt_sposition(sorg, ivec, 1); |
397 |
+ |
ivec[0] = -ivec[0]; ivec[1] = -ivec[1]; /* normal */ |
398 |
|
sprintf(xfargs, "-s %f -t %f %f %f", bsdf_rad, |
399 |
< |
ivec[0]*sph_rad - sph_xoffset, |
388 |
< |
ivec[1]*sph_rad, ivec[2]*sph_rad); |
399 |
> |
sorg[0], sorg[1], sorg[2]); |
400 |
|
nxfa = 6; |
401 |
|
printf("\nvoid colorfunc scale_pat\n"); |
402 |
|
printf("%d bsdf_red bsdf_grn bsdf_blu bsdf2rad.cal\n\t%s\n0\n0\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", |
408 |
< |
scalef, ivec[0]*sph_rad - sph_xoffset, |
398 |
< |
ivec[1]*sph_rad, ivec[2]*sph_rad); |
408 |
> |
scalef, sorg[0], sorg[1], sorg[2]); |
409 |
|
nxfa += 6; |
410 |
|
fname = tfile_name(frpref, dsuffix, i); |
411 |
< |
sprintf(cmdbuf, "gensurf scale_mat %s%d %s %s %s %d %d | xform -mx -my %s", |
411 |
> |
sprintf(cmdbuf, "gensurf scale_mat %s%d %s %s %s %d %d | xform %s", |
412 |
|
frpref, i+1, fname, fname, fname, SAMPRES-1, SAMPRES-1, |
413 |
|
xfargs); |
414 |
|
if (!run_cmd(cmdbuf)) |
415 |
|
return(0); |
416 |
|
} |
417 |
< |
if (front_comp & SDsampT) |
417 |
> |
if (front_comp & SDsampT) /* front transmission */ |
418 |
|
for (i = 0; i < NINCIDENT; i++) { |
419 |
|
get_ivector(ivec, i); |
420 |
|
put_trans_arrow(ivec, 1); |
421 |
+ |
cvt_sposition(sorg, ivec, 1); |
422 |
+ |
ivec[0] = -ivec[0]; ivec[1] = -ivec[1]; /* normal */ |
423 |
|
sprintf(xfargs, "-s %f -t %f %f %f", bsdf_rad, |
424 |
< |
ivec[0]*sph_rad - sph_xoffset, |
413 |
< |
ivec[1]*sph_rad, ivec[2]*sph_rad); |
424 |
> |
sorg[0], sorg[1], sorg[2]); |
425 |
|
nxfa = 6; |
426 |
|
printf("\nvoid colorfunc scale_pat\n"); |
427 |
|
printf("%d bsdf_red bsdf_grn bsdf_blu bsdf2rad.cal\n\t%s\n0\n0\n", |
430 |
|
SDcompXform(vMtx, ivec, Yaxis); |
431 |
|
nxfa = addrot(xfargs, vMtx[0], vMtx[1], vMtx[2]); |
432 |
|
sprintf(xfargs+strlen(xfargs), " -s %f -t %f %f %f", |
433 |
< |
scalef, ivec[0]*sph_rad - sph_xoffset, |
423 |
< |
ivec[1]*sph_rad, ivec[2]*sph_rad); |
433 |
> |
scalef, sorg[0], sorg[1], sorg[2]); |
434 |
|
nxfa += 6; |
435 |
|
fname = tfile_name(ftpref, dsuffix, i); |
436 |
< |
sprintf(cmdbuf, "gensurf scale_mat %s%d %s %s %s %d %d | xform -I -mx -my %s", |
436 |
> |
sprintf(cmdbuf, "gensurf scale_mat %s%d %s %s %s %d %d | xform -I %s", |
437 |
|
ftpref, i+1, fname, fname, fname, SAMPRES-1, SAMPRES-1, |
438 |
|
xfargs); |
439 |
|
if (!run_cmd(cmdbuf)) |
440 |
|
return(0); |
441 |
|
} |
442 |
< |
if (back_comp & SDsampR) |
442 |
> |
if (back_comp & SDsampR) /* rear reflection */ |
443 |
|
for (i = 0; i < NINCIDENT; i++) { |
444 |
|
get_ivector(ivec, i); |
445 |
|
put_mirror_arrow(ivec, -1); |
446 |
< |
fname = tfile_name(brpref, dsuffix, i); |
446 |
> |
cvt_sposition(sorg, ivec, -1); |
447 |
> |
ivec[0] = -ivec[0]; ivec[1] = -ivec[1]; /* normal */ |
448 |
|
sprintf(xfargs, "-s %f -t %f %f %f", bsdf_rad, |
449 |
< |
ivec[0]*sph_rad + sph_xoffset, |
439 |
< |
ivec[1]*sph_rad, ivec[2]*sph_rad); |
449 |
> |
sorg[0], sorg[1], sorg[2]); |
450 |
|
nxfa = 6; |
451 |
|
printf("\nvoid colorfunc scale_pat\n"); |
452 |
|
printf("%d bsdf_red bsdf_grn bsdf_blu bsdf2rad.cal\n\t%s\n0\n0\n", |
455 |
|
SDcompXform(vMtx, ivec, Yaxis); |
456 |
|
nxfa = addrot(xfargs, vMtx[0], vMtx[1], vMtx[2]); |
457 |
|
sprintf(xfargs+strlen(xfargs), " -s %f -t %f %f %f", |
458 |
< |
scalef, ivec[0]*sph_rad + sph_xoffset, |
449 |
< |
ivec[1]*sph_rad, ivec[2]*sph_rad); |
458 |
> |
scalef, sorg[0], sorg[1], sorg[2]); |
459 |
|
nxfa += 6; |
460 |
|
fname = tfile_name(brpref, dsuffix, i); |
461 |
< |
sprintf(cmdbuf, "gensurf scale_mat %s%d %s %s %s %d %d | xform -I -ry 180 -mx -my %s", |
461 |
> |
sprintf(cmdbuf, "gensurf scale_mat %s%d %s %s %s %d %d | xform -I -ry 180 %s", |
462 |
|
brpref, i+1, fname, fname, fname, SAMPRES-1, SAMPRES-1, |
463 |
|
xfargs); |
464 |
|
if (!run_cmd(cmdbuf)) |
465 |
|
return(0); |
466 |
|
} |
467 |
< |
if (back_comp & SDsampT) |
467 |
> |
if (back_comp & SDsampT) /* rear transmission */ |
468 |
|
for (i = 0; i < NINCIDENT; i++) { |
469 |
|
get_ivector(ivec, i); |
470 |
|
put_trans_arrow(ivec, -1); |
471 |
< |
fname = tfile_name(btpref, dsuffix, i); |
471 |
> |
cvt_sposition(sorg, ivec, -1); |
472 |
> |
ivec[0] = -ivec[0]; ivec[1] = -ivec[1]; /* normal */ |
473 |
|
sprintf(xfargs, "-s %f -t %f %f %f", bsdf_rad, |
474 |
< |
ivec[0]*sph_rad + sph_xoffset, |
465 |
< |
ivec[1]*sph_rad, ivec[2]*sph_rad); |
474 |
> |
sorg[0], sorg[1], sorg[2]); |
475 |
|
nxfa = 6; |
476 |
|
printf("\nvoid colorfunc scale_pat\n"); |
477 |
|
printf("%d bsdf_red bsdf_grn bsdf_blu bsdf2rad.cal\n\t%s\n0\n0\n", |
480 |
|
SDcompXform(vMtx, ivec, Yaxis); |
481 |
|
nxfa = addrot(xfargs, vMtx[0], vMtx[1], vMtx[2]); |
482 |
|
sprintf(xfargs+strlen(xfargs), " -s %f -t %f %f %f", |
483 |
< |
scalef, ivec[0]*sph_rad + sph_xoffset, |
475 |
< |
ivec[1]*sph_rad, ivec[2]*sph_rad); |
483 |
> |
scalef, sorg[0], sorg[1], sorg[2]); |
484 |
|
nxfa += 6; |
485 |
|
fname = tfile_name(btpref, dsuffix, i); |
486 |
< |
sprintf(cmdbuf, "gensurf scale_mat %s%d %s %s %s %d %d | xform -ry 180 -mx -my %s", |
486 |
> |
sprintf(cmdbuf, "gensurf scale_mat %s%d %s %s %s %d %d | xform -ry 180 %s", |
487 |
|
btpref, i+1, fname, fname, fname, SAMPRES-1, SAMPRES-1, |
488 |
|
xfargs); |
489 |
|
if (!run_cmd(cmdbuf)) |
532 |
|
put_source(void) |
533 |
|
{ |
534 |
|
printf("\n# Overhead parallel light source\n"); |
535 |
< |
printf("\nvoid light bright\n0\n0\n3 1000 1000 1000\n"); |
535 |
> |
printf("\nvoid light bright\n0\n0\n3 2000 2000 2000\n"); |
536 |
|
printf("\nbright source light\n0\n0\n4 0 0 1 2\n"); |
537 |
|
printf("\n# Material used for labels\n"); |
538 |
|
printf("\nvoid trans vellum\n0\n0\n7 1 1 1 0 0 .5 0\n"); |
552 |
|
-1.01*sph_rad - sph_xoffset, -1.01*sph_rad, -1.01*sph_rad); |
553 |
|
printf("\nvoid brighttext front_text\n3 helvet.fnt . FRONT\n0\n"); |
554 |
|
printf("12\n\t%f %f 0\n\t%f 0 0\n\t0 %f 0\n\t.01 1 -.1\n", |
555 |
< |
-.22*sph_rad - sph_xoffset, -1.2*sph_rad, |
555 |
> |
-.22*sph_rad - sph_xoffset, -1.4*sph_rad, |
556 |
|
.35/5.*sph_rad, -1.6*.35/5.*sph_rad); |
557 |
|
printf("\nfront_text alias front_label_mat vellum\n"); |
558 |
|
printf("\nfront_label_mat polygon front_label\n0\n0\n12"); |
559 |
|
printf("\n\t%f %f 0\n\t%f %f 0\n\t%f %f 0\n\t%f %f 0\n", |
560 |
< |
-.25*sph_rad - sph_xoffset, -1.1*sph_rad, |
561 |
< |
-.25*sph_rad - sph_xoffset, (-1.2-1.6*.35/5.-.1)*sph_rad, |
562 |
< |
.25*sph_rad - sph_xoffset, (-1.2-1.6*.35/5.-.1)*sph_rad, |
563 |
< |
.25*sph_rad - sph_xoffset, -1.1*sph_rad ); |
560 |
> |
-.25*sph_rad - sph_xoffset, -1.3*sph_rad, |
561 |
> |
-.25*sph_rad - sph_xoffset, (-1.4-1.6*.35/5.-.1)*sph_rad, |
562 |
> |
.25*sph_rad - sph_xoffset, (-1.4-1.6*.35/5.-.1)*sph_rad, |
563 |
> |
.25*sph_rad - sph_xoffset, -1.3*sph_rad ); |
564 |
|
} |
565 |
|
if (back_comp) { |
566 |
|
printf("\n%s bubble Back\n0\n0\n4 %f 0 0 %f\n", |
570 |
|
-1.01*sph_rad + sph_xoffset, -1.01*sph_rad, -1.01*sph_rad); |
571 |
|
printf("\nvoid brighttext back_text\n3 helvet.fnt . BACK\n0\n"); |
572 |
|
printf("12\n\t%f %f 0\n\t%f 0 0\n\t0 %f 0\n\t.01 1 -.1\n", |
573 |
< |
-.22*sph_rad + sph_xoffset, -1.2*sph_rad, |
573 |
> |
-.22*sph_rad + sph_xoffset, -1.4*sph_rad, |
574 |
|
.35/4.*sph_rad, -1.6*.35/4.*sph_rad); |
575 |
|
printf("\nback_text alias back_label_mat vellum\n"); |
576 |
|
printf("\nback_label_mat polygon back_label\n0\n0\n12"); |
577 |
|
printf("\n\t%f %f 0\n\t%f %f 0\n\t%f %f 0\n\t%f %f 0\n", |
578 |
< |
-.25*sph_rad + sph_xoffset, -1.1*sph_rad, |
579 |
< |
-.25*sph_rad + sph_xoffset, (-1.2-1.6*.35/4.-.1)*sph_rad, |
580 |
< |
.25*sph_rad + sph_xoffset, (-1.2-1.6*.35/4.-.1)*sph_rad, |
581 |
< |
.25*sph_rad + sph_xoffset, -1.1*sph_rad ); |
578 |
> |
-.25*sph_rad + sph_xoffset, -1.3*sph_rad, |
579 |
> |
-.25*sph_rad + sph_xoffset, (-1.4-1.6*.35/4.-.1)*sph_rad, |
580 |
> |
.25*sph_rad + sph_xoffset, (-1.4-1.6*.35/4.-.1)*sph_rad, |
581 |
> |
.25*sph_rad + sph_xoffset, -1.3*sph_rad ); |
582 |
|
} |
583 |
|
} |
584 |
|
|
590 |
|
const double leg_width = 2.*.75*(sph_xoffset - sph_rad); |
591 |
|
const double leg_height = 2.*sph_rad; |
592 |
|
const int text_lines = 6; |
593 |
< |
const int text_digits = 7; |
593 |
> |
const int text_digits = 8; |
594 |
|
char fmt[16]; |
595 |
|
int i; |
596 |
|
|
597 |
|
printf("\n# BSDF legend with falsecolor scale\n"); |
598 |
|
printf("\nvoid colorfunc lscale\n10 sca_red(Py) sca_grn(Py) sca_blu(Py)"); |
599 |
|
printf("\n\tbsdf2rad.cal -s %f -t 0 %f 0\n0\n0\n", leg_height, -.5*leg_height); |
600 |
< |
sprintf(fmt, "%%%dg", text_digits-2); |
600 |
> |
sprintf(fmt, "%%.%df", text_digits-3); |
601 |
|
for (i = 0; i < text_lines; i++) { |
602 |
|
char vbuf[16]; |
603 |
|
sprintf(vbuf, fmt, pow(10., (i+.5)/text_lines*(max_log10-min_log10)+min_log10)); |
604 |
|
printf("\nlscale brighttext lscale\n"); |
605 |
|
printf("3 helvet.fnt . %s\n0\n12\n", vbuf); |
606 |
|
printf("\t%f %f 0\n", -.45*leg_width, ((i+.9)/text_lines-.5)*leg_height); |
607 |
< |
printf("\t%f 0 0\n", .9*leg_width/strlen(vbuf)); |
607 |
> |
printf("\t%f 0 0\n", .8*leg_width/strlen(vbuf)); |
608 |
|
printf("\t0 %f 0\n", -.9/text_lines*leg_height); |
609 |
|
printf("\t.01 1 -.1\n"); |
610 |
|
} |
631 |
|
return; |
632 |
|
printf("\nvoid brighttext BSDFname\n3 helvet.fnt . \"%s\"\n0\n12\n", bsdf_name); |
633 |
|
printf("\t%f %f 0\n", -.95*leg_width, -.6*leg_height); |
634 |
< |
printf("\t%f 0 0\n", .4/strlen(bsdf_name)*leg_width); |
634 |
> |
printf("\t%f 0 0\n", 1.8/strlen(bsdf_name)*leg_width); |
635 |
|
printf("\t0 %f 0\n", -.1*leg_height); |
636 |
|
printf("\t.01 1 -.1\n"); |
637 |
|
printf("\nBSDFname alias name_mat vellum\n"); |