Scipy.linalg - реализация умножения матрицы, все в C или в Фортране?

Question

Недавно я изучал, как умножение на (в основном разреженных) матрицах реализовано оптимальным образом в python библиотеке scipy.linalg. Из здесь я знаю, что в основном алгоритмы являются обертки некоторых библиотек Fortran или C.

Мой вопрос в том, как я мог узнать, есть ли какие-либо методы умножения в Scipy, которые на самом деле написаны внутри самого python? Я также не очень уверен, как это работает с умножением матриц Numpy.

Answer 1

Умножение матриц, которое требует больших вычислительных ресурсов, обычно не выполняется в Python. Типичное матричное умножение будет O (n ^ 3), поэтому Numpy и Scipy используют умножение матриц на основе C и используют обертки, чтобы привести это к Python.