В MATLAB можно использовать хак, чтобы заставить встроенный алгоритм декомпозиции LU не использовать матрицу перестановок (то есть силу P=I), как показано здесь .
Есть ли подобный хак в python?
А именно, есть ли способ заставить scipy.linalg.lu (или любой другой популярный алгоритм LU) не использовать матрицу перестановок?