(cl:in-package :cl-user)
-(let ((+ifni single-float-positive-infinity)
- (-ifni single-float-negative-infinity))
- (assert (= (* +ifni 1) +ifni))
- (assert (= (* +ifni -0.1) -ifni))
- (assert (= (+ +ifni -0.1) +ifni))
- (assert (= (- +ifni -0.1) +ifni))
- (assert (= (sqrt +ifni) +ifni))
- (assert (= (* -ifni -14) +ifni))
- (assert (= (/ -ifni 0.1) -ifni))
- (assert (= (/ -ifni 100/3) -ifni))
- (assert (< -ifni +ifni))
- (assert (not (< +ifni 100)))
- (assert (not (< +ifni -ifni))))
+(dolist (ifnis (list (cons single-float-positive-infinity
+ single-float-negative-infinity)
+ (cons double-float-positive-infinity
+ double-float-negative-infinity)))
+ (destructuring-bind (+ifni . -ifni) ifnis
+ (assert (= (* +ifni 1) +ifni))
+ (assert (= (* +ifni -0.1) -ifni))
+ (assert (= (+ +ifni -0.1) +ifni))
+ (assert (= (- +ifni -0.1) +ifni))
+ (assert (= (sqrt +ifni) +ifni))
+ (assert (= (* -ifni -14) +ifni))
+ (assert (= (/ -ifni 0.1) -ifni))
+ (assert (= (/ -ifni 100/3) -ifni))
+ (assert (not (= +ifni -ifni)))
+ (assert (= -ifni -ifni))
+ (assert (not (= +ifni 100/3)))
+ (assert (not (= -ifni -1.0 -ifni)))
+ (assert (not (= -ifni -17/02 -ifni)))
+ (assert (< -ifni +ifni))
+ (assert (not (< +ifni 100)))
+ (assert (not (< +ifni 100.0)))
+ (assert (not (< +ifni -ifni)))
+ (assert (< 100 +ifni))
+ (assert (< 100.0 +ifni))
+ (assert (>= 100 -ifni))
+ (assert (not (<= 6/7 (* 3 -ifni))))
+ (assert (not (> +ifni +ifni)))))
+
+;;; ANSI: FLOAT-RADIX should signal an error if its argument is not a
+;;; float.
+;;;
+;;; (Peter Van Eynde's ansi-test suite caught this, and Eric Marsden
+;;; reported a fix for CMU CL, which was ported to sbcl-0.6.12.35.)
+(assert (typep (nth-value 1 (ignore-errors (float-radix "notfloat")))
+ 'type-error))
+
+(assert (typep (nth-value 1 (ignore-errors
+ (funcall (fdefinition 'float-radix) "notfloat")))
+ 'type-error))
+
+;;; Before 0.8.2.14 the cross compiler failed to work with
+;;; denormalized numbers
+(when (subtypep 'single-float 'short-float)
+ (assert (eql least-positive-single-float least-positive-short-float)))
+
+;;; bug found by Paul Dietz: FFLOOR and similar did not work for integers
+(let ((tests '(((ffloor -8 3) (-3.0 1))
+ ((fround -8 3) (-3.0 1))
+ ((ftruncate -8 3) (-2.0 -2))
+ ((fceiling -8 3) (-2.0 -2)))))
+ (loop for (exp res) in tests
+ for real-res = (multiple-value-list (eval exp))
+ do (assert (equal real-res res))))
+
+;;; bug 45b reported by PVE
+(dolist (type '(short single double long))
+ (dolist (sign '(positive negative))
+ (let* ((name (find-symbol (format nil "LEAST-~A-~A-FLOAT"
+ sign type)
+ :cl))
+ (value (symbol-value name)))
+ (assert (zerop (/ value 2))))))
+
+;;; bug found by Paul Dietz: bad rounding on small floats
+(assert (= (fround least-positive-short-float least-positive-short-float) 1.0))
+
+;;; bug found by Peter Seibel: scale-float was only accepting float
+;;; exponents, when it should accept all integers. (also bug #269)
+(assert (= (multiple-value-bind (significand expt sign)
+ (integer-decode-float least-positive-double-float)
+ (* (scale-float (float significand 0.0d0) expt) sign))
+ least-positive-double-float))
+(assert (= (multiple-value-bind (significand expt sign)
+ (decode-float least-positive-double-float)
+ (* (scale-float significand expt) sign))
+ least-positive-double-float))
+(assert (= 0.0 (scale-float 1.0 most-negative-fixnum)))
+(assert (= 0.0d0 (scale-float 1.0d0 (1- most-negative-fixnum))))
+#-(or darwin x86-64) ;; bug 372 / 378
+(progn
+ (assert (raises-error? (scale-float 1.0 most-positive-fixnum)
+ floating-point-overflow))
+ (assert (raises-error? (scale-float 1.0d0 (1+ most-positive-fixnum))
+ floating-point-overflow)))
+
+;;; MISC.564: no out-of-line %ATAN2 for constant folding
+(assert (typep
+ (funcall
+ (compile
+ nil
+ '(lambda (p1)
+ (declare (optimize (speed 3) (safety 2) (debug 3) (space 0))
+ (type complex p1))
+ (phase (the (eql #c(1.0d0 2.0d0)) p1))))
+ #c(1.0d0 2.0d0))
+ 'double-float))