1

我目前正在尝试重构 R 的kappa条件数估计函数的一些函数,该函数通过以下方式估计矩阵 X 的条件数:

  • 计算 X 的 QR 分解。
  • 调用LAPACK'sdtrconLINPACK's dtrco(取决于用户系统的底层依赖项是什么),并计算 R 的条件数,上三角矩阵应该具有与 X 相同的条件数(参见此处)。

我一直在尝试了解LAPACKLINPACK算法的作用,因为它可能对我自己的编码非常有用。

我设法找到了LINPACK使用的算法,这在此处描述,但没有找到LAPACK's 算法的起源。R 函数中的注释kappa表明它们使用不同的算法(请参见此处),但我不确定...

长话短说,我的问题是:

有谁知道LAPACK'dtrconLINPACK'是否使用相同的算法,如果没有, '使用dtrco的是什么算法?LAPACKdtrcon

先感谢您!

4

0 回答 0