1 ------------------------------------------------------------------------------
3 -- GNAT RUN-TIME COMPONENTS --
5 -- G N A T . R A N D O M _ N U M B E R S --
9 -- Copyright (C) 2007-2013, Free Software Foundation, Inc. --
11 -- GNAT is free software; you can redistribute it and/or modify it under --
12 -- terms of the GNU General Public License as published by the Free Soft- --
13 -- ware Foundation; either version 3, or (at your option) any later ver- --
14 -- sion. GNAT is distributed in the hope that it will be useful, but WITH- --
15 -- OUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY --
16 -- or FITNESS FOR A PARTICULAR PURPOSE. --
18 -- As a special exception under Section 7 of GPL version 3, you are granted --
19 -- additional permissions described in the GCC Runtime Library Exception, --
20 -- version 3.1, as published by the Free Software Foundation. --
22 -- You should have received a copy of the GNU General Public License and --
23 -- a copy of the GCC Runtime Library Exception along with this program; --
24 -- see the files COPYING3 and COPYING.RUNTIME respectively. If not, see --
25 -- <http://www.gnu.org/licenses/>. --
27 -- GNAT was originally developed by the GNAT team at New York University. --
28 -- Extensive contributions were provided by Ada Core Technologies Inc. --
30 ------------------------------------------------------------------------------
32 -- Extended pseudo-random number generation
34 -- This package provides a type representing pseudo-random number generators,
35 -- and subprograms to extract various distributions of numbers from them. It
36 -- also provides types for representing initialization values and snapshots of
37 -- internal generator state, which permit reproducible pseudo-random streams.
39 -- The generator currently provided by this package has an extremely long
40 -- period (at least 2**19937-1), and passes the Big Crush test suite, with the
41 -- exception of the two linear complexity tests. Therefore, it is suitable for
42 -- simulations, but should not be used as a cryptographic pseudo-random source
43 -- without additional processing.
45 -- The design of this package effects is simplified compared to the design
46 -- of standard Ada.Numerics packages. There is no separate State type; the
47 -- Generator type itself suffices for this purpose. The parameter modes on
48 -- Reset procedures better reflect the effect of these routines.
50 with System
.Random_Numbers
;
51 with Interfaces
; use Interfaces
;
53 package GNAT
.Random_Numbers
is
55 type Generator
is limited private;
56 subtype Initialization_Vector
is
57 System
.Random_Numbers
.Initialization_Vector
;
59 function Random
(Gen
: Generator
) return Float;
60 function Random
(Gen
: Generator
) return Long_Float;
61 -- Return pseudo-random numbers uniformly distributed on [0 .. 1)
63 function Random
(Gen
: Generator
) return Interfaces
.Integer_32
;
64 function Random
(Gen
: Generator
) return Interfaces
.Unsigned_32
;
65 function Random
(Gen
: Generator
) return Interfaces
.Integer_64
;
66 function Random
(Gen
: Generator
) return Interfaces
.Unsigned_64
;
67 function Random
(Gen
: Generator
) return Integer;
68 function Random
(Gen
: Generator
) return Long_Integer;
69 -- Return pseudo-random numbers uniformly distributed on T'First .. T'Last
70 -- for various builtin integer types.
73 type Result_Subtype
is (<>);
74 Default_Min
: Result_Subtype
:= Result_Subtype
'Val (0);
75 function Random_Discrete
77 Min
: Result_Subtype
:= Default_Min
;
78 Max
: Result_Subtype
:= Result_Subtype
'Last) return Result_Subtype
;
79 -- Returns pseudo-random numbers uniformly distributed on Min .. Max
82 type Result_Subtype
is delta <>;
83 Default_Min
: Result_Subtype
:= 0.0;
84 function Random_Ordinary_Fixed
86 Min
: Result_Subtype
:= Default_Min
;
87 Max
: Result_Subtype
:= Result_Subtype
'Last) return Result_Subtype
;
88 -- Returns pseudo-random numbers uniformly distributed on Min .. Max
91 type Result_Subtype
is delta <> digits <>;
92 Default_Min
: Result_Subtype
:= 0.0;
93 function Random_Decimal_Fixed
95 Min
: Result_Subtype
:= Default_Min
;
96 Max
: Result_Subtype
:= Result_Subtype
'Last) return Result_Subtype
;
97 -- Returns pseudo-random numbers uniformly distributed on Min .. Max
100 type Result_Subtype
is digits <>;
101 function Random_Float
(Gen
: Generator
) return Result_Subtype
;
102 -- Returns pseudo-random numbers uniformly distributed on [0.0 .. 1.0)
104 function Random_Gaussian
(Gen
: Generator
) return Long_Float;
105 function Random_Gaussian
(Gen
: Generator
) return Float;
106 -- Returns pseudo-random numbers normally distributed value with mean 0
107 -- and standard deviation 1.0.
109 procedure Reset
(Gen
: out Generator
);
110 -- Re-initialize the state of Gen from the time of day
113 (Gen
: out Generator
;
114 Initiator
: Initialization_Vector
);
116 (Gen
: out Generator
;
117 Initiator
: Interfaces
.Integer_32
);
119 (Gen
: out Generator
;
120 Initiator
: Interfaces
.Unsigned_32
);
122 (Gen
: out Generator
;
123 Initiator
: Integer);
124 -- Re-initialize Gen based on the Initiator in various ways. Identical
125 -- values of Initiator cause identical sequences of values.
127 procedure Reset
(Gen
: out Generator
; From_State
: Generator
);
128 -- Causes the state of Gen to be identical to that of From_State; Gen
129 -- and From_State will produce identical sequences of values subsequently.
131 procedure Reset
(Gen
: out Generator
; From_Image
: String);
132 function Image
(Gen
: Generator
) return String;
134 -- Reset (Gen2, Image (Gen1))
135 -- has the same effect as Reset (Gen2, Gen1);
137 Max_Image_Width
: constant :=
138 System
.Random_Numbers
.Max_Image_Width
+ 2 + 20 + 5;
139 -- Maximum possible length of result of Image (...)
143 type Generator
is limited record
144 Rep
: System
.Random_Numbers
.Generator
;
146 Have_Gaussian
: Boolean;
147 -- The algorithm used for Random_Gaussian produces deviates in
148 -- pairs. Have_Gaussian is true iff Random_Gaussian has returned one
149 -- member of the pair and Next_Gaussian contains the other.
151 Next_Gaussian
: Long_Float;
152 -- Next random deviate to be produced by Random_Gaussian, if
156 end GNAT
.Random_Numbers
;