[cl-gsl-cvs] CVS update: cl-gsl/test/test-sf.lisp

Update of /project/cl-gsl/cvsroot/cl-gsl/test In directory common-lisp.net:/tmp/cvs-serv10932

Modified Files: test-sf.lisp Log Message: More unit tests added.

Date: Tue Mar 8 04:34:17 2005 Author: edenny

Index: cl-gsl/test/test-sf.lisp diff -u cl-gsl/test/test-sf.lisp:1.2 cl-gsl/test/test-sf.lisp:1.3 --- cl-gsl/test/test-sf.lisp:1.2 Mon Mar 7 03:08:25 2005 +++ cl-gsl/test/test-sf.lisp Tue Mar 8 04:34:16 2005 @@ -26,6 +26,8 @@ (defconstant +tol-1+ (/ +tol0+ 10.0d0)) (defconstant +tol+2+ (* +tol0+ 100.0d0)) (defconstant +tol-2+ (/ +tol0+ 100.0d0)) +(defconstant +tol+3+ (* +tol0+ 1000.0d0)) +(defconstant +tol-3+ (/ +tol0+ 1000.0d0))

(defconstant +sf+ "sf")

@@ -142,103 +144,322 @@ -4.073155089071828216d0 +tol0+ gsl:+success+)

;; bessel-c-j0 +(sf-deftest "bessel-c-j0" (0.1d0) 0.99750156206604003230d0 +tol0+) + ;; bessel-c-j0-e +(sf-deftest-e "bessel-c-j0-e" (0.1d0) + 0.99750156206604003230d0 +tol0+ gsl:+success+) + ;; bessel-c-j1 +(sf-deftest "bessel-c-j1" (0.1d0) 0.04993752603624199756d0 +tol0+) + ;; bessel-c-j1-e +(sf-deftest-e "bessel-c-j1-e" (0.1d0) + 0.04993752603624199756d0 +tol0+ gsl:+success+) + ;; bessel-c-jn +(sf-deftest "bessel-c-jn" (10 20.0d0) 0.18648255802394508321d0 +tol+1+) + ;; bessel-c-jn-e +(sf-deftest-e "bessel-c-jn-e" (10 20.0d0) + 0.18648255802394508321d0 +tol+1+ gsl:+success+) + ;; bessel-c-jn-array ;; bessel-c-y0 +(sf-deftest "bessel-c-y0" (0.1d0) -1.5342386513503668441d0 +tol0+) + ;; bessel-c-y0-e +(sf-deftest-e "bessel-c-y0-e" (0.1d0) + -1.5342386513503668441d0 +tol0+ gsl:+success+) + ;; bessel-c-y1 +(sf-deftest "bessel-c-y1" (0.1d0) -6.45895109470202698800d0 +tol0+) + ;; bessel-c-y1-e +(sf-deftest-e "bessel-c-y1-e" (0.1d0) + -6.45895109470202698800d0 +tol0+ gsl:+success+) + ;; bessel-c-yn +(sf-deftest "bessel-c-jn" (10 20.0d0) 0.18648255802394508321d0 +tol0+) + ;; bessel-c-yn-e +(sf-deftest-e "bessel-c-jn-e" (10 20.0d0) + 0.18648255802394508321d0 +tol0+ gsl:+success+) + ;; bessel-c-yn-array ;; bessel-c-i0 +(sf-deftest "bessel-c-i0" (0.1d0) 1.0025015629340956014d0 +tol0+) + ;; bessel-c-i0-e +(sf-deftest-e "bessel-c-i0-e" (0.1d0) + 1.0025015629340956014d0 +tol0+ gsl:+success+) + ;; bessel-c-i1 +(sf-deftest "bessel-c-i1" (0.1d0) 0.05006252604709269211d0 +tol0+) + ;; bessel-c-i1-e +(sf-deftest-e "bessel-c-i1-e" (0.1d0) + 0.05006252604709269211d0 +tol0+ gsl:+success+) + ;; bessel-c-in +(sf-deftest "bessel-c-in" (5 2.0d0) 0.009825679323131702321d0 +tol0+) + ;; bessel-c-in-e +(sf-deftest-e "bessel-c-in-e" (5 2.0d0) + 0.009825679323131702321d0 +tol0+ gsl:+success+) + ;; bessel-c-in-array ;; bessel-c-i0-scaled +(sf-deftest "bessel-c-i0-scaled" (0.1d0) 0.90710092578230109640d0 +tol0+) + ;; bessel-c-i0-scaled-e +(sf-deftest-e "bessel-c-i0-scaled-e" (0.1d0) + 0.90710092578230109640d0 +tol0+ gsl:+success+) + ;; bessel-c-i1-scaled +(sf-deftest "bessel-c-i1-scaled" (0.1d0) 0.04529844680880932501d0 +tol0+) + ;; bessel-c-i1-scaled-e +(sf-deftest-e "bessel-c-i1-scaled-e" (0.1d0) + 0.04529844680880932501d0 +tol0+ gsl:+success+) + ;; bessel-c-in-scaled +(sf-deftest "bessel-c-in-scaled" (5 2.0d0) 0.0013297610941881578142d0 +tol0+) + ;; bessel-c-in-scaled-e +(sf-deftest-e "bessel-c-in-scaled-e" (5 2.0d0) + 0.0013297610941881578142d0 +tol0+ gsl:+success+) + ;; bessel-c-in-scaled-array ;; bessel-c-k0 +(sf-deftest "bessel-c-k0" (2.0d0) 0.11389387274953343565d0 +tol0+) + ;; bessel-c-k0-e +(sf-deftest-e "bessel-c-k0-e" (2.0d0) + 0.11389387274953343565d0 +tol0+ gsl:+success+) + ;; bessel-c-k1 +(sf-deftest "bessel-c-k1" (2.0d0) 0.13986588181652242728d0 +tol0+) + ;; bessel-c-k1-e +(sf-deftest-e "bessel-c-k1-e" (2.0d0) + 0.13986588181652242728d0 +tol0+ gsl:+success+) + ;; bessel-c-kn +(sf-deftest "bessel-c-kn" (5 2.0d0) 9.431049100596467443d0 +tol+2+) + ;; bessel-c-kn-e +(sf-deftest-e "bessel-c-kn-e" (5 2.0d0) + 9.431049100596467443d0 +tol+2+ gsl:+success+) + ;; bessel-c-kn-array ;; bessel-c-k0-scaled +(sf-deftest "bessel-c-k0-scaled" (2.0d0) 0.8415682150707714179d0 +tol+1+) + ;; bessel-c-k0-scaled-e +(sf-deftest-e "bessel-c-k0-scaled-e" (2.0d0) + 0.8415682150707714179d0 +tol+1+ gsl:+success+) + ;; bessel-c-k1-scaled +(sf-deftest "bessel-c-k1-scaled" (2.0d0) 1.0334768470686885732d0 +tol0+) + ;; bessel-c-k1-scaled-e +(sf-deftest-e "bessel-c-k1-scaled-e" (2.0d0) + 1.0334768470686885732d0 +tol0+ gsl:+success+) + ;; bessel-c-kn-scaled +(sf-deftest "bessel-c-kn-scaled" (5 2.0d0) 69.68655087607675118d0 +tol+3+) + ;; bessel-c-kn-scaled-e +(sf-deftest-e "bessel-c-kn-scaled-e" (5 2.0d0) + 69.68655087607675118d0 +tol+3+ gsl:+success+) + ;; bessel-c-kn-scaled-array ;; bessel-s-j0 +(sf-deftest "bessel-s-j0" (1.0d0) 0.84147098480789650670d0 +tol0+) + ;; bessel-s-j0-e +(sf-deftest-e "bessel-s-j0-e" (1.0d0) + 0.84147098480789650670d0 +tol0+ gsl:+success+) + ;; bessel-s-j1 +(sf-deftest "bessel-s-j1" (1.0d0) 0.30116867893975678925d0 +tol0+) + ;; bessel-s-j1-e +(sf-deftest-e "bessel-s-j1-e" (1.0d0) + 0.30116867893975678925d0 +tol0+ gsl:+success+) + ;; bessel-s-j2 +(sf-deftest "bessel-s-j2" (1.0d0) 0.06203505201137386110d0 +tol0+) + ;; bessel-s-j2-e +(sf-deftest-e "bessel-s-j2-e" (1.0d0) + 0.06203505201137386110d0 +tol0+ gsl:+success+) + ;; bessel-s-jl -;; bessel-s-ji-e +(sf-deftest "bessel-s-jl" (10 10.0d0) 0.06460515449256426427d0 +tol0+) + +;; bessel-s-jl-e +(sf-deftest-e "bessel-s-jl-e" (10 10.0d0) + 0.06460515449256426427d0 +tol0+ gsl:+success+) + ;; bessel-s-jl-array ;; bessel-s-y0 +(sf-deftest "bessel-s-y0" (1.0d0) -0.5403023058681397174d0 +tol0+) + ;; bessel-s-y0-e +(sf-deftest-e "bessel-s-y0-e" (1.0d0) + -0.5403023058681397174d0 +tol0+ gsl:+success+) + ;; bessel-s-y1 +(sf-deftest "bessel-s-y1" (1.0d0) -1.3817732906760362241d0 +tol0+) + ;; bessel-s-y1-e +(sf-deftest-e "bessel-s-y1-e" (1.0d0) + -1.3817732906760362241d0 +tol0+ gsl:+success+) + ;; bessel-s-y2 +(sf-deftest "bessel-s-y2" (1.0d0) -3.605017566159968955d0 +tol0+) + ;; bessel-s-y2-e +(sf-deftest-e "bessel-s-y2-e" (1.0d0) + -3.605017566159968955d0 +tol0+ gsl:+success+) + ;; bessel-s-yl +(sf-deftest "bessel-s-yl" (0 1.0d0) -0.54030230586813972d0 +tol0+) + ;; bessel-s-yl-e +(sf-deftest-e "bessel-s-yl-e" (0 1.0d0) + -0.54030230586813972d0 +tol0+ gsl:+success+) + ;; bessel-s-yl-array ;; bessel-s-i0-scaled +(sf-deftest "bessel-s-i0-scaled" (2.0d0) 0.24542109027781645493d0 +tol0+) + ;; bessel-s-i0-scaled-e +(sf-deftest-e "bessel-s-i0-scaled-e" (2.0d0) + 0.24542109027781645493d0 +tol0+ gsl:+success+) + ;; bessel-s-i1-scaled +(sf-deftest "bessel-s-i1-scaled" (2.0d0) 0.131868364583275317610d0 +tol0+) + ;; bessel-s-i1-scaled-e +(sf-deftest-e "bessel-s-i1-scaled-e" (2.0d0) + 0.131868364583275317610d0 +tol0+ gsl:+success+) + ;; bessel-s-i2-scaled +(sf-deftest "bessel-s-i2-scaled" (2.0d0) 0.0476185434029034785100d0 +tol0+) + ;; bessel-s-i2-scaled-e +(sf-deftest-e "bessel-s-i2-scaled-e" (2.0d0) + 0.0476185434029034785100d0 +tol0+ gsl:+success+) + ;; bessel-s-il-scaled +(sf-deftest "bessel-s-il-scaled" (5 2.0d0) 0.0004851564602127540059d0 +tol0+) + ;; bessel-s-il-scaled-e +(sf-deftest-e "bessel-s-il-scaled-e" (5 2.0d0) + 0.0004851564602127540059d0 +tol0+ gsl:+success+) + ;; bessel-s-il-scaled-array ;; bessel-s-k0-scaled +(sf-deftest "bessel-s-k0-scaled" (2.0d0) 0.7853981633974483096d0 +tol0+) + ;; bessel-s-k0-scaled-e +(sf-deftest-e "bessel-s-k0-scaled-e" (2.0d0) + 0.7853981633974483096d0 +tol0+ gsl:+success+) + ;; bessel-s-k1-scaled +(sf-deftest "bessel-s-k1-scaled" (2.0d0) 1.1780972450961724644d0 +tol0+) + ;; bessel-s-k1-scaled-e +(sf-deftest-e "bessel-s-k1-scaled-e" (2.0d0) + 1.1780972450961724644d0 +tol0+ gsl:+success+) + ;; bessel-s-k2-scaled +(sf-deftest "bessel-s-k2-scaled" (2.0d0) 2.5525440310417070063d0 +tol0+) + ;; bessel-s-k2-scaled-e +(sf-deftest-e "bessel-s-k2-scaled-e" (2.0d0) + 2.5525440310417070063d0 +tol0+ gsl:+success+) + ;; bessel-s-kl-scaled +(sf-deftest "bessel-s-kl-scaled" (5 2.0d0) 138.10735829492005119d0 +tol0+) + ;; bessel-s-kl-scaled-e +(sf-deftest-e "bessel-s-kl-scaled-e" (5 2.0d0) + 138.10735829492005119d0 +tol0+ gsl:+success+) + ;; bessel-s-kl-scaled-array ;; bessel-c-jnu +(sf-deftest "bessel-c-jnu" (1.0d0 1.0d0) 0.4400505857449335160d0 +tol0+) + ;; bessel-c-jnu-e +(sf-deftest-e "bessel-c-jnu-e" (1.0d0 1.0d0) + 0.4400505857449335160d0 +tol0+ gsl:+success+) + ;; bessel-c-ynu +(sf-deftest "bessel-c-ynu" (1.0d0 1.0d0) -0.7812128213002887165d0 +tol+1+) + ;; bessel-c-ynu-e +(sf-deftest-e "bessel-c-ynu-e" (1.0d0 1.0d0) + -0.7812128213002887165d0 +tol+1+ gsl:+success+) + ;; bessel-c-inu +(sf-deftest "bessel-c-inu" (1.0d0 1.0d0) 0.5651591039924850272d0 +tol0+) + ;; bessel-c-inu-e +(sf-deftest-e "bessel-c-inu-e" (1.0d0 1.0d0) + 0.5651591039924850272d0 +tol0+ gsl:+success+) + ;; bessel-c-inu-scaled +(sf-deftest "bessel-c-inu-scaled" (1.0d0 1.0d0) 0.20791041534970844887d0 +tol0+) + ;; bessel-c-inu-scaled-e +(sf-deftest-e "bessel-c-inu-scaled-e" (1.0d0 1.0d0) + 0.20791041534970844887d0 +tol0+ gsl:+success+) + ;; bessel-c-knu +(sf-deftest "bessel-c-knu" (1.0d0 1.0d0) 0.6019072301972345747d0 +tol0+) + ;; bessel-c-knu-e +(sf-deftest-e "bessel-c-knu-e" (1.0d0 1.0d0) + 0.6019072301972345747d0 +tol0+ gsl:+success+) + ;; bessel-c-lnknu +(sf-deftest "bessel-c-lnknu" (1.0d0 1.0d0) -0.5076519482107523309d0 +tol0+) + ;; bessel-c-lnknu-e +(sf-deftest-e "bessel-c-lnknu-e" (1.0d0 1.0d0) + -0.5076519482107523309d0 +tol0+ gsl:+success+) + ;; bessel-c-knu-scaled +(sf-deftest "bessel-c-knu-scaled" (1.0d0 1.0d0) 1.6361534862632582465d0 +tol0+) + ;; bessel-c-knu-scaled-e +(sf-deftest-e "bessel-c-knu-scaled-e" (1.0d0 1.0d0) + 1.6361534862632582465d0 +tol0+ gsl:+success+) + ;; bessel-c-zero-j0 +(sf-deftest "bessel-c-zero-j0" (1) 2.404825557695771d0 +tol+1+) + ;; bessel-c-zero-j0-e +(sf-deftest-e "bessel-c-zero-j0-e" (1) + 2.404825557695771d0 +tol+1+ gsl:+success+) + ;; bessel-c-zero-j1 +(sf-deftest "bessel-c-zero-j1" (1) 3.831705970207512d0 +tol+2+) + ;; bessel-c-zero-j1-e +(sf-deftest-e "bessel-c-zero-j1-e" (1) + 3.831705970207512d0 +tol+2+ gsl:+success+) + ;; bessel-c-zero-jnu +(sf-deftest "bessel-c-zero-jnu" (0.0d0 1) 2.404825557695771d0 +tol+2+) + ;; bessel-c-zero-jnu-e +(sf-deftest-e "bessel-c-zero-jnu-e" (0.0d0 1) + 2.404825557695771d0 +tol+2+ gsl:+success+)

;; clausen (sf-deftest "clausen" ((/ gsl-math:+pi+ 20.0d0)) 0.4478882448133546d0 +tol0+) @@ -634,28 +855,77 @@ (sf-deftest-e "gegenpoly-n-e" (10 1.0d0 1.0d0) 11.0d0 +tol0+ gsl:+success+)

;; gegenpoly-array -;; hyperg-0F1 -;; hyperg-0F1-e -;; hyperg-1F1-int -;; hyperg-1F1-int-e -;; hyperg-1F1 -;; hyperg-1F1-e -;; hyperg-U-int -;; hyperg-U-int-e -;; hyperg-U-int-e10-e -;; hyperg-U -;; hyperg-U-e -;; hyperg-U-e10-e -;; hyperg-2F1 -;; hyperg-2F1-e -;; hyperg-2F1-conj -;; hyperg-2F1-conj-e -;; hyperg-2F1-renorm -;; hyperg-2F1-renorm-e -;; hyperg-2F1-conj-renorm -;; hyperg-2F1-conj-renorm-e -;; hyperg-2F0 -;; hyperg-2F0-e +;; hyperg-0f1 +(sf-deftest "hyperg-0f1" (1.0d0 0.5d0) 1.5660829297563505373d0 +tol+1+) + +;; hyperg-0f1-e +(sf-deftest-e "hyperg-0f1-e" (1.0d0 0.5d0) + 1.5660829297563505373d0 +tol+1+ gsl:+success+) + +;; hyperg-1f1-int +(sf-deftest "hyperg-1f1-int" (1 1 0.5d0) 1.6487212707001281468d0 +tol0+) + +;; hyperg-1f1-int-e +(sf-deftest-e "hyperg-1f1-int-e" (1 1 0.5d0) + 1.6487212707001281468d0 +tol0+ gsl:+success+) + +;; hyperg-1f1 +(sf-deftest "hyperg-1f1" (1.5d0 2.5d0 1.0d0) 1.8834451238277954398d0 +tol0+) + +;; hyperg-1f1-e +(sf-deftest-e "hyperg-1f1-e" (1.5d0 2.5d0 1.0d0) + 1.8834451238277954398d0 +tol0+ gsl:+success+) + +;; hyperg-u-int +(sf-deftest "hyperg-u-int" (1 1 2.0d0) 0.3613286168882225847d0 +tol0+) + +;; hyperg-u-int-e +(sf-deftest-e "hyperg-u-int-e" (1 1 2.0d0) + 0.3613286168882225847d0 +tol0+ gsl:+success+) + +;; hyperg-u-int-e10-e +;; hyperg-u +(sf-deftest "hyperg-u" (0.0001d0 1.0d0 1.0d0) 0.9999999925484179084d0 +tol+2+) + +;; hyperg-u-e +(sf-deftest-e "hyperg-u-e" (0.0001d0 1.0d0 1.0d0) + 0.9999999925484179084d0 +tol+2+ gsl:+success+) + +;; hyperg-u-e10-e +;; hyperg-2f1 +(sf-deftest "hyperg-2f1" (8.0d0 -8.0d0 1.0d0 0.5d0) + 0.13671875d0 +tol0+) + +;; hyperg-2f1-e +(sf-deftest-e "hyperg-2f1-e" (8.0d0 -8.0d0 1.0d0 0.5d0) + 0.13671875d0 +tol0+ gsl:+success+) + +;; hyperg-2f1-conj +(sf-deftest "hyperg-2f1-conj" (1.0d0 1.0d0 1.0d0 0.5d0) + 3.352857095662929028d0 +tol+2+) + +;; hyperg-2f1-conj-e +(sf-deftest-e "hyperg-2f1-conj-e" (1.0d0 1.0d0 1.0d0 0.5d0) + 3.352857095662929028d0 +tol+2+ gsl:+success+) + +;; hyperg-2f1-renorm +(sf-deftest "hyperg-2f1-renorm" (8.0d0 -8.0d0 1.0d0 0.5d0) + 0.13671875d0 +tol0+) + +;; hyperg-2f1-renorm-e +(sf-deftest-e "hyperg-2f1-renorm-e" (8.0d0 -8.0d0 1.0d0 0.5d0) + 0.13671875d0 +tol0+ gsl:+success+) + +;; hyperg-2f1-conj-renorm +(sf-deftest "hyperg-2f1-conj-renorm" (9.0d0 9.0d0 -1.5d0 -0.99d0) + 0.10834020229476124874d0 +tol0+) + +;; hyperg-2f1-conj-renorm-e +(sf-deftest-e "hyperg-2f1-conj-renorm-e" (9.0d0 9.0d0 -1.5d0 -0.99d0) + 0.10834020229476124874d0 +tol0+ gsl:+success+) + +;; hyperg-2f0 +;; hyperg-2f0-e ;; laguerre-1 (sf-deftest "laguerre-1" (0.5d0 -1.0d0) 2.5d0 +tol0+)