| 72 |
|
hlist[0] = r->rno; |
| 73 |
|
hlist[1] = dp->t; |
| 74 |
|
hlist[2] = dp->p; |
| 75 |
< |
peano(spt, 2, urand(ilhash(hlist,3)+dp->n), .01); |
| 75 |
> |
multisamp(spt, 2, urand(ilhash(hlist,3)+dp->n)); |
| 76 |
|
zd = sqrt((dp->t + spt[0])/h->nt); |
| 77 |
|
phi = 2.0*PI * (dp->p + spt[1])/h->np; |
| 78 |
|
xd = cos(phi) * zd; |
| 237 |
|
{ |
| 238 |
|
register int i; |
| 239 |
|
/* set number of divisions */ |
| 240 |
+ |
if (wt < (.25*PI)/ambdiv+FTINY) { |
| 241 |
+ |
hp->nt = hp->np = 0; |
| 242 |
+ |
return; /* zero samples */ |
| 243 |
+ |
} |
| 244 |
|
hp->nt = sqrt(ambdiv * wt / PI) + 0.5; |
| 245 |
< |
hp->np = PI * hp->nt; |
| 245 |
> |
hp->np = PI * hp->nt + 0.5; |
| 246 |
|
/* make axes */ |
| 247 |
|
VCOPY(hp->uz, r->ron); |
| 248 |
|
hp->uy[0] = hp->uy[1] = hp->uy[2] = 0.0; |