Can Sage multiply sparse matrices? If so, does it require some special syntax? I've been trying to run the following code:
M = Matrix(QQ, 1000001, 1000001, {(3,2): 27, (2,98): 71})
M2 = M*M
and it hasn't been working. If M is being converted to a dense matrix for multiplication then of course I wouldn't expect this to run. But it shouldn't be hard to do as sparse matrices, right?
I think the relevant code is here: https://github.com/sagemath/sage/blob...
The algorithm does take into account the sparse nature, but it does have a triple-nested for-loop, enumerating all the indices. It may be that for moderate sparsity this is faster than just enumerating the non-zero entries, but especially for your example it's pretty clear a much more efficient implementation is possible.
Thanks for reporting, this is now trac ticket 26532.
Asked: 2018-10-18 20:49:55 +0100
Seen: 422 times
Last updated: Oct 23 '18
Copyright Sage, 2010. Some rights reserved under creative commons license. Content on this site is licensed under a Creative Commons Attribution Share Alike 3.0 license.
