| 40 |
|
|
| 41 |
|
struct illum *indirect; /* array of indirect illuminances */ |
| 42 |
|
|
| 43 |
+ |
long npixinvw; /* number of pixels in view */ |
| 44 |
+ |
long npixmiss; /* number of pixels missed */ |
| 45 |
|
|
| 46 |
+ |
|
| 47 |
|
main(argc, argv) |
| 48 |
|
int argc; |
| 49 |
|
char *argv[]; |
| 207 |
|
/* set direction vectors */ |
| 208 |
|
for (i = 0; glarang[i] != AEND; i++) |
| 209 |
|
; |
| 210 |
< |
if (i > 0 && glarang[0] <= 0 || glarang[i-1] >= 180) { |
| 210 |
> |
if (i > 0 && (glarang[0] <= 0 || glarang[i-1] >= 180)) { |
| 211 |
|
fprintf(stderr, "%s: glare angles must be between 1 and 179\n", |
| 212 |
|
progname); |
| 213 |
|
exit(1); |
| 226 |
|
indirect = (struct illum *)calloc(nglardirs, sizeof(struct illum)); |
| 227 |
|
if (indirect == NULL) |
| 228 |
|
memerr("indirect illuminances"); |
| 229 |
+ |
npixinvw = npixmiss = 0L; |
| 230 |
|
copystruct(&leftview, &ourview); |
| 231 |
|
copystruct(&rightview, &ourview); |
| 232 |
|
spinvector(leftview.vdir, ourview.vdir, ourview.vup, maxtheta); |
| 235 |
|
setview(&rightview); |
| 236 |
|
indirect[nglarangs].lcos = |
| 237 |
|
indirect[nglarangs].rcos = cos(maxtheta); |
| 238 |
< |
indirect[nglarangs].lsin = |
| 239 |
< |
-(indirect[nglarangs].rsin = sin(maxtheta)); |
| 238 |
> |
indirect[nglarangs].rsin = |
| 239 |
> |
-(indirect[nglarangs].lsin = sin(maxtheta)); |
| 240 |
|
indirect[nglarangs].theta = 0.0; |
| 241 |
|
for (i = 0; i < nglarangs; i++) { |
| 242 |
|
d = (glarang[nglarangs-1] - glarang[i])*(PI/180.); |
| 243 |
|
indirect[nglarangs-i-1].lcos = |
| 244 |
|
indirect[nglarangs+i+1].rcos = cos(d); |
| 245 |
< |
indirect[nglarangs-i-1].lsin = |
| 246 |
< |
-(indirect[nglarangs+i+1].rsin = sin(d)); |
| 245 |
> |
indirect[nglarangs+i+1].rsin = |
| 246 |
> |
-(indirect[nglarangs-i-1].lsin = sin(d)); |
| 247 |
|
d = (glarang[nglarangs-1] + glarang[i])*(PI/180.); |
| 248 |
|
indirect[nglarangs-i-1].rcos = |
| 249 |
|
indirect[nglarangs+i+1].lcos = cos(d); |
| 250 |
< |
indirect[nglarangs+i+1].lsin = |
| 251 |
< |
-(indirect[nglarangs-i-1].rsin = sin(d)); |
| 252 |
< |
indirect[nglarangs-i-1].theta = -(PI/180.)*glarang[i]; |
| 253 |
< |
indirect[nglarangs+i+1].theta = (PI/180.)*glarang[i]; |
| 250 |
> |
indirect[nglarangs-i-1].rsin = |
| 251 |
> |
-(indirect[nglarangs+i+1].lsin = sin(d)); |
| 252 |
> |
indirect[nglarangs-i-1].theta = (PI/180.)*glarang[i]; |
| 253 |
> |
indirect[nglarangs+i+1].theta = -(PI/180.)*glarang[i]; |
| 254 |
|
} |
| 255 |
|
/* open picture */ |
| 256 |
|
if (picture != NULL) { |
| 280 |
|
close_pict(); |
| 281 |
|
if (octree != NULL) |
| 282 |
|
done_rtrace(); |
| 283 |
+ |
if (npixinvw < 100*npixmiss) |
| 284 |
+ |
fprintf(stderr, "%s: warning -- missing %ld%% of samples\n", |
| 285 |
+ |
progname, 100L*npixmiss/npixinvw); |
| 286 |
|
} |
| 287 |
|
|
| 288 |
|
|
| 290 |
|
FVECT vd; |
| 291 |
|
int x, y; |
| 292 |
|
{ |
| 293 |
+ |
static int cury = 10000; |
| 294 |
+ |
static double err, cmpval; |
| 295 |
+ |
long t; |
| 296 |
|
FVECT org; /* dummy variable */ |
| 297 |
|
|
| 298 |
|
if (x <= -hlim) /* left region */ |
| 299 |
|
return(viewray(org, vd, &leftview, |
| 300 |
< |
(x+hlim)/(2.*sampdens)+.5, |
| 301 |
< |
y/(2.*sampdens)+.5)); |
| 300 |
> |
(double)(x+hlim)/(2*sampdens)+.5, |
| 301 |
> |
(double)y/(2*sampdens)+.5)); |
| 302 |
|
if (x >= hlim) /* right region */ |
| 303 |
|
return(viewray(org, vd, &rightview, |
| 304 |
< |
(x-hlim)/(2.*sampdens)+.5, |
| 305 |
< |
y/(2.*sampdens)+.5)); |
| 306 |
< |
/* central region */ |
| 307 |
< |
if (viewray(org, vd, &ourview, .5, y/(2.*sampdens)+.5) < 0) |
| 304 |
> |
(double)(x-hlim)/(2*sampdens)+.5, |
| 305 |
> |
(double)y/(2*sampdens)+.5)); |
| 306 |
> |
/* central region */ |
| 307 |
> |
/* avoid over-counting of poles */ |
| 308 |
> |
if (cury != y) { |
| 309 |
> |
err = 0.0; |
| 310 |
> |
cmpval = sqrt(1.0 - (double)((long)y*y)/((long)vsize*vsize)); |
| 311 |
> |
cury = y; |
| 312 |
> |
} |
| 313 |
> |
err += cmpval; |
| 314 |
> |
if (err <= 0.5) |
| 315 |
|
return(-1); |
| 316 |
+ |
err -= 1.0; |
| 317 |
+ |
if (viewray(org, vd, &ourview, .5, (double)y/(2*sampdens)+.5) < 0) |
| 318 |
+ |
return(-1); |
| 319 |
|
spinvector(vd, vd, ourview.vup, h_theta(x)); |
| 320 |
|
return(0); |
| 321 |
|
} |
| 348 |
|
|
| 349 |
|
printf("BEGIN indirect illuminance\n"); |
| 350 |
|
for (i = 0; i < nglardirs; i++) |
| 351 |
< |
printf("\t%.0f\t%f\n", (180.0/PI)*indirect[i].theta, |
| 352 |
< |
PI * indirect[i].sum / (double)indirect[i].n); |
| 351 |
> |
if (indirect[i].n > FTINY) |
| 352 |
> |
printf("\t%.0f\t%f\n", (180.0/PI)*indirect[i].theta, |
| 353 |
> |
PI * indirect[i].sum / indirect[i].n); |
| 354 |
|
printf("END indirect illuminance\n"); |
| 355 |
|
} |
