15c6c1daeSBarry Smith 2b6bfaf9aSBarry Smith #include <../src/sys/classes/random/randomimpl.h> 35c6c1daeSBarry Smith 45c6c1daeSBarry Smith #define USE_MPI 55c6c1daeSBarry Smith #define SIMPLE_SPRNG 65c6c1daeSBarry Smith EXTERN_C_BEGIN 75c6c1daeSBarry Smith #include <sprng.h> 85c6c1daeSBarry Smith EXTERN_C_END 95c6c1daeSBarry Smith 105c6c1daeSBarry Smith PetscErrorCode PetscRandomSeed_Sprng(PetscRandom r) 115c6c1daeSBarry Smith { 125c6c1daeSBarry Smith PetscFunctionBegin; 135c6c1daeSBarry Smith init_sprng(r->seed,SPRNG_DEFAULT); 145c6c1daeSBarry Smith PetscFunctionReturn(0); 155c6c1daeSBarry Smith } 165c6c1daeSBarry Smith 175c6c1daeSBarry Smith PetscErrorCode PetscRandomGetValue_Sprng(PetscRandom r,PetscScalar *val) 185c6c1daeSBarry Smith { 195c6c1daeSBarry Smith PetscFunctionBegin; 205c6c1daeSBarry Smith #if defined(PETSC_USE_COMPLEX) 215c6c1daeSBarry Smith if (r->iset) { 22b0043f70SBarry Smith *val = PetscRealPart(r->width)*sprng() + PetscRealPart(r->low) + (PetscImaginaryPart(r->width)*sprng() + PetscImaginaryPart(r->low)) * PETSC_i; 235c6c1daeSBarry Smith } else { 245c6c1daeSBarry Smith *val = sprng() + sprng()*PETSC_i; 255c6c1daeSBarry Smith } 265c6c1daeSBarry Smith #else 275c6c1daeSBarry Smith if (r->iset) *val = r->width * sprng() + r->low; 285c6c1daeSBarry Smith else *val = sprng(); 295c6c1daeSBarry Smith #endif 305c6c1daeSBarry Smith PetscFunctionReturn(0); 315c6c1daeSBarry Smith } 325c6c1daeSBarry Smith 33*4d1a973fSToby Isaac PetscErrorCode PetscRandomGetValueReal_Sprng(PetscRandom r,PetscReal *val) 345c6c1daeSBarry Smith { 355c6c1daeSBarry Smith PetscFunctionBegin; 365c6c1daeSBarry Smith #if defined(PETSC_USE_COMPLEX) 375c6c1daeSBarry Smith if (r->iset) *val = PetscRealPart(r->width)*sprng() + PetscRealPart(r->low); 385c6c1daeSBarry Smith else *val = sprng(); 395c6c1daeSBarry Smith #else 405c6c1daeSBarry Smith if (r->iset) *val = r->width * sprng() + r->low; 415c6c1daeSBarry Smith else *val = sprng(); 425c6c1daeSBarry Smith #endif 435c6c1daeSBarry Smith PetscFunctionReturn(0); 445c6c1daeSBarry Smith } 455c6c1daeSBarry Smith 465c6c1daeSBarry Smith static struct _PetscRandomOps PetscRandomOps_Values = { 475c6c1daeSBarry Smith /* 0 */ 485c6c1daeSBarry Smith PetscRandomSeed_Sprng, 495c6c1daeSBarry Smith PetscRandomGetValue_Sprng, 505c6c1daeSBarry Smith PetscRandomGetValueReal_Sprng, 515c6c1daeSBarry Smith 0, 525c6c1daeSBarry Smith /* 5 */ 535c6c1daeSBarry Smith 0 545c6c1daeSBarry Smith }; 555c6c1daeSBarry Smith 565c6c1daeSBarry Smith /*MC 575c6c1daeSBarry Smith PETSCSPRNG- access to the publically available random number generator sprng 585c6c1daeSBarry Smith 595c6c1daeSBarry Smith Options Database Keys: 605c6c1daeSBarry Smith . -random_type <rand,rand48,sprng> 615c6c1daeSBarry Smith 625c6c1daeSBarry Smith Level: beginner 635c6c1daeSBarry Smith 645c6c1daeSBarry Smith PETSc must have been ./configure with the option --download-sprng to use 655c6c1daeSBarry Smith this random number generator. 665c6c1daeSBarry Smith 675c6c1daeSBarry Smith This is NOT currently using a parallel random number generator. Sprng does have 685c6c1daeSBarry Smith an MPI version we should investigate. 695c6c1daeSBarry Smith 705c6c1daeSBarry Smith .seealso: RandomCreate(), RandomSetType(), PETSCRAND, PETSCRAND48 715c6c1daeSBarry Smith M*/ 725c6c1daeSBarry Smith 738cc058d9SJed Brown PETSC_EXTERN PetscErrorCode PetscRandomCreate_Sprng(PetscRandom r) 745c6c1daeSBarry Smith { 755c6c1daeSBarry Smith PetscErrorCode ierr; 765c6c1daeSBarry Smith 775c6c1daeSBarry Smith PetscFunctionBegin; 785c6c1daeSBarry Smith ierr = PetscMemcpy(r->ops,&PetscRandomOps_Values,sizeof(PetscRandomOps_Values));CHKERRQ(ierr); 795c6c1daeSBarry Smith ierr = PetscObjectChangeTypeName((PetscObject)r,PETSCSPRNG);CHKERRQ(ierr); 805c6c1daeSBarry Smith PetscFunctionReturn(0); 815c6c1daeSBarry Smith } 82