Re: [Numpy-discussion] SVD errors

2013-06-04 Thread Amir Mohammadi
Robert Kern robert.kern at gmail.com writes: On Mon, Feb 9, 2009 at 16:25, M Trumpis mtrumpis at berkeley.edu wrote: I played around with a C translation of that test program, and found that dgesvd (but not dgesdd) happens to converge and return all non-negative singular values for both

Re: [Numpy-discussion] SVD errors

2013-06-04 Thread Robert Kern
On Tue, Jun 4, 2013 at 5:18 PM, Amir Mohammadi 183.a...@gmail.com wrote: Hello Everyone, I am trying to write my own svd function to use the dgesvd method from lapack, but my problem is that I cannot find the dgesvd method from this import from numpy.linalg import lapack_lite. Do I have to

Re: [Numpy-discussion] numpy tests errors and failures

2013-06-04 Thread Warren Weckesser
On Tue, Jun 4, 2013 at 7:52 AM, Warren Weckesser warren.weckes...@gmail.com wrote: On Tue, Jun 4, 2013 at 1:20 AM, Tim Burgess tim.burg...@noaa.gov wrote: On Sat, 2013-06-01 at 20:09 -0400, Warren Weckesser wrote: I'm using Ubuntu 12.04, so I suspect I won't be the only one who sees

[Numpy-discussion] suggested change of behavior for interp

2013-06-04 Thread Slavin, Jonathan
Hi, I would like to suggest that the behavior of numpy.interp be changed regarding treatment of situations in which the x-coordinates are not monotonically increasing. Specifically, it seems to me that interp should work correctly when the x-coordinate is decreasing monotonically. Clearly it

[Numpy-discussion] request for SWIG numpy.i users

2013-06-04 Thread Ralf Gommers
Hi, If you're using or are very familiar with SWIG and the numpy.i interface to it, please help to test and/or review https://github.com/numpy/numpy/pull/3148. It's a fairly major update to numpy.i by Bill Spotz, containing the following: - support for 4D arrays and memory managed output

Re: [Numpy-discussion] suggested change of behavior for interp

2013-06-04 Thread Charles R Harris
On Tue, Jun 4, 2013 at 12:07 PM, Slavin, Jonathan jsla...@cfa.harvard.eduwrote: Hi, I would like to suggest that the behavior of numpy.interp be changed regarding treatment of situations in which the x-coordinates are not monotonically increasing. Specifically, it seems to me that interp

Re: [Numpy-discussion] suggested change of behavior for interp

2013-06-04 Thread Eric Firing
On 2013/06/04 2:05 PM, Charles R Harris wrote: On Tue, Jun 4, 2013 at 12:07 PM, Slavin, Jonathan jsla...@cfa.harvard.edu mailto:jsla...@cfa.harvard.edu wrote: Hi, I would like to suggest that the behavior of numpy.interp be changed regarding treatment of situations in which

Re: [Numpy-discussion] suggested change of behavior for interp

2013-06-04 Thread Benjamin Root
Could non-monotonicity be detected as part of the interp process? Perhaps a sign switch in the deltas? I have been bitten by this problem too. Cheers! Ben Root On Jun 4, 2013 9:08 PM, Eric Firing efir...@hawaii.edu wrote: On 2013/06/04 2:05 PM, Charles R Harris wrote: On Tue, Jun 4,

Re: [Numpy-discussion] suggested change of behavior for interp

2013-06-04 Thread Eric Firing
On 2013/06/04 4:15 PM, Benjamin Root wrote: Could non-monotonicity be detected as part of the interp process? Perhaps a sign switch in the deltas? There are two code paths, depending on the number of points to be interpolated. When it is greater than the size of the table, the deltas are