#include <petscsys.h>
#include <petscviewer.h>
int main(int argc,char **argv)
{
  PetscMPIInt rank, size;
  PetscInt seed=100, iters=10;
  PetscRandom rng;
  PetscScalar x;
  PetscInt i;
  
  
  PetscInitialize(&argc,&argv,NULL,NULL);
  MPI_Comm_rank(PETSC_COMM_WORLD,&rank);
  MPI_Comm_size(PETSC_COMM_WORLD,&size);

  PetscOptionsGetInt(NULL, "-num_iters", &iters, NULL);
  PetscOptionsGetInt(NULL, "-seed", &seed, NULL);

  PetscRandomCreate(PETSC_COMM_WORLD,&rng);
  PetscRandomSetType(rng,PETSCSPRNG);
  PetscRandomSetSeed(rng, seed);
  PetscRandomSeed(rng);

  PetscRandomView(rng,PETSC_VIEWER_STDOUT_WORLD);

  for(i=0;i<iters;i++){
    PetscRandomGetValue(rng, &x);
    PetscPrintf(PETSC_COMM_SELF, "[%i]: %g\n", rank, x);
  }

  PetscRandomDestroy(&rng);
  PetscFinalize();
  return 0;
}