From 1d604dc615d00763f303fb035eb8d64e987fc69c Mon Sep 17 00:00:00 2001
From: Raymond Toy <toy.raymond@gmail.com>
Date: Sat, 19 Dec 2015 09:32:37 -0800
Subject: [PATCH] Add tests for asinh, taken from fdlibm-js.

---
 tests/fdlibm.lisp | 34 ++++++++++++++++++++++++++++++++++
 1 file changed, 34 insertions(+)

diff --git a/tests/fdlibm.lisp b/tests/fdlibm.lisp
index 98657efeb..ef1fe5fc3 100644
--- a/tests/fdlibm.lisp
+++ b/tests/fdlibm.lisp
@@ -266,3 +266,37 @@
 		  (kernel:%scalbn most-positive-double-float 2))
     (assert-equal ext:double-float-negative-infinity
 		  (kernel:%scalbn most-negative-double-float 2))))
+
+;;; These tests taken from github.com/rtoy/fdlibm-js
+(define-test asinh-basic-tests
+    (:tag :fdlibm)
+  (assert-eql -0d0 (asinh -0d0))
+  (assert-eql 0d0 (asinh 0d0))
+  (let ((x (scale-float 1d0 -29)))
+    ;; asinh(x) = x for x < 2^-28
+    (assert-eql x (asinh x))
+    (assert-eql (- x) (asinh (- x))))
+  (let ((x (scale-float 1d0 -28)))
+    ;; Case 2 > |x| >= 2^-28
+    (assert-eql 3.725290298461914d-9 (asinh x))
+    (assert-eql -3.725290298461914d-9 (asinh (- x))))
+  (let ((x 1d0))
+    ;; Case 2 > |x| >= 2^-28
+    (assert-eql 0.881373587019543d0 (asinh x))
+    (assert-eql -0.881373587019543d0 (asinh (- x))))
+  (let ((x 5d0))
+    ;; Case 2^28 > |x| > 2
+    (assert-eql 2.3124383412727525d0 (asinh x))
+    (assert-eql -2.3124383412727525d0 (asinh (- x))))
+  (let ((x (scale-float 1d0 28)))
+    ;; Case 2^28 > |x|
+    (assert-eql 20.101268236238415d0 (asinh x))
+    (assert-eql -20.101268236238415d0 (asinh (- x))))
+  (let ((x most-positive-double-float))
+    ;; No overflow for most-positive-double-float
+    (assert-eql 710.4758600739439d0 (asinh x))
+    (assert-eql -710.4758600739439d0 (asinh (- x)))))
+  
+    
+  
+  
\ No newline at end of file
-- 
GitLab