From fe339ae42f89d9fb36f823ef05d4a0922335e044 Mon Sep 17 00:00:00 2001 From: Liam Healy <liam@thinkpad.local> Date: Sun, 14 Sep 2008 22:33:07 -0400 Subject: [PATCH] Restore some tests Multiple fixes to restore tests: change to #'cl-array from #'data, use #'a macro, comment out a few tests do defer repair until later. Two errors in definitions of functions in histogram/updating-accessing.lisp fixed. Some uncommented-out test failures remain. --- histogram/histogram.lisp | 5 ++- histogram/updating-accessing.lisp | 12 +++---- init/generate-tests.lisp | 7 ++-- numerical-integration.lisp | 4 +-- polynomial.lisp | 18 +++++------ random/dirichlet.lisp | 26 +++++++-------- random/discrete.lisp | 16 ++++----- random/multinomial.lisp | 18 +++++------ random/quasi.lisp | 6 ++-- random/shuffling-sampling.lisp | 4 +-- special-functions/bessel.lisp | 54 +++++++++++++++---------------- special-functions/coulomb.lisp | 22 ++++++------- special-functions/gegenbauer.lisp | 6 ++-- special-functions/legendre.lisp | 26 +++++++-------- 14 files changed, 112 insertions(+), 112 deletions(-) diff --git a/histogram/histogram.lisp b/histogram/histogram.lisp index a889a4bd..5a5a98ee 100644 --- a/histogram/histogram.lisp +++ b/histogram/histogram.lisp @@ -1,6 +1,6 @@ ;; The histogram structure ;; Liam Healy, Mon Jan 1 2007 - 11:32 -;; Time-stamp: <2008-08-31 11:03:07EDT histogram.lisp> +;; Time-stamp: <2008-09-14 22:06:19EDT histogram.lisp> ;; $Id$ (in-package :gsl) @@ -32,8 +32,7 @@ pointer :ENOMEM (format nil "for ~a." - (with-output-to-string (stream) - (print-data-object object nil stream)))) + (with-output-to-string (stream) (princ object stream)))) (setf (mpointer object) pointer)) (defmfun alloc-histo-1 (object) diff --git a/histogram/updating-accessing.lisp b/histogram/updating-accessing.lisp index 89164933..1b6cffe2 100644 --- a/histogram/updating-accessing.lisp +++ b/histogram/updating-accessing.lisp @@ -1,11 +1,11 @@ ;; Updating and accessing histogram elements. ;; Liam Healy, Mon Jan 1 2007 - 14:43 -;; Time-stamp: <2008-08-28 21:54:43EDT updating-accessing.lisp> +;; Time-stamp: <2008-09-14 22:04:40EDT updating-accessing.lisp> ;; $Id$ (in-package :gsl) -(defmfun increment-fix (histogram value &optional weight) +(defmfun increment (histogram value &optional weight) ("gsl_histogram_increment" "gsl_histogram_accumulate") ((((mpointer histogram) :pointer) (value :double)) (((mpointer histogram) :pointer) (value :double) (weight :double))) @@ -26,7 +26,7 @@ (defmfun maref ((histogram histogram) i &optional i2 element-type) "gsl_histogram_get" - (((mpointer histogram) :pointer) ((first i) sizet)) + (((mpointer histogram) :pointer) (i sizet)) :definition :method :c-return :double :documentation ; FDL @@ -130,7 +130,7 @@ (set-ranges-uniform histo 0.0d0 10.0d0) (increment histo 2.7d0) (increment histo 6.9d0 2.0d0) - (values (gsl-min-range histo) (gsl-max-range histo))) + (values (min-range histo) (max-range histo))) (letm ((histo (histogram 10))) (set-ranges-uniform histo 0.0d0 10.0d0) (increment histo 2.7d0) @@ -187,8 +187,8 @@ (SET-RANGES-UNIFORM HISTO 0.0d0 10.0d0) (INCREMENT HISTO 2.7d0) (INCREMENT HISTO 6.9d0 2.0d0) - (VALUES (GSL-MIN-RANGE HISTO) - (GSL-MAX-RANGE HISTO))))) + (VALUES (MIN-RANGE HISTO) + (MAX-RANGE HISTO))))) (LISP-UNIT::ASSERT-NUMERICAL-EQUAL (LIST 10) (MULTIPLE-VALUE-LIST diff --git a/init/generate-tests.lisp b/init/generate-tests.lisp index 6ded37c3..5bf4a0af 100644 --- a/init/generate-tests.lisp +++ b/init/generate-tests.lisp @@ -1,6 +1,6 @@ ;; Make tests and examples ;; Liam Healy 2008-09-07 21:00:48EDT generate-tests.lisp -;; Time-stamp: <2008-09-14 18:14:17EDT generate-tests.lisp> +;; Time-stamp: <2008-09-14 18:22:22EDT generate-tests.lisp> ;; $Id: $ (in-package :gsl) @@ -103,6 +103,7 @@ +#| ;;;;;;;;;; ;;; doesn't work for complex because: @@ -111,7 +112,7 @@ '(letm ((v1 (vector-default 3 t)) (v2 (vector-default 3))) (set-all v1 (scalar-default 2)) - (m+ v1 v2)))) + (m+ v1 v2))) ;;; doesn't work for complex because GSL @@ -119,3 +120,5 @@ '(letm ((v1 (vector-default 3)) (v2 (vector-default 3))) (m+ v1 v2))) + +|# diff --git a/numerical-integration.lisp b/numerical-integration.lisp index 7454f738..d3021e71 100644 --- a/numerical-integration.lisp +++ b/numerical-integration.lisp @@ -1,6 +1,6 @@ ;; Numerical integration ;; Liam Healy, Wed Jul 5 2006 - 23:14 -;; Time-stamp: <2008-08-23 19:18:03EDT numerical-integration.lisp> +;; Time-stamp: <2008-09-14 21:55:47EDT numerical-integration.lisp> ;; $Id$ ;;; To do: QAWS, QAWO, QAWF, more tests @@ -233,7 +233,6 @@ (integration-QAG one-sine 0.0d0 pi :gauss15 20 ws)) (letm ((ws (integration-workspace 20))) (integration-QAG one-sine 0.0d0 pi :gauss15 50 ws))) -|# (LISP-UNIT:DEFINE-TEST NUMERICAL-INTEGRATION (LISP-UNIT::ASSERT-NUMERICAL-EQUAL @@ -251,3 +250,4 @@ (LETM ((WS (INTEGRATION-WORKSPACE 20))) (INTEGRATION-QAG ONE-SINE 0.0d0 PI :GAUSS15 50 WS)))) +|# diff --git a/polynomial.lisp b/polynomial.lisp index c04bc3e2..41bde745 100644 --- a/polynomial.lisp +++ b/polynomial.lisp @@ -1,6 +1,6 @@ ;; Polynomials ;; Liam Healy, Tue Mar 21 2006 - 18:33 -;; Time-stamp: <2008-07-06 16:57:52EDT polynomial.lisp> +;; Time-stamp: <2008-09-14 22:01:10EDT polynomial.lisp> ;; $Id$ (in-package :gsl) @@ -173,8 +173,8 @@ #| (make-tests polynomial - (letm ((xa (vector-double-float #(0.0d0 1.0d0 2.0d0 3.0d0))) - (ya (vector-double-float #(2.5d0 7.2d0 32.7d0 91.0d0))) + (letm ((xa (vector-double-float (a 0.0d0 1.0d0 2.0d0 3.0d0))) + (ya (vector-double-float (a 2.5d0 7.2d0 32.7d0 91.0d0))) (dd (vector-double-float 4))) (divided-difference dd xa ya) (list @@ -182,7 +182,7 @@ (polynomial-eval-divided-difference dd xa 1.0d0) (polynomial-eval-divided-difference dd xa 2.0d0) (polynomial-eval-divided-difference dd xa 3.0d0))) - (letm ((vec (vector-double-float #(1.0d0 2.0d0 3.0d0)))) + (letm ((vec (vector-double-float (a 1.0d0 2.0d0 3.0d0)))) (polynomial-eval vec -1.0d0)) (solve-quadratic 1.0d0 0.0d0 1.0d0) (solve-quadratic 1.0d0 -2.0d0 1.0d0) @@ -190,7 +190,7 @@ (solve-cubic -6.0d0 -13.0d0 42.0d0) (solve-cubic-complex -1.0d0 1.0d0 -1.0d0) ;; Example from GSL manual - (polynomial-solve #(-1.0d0 0.0d0 0.0d0 0.0d0 0.0d0 1.0d0))) + (polynomial-solve (a double-float -1.0d0 0.0d0 0.0d0 0.0d0 0.0d0 1.0d0))) |# (LISP-UNIT:DEFINE-TEST POLYNOMIAL @@ -198,8 +198,8 @@ (LIST (LIST 2.5d0 7.2d0 32.7d0 91.0d0)) (MULTIPLE-VALUE-LIST (LETM - ((XA (VECTOR-DOUBLE-FLOAT #(0.0d0 1.0d0 2.0d0 3.0d0))) - (YA (VECTOR-DOUBLE-FLOAT #(2.5d0 7.2d0 32.7d0 91.0d0))) + ((XA (VECTOR-DOUBLE-FLOAT (A 0.0d0 1.0d0 2.0d0 3.0d0))) + (YA (VECTOR-DOUBLE-FLOAT (A 2.5d0 7.2d0 32.7d0 91.0d0))) (DD (VECTOR-DOUBLE-FLOAT 4))) (DIVIDED-DIFFERENCE DD XA YA) (LIST @@ -210,7 +210,7 @@ (LISP-UNIT::ASSERT-NUMERICAL-EQUAL (LIST 2.0d0) (MULTIPLE-VALUE-LIST - (LETM ((VEC (VECTOR-DOUBLE-FLOAT #(1.0d0 2.0d0 3.0d0)))) + (LETM ((VEC (VECTOR-DOUBLE-FLOAT (A 1.0d0 2.0d0 3.0d0)))) (POLYNOMIAL-EVAL VEC -1.0d0)))) (LISP-UNIT::ASSERT-NUMERICAL-EQUAL (LIST (LIST) (LIST)) @@ -242,4 +242,4 @@ #C(0.9999999999999999d0 0.0d0)) (MULTIPLE-VALUE-LIST (POLYNOMIAL-SOLVE - #(-1.0d0 0.0d0 0.0d0 0.0d0 0.0d0 1.0d0))))) + (A DOUBLE-FLOAT -1.0d0 0.0d0 0.0d0 0.0d0 0.0d0 1.0d0))))) diff --git a/random/dirichlet.lisp b/random/dirichlet.lisp index dc314af6..529e1e72 100644 --- a/random/dirichlet.lisp +++ b/random/dirichlet.lisp @@ -1,6 +1,6 @@ ;; Dirichlet distribution ;; Liam Healy, Sun Oct 29 2006 -;; Time-stamp: <2008-08-16 18:45:15EDT dirichlet.lisp> +;; Time-stamp: <2008-09-14 22:08:40EDT dirichlet.lisp> ;; $Id$ (in-package :gsl) @@ -58,15 +58,15 @@ #| (make-tests dirichlet (letm ((rng (random-number-generator *mt19937* 0)) - (alpha (vector-double-float #(1.0d0 2.0d0 3.0d0 4.0d0))) + (alpha (vector-double-float (a 1.0d0 2.0d0 3.0d0 4.0d0))) (theta (vector-double-float 4))) (dirichlet rng alpha theta) - (data theta)) - (letm ((alpha (vector-double-float #(1.0d0 2.0d0 3.0d0 4.0d0))) - (theta (vector-double-float #(0.1d0 0.3d0 0.4d0 0.2d0)))) + (cl-array theta)) + (letm ((alpha (vector-double-float (a 1.0d0 2.0d0 3.0d0 4.0d0))) + (theta (vector-double-float (a 0.1d0 0.3d0 0.4d0 0.2d0)))) (dirichlet-pdf alpha theta)) - (letm ((alpha (vector-double-float #(1.0d0 2.0d0 3.0d0 4.0d0))) - (theta (vector-double-float #(0.1d0 0.3d0 0.4d0 0.2d0)))) + (letm ((alpha (vector-double-float (a 1.0d0 2.0d0 3.0d0 4.0d0))) + (theta (vector-double-float (a 0.1d0 0.3d0 0.4d0 0.2d0)))) (dirichlet-log-pdf alpha theta))) |# @@ -78,21 +78,21 @@ (MULTIPLE-VALUE-LIST (LETM ((RNG (RANDOM-NUMBER-GENERATOR *MT19937* 0)) - (ALPHA (VECTOR-DOUBLE-FLOAT #(1.0d0 2.0d0 3.0d0 4.0d0))) + (ALPHA (VECTOR-DOUBLE-FLOAT (A 1.0d0 2.0d0 3.0d0 4.0d0))) (THETA (VECTOR-DOUBLE-FLOAT 4))) - (DIRICHLET RNG ALPHA THETA) (DATA THETA)))) + (DIRICHLET RNG ALPHA THETA) (CL-ARRAY THETA)))) (LISP-UNIT::ASSERT-NUMERICAL-EQUAL (LIST 2.8800000000000043d0) (MULTIPLE-VALUE-LIST (LETM - ((ALPHA (VECTOR-DOUBLE-FLOAT #(1.0d0 2.0d0 3.0d0 4.0d0))) - (THETA (VECTOR-DOUBLE-FLOAT #(0.1d0 0.3d0 0.4d0 0.2d0)))) + ((ALPHA (VECTOR-DOUBLE-FLOAT (A 1.0d0 2.0d0 3.0d0 4.0d0))) + (THETA (VECTOR-DOUBLE-FLOAT (A 0.1d0 0.3d0 0.4d0 0.2d0)))) (DIRICHLET-PDF ALPHA THETA)))) (LISP-UNIT::ASSERT-NUMERICAL-EQUAL (LIST 1.057790294147856d0) (MULTIPLE-VALUE-LIST (LETM - ((ALPHA (VECTOR-DOUBLE-FLOAT #(1.0d0 2.0d0 3.0d0 4.0d0))) - (THETA (VECTOR-DOUBLE-FLOAT #(0.1d0 0.3d0 0.4d0 0.2d0)))) + ((ALPHA (VECTOR-DOUBLE-FLOAT (A 1.0d0 2.0d0 3.0d0 4.0d0))) + (THETA (VECTOR-DOUBLE-FLOAT (A 0.1d0 0.3d0 0.4d0 0.2d0)))) (DIRICHLET-LOG-PDF ALPHA THETA))))) diff --git a/random/discrete.lisp b/random/discrete.lisp index 07c09dfa..d1a78b51 100644 --- a/random/discrete.lisp +++ b/random/discrete.lisp @@ -1,6 +1,6 @@ ;; Discrete random variables ;; Liam Healy, Sat Nov 11 2006 - 21:51 -;; Time-stamp: <2008-08-16 18:48:55EDT discrete.lisp> +;; Time-stamp: <2008-09-14 21:44:32EDT discrete.lisp> ;; $Id$ (in-package :gsl) @@ -54,13 +54,13 @@ (make-tests discrete ;; Must have two letms because the vector value is not set until ;; the body, but the discrete-random needs that set value. - (letm ((probabilities (vector-double-float #(0.25d0 0.5d0 0.25d0)))) + (letm ((probabilities (vector-double-float (a 0.25d0 0.5d0 0.25d0)))) (letm ((table (discrete-random probabilities)) (rng (random-number-generator *mt19937* 0))) (loop for i from 0 to 10 collect (discrete rng table)))) - (letm ((probabilities (vector-double-float #(0.25d0 0.5d0 0.25d0)))) + (letm ((probabilities (vector-double-float (a 0.25d0 0.5d0 0.25d0)))) (letm ((table (discrete-random probabilities))) (discrete-pdf 1 table)))) |# @@ -69,16 +69,16 @@ (LISP-UNIT::ASSERT-NUMERICAL-EQUAL (LIST (LIST 1 0 1 1 0 1 1 2 1 2 2)) (MULTIPLE-VALUE-LIST - (LETM ((PROBABILITIES - (VECTOR-DOUBLE-FLOAT #(0.25d0 0.5d0 0.25d0)))) + (LETM ((PROBABILITIES (VECTOR-DOUBLE-FLOAT (A 0.25d0 0.5d0 0.25d0)))) (LETM ((TABLE (DISCRETE-RANDOM PROBABILITIES)) - (RNG (RANDOM-NUMBER-GENERATOR *MT19937* 0))) + (RNG (RANDOM-NUMBER-GENERATOR *MT19937* 0))) (LOOP FOR I FROM 0 TO 10 COLLECT - (DISCRETE RNG TABLE)))))) + (DISCRETE RNG TABLE)))))) (LISP-UNIT::ASSERT-NUMERICAL-EQUAL (LIST 0.5d0) (MULTIPLE-VALUE-LIST - (LETM ((PROBABILITIES (VECTOR-DOUBLE-FLOAT #(0.25d0 0.5d0 0.25d0)))) + (LETM ((PROBABILITIES (VECTOR-DOUBLE-FLOAT (A 0.25d0 0.5d0 0.25d0)))) (LETM ((TABLE (DISCRETE-RANDOM PROBABILITIES))) (DISCRETE-PDF 1 TABLE)))))) + diff --git a/random/multinomial.lisp b/random/multinomial.lisp index 2b342ef2..31c614b9 100644 --- a/random/multinomial.lisp +++ b/random/multinomial.lisp @@ -1,6 +1,6 @@ ;; Multinomial distribution ;; Liam Healy, Sat Nov 25 2006 - 16:00 -;; Time-stamp: <2008-08-16 18:52:11EDT multinomial.lisp> +;; Time-stamp: <2008-09-14 21:40:55EDT multinomial.lisp> ;; $Id$ (in-package :gsl) @@ -51,19 +51,16 @@ #| (make-tests multinomial (letm ((rng (random-number-generator *mt19937* 0)) - (p (vector-double-float #(0.1d0 0.2d0 0.3d0 0.4d0))) + (p (vector-double-float (a 0.1d0 0.2d0 0.3d0 0.4d0))) (n (vector-fixnum 4))) (multinomial rng 8 p n) - (data n)) - (letm ((p (vector-double-float #(0.1d0 0.2d0 0.3d0 0.4d0))) - (n (vector-fixnum 4))) - (setf (data n) #(5 0 1 2)) + (cl-array n)) + (letm ((p (vector-double-float (a 0.1d0 0.2d0 0.3d0 0.4d0))) + (n (vector-fixnum 4) (a 5 0 1 2))) (multinomial-pdf p N)) - (letm ((p (vector-double-float #(0.1d0 0.2d0 0.3d0 0.4d0))) - (n (vector-fixnum 4))) - (setf (data n) #(5 0 1 2)) + (letm ((p (vector-double-float (a 0.1d0 0.2d0 0.3d0 0.4d0))) + (n (vector-fixnum 4) (a 5 0 1 2))) (multinomial-log-pdf p n))) -|# (LISP-UNIT:DEFINE-TEST MULTINOMIAL (LISP-UNIT::ASSERT-NUMERICAL-EQUAL @@ -90,3 +87,4 @@ (SETF (DATA N) #(5 0 1 2)) (MULTINOMIAL-LOG-PDF P N))))) +|# diff --git a/random/quasi.lisp b/random/quasi.lisp index 881d3c87..74354f6b 100644 --- a/random/quasi.lisp +++ b/random/quasi.lisp @@ -1,6 +1,6 @@ ;; Quasi-random sequences in arbitrary dimensions. ;; Liam Healy, Sun Jul 16 2006 - 15:54 -;; Time-stamp: <2008-08-16 20:17:54EDT quasi.lisp> +;; Time-stamp: <2008-09-14 22:01:42EDT quasi.lisp> ;; $Id$ (in-package :gsl) @@ -132,7 +132,7 @@ (vec (vector-double-float 2))) (loop repeat 5 do (qrng-get gen vec) - append (coerce (data vec) 'list)))) + append (coerce (cl-array vec) 'list)))) |# (LISP-UNIT:DEFINE-TEST QUASI-RANDOM-NUMBER-GENERATORS @@ -145,5 +145,5 @@ ((GEN (QUASI-RANDOM-NUMBER-GENERATOR 2 *SOBOL*)) (VEC (VECTOR-DOUBLE-FLOAT 2))) (LOOP REPEAT 5 DO (QRNG-GET GEN VEC) APPEND - (COERCE (DATA VEC) 'LIST)))))) + (COERCE (CL-ARRAY VEC) 'LIST)))))) diff --git a/random/shuffling-sampling.lisp b/random/shuffling-sampling.lisp index c5a01a44..4595a248 100644 --- a/random/shuffling-sampling.lisp +++ b/random/shuffling-sampling.lisp @@ -1,6 +1,6 @@ ;; Shuffling and sampling ;; Liam Healy, Sat Dec 2 2006 - 18:40 -;; Time-stamp: <2008-08-16 18:55:38EDT shuffling-sampling.lisp> +;; Time-stamp: <2008-09-14 21:41:43EDT shuffling-sampling.lisp> ;; $Id$ (in-package :gsl) @@ -72,7 +72,6 @@ (v2 (vector-fixnum 10))) (sample rng v2 v1) (data v2))) -|# (LISP-UNIT:DEFINE-TEST SHUFFLING-SAMPLING (LISP-UNIT::ASSERT-NUMERICAL-EQUAL @@ -99,3 +98,4 @@ (V1 (VECTOR-FIXNUM #(1 2 3 4 5 6 7 8))) (V2 (VECTOR-FIXNUM 10))) (SAMPLE RNG V2 V1) (DATA V2))))) +|# diff --git a/special-functions/bessel.lisp b/special-functions/bessel.lisp index 7e2f2c6d..1077598a 100644 --- a/special-functions/bessel.lisp +++ b/special-functions/bessel.lisp @@ -1,6 +1,6 @@ ;; Bessel functions ;; Liam Healy, Fri Mar 17 2006 - 18:42 -;; Time-stamp: <2008-08-10 17:55:45EDT bessel.lisp> +;; Time-stamp: <2008-09-14 22:08:10EDT bessel.lisp> ;; $Id$ (in-package :gsl) @@ -465,61 +465,61 @@ (cylindrical-bessel-Jn 2 4.0d0) (letm ((besarr (vector-double-float 4))) (cylindrical-bessel-Jn-array 2.0d0 besarr 2) - (data besarr)) + (cl-array besarr)) (cylindrical-bessel-Y0 4.0d0) (cylindrical-bessel-Y1 4.0d0) (cylindrical-bessel-Yn 3 4.0d0) (letm ((besarr (vector-double-float 4))) (cylindrical-bessel-Yn-array 2.0d0 besarr 2) - (data besarr)) + (cl-array besarr)) (cylindrical-bessel-I0 4.0d0) (cylindrical-bessel-I1 4.0d0) (cylindrical-bessel-In 3 4.0d0) (letm ((besarr (vector-double-float 4))) (cylindrical-bessel-In-array 2.0d0 besarr 2) - (data besarr)) + (cl-array besarr)) (cylindrical-bessel-I0-scaled 4.0d0) (cylindrical-bessel-I1-scaled 4.0d0) (cylindrical-bessel-In-scaled 3 4.0d0) (letm ((besarr (vector-double-float 4))) (cylindrical-bessel-In-scaled-array 2.0d0 besarr 2) - (data besarr)) + (cl-array besarr)) (cylindrical-bessel-K0 4.0d0) (cylindrical-bessel-K1 4.0d0) (cylindrical-bessel-Kn 2 4.0d0) (letm ((besarr (vector-double-float 4))) (cylindrical-bessel-Kn-array 2.0d0 besarr 2) - (data besarr)) + (cl-array besarr)) (cylindrical-bessel-K0-scaled 4.0d0) (cylindrical-bessel-K1-scaled 4.0d0) (cylindrical-bessel-Kn-scaled 2 4.0d0) (letm ((besarr (vector-double-float 4))) (cylindrical-bessel-Kn-array 2.0d0 besarr 2) - (data besarr)) + (cl-array besarr)) (spherical-bessel-j0 4.0d0) (spherical-bessel-j1 4.0d0) (spherical-bessel-j2 4.0d0) (spherical-bessel-jl 3 4.0d0) (letm ((besarr (vector-double-float 4))) (spherical-bessel-jl-array 4.0d0 besarr) - (data besarr)) + (cl-array besarr)) (letm ((besarr (vector-double-float 4))) (spherical-bessel-jl-steed-array 4.0d0 besarr) - (data besarr)) + (cl-array besarr)) (spherical-bessel-y0 4.0d0) (spherical-bessel-y1 4.0d0) (spherical-bessel-y2 4.0d0) (spherical-bessel-yl 2 4.0d0) (letm ((besarr (vector-double-float 4))) (spherical-bessel-yl-array 4.0d0 besarr) - (data besarr)) + (cl-array besarr)) (spherical-bessel-i0-scaled 4.0d0) (spherical-bessel-i1-scaled 4.0d0) (spherical-bessel-i2-scaled 4.0d0) (spherical-bessel-il-scaled 3 4.0d0) (letm ((besarr (vector-double-float 4))) (spherical-bessel-il-scaled-array 4.0d0 besarr) - (data besarr)) + (cl-array besarr)) (spherical-bessel-k0-scaled 4.0d0) (spherical-bessel-k1-scaled 4.0d0) @@ -527,11 +527,11 @@ (spherical-bessel-kl-scaled 3 4.0d0) (letm ((besarr (vector-double-float 4))) (spherical-bessel-kl-scaled-array 4.0d0 besarr) - (data besarr)) + (cl-array besarr)) (bessel-jnu 3.0d0 4.0d0) - (letm ((besarr (vector-double-float #(1.0d0 2.0d0 3.0d0)))) + (letm ((besarr (vector-double-float (a 1.0d0 2.0d0 3.0d0)))) (spherical-Jnu-array 0.5d0 besarr) - (data besarr)) + (cl-array besarr)) (bessel-Ynu 3.0d0 4.0d0) (bessel-Inu 3.0d0 4.0d0) (bessel-Inu-scaled 3.0d0 4.0d0) @@ -562,7 +562,7 @@ (MULTIPLE-VALUE-LIST (LETM ((BESARR (VECTOR-DOUBLE-FLOAT 4))) (CYLINDRICAL-BESSEL-JN-ARRAY 2.0d0 BESARR 2) - (DATA BESARR)))) + (CL-ARRAY BESARR)))) (LISP-UNIT::ASSERT-NUMERICAL-EQUAL (LIST -0.016940739325064968d0 1.8993556609468549d-16) (MULTIPLE-VALUE-LIST (CYLINDRICAL-BESSEL-Y0 4.0d0))) @@ -579,7 +579,7 @@ (MULTIPLE-VALUE-LIST (LETM ((BESARR (VECTOR-DOUBLE-FLOAT 4))) (CYLINDRICAL-BESSEL-YN-ARRAY 2.0d0 BESARR 2) - (DATA BESARR)))) + (CL-ARRAY BESARR)))) (LISP-UNIT::ASSERT-NUMERICAL-EQUAL (LIST 11.301921952136329d0 2.7297681442535893d-14) (MULTIPLE-VALUE-LIST (CYLINDRICAL-BESSEL-I0 4.0d0))) @@ -596,7 +596,7 @@ (MULTIPLE-VALUE-LIST (LETM ((BESARR (VECTOR-DOUBLE-FLOAT 4))) (CYLINDRICAL-BESSEL-IN-ARRAY 2.0d0 BESARR 2) - (DATA BESARR)))) + (CL-ARRAY BESARR)))) (LISP-UNIT::ASSERT-NUMERICAL-EQUAL (LIST 0.2070019212239867d0 2.241925168997723d-16) (MULTIPLE-VALUE-LIST @@ -617,7 +617,7 @@ (LETM ((BESARR (VECTOR-DOUBLE-FLOAT 4))) (CYLINDRICAL-BESSEL-IN-SCALED-ARRAY 2.0d0 BESARR 2) - (DATA BESARR)))) + (CL-ARRAY BESARR)))) (LISP-UNIT::ASSERT-NUMERICAL-EQUAL (LIST 0.011159676085853023d0 2.0424662435034432d-17) (MULTIPLE-VALUE-LIST (CYLINDRICAL-BESSEL-K0 4.0d0))) @@ -634,7 +634,7 @@ (MULTIPLE-VALUE-LIST (LETM ((BESARR (VECTOR-DOUBLE-FLOAT 4))) (CYLINDRICAL-BESSEL-KN-ARRAY 2.0d0 BESARR 2) - (DATA BESARR)))) + (CL-ARRAY BESARR)))) (LISP-UNIT::ASSERT-NUMERICAL-EQUAL (LIST 0.6092976692566953d0 3.0340122249326356d-16) (MULTIPLE-VALUE-LIST @@ -654,7 +654,7 @@ (MULTIPLE-VALUE-LIST (LETM ((BESARR (VECTOR-DOUBLE-FLOAT 4))) (CYLINDRICAL-BESSEL-KN-ARRAY 2.0d0 BESARR 2) - (DATA BESARR)))) + (CL-ARRAY BESARR)))) (LISP-UNIT::ASSERT-NUMERICAL-EQUAL (LIST -0.18920062382698205d0 1.6804391107692678d-16) (MULTIPLE-VALUE-LIST (SPHERICAL-BESSEL-J0 4.0d0))) @@ -674,7 +674,7 @@ (MULTIPLE-VALUE-LIST (LETM ((BESARR (VECTOR-DOUBLE-FLOAT 4))) (SPHERICAL-BESSEL-JL-ARRAY 4.0d0 BESARR) - (DATA BESARR)))) + (CL-ARRAY BESARR)))) (LISP-UNIT::ASSERT-NUMERICAL-EQUAL (LIST #(-0.18920062382698208d0 0.11611074925915742d0 @@ -682,7 +682,7 @@ (MULTIPLE-VALUE-LIST (LETM ((BESARR (VECTOR-DOUBLE-FLOAT 4))) (SPHERICAL-BESSEL-JL-STEED-ARRAY 4.0d0 BESARR) - (DATA BESARR)))) + (CL-ARRAY BESARR)))) (LISP-UNIT::ASSERT-NUMERICAL-EQUAL (LIST 0.16341090521590299d0 1.4513803955642766d-16) (MULTIPLE-VALUE-LIST (SPHERICAL-BESSEL-Y0 4.0d0))) @@ -702,7 +702,7 @@ (MULTIPLE-VALUE-LIST (LETM ((BESARR (VECTOR-DOUBLE-FLOAT 4))) (SPHERICAL-BESSEL-YL-ARRAY 4.0d0 BESARR) - (DATA BESARR)))) + (CL-ARRAY BESARR)))) (LISP-UNIT::ASSERT-NUMERICAL-EQUAL (LIST 0.12495806717151219d0 5.5492529314587895d-17) (MULTIPLE-VALUE-LIST @@ -726,7 +726,7 @@ (MULTIPLE-VALUE-LIST (LETM ((BESARR (VECTOR-DOUBLE-FLOAT 4))) (SPHERICAL-BESSEL-IL-SCALED-ARRAY 4.0d0 BESARR) - (DATA BESARR)))) + (CL-ARRAY BESARR)))) (LISP-UNIT::ASSERT-NUMERICAL-EQUAL (LIST 0.39269908169872414d0 1.743934249004316d-16) (MULTIPLE-VALUE-LIST @@ -750,7 +750,7 @@ (MULTIPLE-VALUE-LIST (LETM ((BESARR (VECTOR-DOUBLE-FLOAT 4))) (SPHERICAL-BESSEL-KL-SCALED-ARRAY 4.0d0 BESARR) - (DATA BESARR)))) + (CL-ARRAY BESARR)))) (LISP-UNIT::ASSERT-NUMERICAL-EQUAL (LIST 0.43017147387562193d0 7.641380397338472d-16) (MULTIPLE-VALUE-LIST (BESSEL-JNU 3.0d0 4.0d0))) @@ -759,9 +759,9 @@ #(0.6713967071418024d0 0.5130161365618323d0 0.06500818287738516d0)) (MULTIPLE-VALUE-LIST - (LETM ((BESARR (VECTOR-DOUBLE-FLOAT #(1.0d0 2.0d0 3.0d0)))) + (LETM ((BESARR (VECTOR-DOUBLE-FLOAT (a 1.0d0 2.0d0 3.0d0)))) (SPHERICAL-JNU-ARRAY 0.5d0 BESARR) - (DATA BESARR)))) + (CL-ARRAY BESARR)))) (LISP-UNIT::ASSERT-NUMERICAL-EQUAL (LIST -0.1820221159534852d0 2.020851441225493d-15) (MULTIPLE-VALUE-LIST (BESSEL-YNU 3.0d0 4.0d0))) diff --git a/special-functions/coulomb.lisp b/special-functions/coulomb.lisp index 1ac1b570..4da55d2d 100644 --- a/special-functions/coulomb.lisp +++ b/special-functions/coulomb.lisp @@ -1,6 +1,6 @@ ;; Coulumb functions ;; Liam Healy, Sat Mar 18 2006 - 23:23 -;; Time-stamp: <2008-07-08 21:46:07EDT coulomb.lisp> +;; Time-stamp: <2008-09-14 22:29:18EDT coulomb.lisp> ;; $Id$ (in-package :gsl) @@ -133,17 +133,17 @@ (hydrogenicr 3 1 1.0d0 2.5d0) (coulomb-wave-FG 0.0d0 1.0d0 2.0d0 0) (letm ((arr (vector-double-float 3))) - (coulomb-wave-F-array 0.0d0 1.0d0 2.0d0 arr) - (data arr)) + (coulomb-wave-F-array 0.0d0 1.0d0 2.0d0 arr) + (cl-array arr)) (coulomb-wave-fg 1.0d0 2.0d0 2.5d0 1) (letm ((Farr (vector-double-float 3)) (Garr (vector-double-float 3))) (coulomb-wave-FG-array 1.5d0 1.0d0 1.0d0 Farr Garr) - (append (coerce (data Farr) 'list) (coerce (data Garr) 'list))) + (append (coerce (cl-array Farr) 'list) (coerce (cl-array Garr) 'list))) (letm ((arr (vector-double-float 3))) - (coulomb-wave-sphF-array 0.0d0 1.0d0 2.0d0 arr) (data arr)) + (coulomb-wave-sphF-array 0.0d0 1.0d0 2.0d0 arr) (cl-array arr)) (coulomb-cl 1.0d0 2.5d0) (letm ((cl (vector-double-float 3))) - (coulomb-CL-array 0.0d0 1.0d0 cl) (data cl))) + (coulomb-CL-array 0.0d0 1.0d0 cl) (cl-array cl))) |# #| @@ -169,7 +169,7 @@ (MULTIPLE-VALUE-LIST (LETM ((ARR (VECTOR-DOUBLE-FLOAT 3))) (COULOMB-WAVE-F-ARRAY 0.0d0 1.0d0 2.0d0 ARR) - (DATA ARR)))) + (CL-ARRAY ARR)))) (LISP-UNIT::ASSERT-NUMERICAL-EQUAL (LIST 0.07161779967468254d0 0.1278101031568499d0 2.2510703464871114d0 -1.4245543587641651d0 0.0d0 @@ -187,8 +187,8 @@ (LETM ((FARR (VECTOR-DOUBLE-FLOAT 3)) (GARR (VECTOR-DOUBLE-FLOAT 3))) (COULOMB-WAVE-FG-ARRAY 1.5d0 1.0d0 1.0d0 FARR GARR) - (APPEND (COERCE (DATA FARR) 'LIST) - (COERCE (DATA GARR) 'LIST))))) + (APPEND (COERCE (CL-ARRAY FARR) 'LIST) + (COERCE (CL-ARRAY GARR) 'LIST))))) (LISP-UNIT::ASSERT-NUMERICAL-EQUAL (LIST #(0.33089080691634065d0 0.18070642887252675d0 @@ -196,7 +196,7 @@ (MULTIPLE-VALUE-LIST (LETM ((ARR (VECTOR-DOUBLE-FLOAT 3))) (COULOMB-WAVE-SPHF-ARRAY 0.0d0 1.0d0 2.0d0 ARR) - (DATA ARR)))) + (CL-ARRAY ARR)))) (LISP-UNIT::ASSERT-NUMERICAL-EQUAL (LIST 0.0013809146441856027d0 2.759621819430441d-17) (MULTIPLE-VALUE-LIST (COULOMB-CL 1.0d0 2.5d0))) @@ -206,6 +206,6 @@ 0.011428736368066591d0)) (MULTIPLE-VALUE-LIST (LETM ((CL (VECTOR-DOUBLE-FLOAT 3))) - (COULOMB-CL-ARRAY 0.0d0 1.0d0 CL) (DATA CL))))) + (COULOMB-CL-ARRAY 0.0d0 1.0d0 CL) (CL-ARRAY CL))))) |# diff --git a/special-functions/gegenbauer.lisp b/special-functions/gegenbauer.lisp index eccf7299..dcb1a454 100644 --- a/special-functions/gegenbauer.lisp +++ b/special-functions/gegenbauer.lisp @@ -1,6 +1,6 @@ ;; Gegenbauer polynomials ;; Liam Healy, Fri Apr 28 2006 - 20:40 -;; Time-stamp: <2008-07-08 21:53:11EDT gegenbauer.lisp> +;; Time-stamp: <2008-09-14 21:38:13EDT gegenbauer.lisp> ;; $Id$ (in-package :gsl) @@ -50,7 +50,7 @@ (gegenbauer-3 1.0d0 3.0d0) (gegenbauer 4 1.0d0 3.0d0) (letm ((arr (vector-double-float 4))) - (gegenbauer-array 1.0d0 3.0d0 arr) (data arr))) + (gegenbauer-array 1.0d0 3.0d0 arr) (cl-array arr))) |# (LISP-UNIT:DEFINE-TEST GEGENBAUER @@ -70,4 +70,4 @@ (LIST #(1.0d0 6.0d0 35.0d0 204.0d0)) (MULTIPLE-VALUE-LIST (LETM ((ARR (VECTOR-DOUBLE-FLOAT 4))) - (GEGENBAUER-ARRAY 1.0d0 3.0d0 ARR) (DATA ARR))))) + (GEGENBAUER-ARRAY 1.0d0 3.0d0 ARR) (CL-ARRAY ARR))))) diff --git a/special-functions/legendre.lisp b/special-functions/legendre.lisp index b4f1c58d..ae3cd8d7 100644 --- a/special-functions/legendre.lisp +++ b/special-functions/legendre.lisp @@ -1,6 +1,6 @@ ;; Legendre functions ;; Liam Healy, Sat Apr 29 2006 - 19:16 -;; Time-stamp: <2008-07-08 21:55:40EDT legendre.lisp> +;; Time-stamp: <2008-09-14 21:38:13EDT legendre.lisp> ;; $Id$ (in-package :gsl) @@ -261,26 +261,26 @@ (legendre-Pl 4 0.3d0) (letm ((arr (vector-double-float 4))) (legendre-Pl-array 0.5d0 arr) - (data arr)) + (cl-array arr)) (legendre-Q0 3.3d0) (legendre-Q1 3.3d0) (legendre-Ql 2 3.3d0) (legendre-Plm 4 3 0.5d0) (letm ((arr (vector-double-float 4))) (legendre-Plm-array 2 0.5d0 arr) - (data arr)) + (cl-array arr)) (letm ((val (vector-double-float 4)) (deriv (vector-double-float 4))) (legendre-Plm-deriv-array 2 0.5d0 val deriv) - (data deriv)) + (cl-array deriv)) (legendre-sphplm 1200 1100 0.3d0) (letm ((arr (vector-double-float 4))) (legendre-sphPlm-array 4 0.5d0 arr) - (data arr)) + (cl-array arr)) (letm ((val (vector-double-float 4)) (deriv (vector-double-float 4))) (legendre-sphPlm-deriv-array 4 0.5d0 val deriv) - (data deriv)) + (cl-array deriv)) (legendre-conicalp-half 3.5d0 10.0d0) (legendre-conicalp-mhalf 3.5d0 10.0d0) (legendre-conicalp-0 3.5d0 10.0d0) @@ -292,7 +292,7 @@ (legendre-h3d 4 1.0d0 0.5d0) (letm ((arr (vector-double-float 4))) (legendre-h3d-array 1.0d0 0.5d0 arr) - (data arr))) + (cl-array arr))) |# (LISP-UNIT:DEFINE-TEST LEGENDRE @@ -317,7 +317,7 @@ (LIST #(1.0d0 0.5d0 -0.125d0 -0.4375d0)) (MULTIPLE-VALUE-LIST (LETM ((ARR (VECTOR-DOUBLE-FLOAT 4))) - (LEGENDRE-PL-ARRAY 0.5d0 ARR) (DATA ARR)))) + (LEGENDRE-PL-ARRAY 0.5d0 ARR) (CL-ARRAY ARR)))) (LISP-UNIT::ASSERT-NUMERICAL-EQUAL (LIST 0.3128529498822064d0 1.3893461931245028d-16) (MULTIPLE-VALUE-LIST (LEGENDRE-Q0 3.3d0))) @@ -334,14 +334,14 @@ (LIST #(2.25d0 5.625d0 4.21875d0 -4.921875d0)) (MULTIPLE-VALUE-LIST (LETM ((ARR (VECTOR-DOUBLE-FLOAT 4))) - (LEGENDRE-PLM-ARRAY 2 0.5d0 ARR) (DATA ARR)))) + (LEGENDRE-PLM-ARRAY 2 0.5d0 ARR) (CL-ARRAY ARR)))) (LISP-UNIT::ASSERT-NUMERICAL-EQUAL (LIST #(-3.0d0 3.75d0 33.75d0 55.78125d0)) (MULTIPLE-VALUE-LIST (LETM ((VAL (VECTOR-DOUBLE-FLOAT 4)) (DERIV (VECTOR-DOUBLE-FLOAT 4))) (LEGENDRE-PLM-DERIV-ARRAY 2 0.5d0 VAL DERIV) - (DATA DERIV)))) + (CL-ARRAY DERIV)))) (LISP-UNIT::ASSERT-NUMERICAL-EQUAL (LIST 0.30366280894310793d0 3.438761110552081d-15) (MULTIPLE-VALUE-LIST @@ -353,7 +353,7 @@ (MULTIPLE-VALUE-LIST (LETM ((ARR (VECTOR-DOUBLE-FLOAT 4))) (LEGENDRE-SPHPLM-ARRAY 4 0.5d0 ARR) - (DATA ARR)))) + (CL-ARRAY ARR)))) (LISP-UNIT::ASSERT-NUMERICAL-EQUAL (LIST #(-0.6637990386674741d0 -0.2751965434323283d0 @@ -362,7 +362,7 @@ (LETM ((VAL (VECTOR-DOUBLE-FLOAT 4)) (DERIV (VECTOR-DOUBLE-FLOAT 4))) (LEGENDRE-SPHPLM-DERIV-ARRAY 4 0.5d0 VAL DERIV) - (DATA DERIV)))) + (CL-ARRAY DERIV)))) (LISP-UNIT::ASSERT-NUMERICAL-EQUAL (LIST -0.1255299048878925d0 1.3395992025650077d-15) (MULTIPLE-VALUE-LIST @@ -404,4 +404,4 @@ (MULTIPLE-VALUE-LIST (LETM ((ARR (VECTOR-DOUBLE-FLOAT 4))) (LEGENDRE-H3D-ARRAY 1.0d0 0.5d0 ARR) - (DATA ARR))))) + (CL-ARRAY ARR))))) -- GitLab