Mailman 3 February 2004 - SciPy-Dev

Urgent: Bug in linalg.lu_factor linalg.lu_solve
by Nils Wagner 20 Feb '04

20 Feb '04

Hi all, I have observed a very strange behaviour within linalg.lu_factor / linalg.lu_solve. Any idea or suggestion ? Nils from scipy import * def F(A,l): tmp=zeros((n,n),Complex) tmp[:n,:n] = A+l*identity(n)+exp(-l)*identity(n) return tmp def Fs(l): tmp=zeros((n,n),Complex) tmp[:n,:n] = identity(n)-l*exp(-l)*identity(n) return tmp def T(n): tmp=zeros((n,n),Complex) for i in arange(0,n): tmp[i,i] = -2.0 if i < n-1: tmp[i+1,i] = 1.0 tmp[i,i+1] = 1.0 return tmp n = 4 A = T(n) gamma = 2.0+0j rho = 2.0 ns = 180 #for i in arange(0,ns): for i in arange(0,2): l = gamma+rho*exp(2.*pi*i*1j/ns) D1 = F(A,l) Ds = Fs(l) Dinv = linalg.inv(D1) lu, piv = linalg.lu_factor(D1,overwrite_a=0) tr = 0.0 for k in arange(0,n): rhs = Ds[:,k] g1 = linalg.lu_solve((lu,piv),rhs,trans=0,overwrite_b=0) tr = tr + g1[k] # # tr and trace(dot(Dinv,Ds)) should be the same !!!!! # print tr, trace(dot(Dinv,Ds))

1 0

"Building Scipy on Windows" notes available
by Pearu Peterson 19 Feb '04

19 Feb '04

Hi all, Detailed "from scratch" notes about "Building Scipy on Windows" are available at http://cens.ioc.ee/~pearu/scipy/BUILD_WIN32.html It also contains a section on building Win32 port of Atlas libraries that seems to be a major show-stopper for Windows people, and indeed, I found it a rather sensitive subject too. Hopefully these notes will help Windows people to overcome the threshold of building Scipy themself and may enjoy most recent Scipy features from CVS. Best regards, Pearu

3 2

[SciPy-dev] How to avoid arch specific compiler flags in scipy_distutils?
by José Fonseca 15 Feb '04

15 Feb '04

I've finally been doing some work updating the scipy & f2py debian packages and things are slowly progressing. The showstopper ATM, is the fact that the fortran sources are being compiled with architecture specific compiler flags ("-march=pentium4 -msse2 -msse ..."). I know that this is done automatically in scipy_distutils/gnufcompiler.py . Is there any non-intrusive way (i.e., from command line) to prevent that? José Fonseca

2 4