26 return (R)(drand48());
28 return ((R)rand())/((R)RAND_MAX);
32 void Y(srand48)(
long int seed)
37 srand((
unsigned int)seed);
41 void Y(vrand_unit_complex)(C *x,
const INT n)
45 for (k = 0; k < n; k++)
46 x[k] = Y(drand48)() + II * Y(drand48)();
49 void Y(vrand_shifted_unit_double)(R *x,
const INT n)
53 for (k = 0; k < n; k++)
54 x[k] = Y(drand48)() - K(0.5);
57 void Y(vrand_real)(R *x,
const INT n,
const R a,
const R b)
61 for (k = 0; k < n; k++)
62 x[k] = a + Y(drand48)() * (b - a);