21 |
|
#include "ambient.h" |
22 |
|
#include "random.h" |
23 |
|
|
24 |
< |
#ifdef NEWAMB |
24 |
> |
#ifndef OLDAMB |
25 |
|
|
26 |
|
extern void SDsquare2disk(double ds[2], double seedx, double seedy); |
27 |
|
|
68 |
|
setcolor(ar.rcoef, AVGREFL, AVGREFL, AVGREFL); |
69 |
|
else |
70 |
|
copycolor(ar.rcoef, hp->acoef); |
71 |
< |
if (rayorigin(&ar, AMBIENT, hp->rp, ar.rcoef) < 0) { |
72 |
< |
if (!n) memset(ap, 0, sizeof(AMBSAMP)); |
71 |
> |
if (rayorigin(&ar, AMBIENT, hp->rp, ar.rcoef) < 0) |
72 |
|
return(0); |
74 |
– |
} |
73 |
|
if (ambacc > FTINY) { |
74 |
|
multcolor(ar.rcoef, hp->acoef); |
75 |
|
scalecolor(ar.rcoef, 1./AVGREFL); |
78 |
|
hlist[1] = j; |
79 |
|
hlist[2] = i; |
80 |
|
multisamp(spt, 2, urand(ilhash(hlist,3)+n)); |
81 |
< |
if (!n) { /* avoid border samples for n==0 */ |
81 |
> |
/* avoid coincident samples */ |
82 |
> |
if (!n && (0 < i) & (i < hp->ns-1) && |
83 |
> |
(0 < j) & (j < hp->ns-1)) { |
84 |
|
if ((spt[0] < 0.1) | (spt[0] >= 0.9)) |
85 |
|
spt[0] = 0.1 + 0.8*frandom(); |
86 |
|
if ((spt[1] < 0.1) | (spt[1] >= 0.9)) |
99 |
|
if (ar.rt <= FTINY) |
100 |
|
return(0); /* should never happen */ |
101 |
|
multcolor(ar.rcol, ar.rcoef); /* apply coefficient */ |
102 |
< |
if (!n || ar.rt*ap->d < 1.0) /* new/closer distance? */ |
102 |
> |
if (ar.rt*ap->d < 1.0) /* new/closer distance? */ |
103 |
|
ap->d = 1.0/ar.rt; |
104 |
|
if (!n) { /* record first vertex & value */ |
105 |
|
if (ar.rt > 10.0*thescene.cusize) |
195 |
|
if (e2rem <= FTINY) |
196 |
|
goto done; /* nothing left to do */ |
197 |
|
nss = *ep/e2rem*cnt + frandom(); |
198 |
< |
for (n = 1; n <= nss; n++) |
199 |
< |
cnt -= ambsample(hp, i, j, n); |
198 |
> |
for (n = 1; n <= nss && ambsample(hp,i,j,n); n++) |
199 |
> |
--cnt; |
200 |
|
e2rem -= *ep++; /* update remainder */ |
201 |
|
} |
202 |
|
done: |
229 |
|
hp->rp = r; |
230 |
|
hp->ns = n; |
231 |
|
hp->acol[RED] = hp->acol[GRN] = hp->acol[BLU] = 0.0; |
232 |
+ |
memset(hp->sa, 0, sizeof(AMBSAMP)*n*n); |
233 |
|
hp->sampOK = 0; |
234 |
|
/* assign coefficient */ |
235 |
|
copycolor(hp->acoef, rcol); |
615 |
|
double avg_d = 0; |
616 |
|
uint32 flgs = 0; |
617 |
|
FVECT vec; |
618 |
< |
double d, u, v; |
618 |
> |
double u, v; |
619 |
|
double ang, a1; |
620 |
|
int i, j; |
621 |
|
/* don't bother for a few samples */ |
637 |
|
if ((ap->d <= FTINY) | (ap->d >= max_d)) |
638 |
|
continue; /* too far or too near */ |
639 |
|
VSUB(vec, ap->p, hp->rp->rop); |
640 |
< |
d = DOT(vec, hp->rp->ron); |
641 |
< |
d = 1.0/sqrt(DOT(vec,vec) - d*d); |
642 |
< |
u = DOT(vec, uv[0]) * d; |
642 |
< |
v = DOT(vec, uv[1]) * d; |
643 |
< |
if ((r0*r0*u*u + r1*r1*v*v) * ap->d*ap->d <= 1.0) |
640 |
> |
u = DOT(vec, uv[0]); |
641 |
> |
v = DOT(vec, uv[1]); |
642 |
> |
if ((r0*r0*u*u + r1*r1*v*v) * ap->d*ap->d <= u*u + v*v) |
643 |
|
continue; /* occluder outside ellipse */ |
644 |
|
ang = atan2a(v, u); /* else set direction flags */ |
645 |
|
for (a1 = ang-.5*ang_res; a1 <= ang+.5*ang_res; a1 += ang_step) |