| 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)) < |