Skip to content
GitLab
Projects
Groups
Snippets
/
Help
Help
Support
Community forum
Keyboard shortcuts
?
Submit feedback
Contribute to GitLab
Sign in / Register
Toggle navigation
Menu
Open sidebar
antik
gsll
Commits
a286ea0d
Commit
a286ea0d
authored
Mar 06, 2014
by
Liam M. Healy
Browse files
Balance parentheses in test
parent
003c83f8
Changes
1
Hide whitespace changes
Inline
Side-by-side
tests/lu.lisp
View file @
a286ea0d
...
...
@@ -24,10 +24,10 @@
;;; GSL has #define GSL_DBL_EPSILON 2.2204460492503131e-16
;;; which is 2x what double-float-epsilon is.
(
LISP-UNIT:DEFINE-TEST
LU
(
let
((
lisp-unit:*epsilon*
(
*
2
8
double-float-epsilon
)))
(
LISP-UNIT:ASSERT-NUMERICAL-EQUAL
(
LIST
*hilb2-soln*
)
(
MULTIPLE-VALUE-LIST
(
TEST-LU-SOLVE-DIM
*hilb2*
))))
(
let
((
lisp-unit:*epsilon*
(
*
2
8
double-float-epsilon
)))
(
LISP-UNIT:ASSERT-NUMERICAL-EQUAL
(
LIST
*hilb2-soln*
)
(
MULTIPLE-VALUE-LIST
(
TEST-LU-SOLVE-DIM
*hilb2*
))))
(
let
((
lisp-unit:*epsilon*
(
*
2
64
double-float-epsilon
)))
(
LISP-UNIT:ASSERT-NUMERICAL-EQUAL
(
LIST
*hilb3-soln*
)
...
...
@@ -69,22 +69,22 @@
#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
19.679999999999993d0
-5.549999999999997d0
))
(
MULTIPLE-VALUE-LIST
(
LET
((
MATRIX
(
GRID:MAKE-FOREIGN-ARRAY
'DOUBLE-FLOAT
:INITIAL-CONTENTS
'
((
-34.5d0
8.24d0
3.29d0
-8.93d0
)
(
34.12d0
-6.15d0
49.27d0
-13.49d0
)
(
32.5d0
42.73d0
-17.24d0
43.31d0
)
(
-16.12d0
-8.25d0
21.44d0
-49.08d0
))))
(
GRID:MAKE-FOREIGN-ARRAY
'DOUBLE-FLOAT
:INITIAL-CONTENTS
'
((
-34.5d0
8.24d0
3.29d0
-8.93d0
)
(
34.12d0
-6.15d0
49.27d0
-13.49d0
)
(
32.5d0
42.73d0
-17.24d0
43.31d0
)
(
-16.12d0
-8.25d0
21.44d0
-49.08d0
))))
(
VEC
(
GRID:MAKE-FOREIGN-ARRAY
'DOUBLE-FLOAT
:INITIAL-CONTENTS
'
(
-39.66d0
-49.46d0
19.68d0
-5.55d0
))))
(
GRID:MAKE-FOREIGN-ARRAY
'DOUBLE-FLOAT
:INITIAL-CONTENTS
'
(
-39.66d0
-49.46d0
19.68d0
-5.55d0
))))
(
MULTIPLE-VALUE-BIND
(
MATRIX
PERM
)
(
LU-DECOMPOSITION
MATRIX
)
...
...
@@ -107,19 +107,19 @@
#C
(
-30.580000000000002d0
31.67d0
)))
(
MULTIPLE-VALUE-LIST
(
LET
((
MATRIX
(
GRID:MAKE-FOREIGN-ARRAY
'
(
COMPLEX
DOUBLE-FLOAT
)
:INITIAL-CONTENTS
'
((
#C
(
-34.5d0
8.24d0
)
#C
(
3.29d0
-8.93d0
)
#C
(
34.12d0
-6.15d0
)
#C
(
49.27d0
-13.49d0
))
(
#C
(
34.12d0
-6.15d0
)
#C
(
49.27d0
-13.49d0
)
#C
(
32.5d0
42.73d0
)
#C
(
-17.24d0
43.31d0
))
(
#C
(
32.5d0
42.73d0
)
#C
(
-17.24d0
43.31d0
)
#C
(
-16.12d0
-8.25d0
)
#C
(
21.44d0
-49.08d0
))
(
#C
(
-16.12d0
-8.25d0
)
#C
(
21.44d0
-49.08d0
)
#C
(
-39.66d0
-49.46d0
)
#C
(
19.68d0
-5.55d0
)))))
(
GRID:MAKE-FOREIGN-ARRAY
'
(
COMPLEX
DOUBLE-FLOAT
)
:INITIAL-CONTENTS
'
((
#C
(
-34.5d0
8.24d0
)
#C
(
3.29d0
-8.93d0
)
#C
(
34.12d0
-6.15d0
)
#C
(
49.27d0
-13.49d0
))
(
#C
(
34.12d0
-6.15d0
)
#C
(
49.27d0
-13.49d0
)
#C
(
32.5d0
42.73d0
)
#C
(
-17.24d0
43.31d0
))
(
#C
(
32.5d0
42.73d0
)
#C
(
-17.24d0
43.31d0
)
#C
(
-16.12d0
-8.25d0
)
#C
(
21.44d0
-49.08d0
))
(
#C
(
-16.12d0
-8.25d0
)
#C
(
21.44d0
-49.08d0
)
#C
(
-39.66d0
-49.46d0
)
#C
(
19.68d0
-5.55d0
)))))
(
VEC
(
GRID:MAKE-FOREIGN-ARRAY
'
(
COMPLEX
DOUBLE-FLOAT
)
:INITIAL-CONTENTS
'
(
#C
(
-39.66d0
-49.46d0
)
#C
(
19.68d0
-5.55d0
)
#C
(
-8.82d0
25.37d0
)
#C
(
-30.58d0
31.67d0
)))))
(
GRID:MAKE-FOREIGN-ARRAY
'
(
COMPLEX
DOUBLE-FLOAT
)
:INITIAL-CONTENTS
'
(
#C
(
-39.66d0
-49.46d0
)
#C
(
19.68d0
-5.55d0
)
#C
(
-8.82d0
25.37d0
)
#C
(
-30.58d0
31.67d0
)))))
(
MULTIPLE-VALUE-BIND
(
MATRIX
PERM
)
(
LU-DECOMPOSITION
MATRIX
)
...
...
@@ -131,4 +131,4 @@
MATRIX
(
MATRIX-PRODUCT-TRIANGULAR
MATRIX
X
1
:UPPER
:NOTRANS
:NONUNIT
)
1
:LOWER
:NOTRANS
:UNIT
))
'array
'
(
complex
double-float
))))))))
'array
'
(
complex
double-float
))))))))
Write
Preview
Supports
Markdown
0%
Try again
or
attach a new file
.
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Cancel
Please
register
or
sign in
to comment