59 |
|
dir[0] = dir[2] * cos((M_PI/180.)*atof(argv[2*n+3])); |
60 |
|
dir[1] = dir[2] * sin((M_PI/180.)*atof(argv[2*n+3])); |
61 |
|
dir[2] = input_orient * sqrt(1. - dir[2]*dir[2]); |
62 |
+ |
#ifdef DEBUG |
63 |
|
fprintf(stderr, "Computing DSF for incident direction (%.1f,%.1f)\n", |
64 |
|
get_theta180(dir), get_phi360(dir)); |
65 |
+ |
#endif |
66 |
|
rbf = advect_rbf(dir, 15000); |
67 |
+ |
#ifdef DEBUG |
68 |
|
if (rbf == NULL) |
69 |
|
fputs("NULL RBF\n", stderr); |
70 |
|
else |
71 |
|
fprintf(stderr, "Hemispherical reflectance: %.3f\n", rbf->vtotal); |
72 |
+ |
#endif |
73 |
|
printf("void trans tmat\n0\n0\n7 %f %f %f .04 .04 .9 1\n", |
74 |
|
colarr[n][0], colarr[n][1], colarr[n][2]); |
75 |
|
if (showPeaks && rbf != NULL) { |