| 36 |
|
int nglarangs = 0; |
| 37 |
|
double maxtheta; /* maximum angle (in radians) */ |
| 38 |
|
int hsize; /* horizontal size */ |
| 39 |
– |
int hlim; /* central limit of horizontal */ |
| 39 |
|
|
| 40 |
|
struct illum *indirect; /* array of indirect illuminances */ |
| 41 |
|
|
| 190 |
|
fputs(VIEWSTR, stdout); |
| 191 |
|
fprintview(&ourview, stdout); |
| 192 |
|
printf("\n"); |
| 193 |
< |
fputformat("ASCII", stdout); |
| 193 |
> |
fputformat("ascii", stdout); |
| 194 |
|
printf("\n"); |
| 195 |
|
printsources(); /* print glare sources */ |
| 196 |
|
printillum(); /* print illuminances */ |
| 226 |
|
maxtheta = (PI/180.)*glarang[nglarangs-1]; |
| 227 |
|
else |
| 228 |
|
maxtheta = 0.0; |
| 229 |
< |
hlim = sampdens*maxtheta; |
| 231 |
< |
hsize = hlim + sampdens - 1; |
| 229 |
> |
hsize = hlim(0) + sampdens - 1; |
| 230 |
|
if (hsize > (int)(PI*sampdens)) |
| 231 |
|
hsize = PI*sampdens; |
| 232 |
|
indirect = (struct illum *)calloc(nglardirs, sizeof(struct illum)); |
| 292 |
|
} |
| 293 |
|
|
| 294 |
|
|
| 295 |
< |
compdir(vd, x, y, se) /* compute direction for x,y */ |
| 295 |
> |
compdir(vd, x, y) /* compute direction for x,y */ |
| 296 |
|
FVECT vd; |
| 297 |
|
int x, y; |
| 300 |
– |
SPANERR *se; |
| 298 |
|
{ |
| 299 |
+ |
int hl; |
| 300 |
|
FVECT org; /* dummy variable */ |
| 301 |
|
|
| 302 |
< |
if (x <= -hlim) /* left region */ |
| 302 |
> |
hl = hlim(y); |
| 303 |
> |
if (x <= -hl) { /* left region */ |
| 304 |
> |
if (x <= -hl-sampdens) |
| 305 |
> |
return(-1); |
| 306 |
|
return(viewray(org, vd, &leftview, |
| 307 |
< |
(double)(x+hlim)/(2*sampdens)+.5, |
| 307 |
> |
(double)(x+hl)/(2*sampdens)+.5, |
| 308 |
|
(double)y/(2*sampdens)+.5)); |
| 309 |
< |
if (x >= hlim) /* right region */ |
| 309 |
> |
} |
| 310 |
> |
if (x >= hl) { /* right region */ |
| 311 |
> |
if (x >= hl+sampdens) |
| 312 |
> |
return(-1); |
| 313 |
|
return(viewray(org, vd, &rightview, |
| 314 |
< |
(double)(x-hlim)/(2*sampdens)+.5, |
| 314 |
> |
(double)(x-hl)/(2*sampdens)+.5, |
| 315 |
|
(double)y/(2*sampdens)+.5)); |
| 312 |
– |
/* central region */ |
| 313 |
– |
if (se != NULL) { /* avoid over-counting of poles */ |
| 314 |
– |
se->err += se->prob; |
| 315 |
– |
if (se->err <= 0.5) |
| 316 |
– |
return(-1); |
| 317 |
– |
se->err -= 1.0; |
| 316 |
|
} |
| 317 |
+ |
/* central region */ |
| 318 |
|
if (viewray(org, vd, &ourview, .5, (double)y/(2*sampdens)+.5) < 0) |
| 319 |
|
return(-1); |
| 320 |
< |
spinvector(vd, vd, ourview.vup, h_theta(x)); |
| 320 |
> |
spinvector(vd, vd, ourview.vup, h_theta(x,y)); |
| 321 |
|
return(0); |
| 322 |
|
} |
| 323 |
|
|
| 324 |
|
|
| 325 |
< |
setspanerr(se, y) /* initialize span error at y */ |
| 326 |
< |
register SPANERR *se; |
| 327 |
< |
int y; |
| 325 |
> |
double |
| 326 |
> |
pixsize(x, y) /* return the solid angle of pixel at (x,y) */ |
| 327 |
> |
int x, y; |
| 328 |
|
{ |
| 329 |
< |
se->err = 0.0; |
| 330 |
< |
se->prob = sqrt(1.0 - (double)((long)y*y)/((long)vsize*vsize)); |
| 329 |
> |
register int hl, xo; |
| 330 |
> |
double disc; |
| 331 |
> |
|
| 332 |
> |
hl = hlim(y); |
| 333 |
> |
if (x < -hl) |
| 334 |
> |
xo = x+hl; |
| 335 |
> |
else if (x > hl) |
| 336 |
> |
xo = x-hl; |
| 337 |
> |
else |
| 338 |
> |
xo = 0; |
| 339 |
> |
disc = 1. - (double)(xo*xo + y*y)/(sampdens*sampdens); |
| 340 |
> |
if (disc <= FTINY) |
| 341 |
> |
return(0.); |
| 342 |
> |
return(1./(sampdens*sampdens*sqrt(disc))); |
| 343 |
|
} |
| 344 |
|
|
| 345 |
|
|
