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
|
# Tests for the module zsh/mathfunc
%prep
if ! zmodload zsh/mathfunc 2>/dev/null; then
ZTST_unimplemented="The module zsh/mathfunc is not available."
fi
%test
# -g makes pi available in later tests
float -gF 5 pi
(( pi = 4 * atan(1.0) ))
print $pi
0:Basic operation with atan
>3.14159
float -F 5 result
(( result = atan(3,2) ))
print $result
0:atan with two arguments
>0.98279
print $(( atan(1,2,3) ))
1:atan can't take three arguments
?(eval):1: wrong number of arguments: atan(1,2,3)
float r1=$(( rand48() ))
float r2=$(( rand48() ))
float r3=$(( rand48() ))
# Yes, this is a floating point equality test like they tell
# you not to do. As the pseudrandom sequence is deterministic,
# this is the right thing to do in this case.
if (( r1 == r2 )); then
print "Seed not updated correctly the first time"
else
print "First two random numbers differ, OK"
fi
if (( r2 == r3 )); then
print "Seed not updated correctly the second time"
else
print "Second two random numbers differ, OK"
fi
0:rand48 with default initialisation
F:This test fails if your math library doesn't have erand48().
>First two random numbers differ, OK
>Second two random numbers differ, OK
seed=f45677a6cbe4
float r1=$(( rand48(seed) ))
float r2=$(( rand48(seed) ))
seed2=$seed
float r3=$(( rand48(seed) ))
float r4=$(( rand48(seed2) ))
# Yes, this is a floating point equality test like they tell
# you not to do. As the pseudrandom sequence is deterministic,
# this is the right thing to do in this case.
if (( r1 == r2 )); then
print "Seed not updated correctly the first time"
else
print "First two random numbers differ, OK"
fi
if (( r2 == r3 )); then
print "Seed not updated correctly the second time"
else
print "Second two random numbers differ, OK"
fi
if (( r3 == r4 )); then
print "Identical seeds generate identical numbers, OK"
else
print "Indeterminate result from identical seeds"
fi
0:rand48 with pre-generated seed
F:This test fails if your math library doesn't have erand48().
>First two random numbers differ, OK
>Second two random numbers differ, OK
>Identical seeds generate identical numbers, OK
float -F 5 pitest
(( pitest = 4.0 * atan(1) ))
# This is a string test of the output to 5 digits.
if [[ $pi = $pitest ]]; then
print "OK, atan on an integer seemed to work"
else
print "BAD: got $pitest instead of $pi"
fi
0:Conversion of arguments from integer
>OK, atan on an integer seemed to work
float -F 5 result
typeset str
for str in 0 0.0 1 1.5 -1 -1.5; do
(( result = abs($str) ))
print $result
done
0:Use of abs on various numbers
>0.00000
>0.00000
>1.00000
>1.50000
>1.00000
>1.50000
print $(( sqrt(-1) ))
0:Non-negative argument checking for square roots.
>NaN
# Simple test that the pseudorandom number generators are producing
# something that could conceivably be pseudorandom numbers in a
# linear range. Not a detailed quantitative verification.
integer N=10000 isource ok=1
float -F f sum sumsq max max2 av sd
typeset -a randoms
randoms=('f = RANDOM' 'f = rand48()')
for isource in 1 2; do
(( sum = sumsq = max = 0 ))
repeat $N; do
let $randoms[$isource]
(( f > max )) && (( max = f ))
(( sum += f, sumsq += f * f ))
done
(( av = sum / N ))
(( sd = sqrt((sumsq - N * av * av) / (N-1)) ))
(( max2 = 0.5 * max ))
if (( av > max2 * 1.1 )) || (( av < max2 * 0.9 )); then
print "WARNING: average of random numbers is suspicious.
Was testing: $randoms[$isource]"
(( ok = 0 ))
fi
if (( sd < max / 4 )); then
print "WARNING: distribution of random numbers is suspicious.
Was testing: $randoms[$isource]"
(( ok = 0 ))
fi
done
(( ok ))
0:Test random number generator distributions are not grossly broken
float -F 5 g l
(( g = gamma(2), l = lgamma(2) ))
print $g, $l
0:Test Gamma function gamma and lgamma
>1.00000, 0.00000
|
