Commit 003c83f8 by Liam M. Healy

Tests include explicit arguments for copy-to

Test results no longer assume default values of arguments to copy-to,
which can change. This improves the results with fewer execution
errors.
parent d21fa5b5
;; Regression test AXPY for GSLL, automatically generated
;;
;; Copyright 2009, 2011 Liam M. Healy
;; Copyright 2009, 2011, 2014 Liam M. Healy
;; Distributed under the terms of the GNU General Public License
;;
;; This program is free software: you can redistribute it and/or modify
......@@ -33,7 +33,7 @@
'SINGLE-FLOAT :INITIAL-CONTENTS
'(32.5f0 42.73f0 -17.24f0 43.31f0 -16.12f0 -8.25f0 21.44f0 -49.08f0)))
(SCALAR -39.66f0))
(grid:copy-to (AXPY SCALAR V1 V2)))))
(GRID:COPY-TO (AXPY SCALAR V1 V2) 'array 'single-float))))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST
#(1400.77d0 -284.06839999999994d0 -147.7214d0
......@@ -55,7 +55,7 @@
-16.12d0 -8.25d0 21.44d0
-49.08d0)))
(SCALAR -39.66d0))
(GRID:COPY-TO (AXPY SCALAR V1 V2)))))
(GRID:COPY-TO (AXPY SCALAR V1 V2) 'array 'double-float))))
#+fsbv
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST
......@@ -80,7 +80,7 @@
#C(19.68f0 -5.55f0) #C(-8.82f0 25.37f0)
#C(-30.58f0 31.67f0))))
(SCALAR #C(-39.66f0 -49.46f0)))
(GRID:COPY-TO (AXPY SCALAR V1 V2)))))
(GRID:COPY-TO (AXPY SCALAR V1 V2) 'array '(complex single-float)))))
#+fsbv
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST
......@@ -110,4 +110,4 @@
#C(-39.66d0 -49.46d0) #C(19.68d0 -5.55d0)
#C(-8.82d0 25.37d0) #C(-30.58d0 31.67d0))))
(SCALAR #C(-39.66d0 -49.46d0)))
(GRID:COPY-TO (AXPY SCALAR V1 V2))))))
(GRID:COPY-TO (AXPY SCALAR V1 V2) 'array '(complex double-float))))))
;; Regression test BESSEL for GSLL, automatically generated
;;
;; Copyright 2009 Liam M. Healy
;; Copyright 2009, 2014 Liam M. Healy
;; Distributed under the terms of the GNU General Public License
;;
;; This program is free software: you can redistribute it and/or modify
......@@ -31,7 +31,7 @@
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
#(0.35283402861563773d0 0.12894324947440206d0
0.033995719807568436d0 0.007039629755871686d0)
(grid:copy-to (cylindrical-bessel-J-array-order 2.0d0 4 2)))
(grid:copy-to (cylindrical-bessel-J-array-order 2.0d0 4 2) 'array 'double-float))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST -0.016940739325064968d0 1.8993556609468549d-16)
(MULTIPLE-VALUE-LIST (CYLINDRICAL-BESSEL-Y0 4.0d0)))
......@@ -49,7 +49,7 @@
(LET ((BESARR
(GRID:MAKE-FOREIGN-ARRAY 'DOUBLE-FLOAT :DIMENSIONS 4)))
(CYLINDRICAL-BESSEL-YN-ARRAY 2.0d0 BESARR 2)
(GRID:COPY-TO BESARR))))
(GRID:COPY-TO BESARR 'array 'double-float))))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST 11.301921952136329d0 2.7297681442535893d-14)
(MULTIPLE-VALUE-LIST (CYLINDRICAL-BESSEL-I0 4.0d0)))
......@@ -70,7 +70,7 @@
(LET ((BESARR
(GRID:MAKE-FOREIGN-ARRAY 'DOUBLE-FLOAT :DIMENSIONS 4)))
(CYLINDRICAL-BESSEL-IN-ARRAY 2.0d0 BESARR 2)
(GRID:COPY-TO BESARR))))
(GRID:COPY-TO BESARR 'array 'double-float))))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST 0.2070019212239867d0 2.241925168997723d-16)
(MULTIPLE-VALUE-LIST
......@@ -94,7 +94,7 @@
(LET ((BESARR
(GRID:MAKE-FOREIGN-ARRAY 'DOUBLE-FLOAT :DIMENSIONS 4)))
(CYLINDRICAL-BESSEL-IN-SCALED-ARRAY 2.0d0 BESARR 2)
(GRID:COPY-TO BESARR))))
(GRID:COPY-TO BESARR 'array 'double-float))))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST 0.011159676085853023d0 2.0424662435034432d-17)
(MULTIPLE-VALUE-LIST (CYLINDRICAL-BESSEL-K0 4.0d0)))
......@@ -130,7 +130,7 @@
(LET ((BESARR
(GRID:MAKE-FOREIGN-ARRAY 'DOUBLE-FLOAT :DIMENSIONS 4)))
(CYLINDRICAL-BESSEL-KN-ARRAY 2.0d0 BESARR 2)
(GRID:COPY-TO BESARR))))
(GRID:COPY-TO BESARR 'array 'double-float))))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST -0.18920062382698205d0 1.6804391107692678d-16)
(MULTIPLE-VALUE-LIST (SPHERICAL-BESSEL-J0 4.0d0)))
......@@ -151,7 +151,7 @@
(LET ((BESARR
(GRID:MAKE-FOREIGN-ARRAY 'DOUBLE-FLOAT :DIMENSIONS 4)))
(SPHERICAL-BESSEL-JL-ARRAY 4.0d0 BESARR)
(GRID:COPY-TO BESARR))))
(GRID:COPY-TO BESARR 'array 'double-float))))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST
#(-0.18920062382698208d0 0.11611074925915742d0
......@@ -160,7 +160,7 @@
(LET ((BESARR
(GRID:MAKE-FOREIGN-ARRAY 'DOUBLE-FLOAT :DIMENSIONS 4)))
(SPHERICAL-BESSEL-JL-STEED-ARRAY 4.0d0 BESARR)
(GRID:COPY-TO BESARR))))
(GRID:COPY-TO BESARR 'array 'double-float))))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST 0.16341090521590299d0 1.4513803955642766d-16)
(MULTIPLE-VALUE-LIST (SPHERICAL-BESSEL-Y0 4.0d0)))
......@@ -181,7 +181,7 @@
(LET ((BESARR
(GRID:MAKE-FOREIGN-ARRAY 'DOUBLE-FLOAT :DIMENSIONS 4)))
(SPHERICAL-BESSEL-YL-ARRAY 4.0d0 BESARR)
(GRID:COPY-TO BESARR))))
(GRID:COPY-TO BESARR 'array 'double-float))))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST 0.12495806717151219d0 5.5492529314587895d-17)
(MULTIPLE-VALUE-LIST
......@@ -206,7 +206,7 @@
(LET ((BESARR
(GRID:MAKE-FOREIGN-ARRAY 'DOUBLE-FLOAT :DIMENSIONS 4)))
(SPHERICAL-BESSEL-IL-SCALED-ARRAY 4.0d0 BESARR)
(GRID:COPY-TO BESARR))))
(GRID:COPY-TO BESARR 'array 'double-float))))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST 0.39269908169872414d0 1.743934249004316d-16)
(MULTIPLE-VALUE-LIST
......@@ -231,7 +231,7 @@
(LET ((BESARR
(GRID:MAKE-FOREIGN-ARRAY 'DOUBLE-FLOAT :DIMENSIONS 4)))
(SPHERICAL-BESSEL-KL-SCALED-ARRAY 4.0d0 BESARR)
(GRID:COPY-TO BESARR))))
(GRID:COPY-TO BESARR 'array 'double-float))))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST 0.43017147387562193d0 7.641380397338472d-16)
(MULTIPLE-VALUE-LIST (cylindrical-bessel-J 3.0d0 4.0d0)))
......@@ -243,7 +243,8 @@
(grid:copy-to
(cylindrical-bessel-J-array-x
0.5d0
(grid:make-foreign-array 'double-float :initial-contents '(1.0d0 2.0d0 3.0d0))))))
(grid:make-foreign-array 'double-float :initial-contents '(1.0d0 2.0d0 3.0d0)))
'array 'double-float)))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST -0.1820221159534852d0 2.020851441225493d-15)
(MULTIPLE-VALUE-LIST (cylindrical-bessel-Y 3.0d0 4.0d0)))
......
;; Regression test BLAS-COPY for GSLL, automatically generated
;;
;; Copyright 2009, 2010, 2011 Liam M. Healy
;; Copyright 2009, 2010, 2011, 2014 Liam M. Healy
;; Distributed under the terms of the GNU General Public License
;;
;; This program is free software: you can redistribute it and/or modify
......@@ -27,7 +27,7 @@
'SINGLE-FLOAT :INITIAL-CONTENTS '(-34.5f0 8.24f0 3.29f0)))
(V2 (GRID:MAKE-FOREIGN-ARRAY 'SINGLE-FLOAT :DIMENSIONS '3)))
(BLAS-COPY V1 V2)
(GRID:COPY-TO V2))))
(GRID:COPY-TO V2 'array 'single-float))))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST #(-34.5d0 8.24d0 3.29d0))
(MULTIPLE-VALUE-LIST
......@@ -36,7 +36,7 @@
'DOUBLE-FLOAT :INITIAL-CONTENTS '(-34.5d0 8.24d0 3.29d0)))
(V2 (GRID:MAKE-FOREIGN-ARRAY 'DOUBLE-FLOAT :DIMENSIONS '3)))
(BLAS-COPY V1 V2)
(GRID:COPY-TO V2))))
(GRID:COPY-TO V2 'array 'double-float))))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST #(#C(-34.5f0 8.24f0) #C(3.29f0 -8.93f0) #C(34.12f0 -6.15f0)))
(MULTIPLE-VALUE-LIST
......@@ -49,7 +49,7 @@
(GRID:MAKE-FOREIGN-ARRAY
'(COMPLEX SINGLE-FLOAT) :DIMENSIONS '3)))
(BLAS-COPY V1 V2)
(GRID:COPY-TO V2))))
(GRID:COPY-TO V2 'array '(complex single-float)))))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST
#(#C(-34.5d0 8.24d0) #C(3.29d0 -8.93d0)
......@@ -65,4 +65,4 @@
(GRID:MAKE-FOREIGN-ARRAY
'(COMPLEX DOUBLE-FLOAT) :DIMENSIONS '3)))
(BLAS-COPY V1 V2)
(GRID:COPY-TO V2)))))
(GRID:COPY-TO V2 'array '(complex double-float))))))
;; Regression test BLAS-SWAP for GSLL, automatically generated
;;
;; Copyright 2009, 2010, 2011 Liam M. Healy
;; Copyright 2009, 2010, 2011, 2014 Liam M. Healy
;; Distributed under the terms of the GNU General Public License
;;
;; This program is free software: you can redistribute it and/or modify
......@@ -31,7 +31,7 @@
'SINGLE-FLOAT :INITIAL-CONTENTS
'(-8.93f0 34.12f0 -6.15f0))))
(BLAS-SWAP V2 V1)
(LIST (GRID:COPY-TO V1) (GRID:COPY-TO V2)))))
(LIST (GRID:COPY-TO V1 'array 'single-float) (GRID:COPY-TO V2 'array 'single-float)))))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST
(LIST #(-8.93d0 34.12d0 -6.15d0)
......@@ -46,7 +46,7 @@
'DOUBLE-FLOAT :INITIAL-CONTENTS
'(-8.93d0 34.12d0 -6.15d0))))
(BLAS-SWAP V2 V1)
(LIST (GRID:COPY-TO V1) (GRID:COPY-TO V2)))))
(LIST (GRID:COPY-TO V1 'array 'double-float) (GRID:COPY-TO V2 'array 'double-float)))))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST
(LIST
......@@ -64,7 +64,7 @@
:INITIAL-CONTENTS
'(#C(-8.93f0 34.12f0) #C(-6.15f0 49.27f0) #C(-13.49f0 32.5f0)))))
(BLAS-SWAP V2 V1)
(LIST (GRID:COPY-TO V1) (GRID:COPY-TO V2)))))
(LIST (GRID:COPY-TO V1 'array '(complex single-float)) (GRID:COPY-TO V2 'array '(complex single-float))))))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST
(LIST
......@@ -84,4 +84,4 @@
:INITIAL-CONTENTS
'(#C(-8.93d0 34.12d0) #C(-6.15d0 49.27d0) #C(-13.49d0 32.5d0)))))
(BLAS-SWAP V2 V1)
(LIST (GRID:COPY-TO V1) (GRID:COPY-TO V2))))))
(LIST (GRID:COPY-TO V1 'array '(complex double-float)) (GRID:COPY-TO V2 'array '(complex double-float)))))))
;; Regression test COLUMN for GSLL, automatically generated
;;
;; Copyright 2009, 2010, 2011 Liam M. Healy
;; Copyright 2009, 2010, 2011, 2014 Liam M. Healy
;; Distributed under the terms of the GNU General Public License
;;
;; This program is free software: you can redistribute it and/or modify
......@@ -30,7 +30,7 @@
(-8.93f0 34.12f0 -6.15f0)
(49.27f0 -13.49f0 32.5f0))))
(COL (GRID:MAKE-FOREIGN-ARRAY 'SINGLE-FLOAT :DIMENSIONS '3)))
(GRID:COPY-TO (COLUMN M1 1 COL)))))
(GRID:COPY-TO (COLUMN M1 1 COL) 'array 'single-float))))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST #(8.24d0 34.12d0 -13.49d0))
(MULTIPLE-VALUE-LIST
......@@ -41,7 +41,7 @@
(-8.93d0 34.12d0 -6.15d0)
(49.27d0 -13.49d0 32.5d0))))
(COL (GRID:MAKE-FOREIGN-ARRAY 'DOUBLE-FLOAT :DIMENSIONS '3)))
(GRID:COPY-TO (COLUMN M1 1 COL)))))
(GRID:COPY-TO (COLUMN M1 1 COL) 'array 'double-float))))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST #(#C(3.29f0 -8.93f0) #C(-6.15f0 49.27f0) #C(32.5f0 42.73f0)))
(MULTIPLE-VALUE-LIST
......@@ -55,7 +55,7 @@
(COL
(GRID:MAKE-FOREIGN-ARRAY
'(COMPLEX SINGLE-FLOAT) :DIMENSIONS '3)))
(GRID:COPY-TO (COLUMN M1 1 COL)))))
(GRID:COPY-TO (COLUMN M1 1 COL) 'array '(complex single-float)))))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST
#(#C(3.29d0 -8.93d0) #C(-6.15d0 49.27d0)
......@@ -70,7 +70,7 @@
(#C(49.27d0 -13.49d0) #C(32.5d0 42.73d0) #C(-17.24d0 43.31d0)))))
(COL
(GRID:MAKE-FOREIGN-ARRAY '(COMPLEX DOUBLE-FLOAT) :DIMENSIONS '3)))
(GRID:COPY-TO (COLUMN M1 1 COL)))))
(GRID:COPY-TO (COLUMN M1 1 COL) 'array '(complex double-float)))))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST #(-68 52 -5))
(MULTIPLE-VALUE-LIST
......@@ -85,7 +85,7 @@
(GRID:MAKE-FOREIGN-ARRAY
'(SIGNED-BYTE 8)
:DIMENSIONS '3)))
(GRID:COPY-TO (COLUMN M1 1 COL)))))
(GRID:COPY-TO (COLUMN M1 1 COL) 'array '(signed-byte 8)))))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST #(44 163 215))
(MULTIPLE-VALUE-LIST
......@@ -100,7 +100,7 @@
(GRID:MAKE-FOREIGN-ARRAY
'(UNSIGNED-BYTE 8)
:DIMENSIONS '3)))
(GRID:COPY-TO (COLUMN M1 1 COL)))))
(GRID:COPY-TO (COLUMN M1 1 COL) 'array '(unsigned-byte 8)))))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST #(-68 52 -5))
(MULTIPLE-VALUE-LIST
......@@ -114,7 +114,7 @@
(COL
(GRID:MAKE-FOREIGN-ARRAY
'(SIGNED-BYTE 16) :DIMENSIONS '3)))
(GRID:COPY-TO (COLUMN M1 1 COL)))))
(GRID:COPY-TO (COLUMN M1 1 COL) 'array '(signed-byte 16)))))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST #(44 163 215))
(MULTIPLE-VALUE-LIST
......@@ -129,7 +129,7 @@
(GRID:MAKE-FOREIGN-ARRAY
'(UNSIGNED-BYTE 16)
:DIMENSIONS '3)))
(GRID:COPY-TO (COLUMN M1 1 COL)))))
(GRID:COPY-TO (COLUMN M1 1 COL) 'array '(unsigned-byte 16)))))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST #(-68 52 -5))
(MULTIPLE-VALUE-LIST
......@@ -143,7 +143,7 @@
(COL
(GRID:MAKE-FOREIGN-ARRAY
'(SIGNED-BYTE 32) :DIMENSIONS '3)))
(GRID:COPY-TO (COLUMN M1 1 COL)))))
(GRID:COPY-TO (COLUMN M1 1 COL) 'array '(signed-byte 32)))))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST #(44 163 215))
(MULTIPLE-VALUE-LIST
......@@ -157,7 +157,7 @@
(COL
(GRID:MAKE-FOREIGN-ARRAY
'(UNSIGNED-BYTE 32) :DIMENSIONS '3)))
(GRID:COPY-TO (COLUMN M1 1 COL)))))
(GRID:COPY-TO (COLUMN M1 1 COL) 'array '(unsigned-byte 32)))))
#+int64
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST #(-68 52 -5))
......@@ -173,7 +173,7 @@
(GRID:MAKE-FOREIGN-ARRAY
'(SIGNED-BYTE 64)
:DIMENSIONS '3)))
(GRID:COPY-TO (COLUMN M1 1 COL)))))
(GRID:COPY-TO (COLUMN M1 1 COL) 'array '(signed-byte 64)))))
#+int64
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST #(44 163 215))
......@@ -189,5 +189,5 @@
(GRID:MAKE-FOREIGN-ARRAY
'(UNSIGNED-BYTE 64)
:DIMENSIONS '3)))
(GRID:COPY-TO (COLUMN M1 1 COL))))))
(GRID:COPY-TO (COLUMN M1 1 COL) 'array '(unsigned-byte 64))))))
;; Regression test COMBINATION for GSLL, automatically generated
;;
;; Copyright 2009 Liam M. Healy
;; Copyright 2009, 2014 Liam M. Healy
;; Distributed under the terms of the GNU General Public License
;;
;; This program is free software: you can redistribute it and/or modify
......@@ -34,14 +34,14 @@
(MULTIPLE-VALUE-LIST
(LET ((COMB (MAKE-COMBINATION 4 2)))
(INIT-FIRST COMB)
(LOOP COLLECT (COPY-SEQ (GRID:COPY-TO COMB)) WHILE
(LOOP COLLECT (COPY-SEQ (GRID:COPY-TO COMB 'array 'fixnum)) WHILE
(COMBINATION-NEXT COMB)))))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST (LIST #(2 3) #(1 3) #(1 2) #(0 3) #(0 2) #(0 1)))
(MULTIPLE-VALUE-LIST
(LET ((COMB (MAKE-COMBINATION 4 2)))
(INIT-LAST COMB)
(LOOP COLLECT (COPY-SEQ (GRID:COPY-TO COMB)) WHILE
(LOOP COLLECT (COPY-SEQ (GRID:COPY-TO COMB 'array 'fixnum)) WHILE
(COMBINATION-PREVIOUS COMB)))))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST
......
;; Regression test COULOMB for GSLL, automatically generated
;;
;; Copyright 2009 Liam M. Healy
;; Copyright 2009, 2014 Liam M. Healy
;; Distributed under the terms of the GNU General Public License
;;
;; This program is free software: you can redistribute it and/or modify
......@@ -40,7 +40,7 @@
(MULTIPLE-VALUE-LIST
(LET ((ARR (GRID:MAKE-FOREIGN-ARRAY 'DOUBLE-FLOAT :DIMENSIONS 3)))
(COULOMB-WAVE-F-ARRAY 0.0d0 1.0d0 2.0d0 ARR)
(GRID:COPY-TO ARR))))
(GRID:COPY-TO ARR 'array 'double-float))))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST 0.07161779967468254d0 0.1278101031568499d0
2.2510703464871114d0 -1.4245543587641651d0 0.0d0
......@@ -58,8 +58,8 @@
(LET ((FARR (GRID:MAKE-FOREIGN-ARRAY 'DOUBLE-FLOAT :DIMENSIONS 3))
(GARR (GRID:MAKE-FOREIGN-ARRAY 'DOUBLE-FLOAT :DIMENSIONS 3)))
(COULOMB-WAVE-FG-ARRAY 1.5d0 1.0d0 1.0d0 FARR GARR)
(APPEND (COERCE (GRID:COPY-TO FARR) 'LIST)
(COERCE (GRID:COPY-TO GARR) 'LIST)))))
(APPEND (COERCE (GRID:COPY-TO FARR 'array 'double-float) 'LIST)
(COERCE (GRID:COPY-TO GARR 'array 'double-float) 'LIST)))))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST
#(0.33089080691634065d0 0.18070642887252675d0
......@@ -67,7 +67,7 @@
(MULTIPLE-VALUE-LIST
(LET ((ARR (GRID:MAKE-FOREIGN-ARRAY 'DOUBLE-FLOAT :DIMENSIONS 3)))
(COULOMB-WAVE-SPHF-ARRAY 0.0d0 1.0d0 2.0d0 ARR)
(GRID:COPY-TO ARR))))
(GRID:COPY-TO ARR 'array 'double-float))))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST 0.0013809146441856027d0 2.759621819430441d-17)
(MULTIPLE-VALUE-LIST (COULOMB-CL 1.0d0 2.5d0)))
......@@ -78,5 +78,5 @@
(MULTIPLE-VALUE-LIST
(LET ((CL (GRID:MAKE-FOREIGN-ARRAY 'DOUBLE-FLOAT :DIMENSIONS 3)))
(COULOMB-cl-array 0.0d0 1.0d0 CL)
(GRID:COPY-TO CL)))))
(GRID:COPY-TO CL 'array 'double-float)))))
;; Regression test DIRICHLET for GSLL, automatically generated
;;
;; Copyright 2009 Liam M. Healy
;; Copyright 2009, 2014 Liam M. Healy
;; Distributed under the terms of the GNU General Public License
;;
;; This program is free software: you can redistribute it and/or modify
......@@ -29,7 +29,7 @@
(GRID:MAKE-FOREIGN-ARRAY
'DOUBLE-FLOAT :INITIAL-CONTENTS
'(1.0d0 2.0d0 3.0d0 4.0d0))))
(grid:copy-to (sample rng :dirichlet :alpha alpha)))))
(grid:copy-to (sample rng :dirichlet :alpha alpha) 'array 'double-float))))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(list 3.483648000000004d7)
(MULTIPLE-VALUE-LIST
......
;; Regression test GEGENBAUER for GSLL, automatically generated
;;
;; Copyright 2009 Liam M. Healy
;; Copyright 2009, 2014 Liam M. Healy
;; Distributed under the terms of the GNU General Public License
;;
;; This program is free software: you can redistribute it and/or modify
......@@ -36,5 +36,5 @@
(MULTIPLE-VALUE-LIST
(LET ((ARR (GRID:MAKE-FOREIGN-ARRAY 'DOUBLE-FLOAT :DIMENSIONS 4)))
(GEGENBAUER-ARRAY 1.0d0 3.0d0 ARR)
(GRID:COPY-TO ARR)))))
(GRID:COPY-TO ARR 'array 'double-float)))))
;; Regression test GIVENS for GSLL, automatically generated
;;
;; Copyright 2009, 2011 Liam M. Healy
;; Copyright 2009, 2011, 2014 Liam M. Healy
;; Distributed under the terms of the GNU General Public License
;;
;; This program is free software: you can redistribute it and/or modify
......@@ -50,7 +50,7 @@
(SETF (GRID:AREF COSINES I)
(COS (GRID:AREF ANGLES I))))
(GIVENS-ROTATION V1 V2 COSINES SINES)
(LIST (GRID:COPY-TO V1) (GRID:COPY-TO V2)))))
(LIST (GRID:COPY-TO V1 'array 'single-float) (GRID:COPY-TO V2 'array 'single-float)))))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST
(LIST
......@@ -83,5 +83,5 @@
(SETF (GRID:AREF COSINES I)
(COS (GRID:AREF ANGLES I))))
(GIVENS-ROTATION V1 V2 COSINES SINES)
(LIST (GRID:COPY-TO V1) (GRID:COPY-TO V2))))))
(LIST (GRID:COPY-TO V1 'array 'double-float) (GRID:COPY-TO V2 'array 'double-float))))))
;; Regression test INVERSE-MATRIX-PRODUCT for GSLL, automatically generated
;;
;; Copyright 2009, 2010, 2011 Liam M. Healy
;; Copyright 2009, 2010, 2011, 2014 Liam M. Healy
;; Distributed under the terms of the GNU General Public License
;;
;; This program is free software: you can redistribute it and/or modify
......@@ -38,7 +38,7 @@
(-16.12f0 -8.25f0 21.44f0)
(-49.08f0 -39.66f0 -49.46f0))))
(S1 19.68f0))
(GRID:COPY-TO (INVERSE-MATRIX-PRODUCT M1 M2 S1)))))
(GRID:COPY-TO (INVERSE-MATRIX-PRODUCT M1 M2 S1) 'array 'single-float))))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST
#2A((-30.708968222785604d0 5.373707833953741d0
......@@ -61,7 +61,7 @@
(-16.12d0 -8.25d0 21.44d0)
(-49.08d0 -39.66d0 -49.46d0))))
(S1 19.68d0))
(GRID:COPY-TO (INVERSE-MATRIX-PRODUCT M1 M2 S1)))))
(GRID:COPY-TO (INVERSE-MATRIX-PRODUCT M1 M2 S1) 'array 'double-float))))
#+fsbv
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST
......@@ -87,7 +87,7 @@
(#C(-16.12f0 -8.25f0) #C(21.44f0 -49.08f0) #C(-39.66f0 -49.46f0))
(#C(-49.08f0 -39.66f0) #C(-49.46f0 19.68f0) #C(-5.55f0 -8.82f0)))))
(S1 #C(19.68f0 -5.55f0)))
(GRID:COPY-TO (INVERSE-MATRIX-PRODUCT M1 M2 S1)))))
(GRID:COPY-TO (INVERSE-MATRIX-PRODUCT M1 M2 S1) 'array '(complex single-float)))))
#+fsbv
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST
......@@ -116,7 +116,7 @@
(#C(-16.12d0 -8.25d0) #C(21.44d0 -49.08d0) #C(-39.66d0 -49.46d0))
(#C(-49.08d0 -39.66d0) #C(-49.46d0 19.68d0) #C(-5.55d0 -8.82d0)))))
(S1 #C(19.68d0 -5.55d0)))
(GRID:COPY-TO (INVERSE-MATRIX-PRODUCT M1 M2 S1)))))
(GRID:COPY-TO (INVERSE-MATRIX-PRODUCT M1 M2 S1) 'array '(complex double-float)))))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST #(-1.1747805f0 -0.26507667f0 1.3326154f0))
(MULTIPLE-VALUE-LIST
......@@ -131,7 +131,7 @@
'SINGLE-FLOAT :INITIAL-CONTENTS
'(42.73f0 -17.24f0 43.31f0)))
(S1 -16.12f0))
(GRID:COPY-TO (INVERSE-MATRIX-PRODUCT M1 V1 S1)))))
(GRID:COPY-TO (INVERSE-MATRIX-PRODUCT M1 V1 S1) 'array 'single-float))))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST
#(-1.1747804927980594d0 -0.26507665253855167d0
......@@ -147,7 +147,7 @@
(GRID:MAKE-FOREIGN-ARRAY 'DOUBLE-FLOAT :INITIAL-CONTENTS
'(42.73d0 -17.24d0 43.31d0)))
(S1 -16.12d0))
(GRID:COPY-TO (INVERSE-MATRIX-PRODUCT M1 V1 S1)))))
(GRID:COPY-TO (INVERSE-MATRIX-PRODUCT M1 V1 S1) 'array 'double-float))))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST
#(#C(-1.1218005f0 0.24990065f0) #C(-0.7679142f0 -0.9143634f0)
......@@ -166,7 +166,7 @@
:INITIAL-CONTENTS
'(#C(42.73f0 -17.24f0) #C(43.31f0 -16.12f0) #C(-8.25f0 21.44f0))))
(S1 #C(-16.12f0 -8.25f0)))
(GRID:COPY-TO (INVERSE-MATRIX-PRODUCT M1 V1 S1)))))
(GRID:COPY-TO (INVERSE-MATRIX-PRODUCT M1 V1 S1) 'array '(complex single-float)))))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST
#(#C(-1.1218005239477833d0 0.2499006037347623d0)
......@@ -186,5 +186,5 @@
:INITIAL-CONTENTS
'(#C(42.73d0 -17.24d0) #C(43.31d0 -16.12d0) #C(-8.25d0 21.44d0))))
(S1 #C(-16.12d0 -8.25d0)))
(GRID:COPY-TO (INVERSE-MATRIX-PRODUCT M1 V1 S1))))))
(GRID:COPY-TO (INVERSE-MATRIX-PRODUCT M1 V1 S1) 'array '(complex double-float))))))
;; Regression test LEGENDRE for GSLL, automatically generated
;;
;; Copyright 2009 Liam M. Healy
;; Copyright 2009, 2014 Liam M. Healy
;; Distributed under the terms of the GNU General Public License
;;
;; This program is free software: you can redistribute it and/or modify
......@@ -41,7 +41,7 @@
(MULTIPLE-VALUE-LIST
(LET ((ARR (GRID:MAKE-FOREIGN-ARRAY 'DOUBLE-FLOAT :DIMENSIONS 4)))
(LEGENDRE-PL-ARRAY 0.5d0 ARR)
(GRID:COPY-TO ARR))))
(GRID:COPY-TO ARR 'array 'double-float))))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST 0.3128529498822064d0 1.3893461931245028d-16)
(MULTIPLE-VALUE-LIST (LEGENDRE-Q0 3.3d0)))
......@@ -59,7 +59,7 @@
(MULTIPLE-VALUE-LIST
(LET ((ARR (GRID:MAKE-FOREIGN-ARRAY 'DOUBLE-FLOAT :DIMENSIONS 4)))
(LEGENDRE-PLM-ARRAY 2 0.5d0 ARR)
(GRID:COPY-TO ARR))))
(GRID:COPY-TO ARR 'array 'double-float))))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST #(-3.0d0 3.75d0 33.75d0 55.78125d0))
(MULTIPLE-VALUE-LIST
......@@ -67,7 +67,7 @@
(DERIV
(GRID:MAKE-FOREIGN-ARRAY 'DOUBLE-FLOAT :DIMENSIONS 4)))
(LEGENDRE-PLM-DERIV-ARRAY 2 0.5d0 VAL DERIV)
(GRID:COPY-TO DERIV))))
(GRID:COPY-TO DERIV 'array 'double-float))))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST 0.30366280894310793d0 3.5267592419993454d-14)
(MULTIPLE-VALUE-LIST
......@@ -79,7 +79,7 @@
(MULTIPLE-VALUE-LIST
(LET ((ARR (GRID:MAKE-FOREIGN-ARRAY 'DOUBLE-FLOAT :DIMENSIONS 4)))
(LEGENDRE-SPHPLM-ARRAY 4 0.5d0 ARR)
(GRID:COPY-TO ARR))))
(GRID:COPY-TO ARR 'array 'double-float))))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST
#(-0.6637990386674741d0 -0.2751965434323283d0
......@@ -89,7 +89,7 @@
(DERIV
(GRID:MAKE-FOREIGN-ARRAY 'DOUBLE-FLOAT :DIMENSIONS 4)))
(LEGENDRE-SPHPLM-DERIV-ARRAY 4 0.5d0 VAL DERIV)
(GRID:COPY-TO DERIV))))
(GRID:COPY-TO DERIV 'array 'double-float))))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST -0.1255299048878925d0 1.3395992025650077d-15)
(MULTIPLE-VALUE-LIST
......@@ -131,5 +131,5 @@
(MULTIPLE-VALUE-LIST
(LET ((ARR (GRID:MAKE-FOREIGN-ARRAY 'DOUBLE-FLOAT :DIMENSIONS 4)))
(LEGENDRE-H3D-ARRAY 1.0d0 0.5d0 ARR)
(GRID:COPY-TO ARR)))))
(GRID:COPY-TO ARR 'array 'double-float)))))
;; Regression test LU for GSLL, automatically generated
;; with some manual changes to the results
;;
;; Copyright 2009, 2011 Liam M. Healy
;; Copyright 2009, 2011, 2014 Liam M. Healy
;; Distributed under the terms of the GNU General Public License
;;
;; This program is free software: you can redistribute it and/or modify
......@@ -69,7 +69,7 @@
#C(-7.88941207561151d+04 1.95053812987858d+04)
#C(3.95548551241728d+04 -7.76593696255317d+03))))
(MULTIPLE-VALUE-LIST
(TEST-LU-SOLVE-DIM (CREATE-COMPLEX-MATRIX 7)))))
(TEST-LU-SOLVE-DIM (CREATE-COMPLEX-MATRIX 7))))
(LISP-UNIT:ASSERT-NUMERICAL-EQUAL
(LIST
#(-39.65999999999999d0 -49.46000000000001d0
......@@ -97,7 +97,8 @@
MATRIX X 1 :UPPER :NOTRANS
:NONUNIT)
1 :LOWER :NOTRANS
:UNIT))))))))
:UNIT))
'array 'double-float))))))
(LISP-UNIT:ASSERT-NUMERICAL-EQUAL
(LIST
#(#C(-39.65999999999999d0 -49.46000000000001d0)
......@@ -129,4 +130,5 @@
(MATRIX-PRODUCT-TRIANGULAR
MATRIX
(MATRIX-PRODUCT-TRIANGULAR MATRIX X 1 :UPPER :NOTRANS :NONUNIT)
1 :LOWER :NOTRANS :UNIT)))))))))
1 :LOWER :NOTRANS :UNIT))
'array '(complex double-float))))))))
;; Regression test MATHIEU for GSLL, automatically generated
;;
;; Copyright 2009, 2010 Liam M. Healy
;; Copyright 2009, 2010, 2014 Liam M. Healy
;; Distributed under the terms of the GNU General Public License
;;
;; This program is free software: you can redistribute it and/or modify
......@@ -315,7 +315,7 @@
-1.8369701987210297d-16 1.0d0
3.061616997868383d-16))
(MULTIPLE-VALUE-LIST
(GRID:COPY-TO (MATHIEU-CE-ARRAY 0.0d0 (/ dPI 2) 6)))))
(GRID:COPY-TO (MATHIEU-CE-ARRAY 0.0d0 (/ dPI 2) 6) 'array 'double-float))))
(let ((lisp-unit:*epsilon* 1.0d-6))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST
......@@ -328,7 +328,7 @@
1.0767085541408374d0 1.0641859681403019d0
1.054606963632144d0 1.0470843441628828d0))
(MULTIPLE-VALUE-LIST
(GRID:COPY-TO (MATHIEU-CE-ARRAY 20.0d0 0.0d0 16)))))
(GRID:COPY-TO (MATHIEU-CE-ARRAY 20.0d0 0.0d0 16) 'array 'double-float))))
(let ((lisp-unit:*epsilon* 1.0d-6))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST
......@@ -341,5 +341,5 @@
0.9433025704899831d0 8.979725089944524d-16
-0.9570452540612865d0))
(MULTIPLE-VALUE-LIST
(GRID:COPY-TO (MATHIEU-SE-ARRAY 20.0d0 (/ dPI 2) 15 1))))))
(GRID:COPY-TO (MATHIEU-SE-ARRAY 20.0d0 (/ dPI 2) 15 1) 'array 'double-float)))))
;; Regression test MATRIX-ADD-SCALAR for GSLL, automatically generated
;;
;; Copyright 2009 Liam M. Healy
;; Copyright 2009, 2014 Liam M. Healy
;; Distributed under the terms of the GNU General Public License
;;
;; This program is free software: you can redistribute it and/or modify
......@@ -21,16 +21,16 @@
(LISP-UNIT:DEFINE-TEST MATRIX-ADD-SCALAR
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST
#2A((-16.31 26.43 21.48)
(9.259999 52.309998 12.04)
(67.46 4.7000003 50.69)))
#2A((-16.31f0 26.43f0 21.48f0)
(9.259999f0 52.309998f0 12.04f0)
(67.46f0 4.7000003f0 50.69f0)))
(MULTIPLE-VALUE-LIST
(LET ((M1
(GRID:MAKE-FOREIGN-ARRAY 'SINGLE-FLOAT :INITIAL-CONTENTS
'((-34.5 8.24 3.29)
(-8.93 34.12 -6.15)
(49.27 -13.49 32.5)))))
(GRID:COPY-TO (ELT+ M1 18.19d0)))))
(GRID:COPY-TO (ELT+ M1 18.19d0) 'array 'single-float))))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST
#2A((-16.31d0 26.43d0 21.48d0)
......@@ -43,7 +43,7 @@
'((-34.5d0 8.24d0 3.29d0)
(-8.93d0 34.12d0 -6.15d0)
(49.27d0 -13.49d0 32.5d0)))))
(GRID:COPY-TO (ELT+ M1 18.19d0)))))
(GRID:COPY-TO (ELT+ M1 18.19d0) 'array 'double-float))))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST #2A((-45 -49 89) (-72 70 8) (91 13 -115)))
(MULTIPLE-VALUE-LIST
......@@ -51,7 +51,7 @@
(GRID:MAKE-FOREIGN-ARRAY '(SIGNED-BYTE 8) :INITIAL-CONTENTS
'((-64 -68 71) (-91 52 -10)
(73 -5 123)))))
(GRID:COPY-TO (ELT+ M1 18.19d0)))))
(GRID:COPY-TO (ELT+ M1 18.19d0) 'array '(signed-byte 8)))))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST #2A((85 62 207) (134 181 158) (179 233 116)))
(MULTIPLE-VALUE-LIST
......@@ -60,7 +60,7 @@
:INITIAL-CONTENTS
'((67 44 189) (116 163 140)
(161 215 98)))))
(GRID:COPY-TO (ELT+ M1 18.19d0)))))
(GRID:COPY-TO (ELT+ M1 18.19d0) 'array '(unsigned-byte 8)))))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST #2A((-45 -49 89) (-72 70 8) (91 13 141)))
(MULTIPLE-VALUE-LIST
......@@ -69,7 +69,7 @@
:INITIAL-CONTENTS
'((-64 -68 71) (-91 52 -10)
(73 -5 123)))))
(GRID:COPY-TO (ELT+ M1 18.19d0)))))
(GRID:COPY-TO (ELT+ M1 18.19d0) 'array '(signed-byte 16)))))
(LISP-UNIT::ASSERT-NUMERICAL-EQUAL
(LIST #2A((85 62 207) (134 181 158) (179 233 116)))
(MULTIPLE-VALUE-LIST
......@@ -78,7 +78,7 @@
:INITIAL-CONTENTS
'((67 44 189) (116 163 140)
(161 215 98)))))