@@ -58,18 +58,18 @@ contains
58
58
59
59
val = simps(y, 1.0_${k1}$)
60
60
ans = 576.0_${k1}$
61
- call check(error, abs( val - ans) < tol_${k1}$)
61
+ call check(error, val, ans, thr= tol_${k1}$)
62
62
if (allocated(error)) return
63
63
64
64
val = simps(y, 0.5_${k1}$)
65
65
ans = 288.0_${k1}$
66
- call check(error, abs( val - ans) < tol_${k1}$)
66
+ call check(error, val, ans, thr= tol_${k1}$)
67
67
if (allocated(error)) return
68
68
69
69
x = [(0.25_${k1}$*(i-1), i = 1, n)]
70
70
val = simps(y, x)
71
71
ans = 144.0_${k1}$
72
- call check(error, abs( val - ans) < tol_${k1}$)
72
+ call check(error, val, ans, thr= tol_${k1}$)
73
73
end subroutine test_simps_${k1}$
74
74
75
75
subroutine test_simps_weights_${k1}$(error)
@@ -91,7 +91,7 @@ contains
91
91
w = simps_weights(x)
92
92
val = sum(w*y)
93
93
ans = simps(y, x)
94
- call check(error, abs( val - ans) < tol_${k1}$)
94
+ call check(error, val, ans, thr= tol_${k1}$)
95
95
end subroutine test_simps_weights_${k1}$
96
96
97
97
subroutine test_simps_zero_${k1}$(error)
@@ -129,18 +129,18 @@ contains
129
129
130
130
val = simps(y, 1.0_${k1}$)
131
131
ans = 1000.0_${k1}$
132
- call check(error, abs( val - ans) < tol_${k1}$)
132
+ call check(error, val, ans, thr= tol_${k1}$)
133
133
if (allocated(error)) return
134
134
135
135
val = simps(y, 0.5_${k1}$)
136
136
ans = 500.0_${k1}$
137
- call check(error, abs( val - ans) < tol_${k1}$)
137
+ call check(error, val, ans, thr= tol_${k1}$)
138
138
if (allocated(error)) return
139
139
140
140
x = [(0.25_${k1}$*(i-1), i = 1, n)]
141
141
val = simps(y, x)
142
142
ans = 250.0_${k1}$
143
- call check(error, abs( val - ans) < tol_${k1}$)
143
+ call check(error, val, ans, thr= tol_${k1}$)
144
144
if (allocated(error)) return
145
145
end do
146
146
end subroutine test_simps_even_${k1}$
@@ -166,7 +166,7 @@ contains
166
166
w = simps_weights(x)
167
167
val = sum(w*y)
168
168
ans = simps(y, x)
169
- call check(error, abs( val - ans) < tol_${k1}$)
169
+ call check(error, val, ans, thr= tol_${k1}$)
170
170
if (allocated(error)) return
171
171
end do
172
172
end subroutine test_simps_weights_even_${k1}$
@@ -190,18 +190,18 @@ contains
190
190
191
191
val = simps(y, 1.0_${k1}$)
192
192
ans = 125.0_${k1}$
193
- call check(error, abs( val - ans) < tol_${k1}$)
193
+ call check(error, val, ans, thr= tol_${k1}$)
194
194
if (allocated(error)) return
195
195
196
196
val = simps(y, 0.5_${k1}$)
197
197
ans = 62.5_${k1}$
198
- call check(error, abs( val - ans) < tol_${k1}$)
198
+ call check(error, val, ans, thr= tol_${k1}$)
199
199
if (allocated(error)) return
200
200
201
201
x = [(0.25_${k1}$*(i-1), i = 1, n)]
202
202
val = simps(y, x)
203
203
ans = 31.25_${k1}$
204
- call check(error, abs( val - ans) < tol_${k1}$)
204
+ call check(error, val, ans, thr= tol_${k1}$)
205
205
if (allocated(error)) return
206
206
end do
207
207
end subroutine test_simps_six_${k1}$
0 commit comments