Changeset 5566bc397937c2e8979ee6847e8f59f279f1f643

Timestamp:

04/12/12 23:24:03 (14 months ago)

Author:

Raymond Toy <toy.raymond@…>

Children:

78801d6381aaaf4f21967582680f26889582db60

Parents:

2f48eb872f71a8f41187decbc57878cb67fa86de

git-committer:

Raymond Toy <toy.raymond@…> (04/12/12 23:24:03)

Message:

Add more tests for Bessel J. Not all of them pass.

Files:

• rt-tests.lisp (modified) (4 diffs)

rt-tests.lisp

@@ -1674 +1674 @@
 
 ;; 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)))
@@ -1680 +1680 @@
   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)))
@@ -1686 +1686 @@
   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)))
@@ -1692 +1692 @@
   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))
-  nil)
-
+      (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)