- (assert (not (eval ',nform)))
- (assert (not (funcall (lambda () ,nform))))))))
- ;; Source transforms for >= and <= used to be too eager about
- ;; inverting the test, causing NaN issues.
- (test (>= nan 1.0))
- (test (>= 1.0 nan))
- (test (>= 1.0 nan 0.0))
- (test (>= 1.0 0.0 nan))
- (test (>= nan 1.0 0.0))
- (test (<= nan 1.0))
- (test (<= 1.0 nan))
- (test (<= 1.0 nan 2.0))
- (test (<= 1.0 2.0 nan))
- (test (<= nan 1.0 2.0)))))
+ (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)))))