30 #define Rnd64_A 6364136223846793005ULL
32 #define RndF_Mul 5.4210108624275222e-20f
33 #define RndD_Mul 5.4210108624275222e-20
41 static unsigned long long int
42 Rnd64_jump(
unsigned long long int n,
unsigned long long int seed ) {
43 unsigned long long int a_k, c_k, ran;
50 for (i = 0; n; n >>= 1, ++i) {
52 ran = a_k * ran + c_k;
60 #if defined(PLASMA_HAVE_WEAK)
61 #pragma weak CORE_cplrnt = PCORE_cplrnt
62 #define CORE_cplrnt PCORE_cplrnt
65 int bigM,
int m0,
int n0,
unsigned long long int seed )
69 unsigned long long int ran, jump;
71 jump = m0 + n0 * bigM;
73 for (j=0; j<n; ++j ) {
74 ran = Rnd64_jump(
NBELEM*jump, seed );
75 for (i = 0; i < m; ++i) {
94 int bigM,
int m0,
int n0,
unsigned long long int seed )
98 sizeof(
int), &m,
VALUE,
99 sizeof(
int), &n,
VALUE,
101 sizeof(
int), &lda,
VALUE,
102 sizeof(
int), &bigM,
VALUE,
103 sizeof(
int), &m0,
VALUE,
104 sizeof(
int), &n0,
VALUE,
105 sizeof(
unsigned long long int), &seed,
VALUE,
112 #if defined(PLASMA_HAVE_WEAK)
113 #pragma weak CORE_cplrnt_quark = PCORE_cplrnt_quark
114 #define CORE_cplrnt_quark PCORE_cplrnt_quark
125 unsigned long long int seed;