| 40 |
|
double maxarad; /* maximum ambient radius */ |
| 41 |
|
double minarad; /* minimum ambient radius */ |
| 42 |
|
|
| 43 |
+ |
/* |
| 44 |
+ |
* Since we've defined our vectors as float below to save space, |
| 45 |
+ |
* watch out for changes in the definitions of VCOPY() and DOT(). |
| 46 |
+ |
*/ |
| 47 |
|
typedef struct ambval { |
| 48 |
< |
FVECT pos; /* position in space */ |
| 49 |
< |
FVECT dir; /* normal direction */ |
| 48 |
> |
float pos[3]; /* position in space */ |
| 49 |
> |
float dir[3]; /* normal direction */ |
| 50 |
|
int lvl; /* recursion level of parent ray */ |
| 51 |
|
float weight; /* weight of parent ray */ |
| 52 |
|
COLOR val; /* computed ambient value */ |
| 213 |
|
for (j = 0; j < 3; j++) |
| 214 |
|
d += (r->rop[j] - av->pos[j]) * |
| 215 |
|
(av->dir[j] + r->ron[j]); |
| 216 |
< |
if (d < -minarad) |
| 216 |
> |
if (d*0.5 < -minarad*ambacc) |
| 217 |
|
continue; |
| 218 |
|
/* |
| 219 |
|
* Jittering final test reduces image artifacts. |
| 282 |
|
{ |
| 283 |
|
extern int ambcmp(); |
| 284 |
|
extern double sin(), cos(), sqrt(); |
| 285 |
+ |
int hlist[4]; |
| 286 |
|
double phi, xd, yd, zd; |
| 287 |
|
double b, b2; |
| 288 |
|
register AMBSAMP *div; |
| 307 |
|
div = (AMBSAMP *)malloc(ndivs*sizeof(AMBSAMP)); |
| 308 |
|
if (div == NULL) |
| 309 |
|
error(SYSTEM, "out of memory in doambient"); |
| 310 |
< |
} |
| 310 |
> |
} else |
| 311 |
> |
div = NULL; |
| 312 |
|
/* make axes */ |
| 313 |
|
uy[0] = uy[1] = uy[2] = 0.0; |
| 314 |
|
for (k = 0; k < 3; k++) |
| 318 |
|
fcross(ux, r->ron, uy); |
| 319 |
|
normalize(ux); |
| 320 |
|
fcross(uy, ux, r->ron); |
| 321 |
+ |
/* set up urand */ |
| 322 |
+ |
hlist[0] = r->rno; |
| 323 |
|
/* sample divisions */ |
| 324 |
|
arad = 0.0; |
| 325 |
|
ne = 0; |
| 326 |
< |
for (i = 0; i < nt; i++) |
| 326 |
> |
for (i = 0; i < nt; i++) { |
| 327 |
> |
hlist[1] = i; |
| 328 |
|
for (j = 0; j < np; j++) { |
| 329 |
|
rayorigin(&ar, r, AMBIENT, 0.5); /* pretested */ |
| 330 |
< |
zd = sqrt((i+frandom())/nt); |
| 331 |
< |
phi = 2.0*PI * (j+frandom())/np; |
| 330 |
> |
hlist[2] = j; |
| 331 |
> |
hlist[3] = 0; |
| 332 |
> |
zd = sqrt((i+urand(ilhash(hlist,4)))/nt); |
| 333 |
> |
hlist[3] = 1; |
| 334 |
> |
phi = 2.0*PI * (j+urand(ilhash(hlist,4)))/np; |
| 335 |
|
xd = cos(phi) * zd; |
| 336 |
|
yd = sin(phi) * zd; |
| 337 |
|
zd = sqrt(1.0 - zd*zd); |
| 338 |
|
for (k = 0; k < 3; k++) |
| 339 |
|
ar.rdir[k] = xd*ux[k]+yd*uy[k]+zd*r->ron[k]; |
| 340 |
+ |
dimlist[ndims++] = i*np + j + 38813; |
| 341 |
|
rayvalue(&ar); |
| 342 |
+ |
ndims--; |
| 343 |
|
if (ar.rot < FHUGE) |
| 344 |
|
arad += 1.0 / ar.rot; |
| 345 |
< |
if (ns > 0) { /* save division */ |
| 345 |
> |
if (div != NULL) { /* save division */ |
| 346 |
|
div[ne].k = 0.0; |
| 347 |
|
copycolor(div[ne].v, ar.rcol); |
| 348 |
|
div[ne].n = 0; |
| 377 |
|
} else |
| 378 |
|
addcolor(acol, ar.rcol); |
| 379 |
|
} |
| 380 |
+ |
} |
| 381 |
|
for (k = 0; k < ne; k++) { /* compute errors */ |
| 382 |
|
if (div[k].n > 1) |
| 383 |
|
div[k].k /= div[k].n; |
| 393 |
|
/* super-sample */ |
| 394 |
|
for (i = ns; i > 0; i--) { |
| 395 |
|
rayorigin(&ar, r, AMBIENT, 0.5); /* pretested */ |
| 396 |
< |
zd = sqrt((div[0].t+frandom())/nt); |
| 397 |
< |
phi = 2.0*PI * (div[0].p+frandom())/np; |
| 396 |
> |
hlist[1] = div[0].t; |
| 397 |
> |
hlist[2] = div[0].p; |
| 398 |
> |
hlist[3] = 0; |
| 399 |
> |
zd = sqrt((div[0].t+urand(ilhash(hlist,4)+div[0].n))/nt); |
| 400 |
> |
hlist[3] = 1; |
| 401 |
> |
phi = 2.0*PI * (div[0].p+urand(ilhash(hlist,4)+div[0].n))/np; |
| 402 |
|
xd = cos(phi) * zd; |
| 403 |
|
yd = sin(phi) * zd; |
| 404 |
|
zd = sqrt(1.0 - zd*zd); |
| 405 |
|
for (k = 0; k < 3; k++) |
| 406 |
|
ar.rdir[k] = xd*ux[k]+yd*uy[k]+zd*r->ron[k]; |
| 407 |
+ |
dimlist[ndims++] = div[0].t*np + div[0].p + 38813; |
| 408 |
|
rayvalue(&ar); |
| 409 |
+ |
ndims--; |
| 410 |
+ |
rayvalue(&ar); |
| 411 |
|
if (ar.rot < FHUGE) |
| 412 |
|
arad += 1.0 / ar.rot; |
| 413 |
|
/* recompute error */ |
| 425 |
|
|
| 426 |
|
if (ne >= i) { /* extract darkest division */ |
| 427 |
|
ne--; |
| 428 |
< |
if (div[ne].n > 1) |
| 429 |
< |
scalecolor(div[ne].v, 1.0/div[ne].n); |
| 428 |
> |
if (div[ne].n > 1) { |
| 429 |
> |
b = 1.0/div[ne].n; |
| 430 |
> |
scalecolor(div[ne].v, b); |
| 431 |
> |
div[ne].n = 1; |
| 432 |
> |
} |
| 433 |
|
addcolor(acol, div[ne].v); |
| 434 |
|
} |
| 435 |
|
} |
| 442 |
|
arad = maxarad; |
| 443 |
|
else if (arad < minarad) |
| 444 |
|
arad = minarad; |
| 445 |
< |
if (ns > 0) |
| 445 |
> |
if (div != NULL) |
| 446 |
|
free((char *)div); |
| 447 |
|
return(arad); |
| 448 |
|
} |
