1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188 | with Agpl.Interfaces.C.Types;
package body Agpl.Random is
Internal : Object;
procedure Gsl_Init;
pragma Import (C, Gsl_Init, "agpl__random__gsl_init");
procedure Gsl_Reset;
pragma Import (C, Gsl_Reset, "agpl__random__gsl_reset");
procedure Gsl_Reset_Seed (Seed : Integer);
pragma Import (C, Gsl_Reset_Seed, "agpl__random__gsl_reset_seed");
function Gsl_Gaussian (Sigma : Interfaces.C.Types.Double)
return Interfaces.C.Types.Double;
pragma Import (C, Gsl_Gaussian, "agpl__random__gsl_gaussian");
------------------
-- Open_Uniform --
------------------
function Open_Uniform return Open_Uniformly_Distributed is
begin
return Internal.Open_Uniform;
end Open_Uniform;
-------------
-- Uniform --
-------------
function Uniform return Uniformly_Distributed is
begin
return Internal.Uniform;
end Uniform;
--------------
-- Gaussian --
--------------
function Gaussian (Sigma : Float) return Float is
begin
return Internal.Gaussian (Sigma);
end Gaussian;
----------------------
-- Uniform_Discrete --
----------------------
function Uniform_Discrete return Discrete is
Min : constant Integer := Discrete'Pos (Discrete'First);
Max : constant Integer := Discrete'Pos (Discrete'Last);
begin
return Discrete'Val (Get_Integer (Min, Max));
end Uniform_Discrete;
---------------
-- Get_Float --
---------------
function Get_Float (Min, Max : in Float) return Float is
begin
return Uniform * (Max - Min) + Min;
end Get_Float;
---------------
-- Get_Float --
---------------
function Get_Float (This : in Object;
Min, Max : in Float) return Float
is
begin
return This.Uniform * (Max - Min) + Min;
end Get_Float;
-----------------
-- Get_Integer --
-----------------
function Get_Integer (Min, Max : in Integer) return Integer is
begin
return Internal.Get_Integer (Min, Max);
end Get_Integer;
-----------
-- Reset --
-----------
procedure Reset (Initiator : in Integer) is
begin
Internal.Reset (Initiator);
end Reset;
-----------
-- Reset --
-----------
procedure Reset (This : in out Object) is
begin
Ada.Numerics.Float_Random.Reset (This.Gen);
Gsl_Reset;
end Reset;
-----------
-- Reset --
-----------
procedure Reset (This : in out Object; Initiator : in Integer) is
begin
Ada.Numerics.Float_Random.Reset (This.Gen, Initiator);
Gsl_Reset_Seed (Initiator);
end Reset;
-------------
-- Uniform --
-------------
function Uniform (This : in Object) return Uniformly_Distributed is
begin
return Ada.Numerics.Float_Random.Random (This.Gen);
end Uniform;
--------------
-- Gaussian --
--------------
function Gaussian (This : in Object; Sigma : Float) return Float is
pragma Unreferenced (This);
begin
return Float (Gsl_Gaussian (Interfaces.C.Types.Double (Sigma)));
end Gaussian;
-----------------
-- Get_Integer --
-----------------
function Get_Integer
(This : in Object;
Min, Max : in Integer)
return Integer
is
begin
if Max < Min then
return Max;
else
return Integer'Min
(Integer (Float'Floor (This.Get_Float (Float (Min), Float (Max + 1)))),
Max);
end if;
end Get_Integer;
---------------------
-- Discrete_Random --
---------------------
function Discrete_Random (This : in Object) return Discrete is
pragma Unreferenced (This);
Min : constant Integer := Discrete'Pos (Discrete'First);
Max : constant Integer := Discrete'Pos (Discrete'Last);
begin
return Discrete'Val (Get_Integer (Min, Max));
end Discrete_Random;
------------------
-- Open_Uniform --
------------------
function Open_Uniform (This : in Object) return Open_Uniformly_Distributed is
begin
return Float'Min (Ada.Numerics.Float_Random.Random (This.Gen),
Open_Uniformly_Distributed'Last);
end Open_Uniform;
---------------
-- Flip_Coin --
---------------
function Flip_Coin return Boolean is
begin
return Get_Integer (0, 1) = 1;
end Flip_Coin;
begin
Gsl_Init;
Internal.Reset;
end Agpl.Random;
|