[SciPy-dev] matrixmultiply alias
a.schmolck at gmx.net
Fri Jul 25 12:36:15 CDT 2003
Pietro Berkes <p.berkes at biologie.hu-berlin.de> writes:
> 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.
More information about the Scipy-dev