1.0.25.40: fix JOIN-THREAD

[sbcl.git] / tests / float.pure.lisp

diff --git a/tests/float.pure.lisp b/tests/float.pure.lisp
index cebcc4e..1b983b7 100644 (file)
--- a/tests/float.pure.lisp
+++ b/tests/float.pure.lisp
@@ -92,7 +92,8 @@
 (assert (= 0.0 (scale-float 1.0 most-negative-fixnum)))
 (assert (= 0.0d0 (scale-float 1.0d0 (1- most-negative-fixnum))))
 
-(with-test (:fails-on '(or :darwin)) ;; bug 372
+(with-test (:name (:scale-float-overflow :bug-372)
+            :fails-on '(or :ppc :darwin (and :x86 :openbsd))) ;; bug 372
   (progn
     (assert (raises-error? (scale-float 1.0 most-positive-fixnum)
                            floating-point-overflow))
@@ -117,10 +118,116 @@
    #c(1.0d0 2.0d0))
     'double-float))
 
-(assert (typep (nth-value
-                1
-                (ignore-errors
-                  (sb-sys:without-interrupts
-                    (loop repeat 2 summing most-positive-double-float)
-                    (sleep 2))))
-               'floating-point-overflow))
+;;; More out of line functions (%COS, %SIN, %TAN) for constant folding,
+;;; reported by Mika Pihlajamäki
+(funcall (compile nil '(lambda () (cos (tan (round 0))))))
+(funcall (compile nil '(lambda () (sin (tan (round 0))))))
+(funcall (compile nil '(lambda () (tan (tan (round 0))))))
+
+(with-test (:name (:addition-overflow :bug-372)
+            :fails-on '(or :ppc :darwin (and :x86 (or :netbsd :openbsd))))
+  (assert (typep (nth-value
+                  1
+                  (ignore-errors
+                    (sb-sys:without-interrupts
+                     (loop repeat 2 summing most-positive-double-float)
+                     (sleep 2))))
+                 'floating-point-overflow)))
+
+;;; On x86-64 generating complex floats on the stack failed an aver in
+;;; the compiler if the stack slot was the same as the one containing
+;;; the real part of the complex. The following expression was able to
+;;; trigger this in 0.9.5.62.
+(with-test (:name :complex-float-stack)
+  (dolist (type '((complex double-float)
+                  (complex single-float)))
+    (compile nil
+             `(lambda (x0 x1 x2 x3 x4 x5 x6 x7)
+                (declare (type ,type x0 x1 x2 x3 x4 x5 x6 x7))
+                (let ((x0 (+ x0 x0))
+                      (x1 (+ x1 x1))
+                      (x2 (+ x2 x2))
+                      (x3 (+ x3 x3))
+                      (x4 (+ x4 x4))
+                      (x5 (+ x5 x5))
+                      (x6 (+ x6 x6))
+                      (x7 (+ x7 x7)))
+                  (* (+ x0 x1 x2 x3) (+ x4 x5 x6 x7)
+                     (+ x0 x2 x4 x6) (+ x1 x3 x5 x7)
+                     (+ x0 x3 x4 x7) (+ x1 x2 x5 x6)
+                     (+ x0 x1 x6 x7) (+ x2 x3 x4 x5)))))))
+
+
+(with-test (:name :nan-comparisons
+            :fails-on '(or :sparc :mips))
+  (sb-int:with-float-traps-masked (:invalid)
+    (macrolet ((test (form)
+                 (let ((nform (subst '(/ 0.0 0.0) 'nan form)))
+                   `(progn
+                      (assert (eval ',nform))
+                      (assert (eval `(let ((nan (/ 0.0 0.0)))
+                                       ,',form)))
+                      (assert (funcall
+                               (compile nil `(lambda () ,',nform))))
+                      (assert (funcall
+                               (compile nil `(lambda (nan) ,',form))
+                               (/ 0.0 0.0)))))))
+      (test (/= nan nan))
+      (test (/= nan nan nan))
+      (test (/= 1.0 nan 2.0 nan))
+      (test (/= nan 1.0 2.0 nan))
+      (test (not (= nan 1.0)))
+      (test (not (= nan nan)))
+      (test (not (= nan nan nan)))
+      (test (not (= 1.0 nan)))
+      (test (not (= nan 1.0)))
+      (test (not (= 1.0 1.0 nan)))
+      (test (not (= 1.0 nan 1.0)))
+      (test (not (= nan 1.0 1.0)))
+      (test (not (>= nan nan)))
+      (test (not (>= nan 1.0)))
+      (test (not (>= 1.0 nan)))
+      (test (not (>= 1.0 nan 0.0)))
+      (test (not (>= 1.0 0.0 nan)))
+      (test (not (>= nan 1.0 0.0)))
+      (test (not (<= nan nan)))
+      (test (not (<= nan 1.0)))
+      (test (not (<= 1.0 nan)))
+      (test (not (<= 1.0 nan 2.0)))
+      (test (not (<= 1.0 2.0 nan)))
+      (test (not (<= nan 1.0 2.0)))
+      (test (not (< nan nan)))
+      (test (not (< -1.0 nan)))
+      (test (not (< nan 1.0)))
+      (test (not (> nan nan)))
+      (test (not (> -1.0 nan)))
+      (test (not (> nan 1.0))))))
+
+(with-test (:name :log-int/double-accuracy)
+  ;; we used to use single precision for intermediate results
+  (assert (eql 2567.6046442221327d0
+               (log (loop for n from 1 to 1000 for f = 1 then (* f n)
+                          finally (return f))
+                    10d0)))
+  ;; both ways
+  (assert (eql (log 123123123.0d0 10) (log 123123123 10.0d0))))
+
+(with-test (:name :log-base-zero-return-type)
+  (assert (eql 0.0f0 (log 123 (eval 0))))
+  (assert (eql 0.0d0 (log 123.0d0 (eval 0))))
+  (assert (eql 0.0d0 (log 123 (eval 0.0d0))))
+  (let ((f (compile nil '(lambda (x y)
+                          (declare (optimize speed))
+                          (etypecase x
+                            (single-float
+                             (etypecase y
+                               (single-float (log x y))
+                               (double-float (log x y))))
+                            (double-float
+                             (etypecase y
+                               (single-float (log x y))
+                               (double-float (log x y)))))))))
+    (assert (eql 0.0f0 (funcall f 123.0 0.0)))
+    (assert (eql 0.0d0 (funcall f 123.0d0 0.0)))
+    (assert (eql 0.0d0 (funcall f 123.0d0 0.0d0)))
+    (assert (eql 0.0d0 (funcall f 123.0 0.0d0)))))