Re: [Numpy-discussion] 673

2008-04-28 Thread Pearu Peterson
Hi, As far as I am concerned, the issue needs a cosmetic fix of renaming pythonxerbla to python_xerbla and the rest of the issue can be postponed to 1.2. Note that this isn't purely a numpy issue. To fix the issue, system or user provided blas/lapack libraries need to be changed, we can only

[Numpy-discussion] Linear Interpolation Question

2008-04-28 Thread Andrea Gavana
Hi All, I have 2 matrices coming from 2 different simulations: the first column of the matrices is a date (time) at which all the other results in the matrix have been reported (simulation step). In these 2 matrices, very often the simulation steps do not coincide, so I just want to

[Numpy-discussion] Starting to work on runtime plugin system for plugin (automatic sse optimization, etc...)

2008-04-28 Thread David Cournapeau
Hi, I've just started working on a prototype for a plugin system for numpy. The plugin aims at providing a framework for the following user cases: - runtime selection of blas/lapack/etc...: instead of harcoding in the binary one blas/lapack implementation, numpy could choose the SSE

Re: [Numpy-discussion] Distribution Functions Change Behavior

2008-04-28 Thread Hoyt Koepke
I may not understand what you are asking, Rich, but I'm not sure I agree with Alan. A Gaussian fit to data x should fit exactly as well as data fit to ax, a 0, just with a variance a^2 times the original. The only way this would not be true is if: 1. You are not fitting the variance, but only

Re: [Numpy-discussion] Distribution Functions Change Behavior

2008-04-28 Thread Hoyt Koepke
Wait, I think I see what Alan is saying. When you use a gaussian approximation on truncated data, the accuracy of the truncation is very dependent on where in the interval the mean is. If it's near the edges, the results will be worse. The width of the interval, though, is a separate factor.

Re: [Numpy-discussion] leaving matrix indexing discussion

2008-04-28 Thread Christopher Barker
Alan G Isaac wrote: I am starting to feel that I have abused my interlocutors. Not at all -- you have contributed a great deal to the conversation, and your maintenance of: URL:http://www.scipy.org/MatrixIndexing Is particularly valuable. Thank you, -Chris -- Christopher Barker, Ph.D.

Re: [Numpy-discussion] untenable matrix behavior in SVN

2008-04-28 Thread Christopher Barker
oops, typo! Christopher Barker wrote: Gael Varoquaux wrote: I don't know why people are indexing matrices with A[x][y], but they shouldn't. I think there has been a misunderstanding here. I don't think anyone is suggesting that if a coder wants an element of a matrix, that s/he should

Re: [Numpy-discussion] prerelease proposal for matrix behavior

2008-04-28 Thread Alan G Isaac
On Mon, 28 Apr 2008, Timothy Hochberg apparently wrote: Can you clarify what you mean by submatrix extraction? It sounds like you want to be able index into an MxN array and get out a 1xN or Mx1 matrix. If that's the case, wouldn't the natural way to spell that under the

Re: [Numpy-discussion] Matrix Class [was numpy release]

2008-04-28 Thread Bill Spotz
On Apr 24, 2008, at 8:52 PM, Bill Spotz wrote: On Apr 24, 2008, at 5:45 PM, Timothy Hochberg wrote: Bill Spotz wrote: I have generally thought about this in the context of, say, a Krylov-space iterative method, and what that type of interface would lead to the most readable code. Can