Commit 5bebc8d76a83b85cdba6d803d026aae741b6e7c8

Authored by Raymond Toy
1 parent 5f6c6280

Test carlson-rf and carlson-rd.

Showing 1 changed file with 31 additions and 0 deletions
rt-tests.lisp
... ... @@ -833,3 +833,34 @@
833 833 (val (jacobi-dn ek #q.5)))
834 834 (check-accuracy 212 val true))
835 835 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)
... ...