27 |
|
short nt, np; /* number of theta and phi directions */ |
28 |
|
} AMBHEMI; /* ambient sample hemisphere */ |
29 |
|
|
30 |
– |
extern double sin(), cos(), sqrt(); |
30 |
|
|
32 |
– |
|
31 |
|
static int |
32 |
|
ambcmp(d1, d2) /* decreasing order */ |
33 |
|
AMBSAMP *d1, *d2; |
146 |
|
arad += dp->r; |
147 |
|
} |
148 |
|
} |
149 |
< |
if (ns > 0) { /* perform super-sampling */ |
149 |
> |
if (ns > 0 && arad > FTINY && ndivs/arad < minarad) |
150 |
> |
ns = 0; /* close enough */ |
151 |
> |
else if (ns > 0) { /* else perform super-sampling */ |
152 |
|
comperrs(div, &hemi); /* compute errors */ |
153 |
|
qsort(div, ndivs, sizeof(AMBSAMP), ambcmp); /* sort divs */ |
154 |
|
/* super-sample */ |
207 |
|
scalecolor(acol, b); |
208 |
|
if (arad <= FTINY) |
209 |
|
arad = maxarad; |
210 |
< |
else { |
210 |
> |
else |
211 |
|
arad = (ndivs+ns)/arad; |
212 |
– |
if (arad > maxarad) |
213 |
– |
arad = maxarad; |
214 |
– |
} |
212 |
|
if (pg != NULL) { /* reduce radius if gradient large */ |
213 |
|
d = DOT(pg,pg); |
214 |
|
if (d*arad*arad > 1.0) |
222 |
|
pg[i] *= d; |
223 |
|
} |
224 |
|
} |
225 |
< |
return(arad/sqrt(wt)); |
225 |
> |
if ((arad /= sqrt(wt)) > maxarad) |
226 |
> |
arad = maxarad; |
227 |
> |
return(arad); |
228 |
|
oopsy: |
229 |
|
if (div != NULL) |
230 |
|
free((char *)div); |