17 static double _dorander48(PetscRandom_Rander48 *r48)  in _dorander48()  argument
22   accu    = (unsigned long)r48->mult[0] * (unsigned long)r48->seed[0] + (unsigned long)r48->add;  in _dorander48()
25 …accu += (unsigned long)r48->mult[0] * (unsigned long)r48->seed[1] + (unsigned long)r48->mult[1] * …  in _dorander48()
28 …ccu += (unsigned long)r48->mult[0] * r48->seed[2] + (unsigned long)r48->mult[1] * r48->seed[1] + (…  in _dorander48()
29   r48->seed[0] = temp[0];  in _dorander48()
30   r48->seed[1] = temp[1];  in _dorander48()
31   r48->seed[2] = (unsigned short)accu;  in _dorander48()
32 …return ldexp((double)r48->seed[0], -48) + ldexp((double)r48->seed[1], -32) + ldexp((double)r48->se…  in _dorander48()
37   PetscRandom_Rander48 *r48 = (PetscRandom_Rander48 *)r->data;  in PetscRandomSeed_Rander48()  local
40   r48->seed[0] = RANDER48_SEED_0;  in PetscRandomSeed_Rander48()
41   r48->seed[1] = (unsigned short)r->seed;  in PetscRandomSeed_Rander48()
42   r48->seed[2] = (unsigned short)(r->seed >> 16);  in PetscRandomSeed_Rander48()
43   r48->mult[0] = RANDER48_MULT_0;  in PetscRandomSeed_Rander48()
44   r48->mult[1] = RANDER48_MULT_1;  in PetscRandomSeed_Rander48()
45   r48->mult[2] = RANDER48_MULT_2;  in PetscRandomSeed_Rander48()
46   r48->add     = RANDER48_ADD;  in PetscRandomSeed_Rander48()
52   PetscRandom_Rander48 *r48 = (PetscRandom_Rander48 *)r->data;  in PetscRandomGetValue_Rander48()  local
58     if (PetscRealPart(r->width)) *val += PetscRealPart(r->width) * _dorander48(r48);  in PetscRandomGetValue_Rander48()
59 …if (PetscImaginaryPart(r->width)) *val += PetscImaginaryPart(r->width) * _dorander48(r48) * PETSC_…  in PetscRandomGetValue_Rander48()
61     *val = _dorander48(r48) + _dorander48(r48) * PETSC_i;  in PetscRandomGetValue_Rander48()
64   if (r->iset) *val = r->width * _dorander48(r48) + r->low;  in PetscRandomGetValue_Rander48()
65   else *val = _dorander48(r48);  in PetscRandomGetValue_Rander48()
72   PetscRandom_Rander48 *r48 = (PetscRandom_Rander48 *)r->data;  in PetscRandomGetValueReal_Rander48()  local
76   if (r->iset) *val = PetscRealPart(r->width) * _dorander48(r48) + PetscRealPart(r->low);  in PetscRandomGetValueReal_Rander48()
77   else *val = _dorander48(r48);  in PetscRandomGetValueReal_Rander48()
79   if (r->iset) *val = r->width * _dorander48(r48) + r->low;  in PetscRandomGetValueReal_Rander48()
80   else *val = _dorander48(r48);  in PetscRandomGetValueReal_Rander48()
125   PetscRandom_Rander48 *r48;  in PetscRandomCreate_Rander48()  local
128   PetscCall(PetscNew(&r48));  in PetscRandomCreate_Rander48()
130   r->data   = r48;  in PetscRandomCreate_Rander48()