[Numpy-discussion] Anyone with Core i7 and Ubuntu 10.04?

2010-11-08 Thread Ian Goodfellow
I'm wondering if anyone here has successfully built numpy with ATLAS and a Core i7 CPU on Ubuntu 10.04. If so, I could really use your help. I've been trying since August (see my earlier messages to this list) to get numpy running at full speed on my machine with no luck. The Ubuntu packages don't

Re: [Numpy-discussion] Trac unaware of github move

2010-11-08 Thread Sébastien Barthélemy
On Sun, 7 Nov 2010, Ralf Gommers wrote: That will require renaming those files in the source tree from *.txt to *.rst, otherwise there's no way to have github render them properly. Unless I missed something. Would that be fine? I think a *.rst.txt extension would also be recognized by github.

[Numpy-discussion] LapackError:non-native byte order

2010-11-08 Thread LittleBigBrain
Hi everyone, In my system '' is the native byte-order, but unless I change the byte-order label to '=', it won't work in linalg sub-module, but in others works OK. I am not sure whether this is an expected behavior or a bug? import sys sys.byteorder 'little' a.dtype.byteorder ''

Re: [Numpy-discussion] LapackError:non-native byte order

2010-11-08 Thread Pauli Virtanen
ma, 2010-11-08 kello 18:56 +0100, LittleBigBrain kirjoitti: In my system '' is the native byte-order, but unless I change the byte-order label to '=', it won't work in linalg sub-module, but in others works OK. I am not sure whether this is an expected behavior or a bug? import sys

Re: [Numpy-discussion] LapackError:non-native byte order

2010-11-08 Thread Pauli Virtanen
Mon, 08 Nov 2010 19:31:31 +0100, Pauli Virtanen wrote: ma, 2010-11-08 kello 18:56 +0100, LittleBigBrain kirjoitti: In my system '' is the native byte-order, but unless I change the byte-order label to '=', it won't work in linalg sub-module, but in others works OK. I am not sure whether this

[Numpy-discussion] OFFTOPIC: simple databases

2010-11-08 Thread Renato Fabbri
Dear All, i want to find simple databases, like a 5 dimensional with more than 30 samples. i am having difficult times with this. where do you get them? all the best, rf -- GNU/Linux User #479299 skype: fabbri.renato ___ NumPy-Discussion mailing

Re: [Numpy-discussion] LapackError:non-native byte order

2010-11-08 Thread Matthew Brett
Hi, On Mon, Nov 8, 2010 at 10:34 AM, Pauli Virtanen p...@iki.fi wrote: Mon, 08 Nov 2010 19:31:31 +0100, Pauli Virtanen wrote: ma, 2010-11-08 kello 18:56 +0100, LittleBigBrain kirjoitti: In my system '' is the native byte-order, but unless I change the byte-order label to '=', it won't work

Re: [Numpy-discussion] OFFTOPIC: simple databases

2010-11-08 Thread Pauli Virtanen
Mon, 08 Nov 2010 17:00:34 -0200, Renato Fabbri wrote: [clip: offtopic] Please post this on the scipy-user list instead, it's more suitable for misc questions. -- Pauli Virtanen ___ NumPy-Discussion mailing list NumPy-Discussion@scipy.org

[Numpy-discussion] Solving Ax = b: inverse vs cholesky factorization

2010-11-08 Thread Joon
Hi,I was wondering when it is better to store cholesky factor and use it to solve Ax = b, instead of storing the inverse of A. (A is a symmetric, positive-definite matrix.)Even in the repeated case, if I have the inverse of A (invA) stored, then I can solve Ax = b_i, i = 1, ... , n, by x =

[Numpy-discussion] Solving Ax = b: inverse vs cholesky factorization

2010-11-08 Thread Joon
Hi,I was wondering when it is better to store cholesky factor and use it to solve Ax = b, instead of storing the inverse of A. (A is a symmetric, positive-definite matrix.)Even in the repeated case, if I have the inverse of A (invA) stored, then I can solve Ax = b_i, i = 1, ... , n, by x =

Re: [Numpy-discussion] Solving Ax = b: inverse vs cholesky factorization

2010-11-08 Thread Pauli Virtanen
Mon, 08 Nov 2010 13:17:11 -0600, Joon wrote: I was wondering when it is better to store cholesky factor and use it to solve Ax = b, instead of storing the inverse of A. (A is a symmetric, positive-definite matrix.) Even in the repeated case, if I have the inverse of A (invA) stored, then I

[Numpy-discussion] Developmental version numbering with git

2010-11-08 Thread Bruce Southey
Hi, Since the change to git the numpy version in setup.py is '2.0.0.dev' regardless because the prior numbering was determined by svn. Is there a plan to add some numbering system to numpy developmental version? Regardless of the answer, the 'numpy/numpy/version.py' will need to changed

Re: [Numpy-discussion] Solving Ax = b: inverse vs cholesky factorization

2010-11-08 Thread Joon
On Mon, 08 Nov 2010 13:23:46 -0600, Pauli Virtanen p...@iki.fi wrote: Mon, 08 Nov 2010 13:17:11 -0600, Joon wrote: I was wondering when it is better to store cholesky factor and use it to solve Ax = b, instead of storing the inverse of A. (A is a symmetric, positive-definite matrix.) Even in

Re: [Numpy-discussion] Solving Ax = b: inverse vs cholesky factorization

2010-11-08 Thread Joon
On Mon, 08 Nov 2010 13:23:46 -0600, Pauli Virtanen p...@iki.fi wrote: Mon, 08 Nov 2010 13:17:11 -0600, Joon wrote: I was wondering when it is better to store cholesky factor and use it to solve Ax = b, instead of storing the inverse of A. (A is a symmetric, positive-definite matrix.) Even

Re: [Numpy-discussion] Solving Ax = b: inverse vs cholesky factorization

2010-11-08 Thread Bruce Southey
On 11/08/2010 01:38 PM, Joon wrote: On Mon, 08 Nov 2010 13:23:46 -0600, Pauli Virtanen p...@iki.fi wrote: Mon, 08 Nov 2010 13:17:11 -0600, Joon wrote: I was wondering when it is better to store cholesky factor and use it to solve Ax = b, instead of storing the inverse of A. (A is a

[Numpy-discussion] Catching and dealing with floating point errors

2010-11-08 Thread Skipper Seabold
I am doing some optimizations on random samples. In a small number of cases, the objective is not well-defined for a given sample (it's not possible to tell beforehand and hopefully won't happen much in practice). What is the most numpythonic way to handle this? It doesn't look like I can use

Re: [Numpy-discussion] Solving Ax = b: inverse vs cholesky factorization

2010-11-08 Thread Pauli Virtanen
On Mon, 08 Nov 2010 14:06:03 -0600, Bruce Southey wrote: [clip] Numpy uses SVD to get the (pseudo) inverse, which is usually very accurate at getting (pseudo) inverse. numpy.linalg.inv does solve(a, identity(a.shape[0], dtype=a.dtype)) It doesn't use xGETRI since that's not included

Re: [Numpy-discussion] Catching and dealing with floating point errors

2010-11-08 Thread Skipper Seabold
On Mon, Nov 8, 2010 at 3:14 PM, Skipper Seabold jsseab...@gmail.com wrote: I am doing some optimizations on random samples.  In a small number of cases, the objective is not well-defined for a given sample (it's not possible to tell beforehand and hopefully won't happen much in practice).  

Re: [Numpy-discussion] Catching and dealing with floating point errors

2010-11-08 Thread Bruce Southey
On 11/08/2010 02:17 PM, Skipper Seabold wrote: On Mon, Nov 8, 2010 at 3:14 PM, Skipper Seaboldjsseab...@gmail.com wrote: I am doing some optimizations on random samples. In a small number of cases, the objective is not well-defined for a given sample (it's not possible to tell beforehand and

Re: [Numpy-discussion] Catching and dealing with floating point errors

2010-11-08 Thread Warren Weckesser
On Mon, Nov 8, 2010 at 2:17 PM, Skipper Seabold jsseab...@gmail.com wrote: On Mon, Nov 8, 2010 at 3:14 PM, Skipper Seabold jsseab...@gmail.com wrote: I am doing some optimizations on random samples. In a small number of cases, the objective is not well-defined for a given sample (it's not

Re: [Numpy-discussion] Catching and dealing with floating point errors

2010-11-08 Thread Skipper Seabold
On Mon, Nov 8, 2010 at 3:42 PM, Bruce Southey bsout...@gmail.com wrote: On 11/08/2010 02:17 PM, Skipper Seabold wrote: On Mon, Nov 8, 2010 at 3:14 PM, Skipper Seaboldjsseab...@gmail.com  wrote: I am doing some optimizations on random samples.  In a small number of cases, the objective is not

Re: [Numpy-discussion] Catching and dealing with floating point errors

2010-11-08 Thread Skipper Seabold
On Mon, Nov 8, 2010 at 3:45 PM, Warren Weckesser warren.weckes...@enthought.com wrote: On Mon, Nov 8, 2010 at 2:17 PM, Skipper Seabold jsseab...@gmail.com wrote: On Mon, Nov 8, 2010 at 3:14 PM, Skipper Seabold jsseab...@gmail.com wrote: I am doing some optimizations on random samples.  In

Re: [Numpy-discussion] numpy.genfromtxt converters issue

2010-11-08 Thread Damien Moore
Pierre GM pgmdevlist at gmail.com writes: On Nov 6, 2010, at 2:22 PM, Damien Moore wrote: Hi List, I'm trying to import csv data as a numpy array using genfromtxt. [...] Please open a ticket so that I don't forget about it. Thx in advance! The ticket is here:

Re: [Numpy-discussion] Catching and dealing with floating point errors

2010-11-08 Thread Bruce Southey
On 11/08/2010 02:52 PM, Skipper Seabold wrote: On Mon, Nov 8, 2010 at 3:42 PM, Bruce Southeybsout...@gmail.com wrote: On 11/08/2010 02:17 PM, Skipper Seabold wrote: On Mon, Nov 8, 2010 at 3:14 PM, Skipper Seaboldjsseab...@gmail.com wrote: I am doing some optimizations on random samples.

Re: [Numpy-discussion] Catching and dealing with floating point errors

2010-11-08 Thread Skipper Seabold
On Mon, Nov 8, 2010 at 4:04 PM, Bruce Southey bsout...@gmail.com wrote: On 11/08/2010 02:52 PM, Skipper Seabold wrote: On Mon, Nov 8, 2010 at 3:42 PM, Bruce Southeybsout...@gmail.com  wrote: On 11/08/2010 02:17 PM, Skipper Seabold wrote: On Mon, Nov 8, 2010 at 3:14 PM, Skipper

Re: [Numpy-discussion] Developmental version numbering with git

2010-11-08 Thread Matthew Brett
Hi, Since the change to git the numpy version in setup.py is '2.0.0.dev' regardless because the prior numbering was determined by svn. Is there a plan to add some numbering system to numpy developmental version? Regardless of the answer, the 'numpy/numpy/version.py' will need to changed

Re: [Numpy-discussion] Catching and dealing with floating point errors

2010-11-08 Thread Warren Weckesser
On Mon, Nov 8, 2010 at 2:52 PM, Skipper Seabold jsseab...@gmail.com wrote: On Mon, Nov 8, 2010 at 3:45 PM, Warren Weckesser warren.weckes...@enthought.com wrote: On Mon, Nov 8, 2010 at 2:17 PM, Skipper Seabold jsseab...@gmail.com wrote: On Mon, Nov 8, 2010 at 3:14 PM, Skipper Seabold

Re: [Numpy-discussion] Solving Ax = b: inverse vs cholesky factorization

2010-11-08 Thread Nathaniel Smith
On Mon, Nov 8, 2010 at 12:00 PM, Joon groups.and.li...@gmail.com wrote: Another question is, is it better to do cho_solve(cho_factor(A), b) than solve(A, b)? If A is symmetric positive definite, then using the cholesky decomposition should be somewhat faster than using a more general solver.

Re: [Numpy-discussion] Solving Ax = b: inverse vs cholesky factorization

2010-11-08 Thread Anne Archibald
On 8 November 2010 14:38, Joon groups.and.li...@gmail.com wrote: Oh I see. So I guess in invA = solve(Ax, I) and then x = dot(invA, b) case, there are more places where numerical errors occur, than just x = solve(Ax, b) case. That's the heart of the matter, but one can be more specific. You

Re: [Numpy-discussion] updated 1.5.1 release schedule

2010-11-08 Thread Russell E. Owen
In article aanlktimfgckbg8cprygukcvwvqzxqycykgexvx_=8...@mail.gmail.com, Ralf Gommers ralf.gomm...@googlemail.com wrote: On Mon, Nov 8, 2010 at 5:16 AM, Vincent Davis vinc...@vincentdavis.net wrote: On Sun, Nov 7, 2010 at 1:51 AM, Ralf Gommers ralf.gomm...@googlemail.com wrote:

Re: [Numpy-discussion] Solving Ax = b: inverse vs cholesky factorization

2010-11-08 Thread Joon
Thanks, Nathaniel. Your reply was very helpful. -Joon On Mon, 08 Nov 2010 15:47:22 -0600, Nathaniel Smith n...@pobox.com wrote: On Mon, Nov 8, 2010 at 12:00 PM, Joon groups.and.li...@gmail.com wrote: Another question is, is it better to do cho_solve(cho_factor(A), b) than solve(A, b)? If A

Re: [Numpy-discussion] LapackError:non-native byte order

2010-11-08 Thread braingateway
Matthew Brett : Hi, On Mon, Nov 8, 2010 at 10:34 AM, Pauli Virtanen p...@iki.fi wrote: Mon, 08 Nov 2010 19:31:31 +0100, Pauli Virtanen wrote: ma, 2010-11-08 kello 18:56 +0100, LittleBigBrain kirjoitti: In my system '' is the native byte-order, but unless I change the

Re: [Numpy-discussion] Anyone with Core i7 and Ubuntu 10.04?

2010-11-08 Thread David
Hi Ian, On 11/08/2010 11:18 PM, Ian Goodfellow wrote: I'm wondering if anyone here has successfully built numpy with ATLAS and a Core i7 CPU on Ubuntu 10.04. If so, I could really use your help. I've been trying since August (see my earlier messages to this list) to get numpy running at full

Re: [Numpy-discussion] Anyone with Core i7 and Ubuntu 10.04?

2010-11-08 Thread David Warde-Farley
On 2010-11-08, at 8:52 PM, David wrote: Please tell us what error you got - saying that something did not working is really not useful to help you. You need to say exactly what fails, and which steps you followed before that failure. I think what he means is that it's very slow, there's no

Re: [Numpy-discussion] Anyone with Core i7 and Ubuntu 10.04?

2010-11-08 Thread Wes McKinney
On Mon, Nov 8, 2010 at 11:33 PM, David Warde-Farley warde...@iro.umontreal.ca wrote: On 2010-11-08, at 8:52 PM, David wrote: Please tell us what error you got - saying that something did not working is really not useful to help you. You need to say exactly what fails, and which steps you

Re: [Numpy-discussion] Catching and dealing with floating point errors

2010-11-08 Thread Warren Weckesser
On Mon, Nov 8, 2010 at 3:20 PM, Warren Weckesser warren.weckes...@enthought.com wrote: On Mon, Nov 8, 2010 at 2:52 PM, Skipper Seabold jsseab...@gmail.comwrote: On Mon, Nov 8, 2010 at 3:45 PM, Warren Weckesser warren.weckes...@enthought.com wrote: On Mon, Nov 8, 2010 at 2:17 PM,