[cl-gsl-cvs] CVS update: cl-gsl/test/test-sf.lisp
Edgar Denny
edenny at common-lisp.net
Tue Mar 8 03:34:18 UTC 2005
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+)
More information about the Cl-gsl-cvs
mailing list