26 |
|
|
27 |
|
char ourTempDir[TEMPLEN] = ""; /* our temporary directory */ |
28 |
|
|
29 |
< |
const char frpref[] = "frefl"; |
30 |
< |
const char ftpref[] = "ftrans"; |
31 |
< |
const char brpref[] = "brefl"; |
32 |
< |
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_fmat[] = "fBSDFmat"; |
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 |
|
|
283 |
|
refl[2] = 2.*nrm[2]*nrm[2] - 1.; |
284 |
|
|
285 |
|
printf("\n# Mirror arrow\n"); |
286 |
< |
printf("\narrow_mat cylinder inc_dir\n0\n0\n7"); |
286 |
> |
printf("\nshaft_mat cylinder inc_dir\n0\n0\n7"); |
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\n0\n0\n7"); |
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\n0\n0\n8"); |
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], |
315 |
|
int i; |
316 |
|
|
317 |
|
printf("\n# Transmission arrow\n"); |
318 |
< |
printf("\narrow_mat cylinder trans_dir\n0\n0\n7"); |
318 |
> |
printf("\nshaft_mat cylinder trans_dir\n0\n0\n7"); |
319 |
|
printf("\n\t%f %f %f\n\t%f %f %f\n\t%f\n", |
320 |
|
origin[0], origin[1], origin[2], |
321 |
|
origin[0], origin[1], origin[2]-arrow_len, |
322 |
|
arrow_rad); |
323 |
< |
printf("\narrow_mat cone trans_tip\n0\n0\n8"); |
323 |
> |
printf("\ntip_mat cone trans_tip\n0\n0\n8"); |
324 |
|
printf("\n\t%f %f %f\n\t%f %f %f\n\t%f 0\n", |
325 |
|
origin[0], origin[1], origin[2]-arrow_len+.5*tip_len, |
326 |
|
origin[0], origin[1], origin[2]-arrow_len-.5*tip_len, |
379 |
|
printf("\n# Gensurf output corresponding to %d incident directions\n", |
380 |
|
NINCIDENT); |
381 |
|
|
382 |
< |
printf("\nvoid glow arrow_glow\n0\n0\n4 1 0 1 0\n"); |
383 |
< |
printf("\nvoid mixfunc arrow_mat\n4 arrow_glow void 0.25 .\n0\n0\n"); |
382 |
> |
printf("\nvoid glow tip_mat\n0\n0\n4 1 0 1 0\n"); |
383 |
> |
printf("\nvoid mixfunc shaft_mat\n4 tip_mat void 0.25 .\n0\n0\n"); |
384 |
|
|
384 |
– |
if (front_comp & SDsampR) /* front reflection */ |
385 |
|
for (i = 0; i < NINCIDENT; i++) { |
386 |
|
get_ivector(ivec, i); |
387 |
|
nrm[0] = -ivec[0]; nrm[1] = -ivec[1]; nrm[2] = ivec[2]; |
471 |
|
} |
472 |
|
switch (XMLfile[0]) { /* avoid RAYPATH search */ |
473 |
|
case '.': |
474 |
+ |
case '~': |
475 |
|
CASEDIRSEP: |
476 |
|
curdir = ""; |
477 |
|
break; |
751 |
|
strcpy(bsdf_manuf, myBSDF.makr); |
752 |
|
put_matBSDF(argv[1]); |
753 |
|
} else { |
754 |
< |
FILE *fp; |
754 |
> |
FILE *fp[4]; |
755 |
> |
if (argc > 5) { |
756 |
> |
fprintf(stderr, "%s: more than 4 hemispheres!\n", progname); |
757 |
> |
return(1); |
758 |
> |
} |
759 |
|
for (n = 1; n < argc; n++) { |
760 |
< |
fp = fopen(argv[n], "rb"); |
761 |
< |
if (fp == NULL) { |
760 |
> |
fp[n-1] = fopen(argv[n], "rb"); |
761 |
> |
if (fp[n-1] == NULL) { |
762 |
|
fprintf(stderr, "%s: cannot open BSDF interpolant '%s'\n", |
763 |
|
progname, argv[n]); |
764 |
|
return(1); |
765 |
|
} |
766 |
< |
if (getheader(fp, rbf_headline, NULL) < 0) { |
766 |
> |
if (getheader(fp[n-1], rbf_headline, NULL) < 0) { |
767 |
|
fprintf(stderr, "%s: bad BSDF interpolant '%s'\n", |
768 |
|
progname, argv[n]); |
769 |
|
return(1); |
770 |
|
} |
766 |
– |
fclose(fp); |
771 |
|
} |
772 |
|
set_minlog(); |
773 |
|
for (n = 1; n < argc; n++) { |
774 |
< |
fp = fopen(argv[n], "rb"); |
775 |
< |
if (!load_bsdf_rep(fp)) |
774 |
> |
if (fseek(fp[n-1], 0L, SEEK_SET) < 0) { |
775 |
> |
fprintf(stderr, "%s: cannot seek on '%s'\n", |
776 |
> |
progname, argv[n]); |
777 |
|
return(1); |
778 |
< |
fclose(fp); |
778 |
> |
} |
779 |
> |
if (!load_bsdf_rep(fp[n-1])) |
780 |
> |
return(1); |
781 |
> |
fclose(fp[n-1]); |
782 |
|
if (!build_wRBF()) |
783 |
|
return(1); |
784 |
|
} |
789 |
|
put_scale(); |
790 |
|
if (inpXML && myBSDF.mgf) |
791 |
|
convert_mgf(myBSDF.mgf); |
792 |
< |
if (!put_BSDFs()) |
792 |
> |
if (!put_BSDFs()) /* most of the output happens here */ |
793 |
|
return(1); |
794 |
|
cleanup_tmp(); |
795 |
|
return(0); |