| 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 |
|
|
| 385 |
|
for (i = 0; i < NINCIDENT; i++) { |
| 386 |
|
get_ivector(ivec, i); |
| 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 |
|
} |
| 765 |
– |
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); |
