Почему сборка разреженной матрицы выполняется медленнее, чем сборка плотной (полной) матрицы? ⇐ Python

[Anonymous]

У меня есть две функции (i) getKglobal и (ii) getKglobal_Sp. Однако время для разреженной функции примерно в 5 раз медленнее, чем для полноматричной формулировки. Я понимаю, что добиться производительности от разреженности может быть непросто. Но замедление в 5 раз кажется слишком плохим. Какую «большую и очевидную» ошибку я упускаю?
def getKglobal(geom,mat,nel,tdof):
k_global = np.zeros((tdof, tdof))
for i in range(0,nel):
dof_idx = geom.dof_idx
coord = geom.coord
kloc=getkelement(coord,mat)
k_global[np.ix_(dof_idx, dof_idx)]+=kloc
return k_global

def getKglobal_Sp(geom,mat,nel,tdof):
from scipy.sparse import lil_matrix
k_global = lil_matrix((tdof, tdof))
for i in range(0,nel):
dof_idx = geom.dof_idx
coord = geom.coord
kloc=getkelement(coord,mat)
for a, global_a in enumerate(dof_idx):
for b, global_b in enumerate(dof_idx):
k_global[global_a, global_b] += kloc[a, b]
return k_global.tocsr()


Подробнее здесь: https://stackoverflow.com/questions/793 ... x-assembly