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
|
#!/usr/bin/env ruby
require("gsl")
include GSL
module GSL
module Test
Verbose = true
$tests = 0
$passed = 0
$failed = 0
def test(status, desc)
$tests += 1
if !status
$passed += 1
printf("PASS: #{desc}\n")
else
$failed += 1
printf("FAIL: #{desc}\n")
end
end
def test_factor(result, expected, factor, desc)
status = nil
if result == expected
status = false
elsif expected == 0.0
status = (result > expected or result < expected)
else
u = result/expected
status = (u > factor or u < 1.0/factor)
end
$tests += 1
if !status
$passed += 1
printf("PASS: #{desc} (%g observed vs %g expected)\n", result, expected)
else
$failed += 1
printf("FAIL: #{desc} (%.18g observed vs %.18g expected)\n", result, expected)
end
end
def test_factor2(result, expected, factor, desc)
status = nil
if result == expected
status = false
elsif expected == 0.0
status = (result > expected or result < expected)
else
u = result/expected
status = (u > factor or u < 1.0/factor)
end
$tests += 1
if !status
$passed += 1
printf("PASS: #{desc} (%.18g observed vs %g expected)\n", result, expected)
else
$failed += 1
printf("FAIL: #{desc} (%.18g observed vs %.18g expected)\n", result, expected)
end
end
def test_rel(result, expected, relerr, desc)
status = nil
if isnan?(result) or isnan?(expected)
status = isnan?(result) != isnan?(expected)
elsif isinf?(result) or isinf?(expected)
status = isinf?(result) != isinf?(expected)
elsif expected.to_f != 0.0
status = (result - expected).abs/expected.abs > relerr
else
status = result.abs > relerr
end
$tests += 1
if !status
$passed += 1
printf("PASS: #{desc} (%.18g observed vs %g expected)\n", result, expected)
else
$failed += 1
printf("FAIL: #{desc} (%.18g observed vs %.18g expected)\n", result, expected)
end
end
def test_abs(result, expected, abserr, desc)
status = nil
if isnan?(result) or isnan?(expected)
status = isnan?(result) != isnan?(expected)
elsif isinf?(result) or isinf?(expected)
status = isinf?(result) != isinf?(expected)
else
status = (result - expected).abs > abserr
end
$tests += 1
if !status
$passed += 1
# printf("PASS: #{desc} (%g observed vs %g expected)\n", result, expected)
printf("PASS: #{desc} (%.18g observed vs %g expected)\n", result, expected)
else
$failed += 1
printf("FAIL: #{desc} (%.18g observed vs %.18g expected)\n", result, expected)
end
end
def test_int(result, expected, desc)
status = (result != expected)
$tests += 1
if !status
$passed += 1
printf("PASS: #{desc} (%d observed vs %d expected)\n", result, expected)
else
$failed += 1
printf("FAIL: #{desc} (%d observed vs %d expected)\n", result, expected)
end
end
end
end
|
