[SciPy-user] Disable ATLAS
Nils Wagner
dr.nwagner at web.de
Sat Feb 21 12:32:15 EST 2004
SciPy Users List <scipy-user at scipy.net> schrieb am 21.02.04 14:35:12:
On Sat, 21 Feb 2004, Nils Wagner wrote:
> I have disabled ATLAS. The problem with linalg.lu didn't vanish.
> What can we do ?
Are you absolutely sure that ATLAS is not being used (rebuilded linalg
from scratch, removed earlier installations, eg `ldd blas.so` does not
contain references to ATLAS symbols, etc etc)?
Yes.
nwagner at linux:/usr/local/lib/python2.3/site-packages/scipy/linalg> ldd cblas.so
libg2c.so.0 => /usr/lib/libg2c.so.0 (0x40024000)
libm.so.6 => /lib/i686/libm.so.6 (0x40042000)
libgcc_s.so.1 => /lib/libgcc_s.so.1 (0x40065000)
libc.so.6 => /lib/i686/libc.so.6 (0x4006d000)
/lib/ld-linux.so.2 => /lib/ld-linux.so.2 (0x80000000)
Have you tried the two examples of mine from the earlier thread?
Yes. they work fine for me. but lu.py failed again
Please can you run lu.py
from scipy import *
from RandomArray import *
n = 2
A = rand(n,n)+rand(n,n)*1j
X = zeros((n,n),Complex)
Ainv = linalg.inv(A)
R = identity(n)+identity(n)*0j
#
# linalg.solve cannot handle complex rhs !!
#
#R = R.real
#
for i in arange(0,n):
r = R[:,i]
X[:,i] = linalg.solve(A,r)
#
# Should be the inverse of A as well !!!!
#
print
print 'X'
print
print X
#
# Inverse of A
#
print
print 'inv(A)'
print
print Ainv
print
print 'Identity ???'
print
print dot(A,Ainv)
print
print 'Identity ???'
print
print dot(A,X)
Thanks in advance
Nils
Pearu
_______________________________________________
SciPy-user mailing list
SciPy-user at scipy.net
http://www.scipy.net/mailman/listinfo/scipy-user
______________________________________________________________________________
... and the winner is... WEB.DE FreeMail! - Deutschlands beste E-Mail
ist zum 39. Mal Testsieger (PC Praxis 03/04) http://f.web.de/?mc=021191
More information about the SciPy-User
mailing list