Re: [Numpy-discussion] SVD does not converge on "clean" matrix

2011-08-20 Thread Gael Varoquaux
On Sun, Aug 14, 2011 at 09:15:35PM +0200, Charanpal Dhanjal wrote: > Incidentally, I am confused as to why numpy calls the lapack lite > routines - when I call numpy.show_config() it seems to have detected my > ATLAS libraries and I would have expected it to use those. My rule of thumb is to nev

Re: [Numpy-discussion] SVD does not converge on "clean" matrix

2011-08-14 Thread Charanpal Dhanjal
Thanks very much Lou for the information. I tried delving into the C code and found a line in the dlasd4_ routine which reads: for (niter = iter; niter <= MAXITERLOOPS; ++niter) { This is apparently the main loop for this subroutine and the value of MAXITERLOOPS = 100. All I did was increase th

Re: [Numpy-discussion] SVD does not converge on "clean" matrix

2011-08-14 Thread Lou Pecora
Chuck wrote: Fails here also, fedora 15 64 bits AMD 940. There should be a maximum iterations argument somewhere... Chuck ---    ***  Here's the "FIX": Chuck is right.  There is a max iterations.  Here is a re

Re: [Numpy-discussion] SVD does not converge on "clean" matrix

2011-08-14 Thread Charanpal Dhanjal
I had a quick look at the code (https://github.com/numpy/numpy/blob/master/numpy/linalg/linalg.py) and the numpy.linalg.svd function calls lapack_lite.dgesdd (for real matrices) so I guess the non-convergence occurs in this function. As I understood lapack_lite is used by default unless numpy i

Re: [Numpy-discussion] SVD does not converge on "clean" matrix

2011-08-13 Thread Charles R Harris
On Thu, Aug 11, 2011 at 7:23 AM, wrote: > Hi all, > > I get an error message "numpy.linalg.linalg.LinAlgError: SVD did not > converge" when calling numpy.linalg.svd on a "clean" matrix of size (1952, > 895). The matrix is clean in the sense that it contains no NaN or Inf > values. The correspondi

Re: [Numpy-discussion] SVD does not converge on "clean" matrix

2011-08-12 Thread Nadav Horesh
: 12 August 2011 16:33 To: Discussion of Numerical Python Subject: Re: [Numpy-discussion] SVD does not converge on "clean" matrix On Fri, Aug 12, 2011 at 4:03 AM, Charanpal Dhanjal mailto:dhan...@telecom-paristech.fr>> wrote: Thank Nadav for testing out the matrix. I wonder if you

Re: [Numpy-discussion] SVD does not converge on "clean" matrix

2011-08-12 Thread Warren Weckesser
Aug 2011 07:21:09 -0700 > From: Nadav Horesh > Subject: Re: [Numpy-discussion] SVD does not converge on "clean" > matrix > To: Discussion of Numerical Python > Message-ID: > > <26FC23E7C398A64083C980D16001012D246DFC5F90@VA3DIAXVS361.RED001.local> > Conte

Re: [Numpy-discussion] SVD does not converge on "clean" matrix

2011-08-12 Thread Charanpal Dhanjal
umpy 1.5.1 (without ATLAS), and so the reasons for the differences seem to be OS or processor dependent. Any ideas? Charanpal Date: Thu, 11 Aug 2011 07:21:09 -0700 From: Nadav Horesh Subject: Re: [Numpy-discussion] SVD does not converge on "clean" matrix To: Discussion of Numerica

Re: [Numpy-discussion] SVD does not converge on "clean" matrix

2011-08-11 Thread Nadav Horesh
...@telecom-paristech.fr] Sent: 11 August 2011 16:23 To: numpy-discussion@scipy.org Subject: [Numpy-discussion] SVD does not converge on "clean" matrix Hi all, I get an error message "numpy.linalg.linalg.LinAlgError: SVD did not converge" when calling numpy.linalg.svd on a "cl

[Numpy-discussion] SVD does not converge on "clean" matrix

2011-08-11 Thread dhanjal
Hi all, I get an error message "numpy.linalg.linalg.LinAlgError: SVD did not converge" when calling numpy.linalg.svd on a "clean" matrix of size (1952, 895). The matrix is clean in the sense that it contains no NaN or Inf values. The corresponding npz file is available here: https://docs.google.co

Re: [Numpy-discussion] SVD does not converge

2011-06-28 Thread santhu kumar
Pecora) > > > ------ > > Message: 1 > Date: Tue, 28 Jun 2011 10:56:48 -0500 > From: santhu kumar > Subject: [Numpy-discussion] SVD does not converge > To: numpy-discussion@scipy.org > Message-ID: > Co

Re: [Numpy-discussion] SVD does not converge

2011-06-28 Thread Charles R Harris
On Tue, Jun 28, 2011 at 11:36 AM, eat wrote: > Hi, > > On Tue, Jun 28, 2011 at 7:43 PM, Lou Pecora wrote: > >> >> -- >> *From:* santhu kumar >> *To:* numpy-discussion@scipy.org >> *Sent:* Tue, June 28, 2011 11:56:48 AM >

Re: [Numpy-discussion] SVD does not converge

2011-06-28 Thread eat
Hi, On Tue, Jun 28, 2011 at 7:43 PM, Lou Pecora wrote: > > -- > *From:* santhu kumar > *To:* numpy-discussion@scipy.org > *Sent:* Tue, June 28, 2011 11:56:48 AM > *Subject:* [Numpy-discussion] SVD does not converge > > Hello, > > I hav

Re: [Numpy-discussion] SVD does not converge

2011-06-28 Thread Lou Pecora
From: santhu kumar To: numpy-discussion@scipy.org Sent: Tue, June 28, 2011 11:56:48 AM Subject: [Numpy-discussion] SVD does not converge Hello, I have a 380X5 matrix and when I am calculating pseudo-inverse of the matrix using pinv(numpy.linalg) I get the

Re: [Numpy-discussion] SVD does not converge

2011-06-28 Thread Sebastian Berg
Hi, On Tue, 2011-06-28 at 10:56 -0500, santhu kumar wrote: > Hello, > > I have a 380X5 matrix and when I am calculating pseudo-inverse of the > matrix using pinv(numpy.linalg) I get the following error message: > > raise LinAlgError, 'SVD did not converge' > numpy.linalg.linalg.LinAlgError: SVD

Re: [Numpy-discussion] SVD does not converge

2011-06-28 Thread Charles R Harris
On Tue, Jun 28, 2011 at 9:56 AM, santhu kumar wrote: > Hello, > > I have a 380X5 matrix and when I am calculating pseudo-inverse of the > matrix using pinv(numpy.linalg) I get the following error message: > > raise LinAlgError, 'SVD did not converge' > numpy.linalg.linalg.LinAlgError: SVD did not

[Numpy-discussion] SVD does not converge

2011-06-28 Thread santhu kumar
Hello, I have a 380X5 matrix and when I am calculating pseudo-inverse of the matrix using pinv(numpy.linalg) I get the following error message: raise LinAlgError, 'SVD did not converge' numpy.linalg.linalg.LinAlgError: SVD did not converge I have looked in the list that it is a recurring issue b