我的目标是在Python中将矩阵转换为上三角形。我知道该函数scipy.linalg.lu
将进行 LU 分解并计算上三角形和下三角形。我需要重复此过程几次。如果有来自Scipy
、blas
或lapack
的函数只计算上三角矩阵,我的算法的性能将得到显著提高。
你能建议这样的功能吗?非常感谢!
如果您在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
.