[SciPy-User] Sparse matrix multiply

Jaakko Luttinen jaakko.luttinen@aalto...
Thu Mar 22 04:34:52 CDT 2012


Hi!

Why do I get two different results for the code below?

import numpy as np
import scipy.sparse as sp
A = sp.rand(20,20,density=0.1)
B = sp.rand(20,20,density=0.1)
np.multiply(A,B).sum()
# out: 21.058793740984925
A.multiply(B).sum()
# out: 0.76482546226069481

Am I missing something?
I think numpy.multiply should either return the correct answer or an
error that it can't compute the correct answer.

Regards,
Jaakko


More information about the SciPy-User mailing list