17 |
|
|
18 |
|
#include "octree.h" |
19 |
|
|
20 |
+ |
#include "otypes.h" |
21 |
+ |
|
22 |
|
#include "random.h" |
23 |
|
|
24 |
|
#define OCTSCALE 0.5 /* ceil((valid rad.)/(cube size)) */ |
34 |
|
extern char *amblist[]; /* ambient include/exclude list */ |
35 |
|
extern int ambincl; /* include == 1, exclude == 0 */ |
36 |
|
|
37 |
< |
OBJECT ambset[128]; /* ambient include/exclude set */ |
37 |
> |
#define MAXASET 511 /* maximum number of elements in ambient set */ |
38 |
> |
OBJECT ambset[MAXASET+1]={0}; /* ambient include/exclude set */ |
39 |
|
|
40 |
|
double maxarad; /* maximum ambient radius */ |
41 |
|
double minarad; /* minimum ambient radius */ |
77 |
|
char *afile; |
78 |
|
{ |
79 |
|
long ftell(); |
77 |
– |
char **amblp; |
78 |
– |
OBJECT obj; |
80 |
|
AMBVAL amb; |
81 |
< |
/* set up ambient set */ |
81 |
< |
ambset[0] = 0; |
82 |
< |
for (amblp = amblist; *amblp != NULL; amblp++) { |
83 |
< |
if ((obj = modifier(*amblp)) == OVOID) { |
84 |
< |
sprintf(errmsg, "unknown %s modifier \"%s\"", |
85 |
< |
ambincl ? "include" : "exclude", *amblp); |
86 |
< |
error(WARNING, errmsg); |
87 |
< |
continue; |
88 |
< |
} |
89 |
< |
if (!inset(ambset, obj)) |
90 |
< |
insertelem(ambset, obj); |
91 |
< |
} |
81 |
> |
|
82 |
|
maxarad = thescene.cusize / 2.0; /* maximum radius */ |
83 |
|
/* minimum radius */ |
84 |
|
minarad = ambres > 0 ? thescene.cusize/ambres : 0.0; |
86 |
|
/* open ambient file */ |
87 |
|
if (afile != NULL) |
88 |
|
if ((ambfp = fopen(afile, "r+")) != NULL) { |
89 |
< |
while (fread(&amb, sizeof(AMBVAL), 1, ambfp) == 1) |
89 |
> |
while (fread((char *)&amb,sizeof(AMBVAL),1,ambfp) == 1) |
90 |
|
avinsert(&amb, &atrunk, thescene.cuorg, |
91 |
|
thescene.cusize); |
92 |
|
/* align */ |
99 |
|
} |
100 |
|
|
101 |
|
|
102 |
+ |
ambnotify(obj) /* record new modifier */ |
103 |
+ |
OBJECT obj; |
104 |
+ |
{ |
105 |
+ |
static int hitlimit = 0; |
106 |
+ |
register OBJREC *o = objptr(obj); |
107 |
+ |
register char **amblp; |
108 |
+ |
|
109 |
+ |
if (hitlimit || !ismodifier(o->otype)) |
110 |
+ |
return; |
111 |
+ |
for (amblp = amblist; *amblp != NULL; amblp++) |
112 |
+ |
if (!strcmp(o->oname, *amblp)) { |
113 |
+ |
if (ambset[0] >= MAXASET) { |
114 |
+ |
error(WARNING, "too many modifiers in ambient list"); |
115 |
+ |
hitlimit++; |
116 |
+ |
return; /* should this be fatal? */ |
117 |
+ |
} |
118 |
+ |
insertelem(ambset, obj); |
119 |
+ |
return; |
120 |
+ |
} |
121 |
+ |
} |
122 |
+ |
|
123 |
+ |
|
124 |
|
ambient(acol, r) /* compute ambient component for ray */ |
125 |
|
COLOR acol; |
126 |
|
register RAY *r; |
215 |
|
* Jittering final test reduces image artifacts. |
216 |
|
*/ |
217 |
|
wt = sqrt(e1) + sqrt(e2); |
218 |
< |
if (wt > ambacc*(0.9 + 0.2*frandom())) |
218 |
> |
wt *= .9 + .2*frandom(); |
219 |
> |
if (wt > ambacc) |
220 |
|
continue; |
221 |
|
if (wt <= 1e-3) |
222 |
|
wt = 1e3; |
279 |
|
extern int ambcmp(); |
280 |
|
extern double sin(), cos(), sqrt(); |
281 |
|
double phi, xd, yd, zd; |
282 |
+ |
double b, b2; |
283 |
|
register AMBSAMP *div; |
284 |
|
AMBSAMP dnew; |
285 |
|
RAY ar; |
334 |
|
div[ne].n = 0; |
335 |
|
div[ne].t = i; div[ne].p = j; |
336 |
|
/* sum errors */ |
337 |
< |
xd = bright(ar.rcol); |
337 |
> |
b = bright(ar.rcol); |
338 |
|
if (i > 0) { /* from above */ |
339 |
< |
yd = bright(div[ne-np].v) - xd; |
340 |
< |
yd *= yd * 0.25; |
341 |
< |
div[ne].k += yd; |
339 |
> |
b2 = bright(div[ne-np].v) - b; |
340 |
> |
b2 *= b2 * 0.25; |
341 |
> |
div[ne].k += b2; |
342 |
|
div[ne].n++; |
343 |
< |
div[ne-np].k += yd; |
343 |
> |
div[ne-np].k += b2; |
344 |
|
div[ne-np].n++; |
345 |
|
} |
346 |
|
if (j > 0) { /* from behind */ |
347 |
< |
yd = bright(div[ne-1].v) - xd; |
348 |
< |
yd *= yd * 0.25; |
349 |
< |
div[ne].k += yd; |
347 |
> |
b2 = bright(div[ne-1].v) - b; |
348 |
> |
b2 *= b2 * 0.25; |
349 |
> |
div[ne].k += b2; |
350 |
|
div[ne].n++; |
351 |
< |
div[ne-1].k += yd; |
351 |
> |
div[ne-1].k += b2; |
352 |
|
div[ne-1].n++; |
353 |
|
} |
354 |
|
if (j == np-1) { /* around */ |
355 |
< |
yd = bright(div[ne-(np-1)].v) - xd; |
356 |
< |
yd *= yd * 0.25; |
357 |
< |
div[ne].k += yd; |
355 |
> |
b2 = bright(div[ne-(np-1)].v) - b; |
356 |
> |
b2 *= b2 * 0.25; |
357 |
> |
div[ne].k += b2; |
358 |
|
div[ne].n++; |
359 |
< |
div[ne-(np-1)].k += yd; |
359 |
> |
div[ne-(np-1)].k += b2; |
360 |
|
div[ne-(np-1)].n++; |
361 |
|
} |
362 |
|
ne++; |
393 |
|
addcolor(dnew.v, ar.rcol); |
394 |
|
dnew.n = div[0].n + 1; |
395 |
|
dnew.t = div[0].t; dnew.p = div[0].p; |
396 |
< |
yd = bright(dnew.v)/dnew.n - bright(ar.rcol); |
397 |
< |
yd = yd*yd + div[0].k*(div[0].n*div[0].n); |
398 |
< |
dnew.k = yd/(dnew.n*dnew.n); |
396 |
> |
b2 = bright(dnew.v)/dnew.n - bright(ar.rcol); |
397 |
> |
b2 = b2*b2 + div[0].k*(div[0].n*div[0].n); |
398 |
> |
dnew.k = b2/(dnew.n*dnew.n); |
399 |
|
/* reinsert */ |
400 |
|
for (k = 0; k < ne-1 && dnew.k < div[k+1].k; k++) |
401 |
< |
bcopy(&div[k+1], &div[k], sizeof(AMBSAMP)); |
402 |
< |
bcopy(&dnew, &div[k], sizeof(AMBSAMP)); |
401 |
> |
copystruct(&div[k], &div[k+1]); |
402 |
> |
copystruct(&div[k], &dnew); |
403 |
|
|
404 |
|
if (ne >= i) { /* extract darkest division */ |
405 |
|
ne--; |
444 |
|
#endif |
445 |
|
if (ambfp == NULL) |
446 |
|
return; |
447 |
< |
if (fwrite(av, sizeof(AMBVAL), 1, ambfp) != 1) |
447 |
> |
if (fwrite((char *)av, sizeof(AMBVAL), 1, ambfp) != 1) |
448 |
|
goto writerr; |
449 |
|
#ifdef AMBFLUSH |
450 |
|
if (++nunflshed >= AMBFLUSH) { |
473 |
|
|
474 |
|
if ((av = newambval()) == NULL) |
475 |
|
goto memerr; |
476 |
< |
bcopy(aval, av, sizeof(AMBVAL)); |
476 |
> |
copystruct(av, aval); |
477 |
|
VCOPY(ck0, c0); |
478 |
|
while (s*(OCTSCALE/2) > av->rad*ambacc) { |
479 |
|
if (at->kid == NULL) |