numpy: инвертирование верхней треугольной матрицы

В numpy / scipy , каков канонический способ вычисления инверсия верхней треугольной матрицы?

Матрица сохраняется в виде массива 2D numpy с нулевыми субдиагональными элементами, и результат также должен быть сохранен как массив 2D.

edit Лучшее, что я нашел, это scipy.linalg.solve_triangular (A, np.identity (n)) . Это все?