[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[Gcl-devel] clines and dgemm
From: |
Vanuxem Grégory |
Subject: |
[Gcl-devel] clines and dgemm |
Date: |
Tue, 31 May 2005 00:26:28 +0200 |
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
- [Gcl-devel] clines and dgemm,
Vanuxem Grégory <=