Changeset fb7212 for rttests.lisp
 Timestamp:
 04/17/12 23:15:37 (3 years ago)
 Branches:
 master
 Children:
 de65a5
 Parents:
 581a0a08f04a985424c09b6a7b3661d2eb58c3e9
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

rttests.lisp
r5566bc rfb7212 1624 1624 nil) 1625 1625 1626 (defun besselj1/2 (z) 1627 ;; bessel_j(1/2,z) = sin(x)/sqrt(x)*sqrt(2/pi) 1628 (* (/ (sin z) (sqrt z)) 1629 (sqrt (/ 2 (floatpi z))))) 1630 1626 1631 ;; Bessel J for half integer order and real args 1627 1632 (rt:deftest besselj1/2.d.1 … … 1632 1637 for x = (+ 1 (random (/ pi 2))) 1633 1638 for b = (besselj 0.5d0 x) 1634 for true = ( * (/ (sin x) (sqrt x)) (sqrt (/ 2 pi)))1639 for true = (besselj1/2 x) 1635 1640 for result = (checkaccuracy 48.42 b true) 1636 1641 when result … … 1641 1646 (let* ((x 2.3831631289164497d0) 1642 1647 (b (besselj 0.5d0 x)) 1643 (true ( * (/ (sin x) (sqrt x)) (sqrt (/ 2 pi)))))1648 (true (besselj1/2 x))) 1644 1649 (checkaccuracy 48.42 b true)) 1645 1650 nil) … … 1652 1657 for x = (+ 1 (random (/ (floatpi #q1) 2))) 1653 1658 for b = (besselj #q0.5 x) 1654 for true = ( * (/ (sin x) (sqrt x)) (sqrt (/ 2 (floatpi #q1))))1659 for true = (besselj1/2 x) 1655 1660 for result = (checkaccuracy 169.45 b true) 1656 1661 when result … … 1661 1666 (let* ((x #q1.1288834862545916200627583005758663687705443417892789067029865493882q0) 1662 1667 (b (besselj #q0.5 x)) 1663 (true ( * (/ (sin x) (sqrt x)) (sqrt (/ 2 (floatpi #q1))))))1668 (true (besselj1/2 x))) 1664 1669 (checkaccuracy 182.92 b true)) 1665 1670 nil) … … 1668 1673 (let* ((x #q1.1288834862545916200627583005758663687705443417892789067029865493882q0) 1669 1674 (b (besselj #q0.5 x)) 1670 (true ( * (/ (sin x) (sqrt x)) (sqrt (/ 2 (floatpi #q1))))))1675 (true (besselj1/2 x))) 1671 1676 (checkaccuracy 173.28 b true)) 1672 1677 nil) … … 1675 1680 (let* ((x #q1.0360263937639582798798376485114581552570020473846457752365459851056q0) 1676 1681 (b (besselj #q0.5 x)) 1677 (true ( * (/ (sin x) (sqrt x)) (sqrt (/ 2 (floatpi #q1))))))1682 (true (besselj1/2 x))) 1678 1683 (checkaccuracy 169.45 b true)) 1679 1684 nil) … … 1715 1720 (random (/ pi 2))) 1716 1721 for b = (besselj 0.5d0 x) 1717 for true = ( * (/ (sin x) (sqrt x)) (sqrt (/ 2 pi)))1722 for true = (besselj1/2 x) 1718 1723 for result = (checkaccuracy 49.8 b true) 1719 1724 when result … … 1726 1731 (random (/ (floatpi #q1) 2))) 1727 1732 for b = (besselj #q0.5 x) 1728 for true = ( * (/ (sin x) (sqrt x)) (sqrt (/ 2 (floatpi #q1))))1733 for true = (besselj1/2 x) 1729 1734 for result = (checkaccuracy 212 b true) 1730 1735 when result
Note: See TracChangeset
for help on using the changeset viewer.