numpy: инвертирование верхней треугольной матрицы
В numpy / scipy , каков канонический способ вычисления инверсии верхней треугольной матрицы?
Матрица хранится как массив 2D numpy с нулевыми субдиагональными элементами, и результат также должен быть сохранен как 2D-массив.
Лучшее, что я нашел до сих пор, это scipy.linalg.solve_triangular(A, np.identity(n)) . Это оно?
- Итерации по матрице, суммирование по некоторым строкам и добавление результата в другой массив
- Как найти вырожденные строки / столбцы в ковариационной матрице
- nanmean матричных столбцов в numpy
- В Python есть какой-нибудь способ добавить?
- Заменить вспомогательную часть матрицы на другую маленькую матрицу в numpy
Фактически, это не обычная процедура инверсии. scipy.linalg.solve – это канонический способ решения матрично-векторного или матрично-матричного уравнения, и ему может быть дана явная информация о структуре матрицы, которую он будет использовать для выбора правильной подпрограммы (возможно, эквивалент BLAS3 dtrsm в этом случае).
LAPACK включает в себя doptri для этой цели, а scipy.linalg представляет собой открытый интерфейс лапка «C». Если обратная матрица действительно то, что вы хотите, тогда вы можете попробовать ее использовать.
- Переиндексирование уровня MultiIndex в произвольный порядок в Pandas
- glib.GError: Ошибка интерпретации файла изображения JPEG (Неподдерживаемый маркер типа 0x05)
- Как заставить процессы записывать в массив основной программы?
- Как инициализировать массив NumPy с различным значением по умолчанию для каждого столбца?
- Обратная матрица с использованием numpy
- Django – Как отобразить матрицу данных с суммами в конце каждой строки и столбца
- Нумерация, чтобы проверить, существует ли решение, чтобы каждая строка была <0?
- встроенная функция python для уменьшения матрицы
- Получить индексы матрицы из верхнего треугольника
- В numpy, как эффективно перечислить все подматрицы фиксированного размера?
- Numpy / Python: эффективная матрица как умножение декартова произведения входной матрицы