| 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); |
| 513 |
|
if (front_comp) { |
| 514 |
|
printf( |
| 515 |
|
"\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", |
| 516 |
< |
sph_fmat, nsegs, sph_rad, 0.495*PI, sph_xoffset); |
| 516 |
> |
sph_fmat, nsegs, sph_rad, 0.5*PI, sph_xoffset); |
| 517 |
|
printf("\nvoid brighttext front_text\n3 helvet.fnt . FRONT\n0\n"); |
| 518 |
|
printf("12\n\t%f %f 0\n\t%f 0 0\n\t0 %f 0\n\t.01 1 -.1\n", |
| 519 |
|
-.22*sph_rad + sph_xoffset, -1.4*sph_rad, |
| 529 |
|
if (back_comp) { |
| 530 |
|
printf( |
| 531 |
|
"\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", |
| 532 |
< |
sph_bmat, nsegs, sph_rad, 0.495*PI, -sph_xoffset); |
| 532 |
> |
sph_bmat, nsegs, sph_rad, 0.5*PI, -sph_xoffset); |
| 533 |
|
printf("\nvoid brighttext back_text\n3 helvet.fnt . BACK\n0\n"); |
| 534 |
|
printf("12\n\t%f %f 0\n\t%f 0 0\n\t0 %f 0\n\t.01 1 -.1\n", |
| 535 |
|
-.22*sph_rad - sph_xoffset, -1.4*sph_rad, |