Changeset fb7212315ca5f592ef768c2636e7b3cd745473e4 for rttests.lisp
 Timestamp:
 04/17/12 16:15:37 (2 years ago)
 Children:
 de65a54b545af76b6df16f3de36e54763ce97235
 Parents:
 581a0a08f04a985424c09b6a7b3661d2eb58c3e9
 gitcommitter:
 Raymond Toy <rtoy@…> (04/17/12 16:15:37)
 Files:

 1 modified
Legend:
 Unmodified
 Added
 Removed

rttests.lisp
r5566bc rfb7212 1618 1618 nil) 1619 1619 1620 (defun besselj1/2 (z) 1621 ;; bessel_j(1/2,z) = sin(x)/sqrt(x)*sqrt(2/pi) 1622 (* (/ (sin z) (sqrt z)) 1623 (sqrt (/ 2 (floatpi z))))) 1624 1620 1625 ;; Bessel J for half integer order and real args 1621 1626 (rt:deftest besselj1/2.d.1 … … 1626 1631 for x = (+ 1 (random (/ pi 2))) 1627 1632 for b = (besselj 0.5d0 x) 1628 for true = ( * (/ (sin x) (sqrt x)) (sqrt (/ 2 pi)))1633 for true = (besselj1/2 x) 1629 1634 for result = (checkaccuracy 48.42 b true) 1630 1635 when result … … 1635 1640 (let* ((x 2.3831631289164497d0) 1636 1641 (b (besselj 0.5d0 x)) 1637 (true ( * (/ (sin x) (sqrt x)) (sqrt (/ 2 pi)))))1642 (true (besselj1/2 x))) 1638 1643 (checkaccuracy 48.42 b true)) 1639 1644 nil) … … 1646 1651 for x = (+ 1 (random (/ (floatpi #q1) 2))) 1647 1652 for b = (besselj #q0.5 x) 1648 for true = ( * (/ (sin x) (sqrt x)) (sqrt (/ 2 (floatpi #q1))))1653 for true = (besselj1/2 x) 1649 1654 for result = (checkaccuracy 169.45 b true) 1650 1655 when result … … 1655 1660 (let* ((x #q1.1288834862545916200627583005758663687705443417892789067029865493882q0) 1656 1661 (b (besselj #q0.5 x)) 1657 (true ( * (/ (sin x) (sqrt x)) (sqrt (/ 2 (floatpi #q1))))))1662 (true (besselj1/2 x))) 1658 1663 (checkaccuracy 182.92 b true)) 1659 1664 nil) … … 1662 1667 (let* ((x #q1.1288834862545916200627583005758663687705443417892789067029865493882q0) 1663 1668 (b (besselj #q0.5 x)) 1664 (true ( * (/ (sin x) (sqrt x)) (sqrt (/ 2 (floatpi #q1))))))1669 (true (besselj1/2 x))) 1665 1670 (checkaccuracy 173.28 b true)) 1666 1671 nil) … … 1669 1674 (let* ((x #q1.0360263937639582798798376485114581552570020473846457752365459851056q0) 1670 1675 (b (besselj #q0.5 x)) 1671 (true ( * (/ (sin x) (sqrt x)) (sqrt (/ 2 (floatpi #q1))))))1676 (true (besselj1/2 x))) 1672 1677 (checkaccuracy 169.45 b true)) 1673 1678 nil) … … 1709 1714 (random (/ pi 2))) 1710 1715 for b = (besselj 0.5d0 x) 1711 for true = ( * (/ (sin x) (sqrt x)) (sqrt (/ 2 pi)))1716 for true = (besselj1/2 x) 1712 1717 for result = (checkaccuracy 49.8 b true) 1713 1718 when result … … 1720 1725 (random (/ (floatpi #q1) 2))) 1721 1726 for b = (besselj #q0.5 x) 1722 for true = ( * (/ (sin x) (sqrt x)) (sqrt (/ 2 (floatpi #q1))))1727 for true = (besselj1/2 x) 1723 1728 for result = (checkaccuracy 212 b true) 1724 1729 when result