[SciPy-dev] matrixmultiply alias

[Alexander Schmolck]

Pietro Berkes <p.berkes at biologie.hu-berlin.de> writes:

> Hi!
> 
> I would like to propose a function that could be used as alias for 
> scipy.matrixmultiply (code attached in 'matrixmultiply.py'). The function 
> delegate the work to scipy.dot if the arguments it receives aren't 2D, 
> otherwise it uses the BLAS routine 'gemm'. If BLAS fails, an appropriate 
> error message is raised.
> The function should now be backward-compatible and is sensibly faster than the 
> original one for matrices larger than 40x40. Here are the benchmarks (the 
> code to perform the benchmarks is included in 'matmult_bench.py'):

For those who are interested, a similiar patch ("dotblas" submitted by me but
coded mainly by Richard Everson) is now part of Numeric (version 23). For
1000x1000 single float matrices I think the speedup is around factor 80.

'as