X-Git-Url: http://repo.macrolet.net/gitweb/?a=blobdiff_plain;f=tests%2Ffloat.pure.lisp;h=5612a73fbcbb99ab070931236babd78825e7c247;hb=260de2062fca170efdac3e42491d7d866c2d2e56;hp=fee0323bdc6b9a38a7dabc126a85c050049925ab;hpb=5e0af0dad59e01274b0e84b58e5f0904c9890b37;p=sbcl.git

diff --git a/tests/float.pure.lisp b/tests/float.pure.lisp
index fee0323..5612a73 100644
--- a/tests/float.pure.lisp
+++ b/tests/float.pure.lisp
@@ -93,7 +93,7 @@
 (assert (= 0.0d0 (scale-float 1.0d0 (1- most-negative-fixnum))))
 
 (with-test (:name (:scale-float-overflow :bug-372)
-            :fails-on '(or :ppc :darwin)) ;; bug 372
+            :fails-on '(and :darwin :ppc)) ;; bug 372
   (progn
     (assert (raises-error? (scale-float 1.0 most-positive-fixnum)
                            floating-point-overflow))
@@ -125,7 +125,31 @@
 (funcall (compile nil '(lambda () (tan (tan (round 0))))))
 
 (with-test (:name (:addition-overflow :bug-372)
-            :fails-on '(or :ppc :darwin (and :x86 :netbsd)))
+            :fails-on '(or (and :ppc :openbsd)
+                           (and :ppc :darwin)
+                           (and :x86 :netbsd)))
+  (assert (typep (nth-value
+                  1
+                  (ignore-errors
+                    (sb-sys:without-interrupts
+                     (sb-int:set-floating-point-modes :current-exceptions nil
+                                                      :accrued-exceptions nil)
+                     (loop repeat 2 summing most-positive-double-float)
+                     (sleep 2))))
+                 'floating-point-overflow)))
+
+;; This is the same test as above.  Even if the above copy passes,
+;; this copy will fail if SIGFPE handling ends up clearing the FPU
+;; control word, which can happen if the kernel clears the FPU control
+;; (a reasonable thing for it to do) and the runtime fails to
+;; compensate for this (see RESTORE_FP_CONTROL_WORD in interrupt.c).
+;; Note that this only works when running float.pure.lisp alone, as
+;; the preceeding "pure" test files aren't as free of side effects as
+;; we might like.
+(with-test (:name (:addition-overflow :bug-372 :take-2)
+            :fails-on '(or (and :ppc :openbsd)
+                           (and :ppc :darwin)
+                           (and :x86 :netbsd)))
   (assert (typep (nth-value
                   1
                   (ignore-errors
@@ -234,80 +258,136 @@
     (assert (eql 0.0d0 (funcall f 123.0d0 0.0d0)))
     (assert (eql 0.0d0 (funcall f 123.0 0.0d0)))))
 
-;; 1.0.29.44 introduces a ton of changes for complex floats
-;; on x86-64. Huge test of doom to help catch weird corner
-;; cases.
-(with-test (:name :complex-floats)
-  (labels ((equal-enough (x y)
-             (cond ((eql x y))
-                   ((or (complexp x)
-                        (complexp y))
-                    (or (eql (coerce x '(complex double-float))
-                             (coerce y '(complex double-float)))
-                        (and (equal-enough (realpart x) (realpart y))
-                             (equal-enough (imagpart x) (imagpart y)))))
-                   ((numberp x)
-                    (or (eql (coerce x 'double-float) (coerce y 'double-float))
-                        (< (abs (- x y))  1d-5)))))
-           (reflections (x)
-             (values x
-                     (conjugate x)
-                     (complex (- (realpart x)) (imagpart x))
-                     (- x)))
-           (compute (x y r)
-             (list (+ x y) (+ r x) (+ x r)
-                   (- x y) (- r x) (- x r)
-                   (* x y) (* x r) (* r x)
-                   (unless (zerop y)
-                     (/ x y))
-                   (unless (zerop r)
-                     (/ x r))
-                   (unless (zerop x)
-                     (/ r x))
-                   (conjugate x) (conjugate r)
-                   (- x)
-                   (complex r) (complex r r) (complex 0 r)
-                   (= x y) (= r x) (= y r) (= x (complex 0 r))
-                   (eql x y) (eql x (complex r)) (eql y (complex r))
-                   (eql x (complex r r)) (eql y (complex 0 r))))
-           (compute-all (x y r)
-             (multiple-value-bind (x1 x2 x3 x4) (reflections x)
-               (multiple-value-bind (y1 y2 y3 y4) (reflections y)
-                 #.(let ((form '(list)))
-                     (dolist (x '(x1 x2 x3 x4) (reverse form))
-                       (dolist (y '(y1 y2 y3 y4))
-                         (push `(list ,x ,y r
-                                      (append (compute ,x ,y r)
-                                              (compute ,x ,y (- r))))
-                               form))))))))
-    (declare (inline reflections compute compute-all))
-    (let* ((reals     '(0 1 2))
-           (complexes '#.(let ((reals '(0 1 2))
-                               (cpx   '()))
-                           (dolist (x reals (nreverse cpx))
-                             (dolist (y reals)
-                               (push (complex x y) cpx)))))
-           (val       ()))
-      (declare (notinline every))
-      (dolist (r reals (nreverse val))
-        (dolist (x complexes)
-          (dolist (y complexes)
-            (let ((value  (compute-all x y r))
-                  (single (compute-all (coerce x '(complex single-float))
-                                       (coerce y '(complex single-float))
-                                       (coerce r 'single-float)))
-                  (double (compute-all (coerce x '(complex double-float))
-                                       (coerce y '(complex double-float))
-                                       (coerce r 'double-float))))
-              (assert (every (lambda (pos ref single double)
-                               (every (lambda (ref single double)
-                                        (or (and (equal-enough ref single)
-                                                 (equal-enough ref double))
-                                            (and (not (numberp single)) ;; -ve 0s
-                                                 (equal-enough single double))))
-                                      (fourth ref) (fourth single) (fourth double)))
-                             '((0 0) (0 1) (0 2) (0 3)
-                               (1 0) (1 1) (1 2) (1 3)
-                               (2 0) (2 1) (2 2) (2 3)
-                               (3 0) (3 1) (3 2) (3 3))
-                             value single double)))))))))
+;; Bug reported by Eric Marsden on July 15 2009. The compiler
+;; used not to constant fold calls with arguments of type
+;; (EQL foo).
+(with-test (:name :eql-type-constant-fold)
+  (assert (equal '(FUNCTION (T) (VALUES (MEMBER T) &OPTIONAL))
+                 (sb-kernel:%simple-fun-type
+                  (compile nil `(lambda (x)
+                                  (eql #c(1.0 2.0)
+                                       (the (eql #c(1.0 2.0))
+                                         x))))))))
+
+;; Leakage from the host could result in wrong values for truncation.
+(with-test (:name :truncate)
+  (assert (plusp (sb-kernel:%unary-truncate/single-float (expt 2f0 33))))
+  (assert (plusp (sb-kernel:%unary-truncate/double-float (expt 2d0 33))))
+  ;; That'd be one strange host, but just in case
+  (assert (plusp (sb-kernel:%unary-truncate/single-float (expt 2f0 65))))
+  (assert (plusp (sb-kernel:%unary-truncate/double-float (expt 2d0 65)))))
+
+;; On x86-64, we sometimes forgot to clear the higher order bits of the
+;; destination register before using it with an instruction that doesn't
+;; clear the (unused) high order bits. Suspect instructions are operations
+;; with only one operand: for everything else, the destination has already
+;; been loaded with a value, making it safe (by induction).
+;;
+;; The tests are extremely brittle and could be broken by any number of
+;; back- or front-end optimisations. We should just keep the issue above
+;; in mind at all times when working with SSE or similar instruction sets.
+;;
+;; Run only on x86/x86-64m as no other platforms have SB-VM::TOUCH-OBJECT.
+(macrolet ((with-pinned-floats ((count type &rest names) &body body)
+             "Force COUNT float values to be kept live (and hopefully in registers),
+              fill a temporary register with noise, and execute BODY."
+             ;; KLUDGE: SB-VM is locked, and non-x86oids don't have
+             ;; SB-VM::TOUCH-OBJECT.  Don't even READ this body on
+             ;; other platforms.
+             #-(or x86 x86-64)
+             (declare (ignore count type names body))
+             #+(or x86 x86-64)
+             (let ((dummy (loop repeat count
+                                collect (or (pop names)
+                                            (gensym "TEMP")))))
+               `(let ,(loop for i downfrom -1
+                            for var in dummy
+                            for j = (coerce i type)
+                            collect
+                            `(,var ,(complex j j))) ; we don't actually need that, but
+                  (declare (type (complex ,type) ,@dummy)) ; future-proofing can't hurt
+                  ,@(loop for var in dummy
+                          for i upfrom 0
+                          collect `(setf ,var ,(complex i (coerce i type))))
+                  (multiple-value-prog1
+                      (progn
+                        (let ((x ,(complex 1d0 1d0)))
+                          (declare (type (complex double-float) x))
+                          (setf x ,(complex most-positive-fixnum (float most-positive-fixnum 1d0)))
+                          (sb-vm::touch-object x))
+                        (locally ,@body))
+                    ,@(loop for var in dummy
+                            collect `(sb-vm::touch-object ,var)))))))
+  (with-test (:name :clear-sqrtsd :skipped-on '(not (or :x86 :x86-64)))
+    (flet ((test-sqrtsd (float)
+             (declare (optimize speed (safety 1))
+                      (type (double-float (0d0)) float))
+             (with-pinned-floats (14 double-float x0)
+               (let ((x (sqrt float)))
+                 (values (+ x x0) float)))))
+      (declare (notinline test-sqrtsd))
+      (assert (zerop (imagpart (test-sqrtsd 4d0))))))
+
+  (with-test (:name :clear-sqrtsd-single :skipped-on '(not (or :x86 :x86-64)))
+    (flet ((test-sqrtsd-float (float)
+             (declare (optimize speed (safety 1))
+                      (type (single-float (0f0)) float))
+             (with-pinned-floats (14 single-float x0)
+               (let ((x (sqrt float)))
+                 (values (+ x x0) float)))))
+      (declare (notinline test-sqrtsd-float))
+      (assert (zerop (imagpart (test-sqrtsd-float 4f0))))))
+
+  (with-test (:name :clear-cvtss2sd :skipped-on '(not (or :x86 :x86-64)))
+    (flet ((test-cvtss2sd (float)
+             (declare (optimize speed (safety 1))
+                      (type single-float float))
+             (with-pinned-floats (14 double-float x0)
+               (let ((x (float float 0d0)))
+                 (values (+ x x0) (+ 1e0 float))))))
+      (declare (notinline test-cvtss2sd))
+      (assert (zerop (imagpart (test-cvtss2sd 1f0))))))
+
+  (with-test (:name :clear-cvtsd2ss :skipped-on '(not (or :x86 :x86-64)))
+    (flet ((test-cvtsd2ss (float)
+             (declare (optimize speed (safety 1))
+                      (type double-float float))
+             (with-pinned-floats (14 single-float x0)
+               (let ((x (float float 1e0)))
+                 (values (+ x x0) (+ 1d0 float))))))
+      (declare (notinline test-cvtsd2ss))
+      (assert (zerop (imagpart (test-cvtsd2ss 4d0))))))
+
+  (with-test (:name :clear-cvtsi2sd :skipped-on '(not (or :x86 :x86-64)))
+    (flet ((test-cvtsi2sd (int)
+             (declare (optimize speed (safety 0))
+                      (type (unsigned-byte 10) int))
+             (with-pinned-floats (15 double-float x0)
+               (+ (float int 0d0) x0))))
+      (declare (notinline test-cvtsi2sd))
+      (assert (zerop (imagpart (test-cvtsi2sd 4))))))
+
+  (with-test (:name :clear-cvtsi2ss :skipped-on '(not (or :x86 :x86-64)))
+    (flet ((test-cvtsi2ss (int)
+             (declare (optimize speed (safety 0))
+                      (type (unsigned-byte 10) int))
+             (with-pinned-floats (15 single-float x0)
+               (+ (float int 0e0) x0))))
+      (declare (notinline test-cvtsi2ss))
+      (assert (zerop (imagpart (test-cvtsi2ss 4)))))))
+
+(with-test (:name :round-to-bignum)
+  (assert (= (round 1073741822.3d0) 1073741822))
+  (assert (= (round 1073741822.5d0) 1073741822))
+  (assert (= (round 1073741822.7d0) 1073741823))
+  (assert (= (round 1073741823.3d0) 1073741823))
+  (assert (= (round 1073741823.5d0) 1073741824))
+  (assert (= (round 1073741823.7d0) 1073741824)))
+
+(with-test (:name :round-single-to-bignum)
+  (assert (= (round 1e14) 100000000376832))
+  (assert (= (round 1e19) 9999999980506447872)))
+
+(with-test (:name :scaled-%hypot)
+  (assert (<= (abs (complex most-positive-double-float 1d0))
+              (1+ most-positive-double-float))))