Re: [Numpy-discussion] NumPy/SciPy LAPACK version

On 7/16/07, Robert Kern <[EMAIL PROTECTED]> wrote: And we'd certainly appreciate the contribution. I'm tentatively going to say yes, we should start requiring LAPACK 3.0 unless if there is some very important platform that only comes with an older LAPACK. Great! The added benefit is that the

Re: [Numpy-discussion] NumPy/SciPy LAPACK version

Kevin Jacobs <[EMAIL PROTECTED]> wrote: > On 7/16/07, *Kevin Jacobs <[EMAIL PROTECTED] > >* <[EMAIL PROTECTED] > > wrote: > > This is a bit of a SciPy question, but I thought I'd ask here since > I'm already subscribed. I'd like to add s

Re: [Numpy-discussion] NumPy/SciPy LAPACK version

On 7/16/07, Charles R Harris <[EMAIL PROTECTED]> wrote: Hmm, I get a real result for this, although the result is wildly incorrect. Sqrtm isn't part of numpy, where are you getting it from? Mine is coming from pylab and looks remarkably buggy. from scipy.linalg import sqrtm I'm posting on

Re: [Numpy-discussion] NumPy/SciPy LAPACK version

On 7/16/07, Kevin Jacobs <[EMAIL PROTECTED]> <[EMAIL PROTECTED]> wrote: This is a bit of a SciPy question, but I thought I'd ask here since I'm already subscribed. I'd like to add some new LAPACK bindings to SciPy and was wondering if there was a minimum version requirement for LAPACK, since it

Re: [Numpy-discussion] NumPy/SciPy LAPACK version

On 7/16/07, Kevin Jacobs <[EMAIL PROTECTED]> <[EMAIL PROTECTED]> wrote: On 7/16/07, Charles R Harris <[EMAIL PROTECTED]> wrote: > > > > On 7/16/07, Robert Kern < [EMAIL PROTECTED]> wrote: > > > > Kevin Jacobs <[EMAIL PROTECTED]> wrote: > > > Mea culpa on the msqrt example, however I still think

Re: [Numpy-discussion] NumPy/SciPy LAPACK version

On 7/16/07, Charles R Harris <[EMAIL PROTECTED]> wrote: On 7/16/07, Robert Kern <[EMAIL PROTECTED]> wrote: > > Kevin Jacobs <[EMAIL PROTECTED]> wrote: > > Mea culpa on the msqrt example, however I still think it is wrong to > get > > a complex square-root back when a real valued result is expe

Re: [Numpy-discussion] NumPy/SciPy LAPACK version

On 7/16/07, Robert Kern <[EMAIL PROTECTED]> wrote: Kevin Jacobs <[EMAIL PROTECTED]> wrote: > Mea culpa on the msqrt example, however I still think it is wrong to get > a complex square-root back when a real valued result is expected and exists. No, in floating point you accumulate error. Those

Re: [Numpy-discussion] NumPy/SciPy LAPACK version

Kevin Jacobs <[EMAIL PROTECTED]> wrote: > Mea culpa on the msqrt example, however I still think it is wrong to get > a complex square-root back when a real valued result is expected and exists. No, in floating point you accumulate error. Those 1e-22j's are part of the actual result. Some systems l

Re: [Numpy-discussion] NumPy/SciPy LAPACK version

Mea culpa on the msqrt example, however I still think it is wrong to get a complex square-root back when a real valued result is expected and exists. -Kevin On 7/16/07, Hanno Klemm <[EMAIL PROTECTED]> wrote: Kevin, the problem appears to be that sqrtm() gives back an array, rather than a ma

Re: [Numpy-discussion] NumPy/SciPy LAPACK version

Kevin, the problem appears to be that sqrtm() gives back an array, rather than a matrix: >>> import scipy.linalg as sl >>> a = s.matrix([[59, 64, 69],[64, 72, 80],[69, 80, 91]]) >>> type(a) >>> a matrix([[59, 64, 69], [64, 72, 80], [69, 80, 91]]) >>> a*a - N.dot(a,a) matrix([[0

Re: [Numpy-discussion] NumPy/SciPy LAPACK version

Oups, sorry, I missed the 'a=matrix'... 2007/7/16, Matthieu Brucher <[EMAIL PROTECTED]>: Hi, Did you try numpy.dot(sqrtm(a), sqrtm(a)) ? Matthieu 2007/7/16, Kevin Jacobs <[EMAIL PROTECTED]> < [EMAIL PROTECTED]>: > > Hi all, > > This is a bit of a SciPy question, but I thought I'd ask here si

Re: [Numpy-discussion] NumPy/SciPy LAPACK version

Hi, Did you try numpy.dot(sqrtm(a), sqrtm(a)) ? Matthieu 2007/7/16, Kevin Jacobs <[EMAIL PROTECTED]> <[EMAIL PROTECTED]>: Hi all, This is a bit of a SciPy question, but I thought I'd ask here since I'm already subscribed. I'd like to add some new LAPACK bindings to SciPy and was wondering i

[Numpy-discussion] NumPy/SciPy LAPACK version

Hi all, This is a bit of a SciPy question, but I thought I'd ask here since I'm already subscribed. I'd like to add some new LAPACK bindings to SciPy and was wondering if there was a minimum version requirement for LAPACK, since it would be ideal if I could use some of the newer 3.0 features. I