[Gcl-devel] clines and dgemm

[Vanuxem Grégory]

Hi,

I can't use dgemm with clines.
One of the arrays is destroy.
Is it possible to give a 2D long-float array
to an external library ?

Here is my test:
(from
http://lists.gnu.org/archive/html/gcl-devel/2005-05/msg00037.html)


(compile-file "lapack.lisp" :system-p t :c-file t)

# On debian with atlas 
(compiler::link (list "lapack.o") "saved_gcl" "" "-L/usr/lib/atlas
-llapack -lblas -latlas -lg2c" t)

(quit)
--------------------------------------------
./saved_gcl
(setq a (make-array '(500 500) :element-type
'long-float :initial-element 3.2))

(setq a (make-array '(500 500) :element-type
'long-float :initial-element 3.2))

(setq b (make-array '(500 500) :element-type
'long-float :initial-element 7.4))

(setq c (make-array '(500 500) :element-type
'long-float :initial-element 0.0))

(dgemm a b c)
(aref c 0 0)
(dgemm a b c)
(aref c 0 0)
b

Destroyed array with gcl-2.6.5 and gcl-2.6.6
from debian testing (emulated 32bit and 64bits (AMD))

Hope that this can be reproduced easily if it's a bug.

Cheers,

Greg

PJ: lapack.lisp (ddot and dgemm's test
(square row-major matrix))

lapack.lisp
Description: Text document