blob: a066d2e15e2d0a1bfe93d6193c88cbc0c21b18be (
plain)
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
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
|
;;; floatfns-tests.el --- tests for floating point operations -*- lexical-binding: t -*-
;; Copyright 2017-2021 Free Software Foundation, Inc.
;; This file is part of GNU Emacs.
;; GNU Emacs is free software: you can redistribute it and/or modify
;; it under the terms of the GNU General Public License as published by
;; the Free Software Foundation, either version 3 of the License, or
;; (at your option) any later version.
;; GNU Emacs is distributed in the hope that it will be useful,
;; but WITHOUT ANY WARRANTY; without even the implied warranty of
;; MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
;; GNU General Public License for more details.
;; You should have received a copy of the GNU General Public License
;; along with GNU Emacs. If not, see <https://www.gnu.org/licenses/>.
;;; Code:
(require 'ert)
(ert-deftest floatfns-tests-cos ()
(should (= (cos 0) 1.0))
(should (= (cos float-pi) -1.0)))
(ert-deftest floatfns-tests-sin ()
(should (= (sin 0) 0.0)))
(ert-deftest floatfns-tests-tan ()
(should (= (tan 0) 0.0)))
(ert-deftest floatfns-tests-isnan ()
(should (isnan 0.0e+NaN))
(should (isnan -0.0e+NaN))
(should-error (isnan "foo") :type 'wrong-type-argument))
(ert-deftest floatfns-tests-exp ()
(should (= (exp 0) 1.0)))
(ert-deftest floatfns-tests-expt ()
(should (= (expt 2 8) 256)))
(ert-deftest floatfns-tests-log ()
(should (= (log 1000 10) 3.0)))
(ert-deftest floatfns-tests-sqrt ()
(should (= (sqrt 25) 5)))
(ert-deftest floatfns-tests-abs ()
(should (= (abs 10) 10))
(should (= (abs -10) 10)))
(ert-deftest floatfns-tests-logb ()
(should (= (logb 10000) 13)))
(ert-deftest floatfns-tests-ceiling ()
(should (= (ceiling 0.5) 1)))
(ert-deftest floatfns-tests-floor ()
(should (= (floor 1.5) 1)))
(ert-deftest floatfns-tests-round ()
(should (= (round 1.49999999999) 1))
(should (= (round 1.50000000000) 2))
(should (= (round 1.50000000001) 2)))
(ert-deftest floatfns-tests-truncate ()
(should (= (truncate float-pi) 3)))
(ert-deftest floatfns-tests-fceiling ()
(should (= (fceiling 0.5) 1.0)))
(ert-deftest floatfns-tests-ffloor ()
(should (= (ffloor 1.5) 1.0)))
(ert-deftest floatfns-tests-fround ()
(should (= (fround 1.49999999999) 1.0))
(should (= (fround 1.50000000000) 2.0))
(should (= (fround 1.50000000001) 2.0)))
(ert-deftest floatfns-tests-ftruncate ()
(should (= (ftruncate float-pi) 3.0)))
(ert-deftest divide-extreme-sign ()
(should (= (ceiling most-negative-fixnum -1.0) (- most-negative-fixnum)))
(should (= (floor most-negative-fixnum -1.0) (- most-negative-fixnum)))
(should (= (round most-negative-fixnum -1.0) (- most-negative-fixnum)))
(should (= (truncate most-negative-fixnum -1.0) (- most-negative-fixnum))))
(ert-deftest logb-extreme-fixnum ()
(should (= (logb most-negative-fixnum) (1+ (logb most-positive-fixnum)))))
(ert-deftest fround-fixnum ()
(should-error (ffloor 0) :type 'wrong-type-argument)
(should-error (fceiling 0) :type 'wrong-type-argument)
(should-error (ftruncate 0) :type 'wrong-type-argument)
(should-error (fround 0) :type 'wrong-type-argument))
(ert-deftest bignum-to-float ()
;; 122 because we want to go as big as possible to provoke a rounding error,
;; but not too big: 2**122 < 10**37 < 2**123, and the C standard says
;; 10**37 <= DBL_MAX so 2**122 cannot overflow as a double.
(let ((a (1- (ash 1 122))))
(should (or (eql a (1- (floor (float a))))
(eql a (floor (float a))))))
(should (eql (float (+ most-positive-fixnum 1))
(+ (float most-positive-fixnum) 1))))
(ert-deftest bignum-abs ()
(should (= most-positive-fixnum
(- (abs most-negative-fixnum) 1))))
(ert-deftest bignum-expt ()
(dolist (n (list most-positive-fixnum (1+ most-positive-fixnum)
most-negative-fixnum (1- most-negative-fixnum)
(* 5 most-negative-fixnum)
(* 5 (1+ most-positive-fixnum))
-2 -1 0 1 2))
(should (or (<= n 0) (= (expt 0 n) 0)))
(should (= (expt 1 n) 1))
(should (or (< n 0) (= (expt -1 n) (if (zerop (logand n 1)) 1 -1))))
(should (= (expt n 0) 1))
(should (= (expt n 1) n))
(should (= (expt n 2) (* n n)))
(should (= (expt n 3) (* n n n)))))
(ert-deftest bignum-logb ()
(should (= (+ (logb most-positive-fixnum) 1)
(logb (+ most-positive-fixnum 1)))))
(ert-deftest bignum-mod ()
(should (= 0 (mod (1+ most-positive-fixnum) 2.0))))
(ert-deftest bignum-round ()
(let ((ns (list (* most-positive-fixnum most-negative-fixnum)
(1- most-negative-fixnum) most-negative-fixnum
(1+ most-negative-fixnum) -2 1 1 2
(1- most-positive-fixnum) most-positive-fixnum
(1+ most-positive-fixnum)
(* most-positive-fixnum most-positive-fixnum))))
(dolist (n ns)
(should (= n (ceiling n)))
(should (= n (floor n)))
(should (= n (round n)))
(should (= n (truncate n)))
(let ((-n (- n))
(f (float n))
(-f (- (float n))))
(should (= 1 (round n f) (round -n -f) (round f n) (round -f -n)))
(should (= -1 (round -n f) (round n -f) (round f -n) (round -f n))))
(dolist (d ns)
(let ((q (/ n d))
(r (% n d))
(same-sign (eq (< n 0) (< d 0))))
(should (= (ceiling n d)
(+ q (if (and same-sign (not (zerop r))) 1 0))))
(should (= (floor n d)
(- q (if (and (not same-sign) (not (zerop r))) 1 0))))
(should (= (truncate n d) q))
(let ((cdelta (abs (- n (* d (ceiling n d)))))
(fdelta (abs (- n (* d (floor n d)))))
(rdelta (abs (- n (* d (round n d))))))
(should (<= rdelta cdelta))
(should (<= rdelta fdelta))
(should (if (zerop r)
(= 0 cdelta fdelta rdelta)
(or (/= cdelta fdelta)
(zerop (% (round n d) 2)))))))))))
(ert-deftest special-round ()
(dolist (f '(ceiling floor round truncate))
(let ((ns '(-1e+INF 1e+INF -1 -0.0 0.0 0 1 -1e+NaN 1e+NaN)))
(dolist (n ns)
(if (not (<= (abs n) 1))
(should-error (funcall f n))
(should (= n (funcall f n)))
(dolist (d '(-1e+INF 1e+INF))
(should (eq 0 (funcall f n d)))))
(dolist (d ns)
(when (or (zerop d) (= (abs n) 1e+INF) (not (= n n)) (not (= d d)))
(should-error (funcall f n d))))))))
(ert-deftest big-round ()
(should (= (floor 54043195528445955 3)
(floor 54043195528445955 3.0)))
(should (= (floor 1.7976931348623157e+308 5e-324)
(ash (1- (ash 1 53)) 2045))))
(provide 'floatfns-tests)
;;; floatfns-tests.el ends here
|
