Changeset 5bebc8d76a83b85cdba6d803d026aae741b6e7c8

Show
Ignore:
Timestamp:
03/10/11 20:05:08 (3 years ago)
Author:
Raymond Toy <toy.raymond@…>
Children:
8c4f2161315be456cf82774bec04ed43f578ffd8
Parents:
5f6c628014268b683e61baae5470a56b078d1c16
git-committer:
Raymond Toy <toy.raymond@…> (03/10/11 20:05:08)
Message:

Test carlson-rf and carlson-rd.

Files:
1 modified

Legend:

Unmodified
Added
Removed
  • rt-tests.lisp

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