|
2270
|
1 |
"{ Encoding: utf8 }"
|
|
|
2 |
|
|
1606
|
3 |
"{ Package: 'stx:goodies/regression' }"
|
|
|
4 |
|
|
|
5 |
"{ NameSpace: RegressionTests }"
|
|
|
6 |
|
|
|
7 |
TestCase subclass:#QDoubleTests
|
|
|
8 |
instanceVariableNames:''
|
|
|
9 |
classVariableNames:''
|
|
|
10 |
poolDictionaries:''
|
|
1627
|
11 |
category:'tests-Regression-Numbers'
|
|
1606
|
12 |
!
|
|
|
13 |
|
|
|
14 |
!QDoubleTests class methodsFor:'documentation'!
|
|
|
15 |
|
|
|
16 |
documentation
|
|
|
17 |
"
|
|
|
18 |
documentation to be added.
|
|
|
19 |
|
|
|
20 |
[author:]
|
|
|
21 |
cg
|
|
|
22 |
|
|
|
23 |
[instance variables:]
|
|
|
24 |
|
|
|
25 |
[class variables:]
|
|
|
26 |
|
|
|
27 |
[see also:]
|
|
|
28 |
|
|
|
29 |
"
|
|
|
30 |
! !
|
|
|
31 |
|
|
|
32 |
!QDoubleTests methodsFor:'tests'!
|
|
|
33 |
|
|
1620
|
34 |
test99_misc
|
|
2195
|
35 |
self skip:'unfinished work in progress'.
|
|
|
36 |
|
|
1620
|
37 |
self assert:( 1 asQDouble / (3 factorial) - (QDouble invFact at:1) ) = 0.0.
|
|
|
38 |
self assert:( 1 asQDouble / (4 factorial) - (QDouble invFact at:2) ) = 0.0.
|
|
|
39 |
self assert:( 1 asQDouble / (5 factorial) - (QDouble invFact at:3) ) = 0.0.
|
|
|
40 |
self assert:( 1 asQDouble / (6 factorial) - (QDouble invFact at:4) ) = 0.0.
|
|
|
41 |
self assert:( 1 asQDouble / (7 factorial) - (QDouble invFact at:5) ) = 0.0.
|
|
|
42 |
self assert:( 1 asQDouble / (8 factorial) - (QDouble invFact at:6) ) = 0.0.
|
|
|
43 |
self assert:( 1 asQDouble / (9 factorial) - (QDouble invFact at:7) ) = 0.0.
|
|
|
44 |
self assert:( 1 asQDouble / (10 factorial) - (QDouble invFact at:8) ) = 0.0.
|
|
|
45 |
self assert:( 1 asQDouble / (11 factorial) - (QDouble invFact at:9) ) = 0.0.
|
|
|
46 |
self assert:( 1 asQDouble / (12 factorial) - (QDouble invFact at:10) ) = 0.0.
|
|
|
47 |
self assert:( 1 asQDouble / (13 factorial) - (QDouble invFact at:11) ) = 0.0.
|
|
|
48 |
self assert:( 1 asQDouble / (14 factorial) - (QDouble invFact at:12) ) = 0.0.
|
|
|
49 |
self assert:( 1 asQDouble / (15 factorial) - (QDouble invFact at:13) ) = 0.0.
|
|
|
50 |
self assert:( 1 asQDouble / (16 factorial) - (QDouble invFact at:14) ) = 0.0.
|
|
|
51 |
self assert:( 1 asQDouble / (17 factorial) - (QDouble invFact at:15) ) = 0.0.
|
|
|
52 |
|
|
|
53 |
"Created: / 20-06-2017 / 13:06:04 / cg"
|
|
2195
|
54 |
"Modified: / 28-03-2019 / 12:20:51 / Claus Gittinger"
|
|
1620
|
55 |
!
|
|
|
56 |
|
|
1606
|
57 |
test_01_instance_creation
|
|
1617
|
58 |
|q|
|
|
|
59 |
|
|
|
60 |
q := QDouble fromFloat:1.0.
|
|
|
61 |
self assert:(q d0 = 1.0).
|
|
|
62 |
self assert:(q d1 = 0.0).
|
|
|
63 |
self assert:(q d2 = 0.0).
|
|
|
64 |
self assert:(q d3 = 0.0).
|
|
|
65 |
|
|
|
66 |
q := 1.0 asQDouble.
|
|
|
67 |
self assert:(q d0 = 1.0).
|
|
|
68 |
self assert:(q d1 = 0.0).
|
|
|
69 |
self assert:(q d2 = 0.0).
|
|
|
70 |
self assert:(q d3 = 0.0).
|
|
1606
|
71 |
|
|
1617
|
72 |
q := 1e40 asQDouble + 1.0 asQDouble.
|
|
|
73 |
self assert:(q d0 = 1e40).
|
|
|
74 |
self assert:(q d1 = 1.0).
|
|
|
75 |
self assert:(q d2 = 0.0).
|
|
|
76 |
self assert:(q d3 = 0.0).
|
|
|
77 |
|
|
|
78 |
q := 1e40 asQDouble + 1.0.
|
|
|
79 |
self assert:(q d0 = 1e40).
|
|
|
80 |
self assert:(q d1 = 1.0).
|
|
|
81 |
self assert:(q d2 = 0.0).
|
|
|
82 |
self assert:(q d3 = 0.0).
|
|
|
83 |
|
|
|
84 |
q := 1e120 asQDouble + 1e40 + 1.0.
|
|
|
85 |
self assert:(q d0 = 1e120).
|
|
|
86 |
self assert:(q d1 = 1e40).
|
|
|
87 |
self assert:(q d2 = 1.0).
|
|
|
88 |
self assert:(q d3 = 0.0).
|
|
|
89 |
|
|
|
90 |
q := 1e200 asQDouble + 1e120 + 1e40 + 1.0.
|
|
|
91 |
self assert:(q d0 = 1e200).
|
|
|
92 |
self assert:(q d1 = 1e120).
|
|
|
93 |
self assert:(q d2 = 1e40).
|
|
|
94 |
self assert:(q d3 = 1.0).
|
|
1606
|
95 |
|
|
|
96 |
"
|
|
|
97 |
self run:#test_01_instance_creation
|
|
|
98 |
self new test_01_instance_creation
|
|
|
99 |
"
|
|
|
100 |
|
|
|
101 |
"Created: / 12-06-2017 / 17:03:46 / cg"
|
|
1617
|
102 |
"Modified: / 20-06-2017 / 07:48:51 / cg"
|
|
2269
|
103 |
"Modified: / 05-06-2019 / 20:15:30 / Claus Gittinger"
|
|
1606
|
104 |
!
|
|
|
105 |
|
|
|
106 |
test_02_addition
|
|
|
107 |
|d1 d2 l1 l2 sum_l q1 q2 sum_d sum_q|
|
|
|
108 |
|
|
2269
|
109 |
"/ self skip:'unfinished work in progress'.
|
|
2195
|
110 |
|
|
1606
|
111 |
d1 := 1.0.
|
|
|
112 |
d2 := 1.0e-16.
|
|
|
113 |
sum_d := d1 + d2.
|
|
|
114 |
"/ precision is lost with doubles!!
|
|
|
115 |
self assert:( sum_d - 1.0 = 0.0 ).
|
|
|
116 |
|
|
|
117 |
l1 := 1.0 asLongFloat.
|
|
|
118 |
l2 := 1.0e-16 asLongFloat.
|
|
|
119 |
sum_l := l1 + l2.
|
|
|
120 |
"/ no precision is lost with long doubles!!
|
|
|
121 |
self assert:( (sum_l - (1.0 asLongFloat)) > 0.0 ).
|
|
|
122 |
|
|
1607
|
123 |
l1 := 1.0 asLongFloat.
|
|
|
124 |
l2 := 1.0e-30 asLongFloat.
|
|
|
125 |
sum_l := l1 + l2.
|
|
|
126 |
"/ precision is lost!!
|
|
|
127 |
self assert:( (sum_l - (1.0 asLongFloat)) = 0.0 ).
|
|
|
128 |
|
|
1606
|
129 |
q1 := QDouble fromFloat:1.0.
|
|
|
130 |
q2 := QDouble fromFloat:1.0e-16.
|
|
|
131 |
sum_q := q1 + q2.
|
|
|
132 |
|
|
1607
|
133 |
"/ no precision is lost
|
|
|
134 |
self assert:( (sum_q - (QDouble fromFloat:1.0)) > 0.0 ).
|
|
1614
|
135 |
|
|
|
136 |
q1 := QDouble fromFloat:1.0.
|
|
|
137 |
q2 := QDouble fromFloat:2.0.
|
|
|
138 |
sum_q := q1 + q2.
|
|
|
139 |
self assert:( sum_q = 3.0 ).
|
|
1607
|
140 |
|
|
1614
|
141 |
q1 := QDouble fromFloat:1.0.
|
|
|
142 |
q2 := QDouble fromFloat:2.0.
|
|
|
143 |
sum_q := q1 + 2.0.
|
|
|
144 |
self assert:( sum_q = 3.0 ).
|
|
|
145 |
|
|
|
146 |
q1 := QDouble fromFloat:1.0.
|
|
|
147 |
q2 := QDouble fromFloat:2.0.
|
|
|
148 |
sum_q := 1.0 + q2.
|
|
|
149 |
self assert:( sum_q = 3.0 ).
|
|
|
150 |
|
|
1617
|
151 |
q1 := 1e200 asQDouble + 1e120 + 1e40 + 1.0.
|
|
|
152 |
q2 := 1e120.
|
|
|
153 |
sum_q := q1 + q2.
|
|
|
154 |
self assert:(sum_q d0 = 1e200).
|
|
|
155 |
self assert:(sum_q d1 = 2e120).
|
|
|
156 |
self assert:(sum_q d2 = 1e40).
|
|
|
157 |
self assert:(sum_q d3 = 1.0).
|
|
|
158 |
|
|
1606
|
159 |
"
|
|
|
160 |
self run:#test_02_addition
|
|
|
161 |
self new test_02_addition
|
|
|
162 |
"
|
|
|
163 |
|
|
|
164 |
"Created: / 12-06-2017 / 17:05:07 / cg"
|
|
1617
|
165 |
"Modified: / 20-06-2017 / 07:49:32 / cg"
|
|
2269
|
166 |
"Modified: / 05-06-2019 / 20:15:38 / Claus Gittinger"
|
|
1606
|
167 |
!
|
|
|
168 |
|
|
|
169 |
test_03_subtract
|
|
|
170 |
|d1 d2 l1 l2 diff_l q1 q2 diff_d diff_q|
|
|
|
171 |
|
|
2269
|
172 |
"/ self skip:'unfinished work in progress'.
|
|
2195
|
173 |
|
|
1606
|
174 |
d1 := 1.0.
|
|
|
175 |
d2 := 0.99999.
|
|
|
176 |
diff_d := d1 - d2.
|
|
|
177 |
"/ precision is lost with doubles!!
|
|
|
178 |
|
|
|
179 |
l1 := 1.0 asLongFloat.
|
|
|
180 |
l2 := 0.99999 asLongFloat.
|
|
|
181 |
diff_l := l1 - l2.
|
|
|
182 |
"/ no precision is lost with long doubles!!
|
|
|
183 |
|
|
|
184 |
q1 := QDouble fromFloat:1.0.
|
|
|
185 |
q2 := QDouble fromFloat:0.99999.
|
|
|
186 |
|
|
|
187 |
diff_q := q1 - q2.
|
|
|
188 |
|
|
|
189 |
"
|
|
|
190 |
self run:#test_03_subtract
|
|
|
191 |
self new test_03_subtract
|
|
|
192 |
"
|
|
|
193 |
|
|
|
194 |
"Created: / 12-06-2017 / 23:12:00 / cg"
|
|
2269
|
195 |
"Modified: / 05-06-2019 / 20:16:12 / Claus Gittinger"
|
|
1608
|
196 |
!
|
|
|
197 |
|
|
|
198 |
test_04_relops
|
|
2269
|
199 |
"/ self skip:'unfinished work in progress'.
|
|
2195
|
200 |
|
|
1608
|
201 |
self assert:(1.0 = 1.0 asQDouble).
|
|
|
202 |
self assert:(1.0 asQDouble = 1.0).
|
|
|
203 |
self assert:(1.0 asQDouble = 1.0 asQDouble).
|
|
|
204 |
|
|
|
205 |
self assert:(1 = 1 asQDouble).
|
|
|
206 |
self assert:(1 asQDouble = 1).
|
|
|
207 |
self assert:(1 asQDouble = 1 asQDouble).
|
|
|
208 |
|
|
|
209 |
self assert:(1 = 1.0 asQDouble).
|
|
|
210 |
self assert:(1.0 asQDouble = 1).
|
|
|
211 |
self assert:(1.0 asQDouble = 1 asQDouble).
|
|
|
212 |
|
|
|
213 |
self assert:(1.0 = 1 asQDouble).
|
|
|
214 |
self assert:(1 asQDouble = 1.0).
|
|
|
215 |
self assert:(1 asQDouble = 1.0 asQDouble).
|
|
|
216 |
|
|
|
217 |
|
|
|
218 |
self deny:(1.0 = 1.1 asQDouble).
|
|
|
219 |
self deny:(1.0 asQDouble = 1.1).
|
|
|
220 |
self deny:(1.0 asQDouble = 1.1 asQDouble).
|
|
|
221 |
|
|
|
222 |
self deny:(1 = 2 asQDouble).
|
|
|
223 |
self deny:(1 asQDouble = 2).
|
|
|
224 |
self deny:(1 asQDouble = 2 asQDouble).
|
|
|
225 |
|
|
|
226 |
self deny:(1 = 2.0 asQDouble).
|
|
|
227 |
self deny:(1.0 asQDouble = 2).
|
|
|
228 |
self deny:(1.0 asQDouble = 2 asQDouble).
|
|
|
229 |
|
|
|
230 |
self deny:(2.0 = 1 asQDouble).
|
|
|
231 |
self deny:(2 asQDouble = 1.0).
|
|
|
232 |
self deny:(2 asQDouble = 1.0 asQDouble).
|
|
|
233 |
|
|
|
234 |
|
|
|
235 |
|
|
|
236 |
self assert:(1.0 < 1.1 asQDouble).
|
|
|
237 |
self assert:(1.0 asQDouble < 1.1).
|
|
|
238 |
self assert:(1.0 asQDouble < 1.1 asQDouble).
|
|
|
239 |
|
|
|
240 |
self assert:(1 < 2 asQDouble).
|
|
|
241 |
self assert:(1 asQDouble < 2).
|
|
|
242 |
self assert:(1 asQDouble < 2 asQDouble).
|
|
|
243 |
|
|
|
244 |
self assert:(1 < 2.0 asQDouble).
|
|
|
245 |
self assert:(1.0 asQDouble < 2).
|
|
|
246 |
self assert:(1.0 asQDouble < 2 asQDouble).
|
|
|
247 |
|
|
|
248 |
self assert:(1.0 < 2 asQDouble).
|
|
|
249 |
self assert:(1 asQDouble < 2.0).
|
|
|
250 |
self assert:(1 asQDouble < 2.0 asQDouble).
|
|
|
251 |
|
|
|
252 |
|
|
|
253 |
self deny:(1.0 < 1.0 asQDouble).
|
|
|
254 |
self deny:(1.0 asQDouble < 1.0).
|
|
|
255 |
self deny:(1.0 asQDouble < 1.0 asQDouble).
|
|
|
256 |
|
|
|
257 |
self deny:(1 < 1 asQDouble).
|
|
|
258 |
self deny:(1 asQDouble < 1).
|
|
|
259 |
self deny:(1 asQDouble < 1 asQDouble).
|
|
|
260 |
|
|
|
261 |
self deny:(1 < 1.0 asQDouble).
|
|
|
262 |
self deny:(1.0 asQDouble < 1).
|
|
|
263 |
self deny:(1.0 asQDouble < 1 asQDouble).
|
|
|
264 |
|
|
|
265 |
self deny:(1.0 < 1 asQDouble).
|
|
|
266 |
self deny:(1 asQDouble < 1.0).
|
|
|
267 |
self deny:(1 asQDouble < 1.0 asQDouble).
|
|
|
268 |
|
|
|
269 |
|
|
|
270 |
|
|
|
271 |
self assert:(2.0 > 1.1 asQDouble).
|
|
|
272 |
self assert:(2.0 asQDouble > 1.1).
|
|
|
273 |
self assert:(2.0 asQDouble > 1.1 asQDouble).
|
|
|
274 |
|
|
|
275 |
self assert:(3 > 2 asQDouble).
|
|
|
276 |
self assert:(3 asQDouble > 2).
|
|
|
277 |
self assert:(3 asQDouble > 2 asQDouble).
|
|
|
278 |
|
|
|
279 |
self assert:(3 > 2.0 asQDouble).
|
|
|
280 |
self assert:(3.0 asQDouble > 2).
|
|
|
281 |
self assert:(3.0 asQDouble > 2 asQDouble).
|
|
|
282 |
|
|
|
283 |
self assert:(3.0 > 2 asQDouble).
|
|
|
284 |
self assert:(3 asQDouble > 2.0).
|
|
|
285 |
self assert:(3 asQDouble > 2.0 asQDouble).
|
|
|
286 |
|
|
|
287 |
|
|
|
288 |
self deny:(1.0 > 3.0 asQDouble).
|
|
|
289 |
self deny:(1.0 asQDouble > 3.0).
|
|
|
290 |
self deny:(1.0 asQDouble > 3.0 asQDouble).
|
|
|
291 |
|
|
|
292 |
self deny:(1 > 3 asQDouble).
|
|
|
293 |
self deny:(1 asQDouble > 3).
|
|
|
294 |
self deny:(1 asQDouble > 3 asQDouble).
|
|
|
295 |
|
|
|
296 |
self deny:(1 > 3.0 asQDouble).
|
|
|
297 |
self deny:(1.0 asQDouble > 3).
|
|
|
298 |
self deny:(1.0 asQDouble > 3 asQDouble).
|
|
|
299 |
|
|
|
300 |
self deny:(1.0 > 3 asQDouble).
|
|
|
301 |
self deny:(1 asQDouble > 3.0).
|
|
|
302 |
self deny:(1 asQDouble > 3.0 asQDouble).
|
|
|
303 |
|
|
|
304 |
|
|
|
305 |
self assert:(1.0 <= 3.0 asQDouble).
|
|
|
306 |
self assert:(1.0 asQDouble <= 3.0).
|
|
|
307 |
self assert:(1.0 asQDouble <= 3.0 asQDouble).
|
|
|
308 |
|
|
|
309 |
self assert:(1.0 <= 1.0 asQDouble).
|
|
|
310 |
self assert:(1.0 asQDouble <= 1.0).
|
|
|
311 |
self assert:(1.0 asQDouble <= 1.0 asQDouble).
|
|
|
312 |
|
|
|
313 |
self assert:(2.0 >= 1.0 asQDouble).
|
|
|
314 |
self assert:(2.0 asQDouble >= 1.0).
|
|
|
315 |
self assert:(2.0 asQDouble >= 1.0 asQDouble).
|
|
|
316 |
|
|
|
317 |
self assert:(1.0 >= 1.0 asQDouble).
|
|
|
318 |
self assert:(1.0 asQDouble >= 1.0).
|
|
|
319 |
self assert:(1.0 asQDouble >= 1.0 asQDouble).
|
|
|
320 |
|
|
1617
|
321 |
|
|
|
322 |
|
|
|
323 |
self assert: ((1.0 asQDouble + 1e-40) > 1.0 asQDouble).
|
|
|
324 |
self assert: (1.0 asQDouble < (1.0 asQDouble + 1e-40)).
|
|
|
325 |
self assert: (1.0 asQDouble < (1.0 asQDouble + 1e-80)).
|
|
|
326 |
self assert: (1.0 asQDouble < (1.0 asQDouble + 1e-120)).
|
|
|
327 |
self assert: (1.0 asQDouble < (1.0 asQDouble + 1e-300)).
|
|
|
328 |
self assert: (1.0 asQDouble < (1.0 asQDouble + 1e-300)).
|
|
|
329 |
|
|
|
330 |
self assert: ((1.0 asQDouble + 1e-40) > 1.0).
|
|
|
331 |
self assert: (1.0 < (1.0 asQDouble + 1e-40)).
|
|
|
332 |
self assert: (1.0 < (1.0 asQDouble + 1e-80)).
|
|
|
333 |
self assert: (1.0 < (1.0 asQDouble + 1e-120)).
|
|
|
334 |
self assert: (1.0 < (1.0 asQDouble + 1e-300)).
|
|
|
335 |
self assert: (1.0 < (1.0 asQDouble + 1e-300)).
|
|
|
336 |
|
|
1619
|
337 |
self assert: (1e10 asQDouble + 1.0) = 10000000001.0.
|
|
|
338 |
"/ that's the problem with Floats...
|
|
|
339 |
"/ (1e20 asQDouble + 1.0) ~= 100000000000000000001.0.
|
|
|
340 |
"/ (1e20 + 1.0) - 1e20 ~= 1.0.
|
|
|
341 |
self assert: (1e20 asQDouble + 1.0) - 1e20 = 1.0.
|
|
1617
|
342 |
|
|
|
343 |
"/ here the precision of QDOubles is no longer enough:
|
|
|
344 |
"/ self assert: 1.0 asQDouble < (1.0 asQDouble + 1e-310).
|
|
|
345 |
|
|
1608
|
346 |
"
|
|
|
347 |
self run:#test_03_subtract
|
|
|
348 |
self new test_03_subtract
|
|
|
349 |
"
|
|
|
350 |
|
|
|
351 |
"Created: / 13-06-2017 / 17:11:37 / cg"
|
|
1619
|
352 |
"Modified (comment): / 20-06-2017 / 09:18:50 / cg"
|
|
2269
|
353 |
"Modified: / 05-06-2019 / 20:16:21 / Claus Gittinger"
|
|
1614
|
354 |
!
|
|
|
355 |
|
|
|
356 |
test_05_multiply
|
|
|
357 |
|q1 q2 p|
|
|
|
358 |
|
|
2269
|
359 |
"/ self skip:'unfinished work in progress'.
|
|
2195
|
360 |
|
|
1614
|
361 |
q1 := QDouble fromFloat:3.0. "/ DoubleArray(2.0 0.0 0.0 0.0)
|
|
|
362 |
p := 2.0 * q1.
|
|
|
363 |
|
|
|
364 |
self assert:(p d0 = 6.0).
|
|
|
365 |
self assert:(p d1 = 0.0).
|
|
|
366 |
self assert:(p d2 = 0.0).
|
|
|
367 |
self assert:(p d3 = 0.0).
|
|
|
368 |
self assert:(p = 6.0).
|
|
|
369 |
self assert:(6.0 = p).
|
|
|
370 |
|
|
|
371 |
q1 := QDouble fromFloat:2.0. "/ DoubleArray(2.0 0.0 0.0 0.0)
|
|
|
372 |
p := q1 * 3.0.
|
|
|
373 |
|
|
|
374 |
self assert:(p d0 = 6.0).
|
|
|
375 |
self assert:(p d1 = 0.0).
|
|
|
376 |
self assert:(p d2 = 0.0).
|
|
|
377 |
self assert:(p d3 = 0.0).
|
|
|
378 |
self assert:(p = 6.0).
|
|
|
379 |
self assert:(6.0 = p).
|
|
|
380 |
|
|
|
381 |
q1 := QDouble fromFloat:2.0. "/ DoubleArray(2.0 0.0 0.0 0.0)
|
|
|
382 |
q2 := QDouble fromFloat:3.0. "/ DoubleArray(3.0 0.0 0.0 0.0)
|
|
|
383 |
p := q1 * q2.
|
|
|
384 |
|
|
|
385 |
self assert:(p d0 = 6.0).
|
|
|
386 |
self assert:(p d1 = 0.0).
|
|
|
387 |
self assert:(p d2 = 0.0).
|
|
|
388 |
self assert:(p d3 = 0.0).
|
|
|
389 |
self assert:(p = 6.0).
|
|
|
390 |
self assert:(6.0 = p).
|
|
1617
|
391 |
|
|
|
392 |
"/ the 1 is not forgotten!!
|
|
|
393 |
q1 := 1 asQDouble + 1e40.
|
|
|
394 |
q2 := 2.0.
|
|
|
395 |
p := q1 * q2.
|
|
|
396 |
|
|
|
397 |
self assert:(p d0 = 2e40).
|
|
|
398 |
self assert:(p d1 = 2.0).
|
|
|
399 |
self assert:(p d2 = 0.0).
|
|
|
400 |
self assert:(p d3 = 0.0).
|
|
|
401 |
self assert:(p > 2e40).
|
|
1751
|
402 |
self assert:(p - 2e40) = 2.0.
|
|
1617
|
403 |
|
|
|
404 |
"/ nothing is forgotten!!
|
|
|
405 |
q1 := 1 asQDouble + 1e40 + 1e120.
|
|
|
406 |
q2 := 2.0.
|
|
|
407 |
p := q1 * q2.
|
|
|
408 |
|
|
|
409 |
self assert:(p d0 = 2e120).
|
|
|
410 |
self assert:(p d1 = 2e40).
|
|
|
411 |
self assert:(p d2 = 2.0).
|
|
|
412 |
self assert:(p d3 = 0.0).
|
|
|
413 |
self assert:(p > 2e120).
|
|
|
414 |
self assert:(p > (2e120 asQDouble + 2)).
|
|
|
415 |
self assert:(p > (2e120 asQDouble + 1e40)).
|
|
1751
|
416 |
self assert:(p - 2e120) = (2 asQDouble + 2e40).
|
|
|
417 |
self assert:(p - 2e120 - 2e40) = 2.0.
|
|
1617
|
418 |
|
|
1614
|
419 |
"
|
|
|
420 |
self run:#test_05_multiply
|
|
|
421 |
self new test_05_multiply
|
|
|
422 |
"
|
|
|
423 |
|
|
|
424 |
"Created: / 19-06-2017 / 16:58:29 / cg"
|
|
1617
|
425 |
"Modified (comment): / 20-06-2017 / 07:42:52 / cg"
|
|
2269
|
426 |
"Modified: / 05-06-2019 / 20:16:30 / Claus Gittinger"
|
|
1616
|
427 |
!
|
|
|
428 |
|
|
|
429 |
test_06_exp
|
|
|
430 |
|f1 q1|
|
|
|
431 |
|
|
2269
|
432 |
"/ self skip:'unfinished work in progress'.
|
|
2195
|
433 |
|
|
1623
|
434 |
f1 := 2.0 exp.
|
|
|
435 |
q1 := 2.0 asQDouble exp.
|
|
1616
|
436 |
|
|
1623
|
437 |
"/ the first 113 digits are:
|
|
|
438 |
"/ 7.3890560989306502272304274605750078131803155705518473240871278225225737960790577633843124850791217947737531612654...
|
|
|
439 |
"/ 2.0 exp
|
|
|
440 |
"/ '%30.28f' printf:{ 2.0 exp } -> '7.38905609893065040694182243896'
|
|
|
441 |
"/ '%30.28f' printf:{ 2.0 exp asQDouble } -> '7.38905609893065040694182243897'
|
|
|
442 |
"/ '%30.28f' printf:{ 2.0 asQDouble exp } -> '7.38905609893065022723042746058'
|
|
|
443 |
"/ '%40.38f' printf:{ 2.0 asQDouble exp } -> '7.389056098930650227230427460575007822092' (10 iters in exp)
|
|
|
444 |
"/ '%40.38f' printf:{ 2.0 asQDouble exp }
|
|
|
445 |
"/ 2.0 asQDouble exp -> 7.3890560989306502272304274605750078* 220914983921672'
|
|
|
446 |
"/ 1.0 asQDouble exp
|
|
1616
|
447 |
self assert:(q1 d0 = f1).
|
|
|
448 |
|
|
|
449 |
"
|
|
|
450 |
self run:#test_06_exp
|
|
|
451 |
self new test_06_exp
|
|
|
452 |
"
|
|
|
453 |
|
|
|
454 |
"Created: / 19-06-2017 / 17:58:08 / cg"
|
|
1623
|
455 |
"Modified (comment): / 21-06-2017 / 13:45:56 / cg"
|
|
2269
|
456 |
"Modified: / 05-06-2019 / 20:16:35 / Claus Gittinger"
|
|
1617
|
457 |
!
|
|
|
458 |
|
|
1627
|
459 |
test_06b_log10
|
|
|
460 |
|q1 q2 d3a q3a d3b q3b q4 i|
|
|
|
461 |
|
|
2270
|
462 |
"/ self skip:'unfinished work in progress'.
|
|
2195
|
463 |
|
|
1627
|
464 |
q1 := QDouble d0: 5.55111512312578E-18
|
|
|
465 |
d1: -3.08148791101958E-34
|
|
|
466 |
d2: 1.71056941445901E-50
|
|
|
467 |
d3: 0.0.
|
|
|
468 |
|
|
|
469 |
q2 := 10 / q1.
|
|
|
470 |
self assert:(q2 > 1e18).
|
|
|
471 |
|
|
|
472 |
d3a := q2 asFloat ln.
|
|
|
473 |
self assert:(d3a isNaN not).
|
|
|
474 |
|
|
|
475 |
q3a := q2 ln.
|
|
|
476 |
self assert:(q3a isNaN not).
|
|
|
477 |
|
|
|
478 |
d3b := q2 asFloat log10.
|
|
|
479 |
self assert:(d3b isNaN not).
|
|
|
480 |
|
|
|
481 |
q3b := q2 log10.
|
|
|
482 |
self assert:(q3b isNaN not).
|
|
|
483 |
|
|
|
484 |
q4 := q3b floor.
|
|
|
485 |
i := q4 asInteger.
|
|
|
486 |
|
|
|
487 |
|
|
|
488 |
"
|
|
|
489 |
self run:#test_06b_log10
|
|
|
490 |
self new test_06b_log10
|
|
|
491 |
"
|
|
|
492 |
|
|
|
493 |
"Created: / 03-07-2017 / 10:31:53 / cg"
|
|
2270
|
494 |
"Modified: / 05-06-2019 / 20:18:26 / Claus Gittinger"
|
|
1627
|
495 |
!
|
|
|
496 |
|
|
1617
|
497 |
test_07_reading
|
|
|
498 |
|q1 s|
|
|
|
499 |
|
|
2270
|
500 |
"/ self skip:'unfinished work in progress'.
|
|
2195
|
501 |
|
|
1617
|
502 |
q1 := QDouble readFrom: '3.1415926535897'.
|
|
|
503 |
s := '%55.54f' printf:{ q1 }.
|
|
|
504 |
self assert:(s = ' 3.1415926535897').
|
|
|
505 |
|
|
|
506 |
q1 := QDouble readFrom: '3.141592653589793238462643383'.
|
|
|
507 |
s := '%55.54f' printf:{ q1 }.
|
|
|
508 |
self assert:(s = ' 3.141592653589793238462643383').
|
|
|
509 |
|
|
|
510 |
q1 := QDouble readFrom: '3.141592653589793238462643383279502884197169399375105820'.
|
|
|
511 |
s := '%55.54f' printf:{ q1 }.
|
|
|
512 |
self assert:(s = '3.14159265358979323846264338327950288419716939937510582').
|
|
|
513 |
|
|
|
514 |
"
|
|
|
515 |
self run:#test_07_reading
|
|
|
516 |
self new test_07_reading
|
|
|
517 |
"
|
|
|
518 |
|
|
|
519 |
"Created: / 19-06-2017 / 23:29:34 / cg"
|
|
2270
|
520 |
"Modified: / 05-06-2019 / 20:18:51 / Claus Gittinger"
|
|
1617
|
521 |
!
|
|
|
522 |
|
|
|
523 |
test_08_conversion
|
|
|
524 |
|q i|
|
|
|
525 |
|
|
2270
|
526 |
"/ self skip:'unfinished work in progress'.
|
|
2195
|
527 |
|
|
1617
|
528 |
q := 1e10 asQDouble + 1.0.
|
|
|
529 |
i := q asInteger.
|
|
|
530 |
|
|
|
531 |
self assert:(i = 10000000001).
|
|
|
532 |
|
|
|
533 |
"
|
|
|
534 |
self run:#test_08_conversion
|
|
|
535 |
self new test_08_conversion
|
|
|
536 |
"
|
|
|
537 |
|
|
|
538 |
"Created: / 20-06-2017 / 07:51:50 / cg"
|
|
2270
|
539 |
"Modified: / 05-06-2019 / 20:18:48 / Claus Gittinger"
|
|
1606
|
540 |
! !
|
|
|
541 |
|
|
|
542 |
!QDoubleTests class methodsFor:'documentation'!
|
|
|
543 |
|
|
|
544 |
version
|
|
|
545 |
^ '$Header$'
|
|
|
546 |
!
|
|
|
547 |
|
|
|
548 |
version_CVS
|
|
|
549 |
^ '$Header$'
|
|
|
550 |
! !
|
|
|
551 |
|