| 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? */ |