56 |
|
return; |
57 |
|
} |
58 |
|
/* set up sampling */ |
59 |
< |
n = PI * il->sampdens; |
60 |
< |
nalt = sqrt(n/PI) + .5; |
61 |
< |
nazi = PI*nalt + .5; |
59 |
> |
if (il->sampdens <= 0) |
60 |
> |
nalt = nazi = 1; |
61 |
> |
else { |
62 |
> |
n = PI * il->sampdens; |
63 |
> |
nalt = sqrt(n/PI) + .5; |
64 |
> |
nazi = PI*nalt + .5; |
65 |
> |
} |
66 |
|
n = nalt*nazi; |
67 |
|
distarr = (float *)calloc(n, 3*sizeof(float)); |
68 |
|
if (distarr == NULL) |
86 |
|
for (i = 0; i < il->nsamps; i++) { |
87 |
|
/* random direction */ |
88 |
|
h = ilhash(dim, 3) + i; |
89 |
< |
peano(sp, 2, urand(h), .02); |
89 |
> |
multisamp(sp, 2, urand(h)); |
90 |
|
r1 = (dim[1] + sp[0])/nalt; |
91 |
< |
r2 = (dim[2] + sp[1])/nazi; |
91 |
> |
r2 = (dim[2] + sp[1] - .5)/nazi; |
92 |
|
flatdir(dn, r1, r2); |
93 |
|
for (j = 0; j < 3; j++) |
94 |
|
dir[j] = -dn[0]*u[j] - dn[1]*v[j] - dn[2]*fa->norm[j]; |
95 |
|
/* random location */ |
96 |
|
do { |
97 |
< |
peano(sp, 2, urand(h+nmisses), .01); |
97 |
> |
multisamp(sp, 2, urand(h+4862+nmisses)); |
98 |
|
r1 = ur[0] + (ur[1]-ur[0]) * sp[0]; |
99 |
|
r2 = vr[0] + (vr[1]-vr[0]) * sp[1]; |
100 |
|
for (j = 0; j < 3; j++) |
115 |
|
raysamp(distarr+3*(dim[1]*nazi+dim[2]), org, dir, rt); |
116 |
|
} |
117 |
|
rayflush(rt); |
118 |
< |
/* write out the face w/ distribution */ |
119 |
< |
flatout(il, distarr, nalt, nazi, u, v, fa->norm); |
120 |
< |
illumout(il, ob); |
118 |
> |
/* write out the face and its distribution */ |
119 |
> |
if (average(il, distarr, nalt*nazi)) { |
120 |
> |
if (il->sampdens > 0) |
121 |
> |
flatout(il, distarr, nalt, nazi, u, v, fa->norm); |
122 |
> |
illumout(il, ob); |
123 |
> |
} else if (!(il->flags & IL_LIGHT)) |
124 |
> |
printobj(il->altmat, ob); |
125 |
|
/* clean up */ |
126 |
|
freeface(ob); |
127 |
|
free((char *)distarr); |
146 |
|
if (ob->oargs.nfargs != 4) |
147 |
|
objerror(ob, USER, "bad # of arguments"); |
148 |
|
/* set up sampling */ |
149 |
< |
n = 4.*PI * il->sampdens; |
150 |
< |
nalt = sqrt(n/PI) + .5; |
151 |
< |
nazi = PI*nalt + .5; |
149 |
> |
if (il->sampdens <= 0) |
150 |
> |
nalt = nazi = 1; |
151 |
> |
else { |
152 |
> |
n = 4.*PI * il->sampdens; |
153 |
> |
nalt = sqrt(n/PI) + .5; |
154 |
> |
nazi = PI*nalt + .5; |
155 |
> |
} |
156 |
|
n = nalt*nazi; |
157 |
|
distarr = (float *)calloc(n, 3*sizeof(float)); |
158 |
|
if (distarr == NULL) |
163 |
|
for (dim[2] = 0; dim[2] < nazi; dim[2]++) |
164 |
|
for (i = 0; i < il->nsamps; i++) { |
165 |
|
/* next sample point */ |
166 |
< |
peano(sp, 4, urand(ilhash(dim,3)+i), .02); |
166 |
> |
multisamp(sp, 4, urand(ilhash(dim,3)+i)); |
167 |
|
/* random direction */ |
168 |
|
r1 = (dim[1] + sp[0])/nalt; |
169 |
< |
r2 = (dim[2] + sp[1])/nazi; |
169 |
> |
r2 = (dim[2] + sp[1] - .5)/nazi; |
170 |
|
rounddir(dir, r1, r2); |
171 |
|
/* random location */ |
172 |
|
mkaxes(u, v, dir); /* yuck! */ |
184 |
|
raysamp(distarr+3*(dim[1]*nazi+dim[2]), org, dir, rt); |
185 |
|
} |
186 |
|
rayflush(rt); |
187 |
< |
/* write out the sphere w/ distribution */ |
188 |
< |
roundout(il, distarr, nalt, nazi); |
189 |
< |
illumout(il, ob); |
187 |
> |
/* write out the sphere and its distribution */ |
188 |
> |
if (average(il, distarr, nalt*nazi)) { |
189 |
> |
if (il->sampdens > 0) |
190 |
> |
roundout(il, distarr, nalt, nazi); |
191 |
> |
else |
192 |
> |
objerror(ob, WARNING, "diffuse distribution"); |
193 |
> |
illumout(il, ob); |
194 |
> |
} else if (!(il->flags & IL_LIGHT)) |
195 |
> |
printobj(il->altmat, ob); |
196 |
|
/* clean up */ |
197 |
|
free((char *)distarr); |
198 |
|
} |
215 |
|
/* get/check arguments */ |
216 |
|
co = getcone(ob, 0); |
217 |
|
/* set up sampling */ |
218 |
< |
n = PI * il->sampdens; |
219 |
< |
nalt = sqrt(n/PI) + .5; |
220 |
< |
nazi = PI*nalt + .5; |
218 |
> |
if (il->sampdens <= 0) |
219 |
> |
nalt = nazi = 1; |
220 |
> |
else { |
221 |
> |
n = PI * il->sampdens; |
222 |
> |
nalt = sqrt(n/PI) + .5; |
223 |
> |
nazi = PI*nalt + .5; |
224 |
> |
} |
225 |
|
n = nalt*nazi; |
226 |
|
distarr = (float *)calloc(n, 3*sizeof(float)); |
227 |
|
if (distarr == NULL) |
233 |
|
for (dim[2] = 0; dim[2] < nazi; dim[2]++) |
234 |
|
for (i = 0; i < il->nsamps; i++) { |
235 |
|
/* next sample point */ |
236 |
< |
peano(sp, 4, urand(ilhash(dim,3)+i), .02); |
236 |
> |
multisamp(sp, 4, urand(ilhash(dim,3)+i)); |
237 |
|
/* random direction */ |
238 |
|
r1 = (dim[1] + sp[0])/nalt; |
239 |
< |
r2 = (dim[2] + sp[1])/nalt; |
239 |
> |
r2 = (dim[2] + sp[1] - .5)/nazi; |
240 |
|
flatdir(dn, r1, r2); |
241 |
|
for (j = 0; j < 3; j++) |
242 |
|
dir[j] = -dn[0]*u[j] - dn[1]*v[j] - dn[2]*co->ad[j]; |
254 |
|
raysamp(distarr+3*(dim[1]*nazi+dim[2]), org, dir, rt); |
255 |
|
} |
256 |
|
rayflush(rt); |
257 |
< |
/* write out the ring w/ distribution */ |
258 |
< |
flatout(il, distarr, nalt, nazi, u, v, co->ad); |
259 |
< |
illumout(il, ob); |
257 |
> |
/* write out the ring and its distribution */ |
258 |
> |
if (average(il, distarr, nalt*nazi)) { |
259 |
> |
if (il->sampdens > 0) |
260 |
> |
flatout(il, distarr, nalt, nazi, u, v, co->ad); |
261 |
> |
illumout(il, ob); |
262 |
> |
} else if (!(il->flags & IL_LIGHT)) |
263 |
> |
printobj(il->altmat, ob); |
264 |
|
/* clean up */ |
265 |
|
freecone(ob); |
266 |
|
free((char *)distarr); |
291 |
|
if (rt->nrays <= 0) |
292 |
|
return; |
293 |
|
bzero(rt->buf+6*rt->nrays, 6*sizeof(float)); |
294 |
+ |
errno = 0; |
295 |
|
if ( process(rt->pd, (char *)rt->buf, (char *)rt->buf, |
296 |
|
3*sizeof(float)*rt->nrays, |
297 |
|
6*sizeof(float)*(rt->nrays+1)) < |