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
|
#! /usr/bin/env python
import openturns as ot
ot.TESTPREAMBLE()
# Test the initialization method
ot.RandomGenerator.SetSeed(0)
# Test the state accessor for reading
initialState = ot.RandomGenerator.GetState()
print("initialState=", repr(initialState))
# Test the generation of one realization
realization = ot.RandomGenerator.Generate()
print("one realization=%.6f" % realization)
state = ot.RandomGenerator.GetState()
print("state=", repr(state))
print("new state index=", state.getIndex())
# Test the generation of several realizations
size = 10000
several = ot.RandomGenerator.Generate(size)
print(size, " realizations, first=%.6f" % several[0], ", last=%.6f" % several[size - 1])
# Test the state accessor for writing
ot.RandomGenerator.SetState(initialState)
print("restoring initial state, one realization=%.6f" % ot.RandomGenerator.Generate())
# Test the uniformity of the random generator
size = 100000
manyPoints = ot.RandomGenerator.Generate(size)
slice = 10
frequencies = ot.Point(slice, 0.0)
mean = 0.0
for i in range(size):
frequencies[int(slice * manyPoints[i])] = (
frequencies[int(slice * manyPoints[i])] + 1
)
mean = mean + manyPoints[i]
mean = mean / size
print("mean of", size, "realizations= %.12f" % mean)
for i in range(slice):
frequencies[i] = frequencies[i] / size
print(
"frequency in slice [",
(1.0 * i) / slice,
",",
(i + 1.0) / slice,
"]=",
frequencies[i],
)
print(
"One integer generation in [[0, ",
slice - 1,
"]]=",
ot.RandomGenerator.IntegerGenerate(slice),
)
print(
"Another integer generation in [[0, ",
slice - 1,
"]]=",
ot.RandomGenerator.IntegerGenerate(slice),
)
manyUnsignedInteger = ot.RandomGenerator.IntegerGenerate(size, slice)
frequencies = ot.Point(slice, 0.0)
mean = 0.0
for i in range(size):
frequencies[manyUnsignedInteger[i]] = frequencies[manyUnsignedInteger[i]] + 1
mean = mean + manyUnsignedInteger[i]
mean = mean / size
print("mean of ", size, " integer realizations=", mean)
for i in range(slice):
frequencies[i] = frequencies[i] / size
print("frequency for value ", i, "=", frequencies[i])
|
