[Numpy-discussion] Matlab/Numeric/numarray benchmarks

[Simon Burton]

Hi all,

Here are some benchmarks, measured in seconds on a 1.5GHz celeron.
Each test does a matrix add (1000x1000), mul (1000x1000) and eigenvalue 
find (500x500).

Matlab:
0.0562
1.5180
3.7630

Numeric:
0.0962309813499
1.73247330189
3.72153270245

numarray:
7.17220497131
19.3960719109
5.72376401424

I have attached the code. Looks like numarray is way behind on the basic 
linear algebra stuff.

We have (so far) chosen to go with numarray for our scientific 
computations, but will be needing fast add/multiply.
I am surmising that these methods just have not been pluged into the 
native BLAS/ATLAS lib.

We will also be needing other solvers from LAPACK such as dpotrf/dposv, 
and some of the special functions (bessel) already implemented
in scipy.

I understand that Todd is working on scipy/numarray compatability. How 
is that progressing, and what should we be doing to get the above 
functionality into numarray ? I have been pokeing around the code 
already, and am able to help out with that.

bye for now,

Simon.

-- 
Simon Burton, B.Sc.
Licensed PO Box 8066
ANU Canberra 2601
Australia
Ph. 61 02 6249 6940
http://arrowtheory.com 

-------------- next part --------------
A non-text attachment was scrubbed...
Name: benchmark.py
Type: text/x-python
Size: 1786 bytes
Desc: not available
URL: <http://mail.python.org/pipermail/numpy-discussion/attachments/20050103/dbc2340b/attachment.py>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: bench.m
Type: text/x-objcsrc
Size: 272 bytes
Desc: not available
URL: <http://mail.python.org/pipermail/numpy-discussion/attachments/20050103/dbc2340b/attachment.bin>