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; |
753 |
|
} else { |
754 |
|
FILE *fp[4]; |
755 |
|
if (argc > 5) { |
756 |
< |
fprintf(stderr, "%s: too many input files\n", progname); |
756 |
> |
fprintf(stderr, "%s: more than 4 hemispheres!\n", progname); |
757 |
|
return(1); |
758 |
|
} |
759 |
|
for (n = 1; n < argc; n++) { |
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); |