Эффективно рассчитывать сумму редкого плотного умножения в Python ⇐ Python

[Anonymous]

У меня есть проблема, когда у меня есть список редких квадратных массивов (матрицы CSC) a = [a1, a2, ..., an], а затем список плотных квадратных массивов
b = [b1, b2, ..., bn] < /p>
Я хочу рассчитать сумму C = A1 B1 + a2 B2 + ... An Bn. Разреженные массивы очень скудны, а размеры вовлеченных матриц могут составлять до миллиона миллионов. Я хотел бы сделать это параллельно. Тем не менее, все, что я пробовал до сих пор, не смогло обеспечить какого-либо ускорения. Я знаю, что NUMBA трудно использовать, учитывая, что мы умножаем скудность с плотными матрицами. Я хотел бы несколько предложений написать эффективный параллельный код для вышеизложенного.

Подробнее здесь: https://stackoverflow.com/questions/796 ... -in-python