blob: d887939c999542e41c1f42ad23b17491e8760ea3 (
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
|
;;; floatfns-tests.el --- tests for floating point operations -*- lexical-binding: t -*-
;; Copyright 2017-2022 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 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
|
