Changeset 5566bc397937c2e8979ee6847e8f59f279f1f643

Show
Ignore:
Timestamp:
04/12/12 23:24:03 (2 years 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:
1 modified

Legend:

Unmodified
Added
Removed
  • rt-tests.lisp

    r2f48eb r5566bc  
    16741674 
    16751675;; Bessel J for complex args 
    1676 (rt:deftest bessel-j-complex.pos-order.d.1 
     1676(rt:deftest bessel-j-complex-arg.d.1 
    16771677    (let ((b (bessel-j 0d0 #c(1d0 1))) 
    16781678          (true #c(0.9376084768060293d0 -0.4965299476091221d0))) 
     
    16801680  nil) 
    16811681 
    1682 (rt:deftest bessel-j-complex.pos-order.d.2 
     1682(rt:deftest bessel-j-complex-arg.d.2 
    16831683    (let ((b (bessel-j 1d0 #c(1d0 1))) 
    16841684          (true #c(0.6141603349229036d0 0.3650280288270878d0))) 
     
    16861686  nil) 
    16871687 
    1688 (rt:deftest bessel-j-complex.pos-order.d.3 
     1688(rt:deftest bessel-j-complex-arg.d.3 
    16891689    (let ((b (bessel-j 2d0 #c(1d0 1))) 
    16901690          (true #c(0.0415798869439621d0 0.2473976415133063d0))) 
     
    16921692  nil) 
    16931693 
    1694 #+nil 
    1695 (rt:deftest bessel-j-complex.pos-order.d.4 
     1694(rt:deftest bessel-j-complex-arg.d.4 
    16961695    (let ((b (bessel-j 2.3d0 #c(1d0 1))) 
     1696          (true #c(-0.0141615213034667d0 0.1677798241687935d0))) 
     1697      (check-accuracy 48.56 b true)) 
     1698  nil) 
     1699 
     1700(rt:deftest bessel-j-complex-arg.d.5 
     1701    (let ((b (bessel-j -2.3d0 #c(1d0 1))) 
    16971702          (true #c(0.1920598664138632d0 -0.5158676904105332d0))) 
    1698       (check-accuracy 53 b true)) 
    1699   nil) 
    1700  
     1703      (check-accuracy 50.97 b true)) 
     1704  nil) 
     1705 
     1706(rt:deftest bessel-j-1/2-complex.d.1 
     1707    (loop for k from 0 below 10 
     1708          for x = (complex (random (/ pi 2)) 
     1709                           (random (/ pi 2))) 
     1710          for b = (bessel-j 0.5d0 x) 
     1711          for true = (* (/ (sin x) (sqrt x)) (sqrt (/ 2 pi))) 
     1712          for result = (check-accuracy 49.8 b true) 
     1713          when result 
     1714            append (list (list (list k x) result))) 
     1715  nil) 
     1716 
     1717(rt:deftest bessel-j-1/2-complex.q.1 
     1718    (loop for k from 0 below 10 
     1719          for x = (complex (random (/ (float-pi #q1) 2)) 
     1720                           (random (/ (float-pi #q1) 2))) 
     1721          for b = (bessel-j #q0.5 x) 
     1722          for true = (* (/ (sin x) (sqrt x)) (sqrt (/ 2 (float-pi #q1)))) 
     1723          for result = (check-accuracy 212 b true) 
     1724          when result 
     1725            append (list (list (list k x) result))) 
     1726  nil)