24 return (R)(drand48());
26 return ((R)rand())/((R)RAND_MAX);
30 void Y(srand48)(
long int seed)
35 srand((
unsigned int)seed);
39 void Y(vrand_unit_complex)(C *x,
const INT n)
43 for (k = 0; k < n; k++)
44 x[k] = Y(drand48)() + II * Y(drand48)();
47 void Y(vrand_shifted_unit_double)(R *x,
const INT n)
51 for (k = 0; k < n; k++)
52 x[k] = Y(drand48)() - K(0.5);
55 void Y(vrand_real)(R *x,
const INT n,
const R a,
const R b)
59 for (k = 0; k < n; k++)
60 x[k] = a + Y(drand48)() * (b - a);