1 |
< |
/* Copyright (c) 1991 Regents of the University of California */ |
1 |
> |
/* Copyright (c) 1995 Regents of the University of California */ |
2 |
|
|
3 |
|
#ifndef lint |
4 |
|
static char SCCSid[] = "$SunId$ LBL"; |
163 |
|
nalt = nazi = 1; |
164 |
|
else { |
165 |
|
n = 4.*PI * il->sampdens; |
166 |
< |
nalt = sqrt(n/PI) + .5; |
167 |
< |
nazi = PI*nalt + .5; |
166 |
> |
nalt = sqrt(2./PI*n) + .5; |
167 |
> |
nazi = PI/2.*nalt + .5; |
168 |
|
} |
169 |
|
n = nalt*nazi; |
170 |
|
distarr = (float *)calloc(n, 3*sizeof(float)); |
260 |
|
r1 = r3*cos(r2); |
261 |
|
r2 = r3*sin(r2); |
262 |
|
for (j = 0; j < 3; j++) |
263 |
< |
org[j] = CO_P0(co)[j] + r1*u[j] + r1*v[j] + |
263 |
> |
org[j] = CO_P0(co)[j] + r1*u[j] + r2*v[j] + |
264 |
|
.001*co->ad[j]; |
265 |
|
|
266 |
|
/* send sample */ |