是否有已发布的 RQ 分解列主要算法?

计算科学 线性代数
2021-12-15 11:10:23

我正在重构一个现有算法,其中 RQ 分解(而不是更常见的 QR)将非常有用。

关于该主题的大多数常见书籍(例如Golub 和 Van Loan)广泛讨论 QR 并提供常见实现的算法(例如,针对 Householder 的算法 5.2.1 或针对 Modified Gram Schmidt 的算法 5.2.5)。

不幸的是,它没有讨论 RQ,除了关于 LAPACK 例程的小注释。维基百科暗示差异是微不足道的,其中“RQ 分解是 A 的行的 Gram-Schmidt 正交化,从最后一行开始”。

由于我正在处理的矩阵是列优先的,所以这并不是特别有吸引力。

是否有著名的 RQ 列主要实现,尤其是使用 MGS 或 Householder?或者即使有一些对列数据有效的阻塞算法?

其他细节:这个特定的问题不需要旋转,我只需要 R - 我不需要 Q。

0个回答
没有发现任何回复~