13 long long i246m1 = 0x00003FFFFFFFFFFFLL;
14 double d2m46 = pow(0.5,46.0);
16 Lx = (
long long) *X; La = (
long long) A;
18 Lx = (Lx * La) & i246m1;
21 return (d2m46 * (
double) Lx);
34 long long i246m1 = 0x00003FFFFFFFFFFFLL;
35 double d2m46 = pow(0.5,46.0);
37 Lx = (
long long) *X; La = (
long long) A;
41 Lx = (Lx * La) & i246m1;
42 Y[i] = d2m46 * (double) Lx;
void vranlc(int N, double *X, double A, Y)
double randlc(double *X, double *A)