Fix make-array transforms.
[sbcl.git] / tests / float.pure.lisp
index ff01c2f..5612a73 100644 (file)
@@ -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 :darwin) ;; bug 372
+            :fails-on '(and :darwin :ppc)) ;; bug 372
   (progn
     (assert (raises-error? (scale-float 1.0 most-positive-fixnum)
                            floating-point-overflow))
 (funcall (compile nil '(lambda () (tan (tan (round 0))))))
 
 (with-test (:name (:addition-overflow :bug-372)
-            :fails-on '(or (and :ppc :openbsd) :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
                                        (the (eql #c(1.0 2.0))
                                          x))))))))
 
-;; This was previously x86-only, with note:
-;;   The x86 port used not to reduce the arguments of transcendentals
-;;   correctly. On other platforms, we trust libm to DTRT.
-;; but it doesn't cost any real amount to just test them all
-(with-test (:name :range-reduction
-            :fails-on '(and :x86-64 (or :linux :darwin)))
-  (flet ((almost= (x y)
-           (< (abs (- x y)) 1d-5)))
-    (macrolet ((foo (op value)
-                 `(assert (almost= (,op (mod ,value (* 2 pi)))
-                                   (,op ,value)))))
-      (let ((big (* pi (expt 2d0 70)))
-            (mid (coerce most-positive-fixnum 'double-float))
-            (odd (* pi most-positive-fixnum)))
-        (foo sin big)
-        (foo sin mid)
-        (foo sin odd)
-        (foo sin (/ odd 2d0))
-
-        (foo cos big)
-        (foo cos mid)
-        (foo cos odd)
-        (foo cos (/ odd 2d0))
-
-        (foo tan big)
-        (foo tan mid)
-        (foo tan odd)))))
-
 ;; 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))))
 (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")))))
                (+ (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))))