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. */ |