【发布时间】:2020-09-26 05:50:03
【问题描述】:
我的目标是在 Python 中将矩阵转换为 上三角 形式。我知道函数 scipy.linalg.lu 将进行 LU 分解并计算 both 上三角和下三角。我需要多次重复这个过程。如果有来自Scipy、blas 或lapack 的函数只计算上三角矩阵,我的算法的性能将得到显着提高。
你能推荐这样的功能吗?非常感谢!
【问题讨论】:
标签: python-3.x matrix scipy lapack blas
我的目标是在 Python 中将矩阵转换为 上三角 形式。我知道函数 scipy.linalg.lu 将进行 LU 分解并计算 both 上三角和下三角。我需要多次重复这个过程。如果有来自Scipy、blas 或lapack 的函数只计算上三角矩阵,我的算法的性能将得到显着提高。
你能推荐这样的功能吗?非常感谢!
【问题讨论】:
标签: python-3.x matrix scipy lapack blas
如果您点击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。
【讨论】:
getrf 同时计算 U 和 L。