Scipy，Blas或Lapack中仅计算上三角矩阵的函数

Question

我的目标是在Python中将矩阵转换为上三角形式。我知道函数scipy.linalg.lu将进行LU分解并计算both上和下三角。我需要重复此过程几次。如果Scipy，blas或lapack中的函数仅计算上三角矩阵，则我的算法的性能将得到显着改善。

您能建议这样的功能吗？非常感谢！

Answer 1

如果您在scipy.linalg.lu_factor页面上单击“源”，则https://docs.scipy.org/doc/scipy-0.14.0/reference/generated/scipy.linalg.lu_factor.html会将您带到https://github.com/scipy/scipy/blob/v0.14.0/scipy/linalg/decomp_lu.py#L17

这表明您需要的LAPACK函数为*getrf。