<html lang='en'>
<head>
<meta content='text/html; charset=utf-8' http-equiv='Content-Type'>
<title>
GitLab
</title>
</meta>
</head>
<style>
img {
max-width: 100%;
height: auto;
}
p.details {
font-style:italic;
color:#777
}
.footer p {
font-size:small;
color:#777
}
pre.commit-message {
white-space: pre-wrap;
}
.file-stats a {
text-decoration: none;
}
.file-stats .new-file {
color: #090;
}
.file-stats .deleted-file {
color: #B00;
}
</style>
<body>
<div class='content'>
<h3>Raymond Toy pushed to branch master at <a href="https://gitlab.common-lisp.net/cmucl/cmucl">cmucl / cmucl</a></h3>
<h4>
Commits:
</h4>
<ul>
<li>
<strong><a href="https://gitlab.common-lisp.net/cmucl/cmucl/commit/106b9498b22385eeea7968f8b7f04ded5f5325c2">106b9498</a></strong>
<div>
<span>by Raymond Toy</span>
<i>at 2015-12-20T12:18:26Z</i>
</div>
<pre class='commit-message'>Add log tests.</pre>
</li>
<li>
<strong><a href="https://gitlab.common-lisp.net/cmucl/cmucl/commit/d1ef807a36b6f93b9d9ce8ae89e6d131f1be49d7">d1ef807a</a></strong>
<div>
<span>by Raymond Toy</span>
<i>at 2015-12-20T12:33:05Z</i>
</div>
<pre class='commit-message'>Add sinh tests.</pre>
</li>
</ul>
<h4>1 changed file:</h4>
<ul>
<li class='file-stats'>
<a href='#diff-0'>
tests/fdlibm.lisp
</a>
</li>
</ul>
<h4>Changes:</h4>
<li id='diff-0'>
<a href='https://gitlab.common-lisp.net/cmucl/cmucl/compare/4c29796b2ba54474fff55a484e0647ab9ca855f1...d1ef807a36b6f93b9d9ce8ae89e6d131f1be49d7#diff-0'>
<strong>
tests/fdlibm.lisp
</strong>
</a>
<hr>
<pre class="highlight"><code><span style="color: #000000;background-color: #ffdddd">--- a/tests/fdlibm.lisp
</span><span style="color: #000000;background-color: #ddffdd">+++ b/tests/fdlibm.lisp
</span><span style="color: #aaaaaa">@@ -403,3 +403,115 @@
</span> (assert-error 'floating-point-overflow (exp 1000d0))
;; exp(-1000) underflow
(assert-eql 0d0 (exp -1000d0)))
<span style="color: #000000;background-color: #ddffdd">+
+(define-test log-basic-tests
+ (:tag :fdlibm)
+ (assert-eql 0d0 (log 1d0))
+ (assert-eql 1d0 (log (exp 1d0)))
+ (assert-eql -1d0 (log (exp -1d0)))
+ (assert-eql 0.5d0 (log (sqrt (exp 1d0))))
+ (assert-eql -0.5d0 (log (sqrt (exp -1d0))))
+ ;; Test a denormal arg
+ (assert-eql -709.08956571282410d0 (log (scale-float 1d0 -1023)))
+ ;; Largest double value
+ (assert-eql 709.7827128933840d0 (log most-positive-double-float))
+ ;; Tests case 0 < f < 2^-20, k = 0
+ ;; log(1+2^-21)
+ (assert-eql 4.7683704451632344d-7 (log (+ 1 (scale-float 1d0 -21))))
+ ;; Tests case 0 < f < 2^-20, k = 1
+ ;; log(2 + 2^-20)
+ (assert-eql 0.6931476573969898d0 (log (+ 2(scale-float 1d0 -20))))
+ (assert-eql 1.3862943611198906d0 (log 4d0))
+ ;; Tests main path, i > 0, k = 0
+ (assert-eql 0.3220828910287846d0
+ (log (kernel:make-double-float (+ #x3ff00000 #x6147a) 0)))
+ ;; Tests main path, i > 0, k = 1
+ (assert-eql 0.35065625373947773d0
+ (log (kernel:make-double-float (+ #x3ff00000 #x6b851) 0)))
+ ;; Tests main path, i > 0, k = -1
+ (assert-eql -0.3710642895311607d0
+ (log (kernel:make-double-float (+ #x3fe00000 #x6147a) 0)))
+ ;; Tests main path, i < 0, k = 0
+ (assert-eql 0.3220821999597803d0
+ (log (kernel:make-double-float (+ #x3ff00000 #x61479) 0)))
+ ;; Tests main path, i < 0, k = 1
+ (assert-eql 1.0152293805197257d0
+ (log (kernel:make-double-float (+ #x40000000 #x61479) 0)))
+ ;; Tests main path, i < 0, k = -1
+ (assert-eql -0.37106498060016496d0
+ (log (kernel:make-double-float (+ #x3fe00000 #x61479) 0))))
+
+(define-test log-consistency
+ (:tag :fdlibm)
+ ;; |log(x) + log(1/x)| < 1.77635684e-15, x = 1.2^k, 0 <= k < 2000
+ ;; The threshold is experimentally determined
+ (let ((x 1d0)
+ (max-value -1d0))
+ (declare (double-float max-value)
+ (type (double-float 1d0) x))
+ (dotimes (k 2000)
+ (let ((y (abs (+ (log x) (log (/ x))))))
+ (setf max-value (max max-value y))
+ (setf x (* x 1.4d0))))
+ (assert-true (< max-value 1.77635684d-15)))
+ ;; |exp(log(x)) - x|/x < 5.6766649d-14, x = 1.4^k, 0 <= k < 2000
+ (let ((x 1d0)
+ (max-error 0d0))
+ (declare (double-float max-error)
+ (type (double-float 1d0) x))
+ (dotimes (k 2000)
+ (let ((y (abs (/ (- (exp (log x)) x) x))))
+ (setf max-error (max max-error y))
+ (setf x (* x 1.4d0))))
+ (assert-true (< max-error 5.6766649d-14)))
+ ;; |exp(log(x)) - x|/x < 5.68410245d-14, x = 1.4^(-k), 0 <= k < 2000
+ (let ((x 1d0)
+ (max-error 0d0))
+ (declare (double-float max-error)
+ (type (double-float (0d0)) x))
+ (dotimes (k 2000)
+ (let ((y (abs (/ (- (exp (log x)) x) x))))
+ (setf max-error (max max-error y))
+ (setf x (/ x 1.4d0))))
+ (assert-true (< max-error 5.68410245d-14))))
+
+(define-test sinh-basic-tests
+ (:tag :fdlibm)
+ (assert-eql +0d0 (sinh 0d0))
+ (assert-eql -0d0 (sinh -0d0))
+ ;; sinh(x) = x, |x| < 2^-28
+ (let ((x (scale-float 1d0 -29)))
+ (assert-eql x (sinh x))
+ (assert-eql (- x) (sinh (- x))))
+ ;; case |x| < 1
+ (assert-eql 0.5210953054937474d0 (sinh 0.5d0))
+ (assert-eql -0.5210953054937474d0 (sinh -0.5d0))
+ ;; sinh(10*log(2)) = 1048575/2048, case |x| < 22
+ (let ((x (* 10 (log 2d0)))
+ (y (float 1048575/2048 1d0)))
+ (assert-eql y (sinh x))
+ (assert-eql (- y) (sinh (- x))))
+ ;; sinh(10), case |x| < 22
+ (let ((y 11013.232874703393d0))
+ (assert-eql y (sinh 10d0))
+ (assert-eql (- y) (sinh -10d0)))
+ ;; sinh(32*log(2)), case |x| in [22, log(maxdouble)]
+ (let ((x (* 32 (log 2d0)))
+ (y 2.1474836479999983d9))
+ (assert-eql y (sinh x))
+ (assert-eql (- y) (sinh (- x))))
+ ;; sinh(100), case |x| in [22, log(maxdouble)]
+ (let ((y 1.3440585709080678d43))
+ (assert-eql y (sinh 100d0))
+ (assert-eql (- y) (sinh -100d0)))
+ ;; sinh(710....), no overflow, case |x| in [log(maxdouble), overflowthreshold]
+ (let ((x 710.4758600739439d0)
+ (y 1.7976931348621744d308))
+ (assert-eql y (sinh x))
+ (assert-eql (- y) (sinh (- x))))
+ ;; sinh(710.475860073944), overflow, case |x| > ovfthreshold]
+ (let ((x 710.475860073944d0))
+ (assert-error 'floating-point-overflow (sinh x))
+ (assert-error 'floating-point-overflow (sinh (- x))))
+ (assert-error 'floating-point-overflow (sinh 1000d0))
+ (assert-error 'floating-point-overflow (sinh -1000d0)))
</span></code></pre>
<br>
</li>
</div>
<div class='footer' style='margin-top: 10px;'>
<p>
—
<br>
<a href="https://gitlab.common-lisp.net/cmucl/cmucl/compare/4c29796b2ba54474fff55a484e0647ab9ca855f1...d1ef807a36b6f93b9d9ce8ae89e6d131f1be49d7">View it on GitLab</a>.
<br>
You're receiving this email because of your account on gitlab.common-lisp.net.
If you'd like to receive fewer emails, you can adjust your notification settings.
</p>
</div>
</body>
</html>