tebdol

Simulation of ultracold atoms in optical lattices
git clone https://miroslavurbanek.com/git/tebdol.git
Log | Files | Refs | README

commit 0de51d6170e4379061104c85481decb7219f7f97
parent ea1b5401f5594b2ccdd875c860fd3859adc7da16
Author: Miroslav Urbanek <mu@miroslavurbanek.com>
Date:   Wed, 28 Jun 2017 15:17:48 +0200

Clean up LAPACK functions

Fix definitions of LAPACK functions and types of arrays passed to
them.

Diffstat:
tebdol/blas.lisp | 8++++----
tebdol/exp.lisp | 4++--
2 files changed, 6 insertions(+), 6 deletions(-)

diff --git a/tebdol/blas.lisp b/tebdol/blas.lisp @@ -62,10 +62,10 @@ (define-alien-routine ("zaxpy_" zaxpy) void (n int :copy) - (a (* double)) - (x (* double)) + (za (* double)) + (zx (* double)) (incx int :copy) - (y (* double)) + (zy (* double)) (incy int :copy)) (define-alien-routine ("zgemm_" zgemm) void @@ -114,7 +114,7 @@ (work (* double)) (lwork int :copy) (rwork (* double)) - (iwork (* double)) + (iwork (* int)) (info int :out)) (define-alien-routine ("zgesvdx_" zgesvdx) void diff --git a/tebdol/exp.lisp b/tebdol/exp.lisp @@ -45,8 +45,8 @@ (m) (setf lwork (floor (realpart (aref work 0)))) (setf work (make-blas-array lwork)) - (setf lrwork (floor (realpart (aref rwork 0)))) - (setf rwork (make-blas-array lrwork)) + (setf lrwork (floor (aref rwork 0))) + (setf rwork (make-double-array lrwork)) (setf liwork (aref iwork 0)) (setf iwork (make-integer-array liwork)) (m))))