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 5566bc397937c2e8979ee6847e8f59f279f1f643 (commit) via 2f48eb872f71a8f41187decbc57878cb67fa86de (commit) from 5daefcc43f8b88741634f7a982bdb8b976739a1c (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 5566bc397937c2e8979ee6847e8f59f279f1f643 Author: Raymond Toy <toy.raymond@gmail.com> Date: Thu Apr 12 23:24:03 2012 -0700 Add more tests for Bessel J. Not all of them pass. diff --git a/rt-tests.lisp b/rt-tests.lisp index a4dac02..5b9f5ad 100644 --- a/rt-tests.lisp +++ b/rt-tests.lisp @@ -1673,28 +1673,54 @@ nil) ;; Bessel J for complex args -(rt:deftest bessel-j-complex.pos-order.d.1 +(rt:deftest bessel-j-complex-arg.d.1 (let ((b (bessel-j 0d0 #c(1d0 1))) (true #c(0.9376084768060293d0 -0.4965299476091221d0))) (check-accuracy 50.73 b true)) nil) -(rt:deftest bessel-j-complex.pos-order.d.2 +(rt:deftest bessel-j-complex-arg.d.2 (let ((b (bessel-j 1d0 #c(1d0 1))) (true #c(0.6141603349229036d0 0.3650280288270878d0))) (check-accuracy 52.51 b true)) nil) -(rt:deftest bessel-j-complex.pos-order.d.3 +(rt:deftest bessel-j-complex-arg.d.3 (let ((b (bessel-j 2d0 #c(1d0 1))) (true #c(0.0415798869439621d0 0.2473976415133063d0))) (check-accuracy 50.41 b true)) nil) -#+nil -(rt:deftest bessel-j-complex.pos-order.d.4 +(rt:deftest bessel-j-complex-arg.d.4 (let ((b (bessel-j 2.3d0 #c(1d0 1))) + (true #c(-0.0141615213034667d0 0.1677798241687935d0))) + (check-accuracy 48.56 b true)) + nil) + +(rt:deftest bessel-j-complex-arg.d.5 + (let ((b (bessel-j -2.3d0 #c(1d0 1))) (true #c(0.1920598664138632d0 -0.5158676904105332d0))) - (check-accuracy 53 b true)) + (check-accuracy 50.97 b true)) nil) +(rt:deftest bessel-j-1/2-complex.d.1 + (loop for k from 0 below 10 + for x = (complex (random (/ pi 2)) + (random (/ pi 2))) + for b = (bessel-j 0.5d0 x) + for true = (* (/ (sin x) (sqrt x)) (sqrt (/ 2 pi))) + for result = (check-accuracy 49.8 b true) + when result + append (list (list (list k x) result))) + nil) + +(rt:deftest bessel-j-1/2-complex.q.1 + (loop for k from 0 below 10 + for x = (complex (random (/ (float-pi #q1) 2)) + (random (/ (float-pi #q1) 2))) + for b = (bessel-j #q0.5 x) + for true = (* (/ (sin x) (sqrt x)) (sqrt (/ 2 (float-pi #q1)))) + for result = (check-accuracy 212 b true) + when result + append (list (list (list k x) result))) + nil) commit 2f48eb872f71a8f41187decbc57878cb67fa86de Author: Raymond Toy <toy.raymond@gmail.com> Date: Thu Apr 12 18:27:08 2012 -0700 Add more tests. diff --git a/rt-tests.lisp b/rt-tests.lisp index 712b96c..a4dac02 100644 --- a/rt-tests.lisp +++ b/rt-tests.lisp @@ -1679,3 +1679,22 @@ (check-accuracy 50.73 b true)) nil) +(rt:deftest bessel-j-complex.pos-order.d.2 + (let ((b (bessel-j 1d0 #c(1d0 1))) + (true #c(0.6141603349229036d0 0.3650280288270878d0))) + (check-accuracy 52.51 b true)) + nil) + +(rt:deftest bessel-j-complex.pos-order.d.3 + (let ((b (bessel-j 2d0 #c(1d0 1))) + (true #c(0.0415798869439621d0 0.2473976415133063d0))) + (check-accuracy 50.41 b true)) + nil) + +#+nil +(rt:deftest bessel-j-complex.pos-order.d.4 + (let ((b (bessel-j 2.3d0 #c(1d0 1))) + (true #c(0.1920598664138632d0 -0.5158676904105332d0))) + (check-accuracy 53 b true)) + nil) + ----------------------------------------------------------------------- Summary of changes: rt-tests.lisp | 47 ++++++++++++++++++++++++++++++++++++++++++++++- 1 files changed, 46 insertions(+), 1 deletions(-) hooks/post-receive -- OCT: A portable Lisp implementation for quad-double precision floats