530 |
|
|
531 |
|
/* Convert 1-dimensional random variable to N-dimensional */ |
532 |
|
void |
533 |
< |
SDmultiSamp(double t[], int n, double randX) |
533 |
> |
SDmultiSamp(RREAL t[], int n, double randX) |
534 |
|
{ |
535 |
|
unsigned nBits; |
536 |
|
double scale; |
564 |
|
{ |
565 |
|
/* convert to position on hemisphere */ |
566 |
|
SDmultiSamp(ioVec, 2, randX); |
567 |
< |
SDsquare2disk(ioVec, ioVec[0], ioVec[1]); |
567 |
> |
square2disk(ioVec, ioVec[0], ioVec[1]); |
568 |
|
ioVec[2] = 1. - ioVec[0]*ioVec[0] - ioVec[1]*ioVec[1]; |
569 |
|
ioVec[2] = sqrt(ioVec[2]*(ioVec[2]>0)); |
570 |
|
if (!outFront) /* going out back? */ |