140 |
|
} |
141 |
|
if (dstrsrc > FTINY) { |
142 |
|
/* distribute source direction */ |
143 |
< |
dimlist[ndims++] = sn; |
144 |
< |
for (i = 0; i < 3; i++) { |
145 |
< |
dimlist[ndims] = i + 8831; |
146 |
< |
vd[i] = dstrsrc * source[sn].ss * |
147 |
< |
(1.0 - 2.0*urand(urind(ilhash(dimlist,ndims+1),samplendx))); |
148 |
< |
} |
149 |
< |
ndims--; |
143 |
> |
dimlist[ndims] = sn + 8831; |
144 |
> |
peano(vd, 3, urand(ilhash(dimlist,ndims+1)+samplendx), .01); |
145 |
> |
for (i = 0; i < 3; i++) |
146 |
> |
vd[i] = dstrsrc * source[sn].ss * (1. - 2.*vd[i]); |
147 |
|
if (source[sn].sflags & SFLAT) { /* project offset */ |
148 |
|
d = DOT(vd, source[sn].snorm); |
149 |
|
for (i = 0; i < 3; i++) |