Можно ли использовать параллельную ARPACK в scipy?

Question

Я использую scipy.sparse.linalg.eigs на некоторых больших матрицах, и не удивительно, что это требует времени. Итак, я искал способы ускорить это. Я понимаю, что под капотом scipy-код использует ARPACK, и есть параллельная версия ARPACK, которая использует MPI. Возможно ли использовать scipy параллельную версию ARPACK без особых проблем? Если да, то как?

(следует отметить, что эквивалентно MATLAB по eigs, действительно кажется, многопоточный, так что может быть наименее болезненным вариантом.)

Answer 1

Казалось бы, что (MPI-) параллельная версия ARPACK совершенно другой проект под названием PARPACK:.

«параллельная версия библиотеки ARPACK в настоящее время воевавших в передачи сообщений слои в настоящее время поддерживаются следующие BLACS и MPI параллельный ARPACK (PARPACK) предоставляется в качестве дополнения к текущему ARPACK. библиотеки (версия 2.1) ».

Вы посмотрели PETsc4py? Или, может быть, даже

«исследовать вызов параллельного разреженный линейной алгебры библиотеки как касп или cuSPARSE из Python, если скорость ваше беспокойство и у Вас есть NVIDIA GPU.»

(см this answer)