Re: [Numpy-discussion] (no subject)

2006-06-27 Thread kanako
:―― INFORMATION ―: 不正・悪質なサイトを一切排除しておりますので 安心してご利用ください。 http://love-match.bz/pc/?03 :――: *・゜゜・*:.。. .。.:*・゜゜・*:.。..。:*・゜゜・*:.。..。:**・゜゜・* ━━━  □■ 不倫・ワリキリ専門の無料出会いサイト『Love☆Match』  --

Re: [Numpy-discussion] (no subject)

2006-06-27 Thread emi
:―― INFORMATION ―: 不正・悪質なサイトを一切排除しておりますので 安心してご利用ください。 http://love-match.bz/pc/?010 :――: *・゜゜・*:.。. .。.:*・゜゜・*:.。..。:*・゜゜・*:.。..。:**・゜゜・* お金と時間を持て余している人妻の間で、噂になってるあのサイト   [登録・利用料全て無料] http://love-match.bz/pc/?010 ━━━

Re: [Numpy-discussion] (no subject)

2006-06-27 Thread miku
:―― INFORMATION ―: 不正・悪質なサイトを一切排除しておりますので 安心してご利用ください。 http://love-match.bz/pc/?07 :――: *・゜゜・*:.。. .。.:*・゜゜・*:.。..。:*・゜゜・*:.。..。:**・゜゜・* ━━━  □■ 不倫・ワリキリ専門の無料出会いサイト『Love☆Match』  --

Re: [Numpy-discussion] Should cholesky return upper or lower triangular matrix?

2006-06-27 Thread Keith Goodman
On 6/27/06, Robert Kern <[EMAIL PROTECTED]> wrote: > Keith Goodman wrote: > > Isn't the Cholesky decomposition by convention an upper triangular > > matrix? I noticed, by porting Octave code, that linalg.cholesky > > returns the lower triangular matrix. > > > > References: > > > > http://mathworld.

Re: [Numpy-discussion] Should cholesky return upper or lower triangular matrix?

2006-06-27 Thread Robert Kern
Keith Goodman wrote: > Isn't the Cholesky decomposition by convention an upper triangular > matrix? I noticed, by porting Octave code, that linalg.cholesky > returns the lower triangular matrix. > > References: > > http://mathworld.wolfram.com/CholeskyDecomposition.html > http://www.mathworks.com

[Numpy-discussion] Should cholesky return upper or lower triangular matrix?

2006-06-27 Thread Keith Goodman
Isn't the Cholesky decomposition by convention an upper triangular matrix? I noticed, by porting Octave code, that linalg.cholesky returns the lower triangular matrix. References: http://mathworld.wolfram.com/CholeskyDecomposition.html http://www.mathworks.com/access/helpdesk/help/techdoc/ref/cho

Re: [Numpy-discussion] Numpy Benchmarking

2006-06-27 Thread Robert Kern
Keith Goodman wrote: > Scipy computes the inverse of a matrix faster than numpy (except if > the dimensions of x are small). But scipy is slower than numpy for > eigh (I only checked for symmetric positive definite matrices): Looks like scipy uses *SYEV and numpy uses the better *SYEVD (the D stan

Re: [Numpy-discussion] Numpy Benchmarking

2006-06-27 Thread Keith Goodman
On 6/27/06, Keith Goodman <[EMAIL PROTECTED]> wrote: > On 6/27/06, Travis Oliphant <[EMAIL PROTECTED]> wrote: > > > The numpy.dual library exists so you can use the SciPy calls if the > > person has SciPy installed or the NumPy ones otherwise. It exists > > precisely for the purpose of seamlessly

Re: [Numpy-discussion] Numpy Benchmarking

2006-06-27 Thread Keith Goodman
On 6/27/06, Travis Oliphant <[EMAIL PROTECTED]> wrote: > The numpy.dual library exists so you can use the SciPy calls if the > person has SciPy installed or the NumPy ones otherwise. It exists > precisely for the purpose of seamlessly taking advantage of > algorithms/interfaces that exist in NumP

Re: [Numpy-discussion] Numpy Benchmarking

2006-06-27 Thread Travis Oliphant
Keith Goodman wrote: > On 6/27/06, Dennis V. Perepelitsa <[EMAIL PROTECTED]> wrote: > > >> I've run some benchmarks comparing the performance of scipy, numpy, >> Numeric and numarray vs. MATLAB. >> > > I enjoyed looking at the results. > > The most interesting result, for me, was that inver

Re: [Numpy-discussion] Numpy Benchmarking

2006-06-27 Thread Travis Oliphant
Dennis V. Perepelitsa wrote: > Hi, all. > > I've run some benchmarks comparing the performance of scipy, numpy, > Numeric and numarray vs. MATLAB. There's also the beginnings of a > benchmark framework included. The results are online at: > > http://web.mit.edu/jonas/www/bench/ > > They were

Re: [Numpy-discussion] Numpy Benchmarking

2006-06-27 Thread Keith Goodman
On 6/27/06, Dennis V. Perepelitsa <[EMAIL PROTECTED]> wrote: > I've run some benchmarks comparing the performance of scipy, numpy, > Numeric and numarray vs. MATLAB. I enjoyed looking at the results. The most interesting result, for me, was that inverting a matrix is much faster in scipy than nu

Re: [Numpy-discussion] Numpy Benchmarking

2006-06-27 Thread Robert Kern
Dennis V. Perepelitsa wrote: > Hi, all. > > I've run some benchmarks comparing the performance of scipy, numpy, > Numeric and numarray vs. MATLAB. There's also the beginnings of a > benchmark framework included. The results are online at: > > http://web.mit.edu/jonas/www/bench/ > > They w

[Numpy-discussion] Numpy Benchmarking

2006-06-27 Thread Dennis V. Perepelitsa
Hi, all. I've run some benchmarks comparing the performance of scipy, numpy, Numeric and numarray vs. MATLAB. There's also the beginnings of a benchmark framework included. The results are online at: http://web.mit.edu/jonas/www/bench/ They were produced on a Thinkpad T42 with an Intel Pe

Re: [Numpy-discussion] int64 wierdness

2006-06-27 Thread Andrew Straw
An SVN checkout from a week or two ago looks OK on my amd64 machine: [EMAIL PROTECTED]:~$ python Python 2.4.3 (#2, Apr 27 2006, 14:43:32) [GCC 4.0.3 (Ubuntu 4.0.3-1ubuntu5)] on linux2 Type "help", "copyright", "credits" or "license" for more information. >>> import numpy >>> numpy.__version__ '0

Re: [Numpy-discussion] int64 wierdness

2006-06-27 Thread EI
numpy.__version__ says 0.9.8. Python 2.4.2, GCC 4.1, OpenSuSE 10.1 (x86_64).Thanks Travis,EugeneOn 6/27/06, Travis Oliphant < [EMAIL PROTECTED]> wrote:EI wrote:> Hi,>> I'm running python 2.4 on a 64bit linux and get strange results:> (int(9))**2 is equal to 81, as it should, but> (int64(9))**2 is

[Numpy-discussion] numexpr does sum.

2006-06-27 Thread Tim Hochberg
I managed to get basic support for sum and prod into numpexpr. I need to tie up some loose ends, for instance only floats are currently supported, but these should be easy. To return to the recently posted multidimensional distance program, this now works: expr = numexpr("sum((a - b)**2, a

Re: [Numpy-discussion] int64 wierdness

2006-06-27 Thread Travis Oliphant
EI wrote: > Hi, > > I'm running python 2.4 on a 64bit linux and get strange results: > (int(9))**2 is equal to 81, as it should, but > (int64(9))**2 is equal to 0 Thanks for the bug-report. Please provide the version of NumPy you are using so we can track it down, or suggest an upgrade. -Travis

[Numpy-discussion] int64 wierdness

2006-06-27 Thread EI
Hi,I'm running python 2.4 on a 64bit linux and get strange results:(int(9))**2 is equal to 81, as it should, but(int64(9))**2 is equal to 0Is it a bug or a feature?Eugene Using Tomcat but need to do more? Need to support web services, security? Get stuff done quickly with pre-integrated technology

Re: [Numpy-discussion] Upgrading from numpy 0.9.7.2416 to 0.9.9.2683

2006-06-27 Thread Travis Oliphant
Stefan van der Walt wrote: > On Tue, Jun 27, 2006 at 09:45:57AM -0700, Keith Goodman wrote: > > With r2691 I see > > In [7]: x = N.asmatrix(N.zeros((3,2)),float) > > In [8]: y = N.asmatrix(N.rand(3,1)) > > In [12]: x[:,1] = y > 0.5 > > In [13]: x > Out[13]: > matrix([[ 0., 1.], >[ 0.,

Re: [Numpy-discussion] Upgrading from numpy 0.9.7.2416 to 0.9.9.2683

2006-06-27 Thread Keith Goodman
On 6/27/06, Travis Oliphant <[EMAIL PROTECTED]> wrote: > Keith Goodman wrote: > > This works in numpy 0.9.7.2416 but doesn't work in numpy 0.9.9.2683: > > > > Numpy 0.9.9.2683 > > > > x = asmatrix(zeros((3,2), float)) > > y = asmatrix(rand(3,1)) > > y > > > > matrix([[ 0.49865026], > >[ 0.

Re: [Numpy-discussion] Upgrading from numpy 0.9.7.2416 to 0.9.9.2683

2006-06-27 Thread Travis Oliphant
Keith Goodman wrote: > This works in numpy 0.9.7.2416 but doesn't work in numpy 0.9.9.2683: > > Numpy 0.9.9.2683 > > x = asmatrix(zeros((3,2), float)) > y = asmatrix(rand(3,1)) > y > > matrix([[ 0.49865026], >[ 0.82675808], >[ 0.30285247]]) > > x[:,1] = y > 0.5 > x > > matrix([[ 0.,

Re: [Numpy-discussion] Upgrading from numpy 0.9.7.2416 to 0.9.9.2683

2006-06-27 Thread Stefan van der Walt
On Tue, Jun 27, 2006 at 09:45:57AM -0700, Keith Goodman wrote: > This works in numpy 0.9.7.2416 but doesn't work in numpy 0.9.9.2683: > > Numpy 0.9.9.2683 > > x = asmatrix(zeros((3,2), float)) > y = asmatrix(rand(3,1)) > y > > matrix([[ 0.49865026], >[ 0.82675808], >[ 0.30285247]

[Numpy-discussion] Upgrading from numpy 0.9.7.2416 to 0.9.9.2683

2006-06-27 Thread Keith Goodman
This works in numpy 0.9.7.2416 but doesn't work in numpy 0.9.9.2683: Numpy 0.9.9.2683 x = asmatrix(zeros((3,2), float)) y = asmatrix(rand(3,1)) y matrix([[ 0.49865026], [ 0.82675808], [ 0.30285247]]) x[:,1] = y > 0.5 x matrix([[ 0., 0.], [ 0., 0.], <--- this should be o

Re: [Numpy-discussion] (no subject)

2006-06-27 Thread yukihana
:―― INFORMATION ―: 不正・悪質なサイトを一切排除しておりますので 安心してご利用ください。 http://love-match.bz/pc/?02 :――: *・゜゜・*:.。. .。.:*・゜゜・*:.。..。:*・゜゜・*:.。..。:**・゜゜・* ━━━  □■ 不倫・ワリキリ専門の無料出会いサイト『Love☆Match』  --

Re: [Numpy-discussion] numpy.linalg.pinv has no docstring

2006-06-27 Thread joris
On Tuesday 27 June 2006 09:37, Nils Wagner wrote: [NW]: Please can you add a docstring to numpy.linalg.pinv. In case it might help, I added an example to the Numpy Example List (http://www.scipy.org/Numpy_Example_List) which illustrates the use of pinv(). J. Disclaimer: http://www.kuleu

Re: [Numpy-discussion] (no subject)

2006-06-27 Thread hitomi
:―― INFORMATION ―: 不正・悪質なサイトを一切排除しておりますので 安心してご利用ください。 http://love-match.bz/pc/?09 :――: *・゜゜・*:.。. .。.:*・゜゜・*:.。..。:*・゜゜・*:.。..。:**・゜゜・* お金と時間を持て余している人妻の間で、噂になってるあのサイト   [登録・利用料全て無料] http://love-match.bz/pc/?09 ━

[Numpy-discussion] numpy.linalg.pinv has no docstring

2006-06-27 Thread Nils Wagner
Hi Travis, Just now I saw that you have fixed the failing test. You have used pinv (pseudo inverse). Please can you add a docstring to numpy.linalg.pinv. Thanks in advance. Nils In [4]: numpy.linalg.pinv? Type: function Base Class: String Form: Namespace: Interactive Fil