X-Git-Url: http://repo.macrolet.net/gitweb/?a=blobdiff_plain;f=tests%2Ffloat.pure.lisp;h=62a54b09e54adf8ae749500a5cedbba33cbd7fc8;hb=d7c9385e5a7b27d38884ebed9df6e9e5c3ac5cbf;hp=6c4815c204ee10a64998a8f33da0aac05796bfd4;hpb=11cfbfd7e07650aa5dae84ed6170a0d041e25acb;p=sbcl.git diff --git a/tests/float.pure.lisp b/tests/float.pure.lisp index 6c4815c..62a54b0 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 :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 (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 @@ -245,12 +269,16 @@ (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 :sunos))) +;; The x86 port used not to reduce the arguments of transcendentals +;; correctly. +;; This test is valid only for x86: The x86 port uses the builtin x87 +;; FPU instructions to implement the trigonometric functions; other +;; ports rely on the system's math library. These two differ in the +;; precision of pi used for the range reduction and so yield results +;; that can differ by arbitrarily large amounts for large inputs. +;; The test expects the x87 results. +(with-test (:name (:range-reduction :x87) + :skipped-on '(not :x86)) (flet ((almost= (x y) (< (abs (- x y)) 1d-5))) (macrolet ((foo (op value) @@ -295,6 +323,12 @@ (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"))))) @@ -373,3 +407,19 @@ (+ (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))))