35 static PetscErrorCode PetscRandomSeed_Rander48(PetscRandom r)  in PetscRandomSeed_Rander48()  argument
37   PetscRandom_Rander48 *r48 = (PetscRandom_Rander48 *)r->data;  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()
50 static PetscErrorCode PetscRandomGetValue_Rander48(PetscRandom r, PetscScalar *val)  in PetscRandomGetValue_Rander48()  argument
52   PetscRandom_Rander48 *r48 = (PetscRandom_Rander48 *)r->data;  in PetscRandomGetValue_Rander48()
56   if (r->iset) {  in PetscRandomGetValue_Rander48()
57     *val = PetscRealPart(r->low) + PetscImaginaryPart(r->low) * PETSC_i;  in PetscRandomGetValue_Rander48()
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()
64   if (r->iset) *val = r->width * _dorander48(r48) + r->low;  in PetscRandomGetValue_Rander48()
70 static PetscErrorCode PetscRandomGetValueReal_Rander48(PetscRandom r, PetscReal *val)  in PetscRandomGetValueReal_Rander48()  argument
72   PetscRandom_Rander48 *r48 = (PetscRandom_Rander48 *)r->data;  in PetscRandomGetValueReal_Rander48()
76   if (r->iset) *val = PetscRealPart(r->width) * _dorander48(r48) + PetscRealPart(r->low);  in PetscRandomGetValueReal_Rander48()
79   if (r->iset) *val = r->width * _dorander48(r48) + r->low;  in PetscRandomGetValueReal_Rander48()
85 static PetscErrorCode PetscRandomDestroy_Rander48(PetscRandom r)  in PetscRandomDestroy_Rander48()  argument
88   PetscCall(PetscFree(r->data));  in PetscRandomDestroy_Rander48()
123 PETSC_EXTERN PetscErrorCode PetscRandomCreate_Rander48(PetscRandom r)  in PetscRandomCreate_Rander48()  argument
130   r->data   = r48;  in PetscRandomCreate_Rander48()
131   r->ops[0] = PetscRandomOps_Values;  in PetscRandomCreate_Rander48()
132   PetscCall(PetscObjectChangeTypeName((PetscObject)r, PETSCRANDER48));  in PetscRandomCreate_Rander48()