我正在重构一个现有算法,其中 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。