#include <PhotosRandom.h>

Inheritance diagram for Photospp::PhotosRandom:

Static Public Member Functions
static void	setSeed (int s1, int s2)

static void	initialize ()

static double	randomReal ()

static void	setSeed (int s1, int s2)

static void	initialize ()

static double	randomReal ()

Static Protected Attributes
static bool	init = false

static int	iseed [2] = { 1802, 9373 }

static int	i97 = 96

static int	j97 = 32

static double	uran [97] = { 0.0 }

static double	cran = 362436.0 /16777216.0

static const double	cdran = 7654321.0 /16777216.0

static const double	cmran = 16777213.0/16777216.0

Detailed Description

Definition at line 20 of file PhotosRandom.h.

Member Function Documentation

void PhotosRandom::initialize ( )

static

Definition at line 45 of file PhotosRandom.cxx.

References Photospp::Log::Debug(), i, init, iseed, T, and uran.

Referenced by Photospp::Photos::initialize().

 {
         long IS1,IS2,IS3,IS4,IS5;
         double S,T;
 
 // Calculate Marsaglia and Zaman seeds (by F. James)
         IS1=(iseed[0]/177)%177+2;
         IS2= iseed[0]%177+2;
         IS3=(iseed[1]/169)%178+1;
         IS4= iseed[1]%169;
         for(int i=0;i<97;i++)
         {
                 S=0.0;
                 T=0.5;
                 for(int j=0;j<24;j++)
                 {
                         IS5=( ((IS1*IS2)%179)*IS3 )%179;
                         IS1=IS2;
                         IS2=IS3;
                         IS3=IS5;
                         IS4=(53*IS4+1)%169;
                         if( (IS4*IS5)%64>=32) S=S+T;
                         T=0.5*T;
                 }
                 uran[i]=S;
         }
         init=true;
         Log::Debug(0)<<"PhotosRandom::inititalize(): seed: "<<iseed[0]<<", "<<iseed[1]<<std::endl;
 }

static void Photospp::PhotosRandom::initialize ( )

static

static double Photospp::PhotosRandom::randomReal ( )

static

double PhotosRandom::randomReal ( )

static

Definition at line 91 of file PhotosRandom.cxx.

References cdran, cmran, cran, Photospp::Log::Fatal(), i97, init, j97, and uran.

Referenced by Photospp::Photos::setRandomGenerator().

 {
         if(!init) Log::Fatal("PhotosRandom::randomReal(): generator not initialized",1);
         double ret=0.0;
         while(true)
         {
                 ret = uran[i97]-uran[j97];
                 if(ret<0.0) ret+=1.;
                 uran[i97]=ret;
                 i97--;
                 if(i97<0) i97=96;
                 j97--;
                 if(j97<0) j97=96;
                 cran-=cdran;
                 if(cran<0.0) cran+=cmran;
                 ret-=cran;
                 if(ret<0.0) ret+=1.0;
                 if(ret>0.0) break;
         }
         return ret;
 }

static void Photospp::PhotosRandom::setSeed	(	int	s1,
		int	s2
	)

static

void PhotosRandom::setSeed	(	int	s1,
		int	s2
	)

static

Definition at line 25 of file PhotosRandom.cxx.

References Photospp::Log::Fatal(), and iseed.

Referenced by Photospp::Photos::setSeed().

 {
         if(s1<0 || s1>31327) Log::Fatal("PhotosRandom::setSeed(): Seed(1) out of range [0,31327]",8);
         if(s2<0 || s2>30080) Log::Fatal("PhotosRandom::setSeed(): Seed(2) out of range [0,30080]",9);
         iseed[0]=s1;
         iseed[1]=s2;
 }