Changeset 5bebc8


Ignore:
Timestamp:
03/11/11 04:05:08 (4 years ago)
Author:
Raymond Toy <toy.raymond@…>
Branches:
master, cvs-ids
Children:
8c4f21
Parents:
5f6c62
Message:

Test carlson-rf and carlson-rd.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • rt-tests.lisp

    r6c9124 r5bebc8  
    835835      (check-accuracy 212 val true))
    836836  nil)
     837
     838(rt:deftest oct.carlson-rf.1d
     839    ;; Rf(0,2,1) = integrate(1/sqrt(1-s^4), s, 0 ,1)
     840    ;;           = 1/4*beta(1/2,1/2)
     841    ;;           = sqrt(%pi)/4*gamma(1/4)/gamma(3/4)
     842    (let ((rf (carlson-rf 0d0 2d0 1d0))
     843          (true 1.31102877714605990523241979494d0))
     844      (check-accuracy 53 rf true))
     845  nil)
     846
     847(rt:deftest oct.carlson-rf.1q
     848    ;; Rf(0,2,1) = integrate(1/sqrt(1-s^4), s, 0 ,1)
     849    (let ((rf (carlson-rf #q0 #q2 #q1))
     850          (true #q1.311028777146059905232419794945559706841377475715811581408410851900395q0))
     851      (check-accuracy 212 rf true))
     852  nil)
     853
     854(rt:deftest oct.carlson-rd.1d
     855    ;; Rd(0,2,1) = 3*integrate(s^2/sqrt(1-s^4), s, 0 ,1)
     856    ;;             = 3*beta(3/4,1/2)/4
     857    ;;             = 3*sqrt(%pi)*gamma(3/4)/gamma(1/4)
     858    (let ((rd (carlson-rd 0d0 2d0 1d0))
     859          (true 1.7972103521033883d0))
     860      (check-accuracy 51 rd true))
     861  nil)
     862
     863(rt:deftest oct.carlson-rd.1q
     864    (let ((rd (carlson-rd #q0 #q2 #q1))
     865          (true #q1.797210352103388311159883738420485817340818994823477337395512429419599q0))
     866      (check-accuracy 212 rd true))
     867  nil)
Note: See TracChangeset for help on using the changeset viewer.