123 |
|
MAT4 ixfm |
124 |
|
) |
125 |
|
{ |
126 |
< |
int nalt, nazi; |
126 |
> |
int nalt=1, nazi=1; |
127 |
|
SRCINDEX si; |
128 |
|
RAY sr; |
129 |
|
FVECT v; |
130 |
|
double d; |
131 |
|
int i, j; |
132 |
|
/* get sampling density */ |
133 |
< |
if (il->sampdens <= 0) { |
134 |
< |
nalt = nazi = 1; |
135 |
< |
} else { |
133 |
> |
if (il->sd == NULL && il->sampdens > 0) { |
134 |
|
i = PI * il->sampdens; |
135 |
|
nalt = sqrt(i/PI) + .5; |
136 |
|
nazi = PI*nalt + .5; |
154 |
|
d = -1.0001*il->thick - 5.*FTINY; |
155 |
|
else |
156 |
|
d = 5.*FTINY; |
157 |
< |
for (i = 3; i--; ) |
160 |
< |
sr.rorg[i] += d*nrm[i]; |
157 |
> |
VSUM(sr.rorg, sr.rorg, nrm, d); |
158 |
|
samplendx++; /* increment sample counter */ |
159 |
|
if (!srcray(&sr, NULL, &si)) |
160 |
|
break; /* end of sources */ |