76 |
|
double mult; |
77 |
|
int i; |
78 |
|
|
79 |
< |
getperpendicular(src->ss[SU], src->snorm); |
79 |
> |
getperpendicular(src->ss[SU], src->snorm, rand_samp); |
80 |
|
mult = .5 * sqrt( src->ss2 ); |
81 |
|
for (i = 0; i < 3; i++) |
82 |
|
src->ss[SU][i] *= mult; |
234 |
|
/* set sampling vectors */ |
235 |
|
for (i = 0; i < 3; i++) |
236 |
|
src->ss[SU][i] = .5 * co->al * co->ad[i]; |
237 |
< |
getperpendicular(src->ss[SW], co->ad); |
237 |
> |
getperpendicular(src->ss[SW], co->ad, rand_samp); |
238 |
|
for (i = 0; i < 3; i++) |
239 |
|
src->ss[SW][i] *= .8559 * CO_R0(co); |
240 |
|
fcross(src->ss[SV], src->ss[SW], co->ad); |