r -the random number generator context
Notes
The random type has to be set by PetscRandomSetType().
This is only a primative "parallel" random number generator, it should NOT
be used for sophisticated parallel Monte Carlo methods since it will very likely
not have the correct statistics across processors. You can provide your own
parallel generator using PetscRandomRegister();
If you create a PetscRandom() using PETSC_COMM_SELF on several processors then
the SAME random numbers will be generated on all those processors. Use PETSC_COMM_WORLD
or the appropriate parallel communicator to eliminate this issue.
Use VecSetRandom() to set the elements of a vector to random numbers.
Example of Usage
PetscRandomCreate(PETSC_COMM_SELF,&r);
PetscRandomSetType(r,PETSCRAND48);
PetscRandomGetValue(r,&value1);
PetscRandomGetValueReal(r,&value2);
PetscRandomDestroy(&r);
See Also
PetscRandomSetType(), PetscRandomGetValue(), PetscRandomGetValueReal(), PetscRandomSetInterval(),
PetscRandomDestroy(), VecSetRandom(), PetscRandomType
Level
intermediate
Location
src/sys/classes/random/interface/randomc.c
Examples
src/sys/classes/random/tutorials/ex1.c.html
src/sys/classes/random/tutorials/ex2.c.html
src/vec/vec/utils/tagger/tutorials/ex1.c.html
src/vec/vec/tutorials/performance.c.html
src/mat/tutorials/ex2.c.html
src/mat/tutorials/ex9.c.html
src/dm/tutorials/swarm_ex3.c.html
src/dm/field/tutorials/ex1.c.html
src/ksp/pc/tutorials/ex3.c.html
src/ksp/ksp/tutorials/ex5.c.html
src/ksp/ksp/tutorials/ex11.c.html
Implementations
PetscRandomCreate_Rand in src/sys/classes/random/impls/rand/rand.c
PetscRandomCreate_Rand48 in src/sys/classes/random/impls/rand48/rand48.c
PetscRandomCreate_Rander48 in src/sys/classes/random/impls/rander48/rander48.c
PetscRandomCreate_Random123 in src/sys/classes/random/impls/random123/random123.c
PetscRandomCreate_Sprng in src/sys/classes/random/impls/sprng/sprng.c
Index of all Sys routines
Table of Contents for all manual pages
Index of all manual pages