[SciPy-User] lstsq error under Windows?

[bruce.labitt at autoliv.com]

scipy-user-bounces at scipy.org wrote on 05/04/2010 04:10:04 PM:

> On Tue, May 4, 2010 at 3:48 PM,  <bruce.labitt at autoliv.com> wrote:
> > I have found an issue with scipy.linalg.lstsq, I think.  The following
> > code works in Ubuntu 10.04 x86-64, but not in WinXP-32.  I think it 
should
> > work in WinXP.  Here is a minimum example:
> >
> > """==== program testlstsq.py ======================"""
> > from scipy.linalg import eig, lstsq
> > from numpy import angle, zeros, pi, arcsin
> >
> > A = zeros((4,1), dtype='complex')
> > B = zeros((4,1), dtype='complex')
> >
> > A[0,0] = -0.535412460549-2.65798938848e-17j
> > A[1,0] = -0.369432866546-0.131765700574j
> > A[2,0] = -0.222906796932-0.263237285725j
> > A[3,0] = -0.069087096386-0.38609560454j
> >
> > B[0,0] = -0.369432866546-0.131765700574j
> > B[1,0] = -0.222906796932-0.263237285725j
> > B[2,0] = -0.069087096386-0.38609560454j
> > B[3,0] =  0.0882283631514-0.528093039953j
> >
> > try:
> >    print 'Got here'
> >    phi = lstsq(A,B)
> >    print 'Finished lstsq'
> > except:
> >    print 'Exception Occurred'
> > else:
> >    for a in range(len(phi)):
> >        print 'phi[',a,']=',phi[a]
> >
> >    w = -angle( eig(phi[0])[0][:] )
> >    d = 0.5
> >    aa = arcsin( w / (2.0*pi) )*180./pi     # in degrees
> >    print 'aa unsorted =', aa
> > """ ==== end of testlstsq.py =========================="""
> >
> > If this is run under ipython, or python on Ubuntu, the answer is:
> >
> > In [3]: run testlstsq.py
> > Got here
> > Finished lstsq
> > phi[ 0 ]= [[ 0.88271111+0.38811028j]]
> > phi[ 1 ]= [-0.04649206-0.01274722j]
> > phi[ 2 ]= 1
> > phi[ 3 ]= [ 0.84459073]
> > aoa unsorted = [-3.780145]
> >
> > If testlstsq.py is run under WinXP-32 one gets the following result:
> >> python testlstsq.py
> > Got here
> >  ** On entry to ZGELSS parameter number 12 had an illegal value
> >
> > I think that ZGELSS is in LAPACK.  After that, I am in over my head. 
 Can
> > someone help with this?
> >
> > Under Windows I am running:
> > WinXP-x86-32
> > Python(x,y)2.6.2.0 --> Python 2.6.2
> > numpy 1.3.0
> > scipy 0.7.1
> > ipython 0.10
> 
> no problem here, windowsXP, python 2.5, numpy 1.4.0, scipy 
0.8dev_something

Is this WinXP-32?

I notice you have an earlier python, and a later version of numpy and 
scipy.  I thought numpy 1.4.0 was "recalled" a while back.  I'll try 
numpy.linalg to see if there is a difference.

> 
> Got here
> Finished lstsq
> phi[ 0 ]= [[ 0.88271111+0.38811028j]]
> phi[ 1 ]= [-0.04649206-0.01274722j]
> phi[ 2 ]= 1
> phi[ 3 ]= [ 0.84459073]
> aa unsorted = [-3.780145]
> 
> after replacing scipy.linalg with numpy.linalg, same result
> 
> Got here
> Finished lstsq
> phi[ 0 ]= [[ 0.88271111+0.38811028j]]
> phi[ 1 ]= [-0.04649206-0.01274722j]
> phi[ 2 ]= 1
> phi[ 3 ]= [ 0.84459073]
> aa unsorted = [-3.780145]
> 
> Josef

Thanks for trying this!
-Bruce

> 
> >
> > For Linux, Ubuntu 10.04 x86-64
> > Python 2.6.5
> > numpy 1.30
> > scipy 0.7.0-2
> > ipython 0.10-1
> >
> > Thanks for any and all help!
> > -Bruce
> >
> >
> >



******************************
Neither the footer nor anything else in this E-mail is intended to or constitutes an <br>electronic signature and/or legally binding agreement in the absence of an <br>express statement or Autoliv policy and/or procedure to the contrary.<br>This E-mail and any attachments hereto are Autoliv property and may contain legally <br>privileged, confidential and/or proprietary information.<br>The recipient of this E-mail is prohibited from distributing, copying, forwarding or in any way <br>disseminating any material contained within this E-mail without prior written <br>permission from the author. If you receive this E-mail in error, please <br>immediately notify the author and delete this E-mail.  Autoliv disclaims all <br>responsibility and liability for the consequences of any person who fails to <br>abide by the terms herein. <br>
******************************