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] [AstroPy] Porting "IDL Astronomy User's Library" to numpy

On Jul 16, 2007, at 1:19 PM, W.T. Bridgman wrote: > Perry, > > I believe some of those documents are getting a bit dated. They > still refer to only supporting numarray vs Numeric. Don't those need > to be updated to specify numpy? > Yes, that's certainly true. Having said that, it's probably g

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] [AstroPy] Porting "IDL Astronomy User's Library" to numpy

On Jul 16, 2007, at 9:13 AM, Yannick Copin wrote: > Hi, > > I'd be interested in some astronomical utilities from the IDL > Astronomy User's > Library (http://idlastro.gsfc.nasa.gov/contents.html) converted to > python/numpy. I had a look to idl2python > (http://software.pseudogreen.org/i2py/),

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] Porting "IDL Astronomy User's Library" to numpy

On 7/16/07, Yannick Copin <[EMAIL PROTECTED]> wrote: Hi, I'd be interested in some astronomical utilities from the IDL Astronomy User's Library (http://idlastro.gsfc.nasa.gov/contents.html) converted to python/numpy. I had a look to idl2python (http://software.pseudogreen.org/i2py/), but the au

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

[Numpy-discussion] Porting "IDL Astronomy User's Library" to numpy

Hi, I'd be interested in some astronomical utilities from the IDL Astronomy User's Library (http://idlastro.gsfc.nasa.gov/contents.html) converted to python/numpy. I had a look to idl2python (http://software.pseudogreen.org/i2py/), but the automatic translation fails, mostly because (I think)

Re: [Numpy-discussion] Problems building numpy

As for my unrelated question, I was still wondering if anyone has any information about the relative merits of MKL vs ATLAS etc. MKL is parallized, so until ATLAS is as well, MKL has the upper hand. Matthieu ___ Numpy-discussion mailing list Numpy-d

Re: [Numpy-discussion] Problems building numpy

Hi, Thanks for the responses. On Fri, 13 Jul 2007 15:37:52 -0700, "Christopher Barker" <[EMAIL PROTECTED]> said: > CygWin is NOT the same as MinGW. They both are gcc, but IIUC, cygwin > builds against the cygwin unix-like standard libs, and MinGW builds > against the system libs -- so only MinGW