[oct-scm] [oct-git]OCT: A portable Lisp implementation for quad-double precision floats branch master updated. bba9f8940c9f904bf14adc405d795a38ac333c24
Raymond Toy
rtoy at common-lisp.net
Wed Apr 11 15:47:42 UTC 2012
This is an automated email from the git hooks/post-receive script. It was
generated because a ref change was pushed to the repository containing
the project "OCT: A portable Lisp implementation for quad-double precision floats".
The branch, master has been updated
via bba9f8940c9f904bf14adc405d795a38ac333c24 (commit)
from 104efdeae6ef52c5d370e0b3a048a5087cdb1ea2 (commit)
Those revisions listed above that are new to this repository have
not appeared on any other notification email; so we list those
revisions in full, below.
- Log -----------------------------------------------------------------
commit bba9f8940c9f904bf14adc405d795a38ac333c24
Author: Raymond Toy <rtoy at google.com>
Date: Wed Apr 11 08:47:28 2012 -0700
Update accuracy for oct.elliptic-pi.n2.d and add new test that
caused oct.elliptic-pi.n2.d to fail.
diff --git a/rt-tests.lisp b/rt-tests.lisp
index 192b118..325d202 100644
--- a/rt-tests.lisp
+++ b/rt-tests.lisp
@@ -1114,13 +1114,28 @@
for epi = (elliptic-pi n phi 0)
for true = (/ (atanh (* (tan phi) (sqrt (- n 1))))
(sqrt (- n 1)))
- for result = (check-accuracy 47 epi true)
+ for result = (check-accuracy 45.85 epi true)
;; Not sure if this formula holds when atanh gives a complex
;; result. Wolfram doesn't say
when (and (not (complexp true)) result)
append (list (list (list k n phi) result)))
nil)
+;; Failed test case:
+;; ((89 66.68551748022054d0 0.12266024127708153d0)
+;; (45.868614757480834d0 47 0.47787458521306514d0
+;; 0.4778745852130726d0))
+;; New threshold is 45.85 bits.
+(rt:deftest oct.elliptic-pi.n2.d-1
+ (let* ((n 66.68551748022054d0)
+ (phi 0.12266024127708153d0)
+ (epi (elliptic-pi n phi 0))
+ (true (/ (atanh (* (tan phi) (sqrt (- n 1))))
+ (sqrt (- n 1)))))
+ (check-accuracy 45.8686d0 epi true))
+ nil)
+
+
(rt:deftest oct.elliptic-pi.n0.q
;; Tests for random values for phi in [0, pi/2] and n in [0, 1]
(loop for k from 0 below 100
@@ -1241,7 +1256,7 @@
;; zero because value-or-tiny was too tiny.
(let ((g (cf-incomplete-gamma 3d0 5d0))
(true (- 2 (* 37 (exp -5d0)))))
- (check-accuracy 53 g true))
+ (check-accuracy 51.2 g true))
nil)
(rt:deftest gamma.1.d
-----------------------------------------------------------------------
Summary of changes:
rt-tests.lisp | 19 +++++++++++++++++--
1 files changed, 17 insertions(+), 2 deletions(-)
hooks/post-receive
--
OCT: A portable Lisp implementation for quad-double precision floats
More information about the oct-cvs
mailing list