| 275 |
|
{ |
| 276 |
|
const double arrow_len = 1.2*bsdf_rad; |
| 277 |
|
const double tip_len = 0.2*bsdf_rad; |
| 278 |
+ |
static int cnt = 1; |
| 279 |
|
FVECT refl; |
| 280 |
|
int i; |
| 281 |
|
|
| 283 |
|
refl[1] = 2.*nrm[2]*nrm[1]; |
| 284 |
|
refl[2] = 2.*nrm[2]*nrm[2] - 1.; |
| 285 |
|
|
| 286 |
< |
printf("\n# Mirror arrow\n"); |
| 287 |
< |
printf("\nshaft_mat cylinder inc_dir\n0\n0\n7"); |
| 286 |
> |
printf("\n# Mirror arrow #%d\n", cnt); |
| 287 |
> |
printf("\nshaft_mat cylinder inc_dir%d\n0\n0\n7", cnt); |
| 288 |
|
printf("\n\t%f %f %f\n\t%f %f %f\n\t%f\n", |
| 289 |
|
origin[0], origin[1], origin[2]+arrow_len, |
| 290 |
|
origin[0], origin[1], origin[2], |
| 291 |
|
arrow_rad); |
| 292 |
< |
printf("\nshaft_mat cylinder mir_dir\n0\n0\n7"); |
| 292 |
> |
printf("\nshaft_mat cylinder mir_dir%d\n0\n0\n7", cnt); |
| 293 |
|
printf("\n\t%f %f %f\n\t%f %f %f\n\t%f\n", |
| 294 |
|
origin[0], origin[1], origin[2], |
| 295 |
|
origin[0] + arrow_len*refl[0], |
| 296 |
|
origin[1] + arrow_len*refl[1], |
| 297 |
|
origin[2] + arrow_len*refl[2], |
| 298 |
|
arrow_rad); |
| 299 |
< |
printf("\ntip_mat cone mir_tip\n0\n0\n8"); |
| 299 |
> |
printf("\ntip_mat cone mir_tip%d\n0\n0\n8", cnt); |
| 300 |
|
printf("\n\t%f %f %f\n\t%f %f %f\n\t%f 0\n", |
| 301 |
|
origin[0] + (arrow_len-.5*tip_len)*refl[0], |
| 302 |
|
origin[1] + (arrow_len-.5*tip_len)*refl[1], |
| 305 |
|
origin[1] + (arrow_len+.5*tip_len)*refl[1], |
| 306 |
|
origin[2] + (arrow_len+.5*tip_len)*refl[2], |
| 307 |
|
2.*arrow_rad); |
| 308 |
+ |
++cnt; |
| 309 |
|
} |
| 310 |
|
|
| 311 |
|
/* Put out transmitted direction arrow for the given incident vector */ |
| 314 |
|
{ |
| 315 |
|
const double arrow_len = 1.2*bsdf_rad; |
| 316 |
|
const double tip_len = 0.2*bsdf_rad; |
| 317 |
+ |
static int cnt = 1; |
| 318 |
|
int i; |
| 319 |
|
|
| 320 |
< |
printf("\n# Transmission arrow\n"); |
| 321 |
< |
printf("\nshaft_mat cylinder trans_dir\n0\n0\n7"); |
| 320 |
> |
printf("\n# Transmission arrow #%d\n", cnt); |
| 321 |
> |
printf("\nshaft_mat cylinder trans_dir%d\n0\n0\n7", cnt); |
| 322 |
|
printf("\n\t%f %f %f\n\t%f %f %f\n\t%f\n", |
| 323 |
|
origin[0], origin[1], origin[2], |
| 324 |
|
origin[0], origin[1], origin[2]-arrow_len, |
| 325 |
|
arrow_rad); |
| 326 |
< |
printf("\ntip_mat cone trans_tip\n0\n0\n8"); |
| 326 |
> |
printf("\ntip_mat cone trans_tip%d\n0\n0\n8", cnt); |
| 327 |
|
printf("\n\t%f %f %f\n\t%f %f %f\n\t%f 0\n", |
| 328 |
|
origin[0], origin[1], origin[2]-arrow_len+.5*tip_len, |
| 329 |
|
origin[0], origin[1], origin[2]-arrow_len-.5*tip_len, |
| 330 |
< |
2.*arrow_rad); |
| 330 |
> |
2.*arrow_rad); |
| 331 |
> |
++cnt; |
| 332 |
|
} |
| 333 |
|
|
| 334 |
|
/* Compute rotation (x,y,z) => (xp,yp,zp) */ |