#include "petscsys.h" PetscErrorCode PetscRandomCreate(MPI_Comm comm,PetscRandom *r)Collective on MPI_Comm
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.
PetscRandomCreate(PETSC_COMM_SELF,&r); PetscRandomSetType(r,PETSCRAND48); PetscRandomGetValue(r,&value1); PetscRandomGetValueReal(r,&value2); PetscRandomDestroy(&r);
Level:intermediate
Location:src/sys/classes/random/interface/randomc.c
Index of all Sys routines
Table of Contents for all manual pages
Index of all manual pages
comm | - MPI communicator | |
r | - the random number generator context |