1! 2! Program to test random number generation routines from fortran. 3! 4 program main 5 6#include <petsc/finclude/petscsys.h> 7 use petscsys 8 implicit none 9 10 PetscErrorCode ierr 11 PetscRandom r 12 PetscScalar rand 13 14 PetscCallA(PetscInitialize(ierr)) 15 16 PetscCallA(PetscRandomCreate(PETSC_COMM_WORLD,r,ierr)) 17 PetscCallA(PetscRandomSetFromOptions(r,ierr)) 18 PetscCallA(PetscRandomGetValue(r,rand,ierr)) 19 print*, 'Random value:',rand 20 PetscCallA(PetscRandomDestroy(r,ierr)) 21 PetscCallA(PetscFinalize(ierr)) 22 end 23 24! 25!/*TEST 26! 27! test: 28! requires: !complex 29! 30!TEST*/ 31