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