我目前正在尝试重构 R 的kappa
条件数估计函数的一些函数,该函数通过以下方式估计矩阵 X 的条件数:
- 计算 X 的 QR 分解。
- 调用
LAPACK
'sdtrcon
或LINPACK
'sdtrco
(取决于用户系统的底层依赖项是什么),并计算 R 的条件数,上三角矩阵应该具有与 X 相同的条件数(参见此处)。
我一直在尝试了解LAPACK
和LINPACK
算法的作用,因为它可能对我自己的编码非常有用。
我设法找到了LINPACK
使用的算法,这在此处描述,但没有找到LAPACK
's 算法的起源。R 函数中的注释kappa
表明它们使用不同的算法(请参见此处),但我不确定...
长话短说,我的问题是:
有谁知道LAPACK
'dtrcon
和LINPACK
'是否使用相同的算法,如果没有, '使用dtrco
的是什么算法?LAPACK
dtrcon
先感谢您!