| 128 |
|
dv[0] = cos(M_PI/180.*inpfile[i].phi)*dv[2]; |
| 129 |
|
dv[1] = sin(M_PI/180.*inpfile[i].phi)*dv[2]; |
| 130 |
|
dv[2] = sqrt(1. - dv[2]*dv[2]); |
| 131 |
< |
pos_from_vec(inpfile[i].igp, dv); |
| 131 |
> |
if (inpfile[i].theta <= FTINY) |
| 132 |
> |
inpfile[i].igp[0] = inpfile[i].igp[1] = grid_res/2 - 1; |
| 133 |
> |
else |
| 134 |
> |
pos_from_vec(inpfile[i].igp, dv); |
| 135 |
|
return(1); |
| 136 |
|
} |
| 137 |
|
|
| 269 |
|
} |
| 270 |
|
#ifdef DEBUG |
| 271 |
|
fprintf(stderr, "Minimum BSDF = %.4f\n", bsdf_min); |
| 272 |
+ |
fprintf(stderr, "Integrated hemispherical value = %.4f\n", dsf_list->vtotal); |
| 273 |
|
#endif |
| 274 |
|
min_log = log(bsdf_min*.5 + 1e-5); |
| 275 |
|
#if 1 /* produce spheres at meas. */ |