(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 '(or :ppc :darwin (and :x86 :openbsd))) ;; 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 :ppc :darwin (and :x86 :netbsd)))
+ :fails-on '(or :ppc :darwin (and (or :x86 :x86-64)
+ (or :netbsd :openbsd))))
(assert (typep (nth-value
1
(ignore-errors
(with-test (:name :nan-comparisons
- :fails-on '(or :x86-64 :sparc :mips))
+ :fails-on '(or :sparc :mips))
(sb-int:with-float-traps-masked (:invalid)
(macrolet ((test (form)
(let ((nform (subst '(/ 0.0 0.0) 'nan form)))
(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)))))