tebdol

commit 1f7921ed6268a66b5ba9b71cc096f723ccfa2cdb
parent 16d0947435afaab8ddaf9d66078e108e0d96e483
Author: Miroslav Urbanek <mu@miroslavurbanek.com>
Date:   Wed, 28 Jun 2017 09:24:36 +0200

Update the ZGESDD arguments

Required dimension of RWORK has changed in the LAPACK 3.6.1
documentation. Also call make-double-array instead of make-array to
simplify the code.

Diffstat:
tebdol/array.lisp
 | 
7
++++---

1 file changed, 4 insertions(+), 3 deletions(-)
diff --git a/tebdol/array.lisp b/tebdol/array.lisp
@@ -174,13 +174,14 @@
 	   (n (reduce #'* dl))
 	   (min (min m n))
 	   (max (max m n))
-	   (s (make-array min :element-type 'double-float))
+	   (s (make-double-array min))
 	   (u (make-blas-array (append (list min) dr)))
 	   (vt (make-blas-array (append dl (list min))))
 	   (work (make-blas-array 1))
 	   (lwork -1)
-	   (rwork (make-array (* min (max (+ 7 (* 5 min)) (+ (* 2 max) (* 2 min) 1)))
-			      :element-type 'double-float))
+	   (rwork (make-double-array
+		   (* min (max (+ (* 5 min) 5)
+			       (+ (* 2 max) (* 2 min) 1)))))
 	   (iwork (make-integer-array (* 8 min))))
 
       ;; svd