X-Git-Url: http://info.iut-bm.univ-fcomte.fr/pub/gitweb/simgrid.git/blobdiff_plain/6de03ecc4e630732984a0673512a5d15fd75e270..4bd1f48f0bf1ad1703be680ec2a38d626c6a2668:/src/xbt/RngStream.c diff --git a/src/xbt/RngStream.c b/src/xbt/RngStream.c index e4dbbb5a74..816d5d6497 100644 --- a/src/xbt/RngStream.c +++ b/src/xbt/RngStream.c @@ -17,19 +17,16 @@ * \***********************************************************************/ - #include "xbt/RngStream.h" #include "xbt/sysdep.h" #include #include #include - /*---------------------------------------------------------------------*/ /* Private part. */ /*---------------------------------------------------------------------*/ - #define norm 2.328306549295727688e-10 #define m1 4294967087.0 #define m2 4294944443.0 @@ -42,14 +39,9 @@ #define two53 9007199254740992.0 #define fact 5.9604644775390625e-8 /* 1 / 2^24 */ - - - -/* Default initial seed of the package. Will be updated to become - the seed of the next created stream. */ +/* Default initial seed of the package. Will be updated to become the seed of the next created stream. */ static double nextSeed[6] = { 12345, 12345, 12345, 12345, 12345, 12345 }; - /* The following are the transition matrices of the two MRG components */ /* (in matrix form), raised to the powers -1, 1, 2^76, and 2^127, resp.*/ static double InvA1[3][3] = { /* Inverse of A1p0 */ @@ -100,12 +92,6 @@ static double A2p127[3][3] = { { 2824425944.0, 32183930.0, 2093834863.0 } }; - - - - -/*-------------------------------------------------------------------------*/ - static double MultModM (double a, double s, double c, double m) /* Compute (a*s + c) % m. m must be < 2^35. Works also for s, c < 0 */ { @@ -127,9 +113,6 @@ static double MultModM (double a, double s, double c, double m) return v; } - -/*-------------------------------------------------------------------------*/ - static void MatVecModM (double A[3][3], double s[3], double v[3], double m) /* Returns v = A*s % m. Assumes that -m < s[i] < m. */ /* Works even if v = s. */ @@ -145,11 +128,7 @@ static void MatVecModM (double A[3][3], double s[3], double v[3], double m) v[i] = x[i]; } - -/*-------------------------------------------------------------------------*/ - -static void MatMatModM (double A[3][3], double B[3][3], double C[3][3], - double m) +static void MatMatModM (double A[3][3], double B[3][3], double C[3][3], double m) /* Returns C = A*B % m. Work even if A = C or B = C or A = B = C. */ { int i, j; @@ -167,9 +146,6 @@ static void MatMatModM (double A[3][3], double B[3][3], double C[3][3], } } - -/*-------------------------------------------------------------------------*/ - static void MatTwoPowModM (double A[3][3], double B[3][3], double m, long e) /* Compute matrix B = (A^(2^e) % m); works even if A = B */ { @@ -187,9 +163,6 @@ static void MatTwoPowModM (double A[3][3], double B[3][3], double m, long e) MatMatModM (B, B, B, m); } - -/*-------------------------------------------------------------------------*/ - static void MatPowModM (double A[3][3], double B[3][3], double m, long n) /* Compute matrix B = A^n % m ; works even if A = B */ { @@ -215,9 +188,6 @@ static void MatPowModM (double A[3][3], double B[3][3], double m, long n) } } - -/*-------------------------------------------------------------------------*/ - static double U01 (RngStream g) { long k; @@ -248,9 +218,6 @@ static double U01 (RngStream g) return (g->Anti) ? (1 - u) : u; } - -/*-------------------------------------------------------------------------*/ - static double U01d (RngStream g) { double u; @@ -265,13 +232,9 @@ static double U01d (RngStream g) } } - -/*-------------------------------------------------------------------------*/ - static int CheckSeed (unsigned long seed[6]) { - /* Check that the seeds are legitimate values. Returns 0 if legal seeds, - -1 otherwise */ + /* Check that the seeds are legitimate values. Returns 0 if legal seeds, -1 otherwise */ int i; for (i = 0; i < 3; ++i) { @@ -306,12 +269,10 @@ static int CheckSeed (unsigned long seed[6]) return 0; } - /*---------------------------------------------------------------------*/ /* Public part. */ /*---------------------------------------------------------------------*/ - RngStream RngStream_CreateStream (const char name[]) { int i; @@ -340,8 +301,6 @@ RngStream RngStream_CreateStream (const char name[]) return g; } -/*-------------------------------------------------------------------------*/ - void RngStream_DeleteStream (RngStream * p) { if (*p == NULL) @@ -351,8 +310,6 @@ void RngStream_DeleteStream (RngStream * p) *p = NULL; } -/*-------------------------------------------------------------------------*/ - RngStream RngStream_CopyStream (const RngStream src) { RngStream g; @@ -372,8 +329,6 @@ RngStream RngStream_CopyStream (const RngStream src) return g; } -/*-------------------------------------------------------------------------*/ - void RngStream_ResetStartStream (RngStream g) { int i; @@ -381,8 +336,6 @@ void RngStream_ResetStartStream (RngStream g) g->Cg[i] = g->Bg[i] = g->Ig[i]; } -/*-------------------------------------------------------------------------*/ - void RngStream_ResetNextSubstream (RngStream g) { int i; @@ -392,8 +345,6 @@ void RngStream_ResetNextSubstream (RngStream g) g->Cg[i] = g->Bg[i]; } -/*-------------------------------------------------------------------------*/ - void RngStream_ResetStartSubstream (RngStream g) { int i; @@ -401,8 +352,6 @@ void RngStream_ResetStartSubstream (RngStream g) g->Cg[i] = g->Bg[i]; } -/*-------------------------------------------------------------------------*/ - int RngStream_SetPackageSeed (unsigned long seed[6]) { int i; @@ -413,8 +362,6 @@ int RngStream_SetPackageSeed (unsigned long seed[6]) return 0; /* SUCCESS */ } -/*-------------------------------------------------------------------------*/ - int RngStream_SetSeed (RngStream g, unsigned long seed[6]) { int i; @@ -425,8 +372,6 @@ int RngStream_SetSeed (RngStream g, unsigned long seed[6]) return 0; /* SUCCESS */ } -/*-------------------------------------------------------------------------*/ - void RngStream_AdvanceState (RngStream g, long e, long c) { double B1[3][3], C1[3][3], B2[3][3], C2[3][3]; @@ -456,8 +401,6 @@ void RngStream_AdvanceState (RngStream g, long e, long c) MatVecModM (C2, &g->Cg[3], &g->Cg[3], m2); } -/*-------------------------------------------------------------------------*/ - void RngStream_GetState (RngStream g, unsigned long seed[6]) { int i; @@ -465,8 +408,6 @@ void RngStream_GetState (RngStream g, unsigned long seed[6]) seed[i] = g->Cg[i]; } -/*-------------------------------------------------------------------------*/ - void RngStream_WriteState (RngStream g) { int i; @@ -483,8 +424,6 @@ void RngStream_WriteState (RngStream g) printf ("%lu }\n\n", (unsigned long) g->Cg[5]); } -/*-------------------------------------------------------------------------*/ - void RngStream_WriteStateFull (RngStream g) { int i; @@ -515,22 +454,16 @@ void RngStream_WriteStateFull (RngStream g) printf ("%lu }\n\n", (unsigned long) g->Cg[5]); } -/*-------------------------------------------------------------------------*/ - void RngStream_IncreasedPrecis (RngStream g, int incp) { g->IncPrec = incp; } -/*-------------------------------------------------------------------------*/ - void RngStream_SetAntithetic (RngStream g, int a) { g->Anti = a; } -/*-------------------------------------------------------------------------*/ - double RngStream_RandU01 (RngStream g) { if (g->IncPrec) @@ -539,8 +472,6 @@ double RngStream_RandU01 (RngStream g) return U01 (g); } -/*-------------------------------------------------------------------------*/ - int RngStream_RandInt (RngStream g, int i, int j) { return i + (int) ((j - i + 1.0) * RngStream_RandU01 (g));