88 |
|
h = ilhash(dim, 3) + i; |
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]; |
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! */ |
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]; |
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)) < |