Changeset fb7212315ca5f592ef768c2636e7b3cd745473e4 for rt-tests.lisp

Timestamp:

04/17/12 16:15:37 (13 months ago)

Author:

Raymond Toy <rtoy@…>

Children:

de65a54b545af76b6df16f3de36e54763ce97235

Parents:

581a0a08f04a985424c09b6a7b3661d2eb58c3e9

git-committer:

Raymond Toy <rtoy@…> (04/17/12 16:15:37)

Message:

Add function to compute J(1/2,z) and use it in the tests.

Files:

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

rt-tests.lisp

@@ -1618 +1618 @@
   nil)
 
+(defun bessel-j-1/2 (z)
+  ;; bessel_j(1/2,z) = sin(x)/sqrt(x)*sqrt(2/pi)
+  (* (/ (sin z) (sqrt z))
+     (sqrt (/ 2 (float-pi z)))))
+
 ;; Bessel J for half integer order and real args
 (rt:deftest bessel-j-1/2.d.1
@@ -1626 +1631 @@
           for x = (+ 1 (random (/ pi 2)))
           for b = (bessel-j 0.5d0 x)
-          for true = (* (/ (sin x) (sqrt x)) (sqrt (/ 2 pi)))
+          for true = (bessel-j-1/2 x)
           for result = (check-accuracy 48.42 b true)
           when result
@@ -1635 +1640 @@
     (let* ((x 2.3831631289164497d0)
            (b (bessel-j 0.5d0 x))
-           (true (* (/ (sin x) (sqrt x)) (sqrt (/ 2 pi)))))
+           (true (bessel-j-1/2 x)))
       (check-accuracy 48.42 b true))
   nil)
@@ -1646 +1651 @@
           for x = (+ 1 (random (/ (float-pi #q1) 2)))
           for b = (bessel-j #q0.5 x)
-          for true = (* (/ (sin x) (sqrt x)) (sqrt (/ 2 (float-pi #q1))))
+          for true = (bessel-j-1/2 x)
           for result = (check-accuracy 169.45 b true)
           when result
@@ -1655 +1660 @@
     (let* ((x #q1.1288834862545916200627583005758663687705443417892789067029865493882q0)
            (b (bessel-j #q0.5 x))
-           (true (* (/ (sin x) (sqrt x)) (sqrt (/ 2 (float-pi #q1))))))
+           (true (bessel-j-1/2 x)))
       (check-accuracy 182.92 b true))
   nil)
@@ -1662 +1667 @@
     (let* ((x #q1.1288834862545916200627583005758663687705443417892789067029865493882q0)
            (b (bessel-j #q0.5 x))
-           (true (* (/ (sin x) (sqrt x)) (sqrt (/ 2 (float-pi #q1))))))
+           (true (bessel-j-1/2 x)))
       (check-accuracy 173.28 b true))
   nil)
@@ -1669 +1674 @@
     (let* ((x #q1.0360263937639582798798376485114581552570020473846457752365459851056q0)
            (b (bessel-j #q0.5 x))
-           (true (* (/ (sin x) (sqrt x)) (sqrt (/ 2 (float-pi #q1))))))
+           (true (bessel-j-1/2 x)))
       (check-accuracy 169.45 b true))
   nil)
@@ -1709 +1714 @@
                            (random (/ pi 2)))
           for b = (bessel-j 0.5d0 x)
-          for true = (* (/ (sin x) (sqrt x)) (sqrt (/ 2 pi)))
+          for true = (bessel-j-1/2 x)
           for result = (check-accuracy 49.8 b true)
           when result
@@ -1720 +1725 @@
                            (random (/ (float-pi #q1) 2)))
           for b = (bessel-j #q0.5 x)
-          for true = (* (/ (sin x) (sqrt x)) (sqrt (/ 2 (float-pi #q1))))
+          for true = (bessel-j-1/2 x)
           for result = (check-accuracy 212 b true)
           when result