Я делаю простой разреженный рост матрицы, a**16
, используя scipy-0.17. (Обратите внимание, что не умножение по элементам). Однако на моих машинах (работает стабилизатор Debian и Ubuntu LTS) это в десять раз медленнее, чем использование цикла for или что-то глупое, как a*a*a*a*a*a*a*a*a*a*a*a*a*a*a*a
. Это не имеет смысла, поэтому я предполагаю, что я делаю что-то неправильно, но что?
import scipy.sparse
from time import time
a=scipy.sparse.rand(2049,2049,.002)
print ("Trying exponentiation (a**16)")
t=time()
x=a**16
print (repr(x))
print ("Exponentiation took %f seconds\n" % (time()-t))
print ("Trying expansion (a*a*a*...*a*a)")
t=time()
y=a*a*a*a*a*a*a*a*a*a*a*a*a*a*a*a
print (repr(y))
print ("Expansion took %f seconds\n" % (time()-t))
print ("Trying a for loop (z=z*a)")
t=time()
z=scipy.sparse.eye(2049)
for i in range(16):
z=z*a
print (repr(z))
print ("Looping took %f seconds\n" % (time()-t))
# Sanity check, all approximately the same answer, right?
assert (abs(x-z)>=1e-9).nnz==0
assert (abs(x-y)>=1e-9).nnz==0