#include <petscsys.h>

#define SIMPLE_SPRNG
#define USE_MPI 
#include <sprng.h>

int main(int argc,char **argv)
{
  PetscMPIInt rank, size;
  PetscInt seed=100, iters=10;
  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);

  init_sprng(seed, SPRNG_DEFAULT);

  print_sprng();

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