[SciPy-dev] Linalg2 benchmarks

Jochen Küpper jochen at jochen-kuepper.de
Fri Apr 5 00:19:07 CST 2002


On Thu, 4 Apr 2002 22:01:21 -0700 Travis Oliphant wrote:

Travis> The other question is are you using ATLAS for Numeric as well?  

Yes, sure.  Anybody on this list not linking Numeric against
lapack/blas?

Travis> Which version of Numeric do you have?

usually latest cvs

Travis> Those numbers mean your Numeric must be using a more optimized
Travis> version of lapack, anyway.

Numeric is linked against lapack and blas. And surely I don't use the
Fortran reference implementation of BLAS...

For the PIII results I am not sure whether Numeric might be linked
against an older ATLAS version.  I'll run the test tomorrow, making
sure I use the same lapack and blas for Numeric and scipy.  

For the Duron results below scipy and Numeric use the same
ATLAS-3.3.13 (exactly the same, gcc-3 compiled IIRC).

,----
| >>> scipy.linalg2.basic.test()
| ................................
|            Finding matrix determinant
|       ==================================
|       |    contiguous     |   non-contiguous
| ----------------------------------------------
|  size |  scipy  | Numeric |  scipy  | Numeric
|    20 |   0.39  |   0.78  |   0.38  |   0.90     (secs for 2000 calls)
|   100 |   0.83  |   1.59  |   0.80  |   2.12     (secs for 300 calls)
|   500 |   0.92  |   1.14  |   0.93  |   1.34     (secs for 4 calls)
| .
|       Solving system of linear equations
|       ==================================
|       |    contiguous     |   non-contiguous
| ----------------------------------------------
|  size |  scipy  | Numeric |  scipy  | Numeric
|    20 |   0.56  |   0.67  |   0.54  |   0.75     (secs for 2000 calls)
|   100 |   0.88  |   1.16  |   0.88  |   1.88     (secs for 300 calls)
|   500 |   0.99  |   1.06  |   0.99  |   1.27     (secs for 4 calls)
| .
|            Finding matrix inverse
|       ==================================
|       |    contiguous     |   non-contiguous
| ----------------------------------------------
|  size |  scipy  | Numeric |  scipy  | Numeric
|    20 |   0.90  |   1.23  |   0.90  |   1.35     (secs for 2000 calls)
|   100 |   2.32  |   3.46  |   2.30  |   3.95     (secs for 300 calls)
|   500 |   2.68  |   3.52  |   2.67  |   3.72     (secs for 4 calls)
| .
| ----------------------------------------------------------------------
| Ran 35 tests in 53.275s
|  
| OK
`----


Greetings,
Jochen
-- 
Einigkeit und Recht und Freiheit                http://www.Jochen-Kuepper.de
    Liberté, Égalité, Fraternité                GnuPG key: 44BCCD8E
        Sex, drugs and rock-n-roll




More information about the Scipy-dev mailing list