c++ - c++ 中的薄 QR 分解

Question

是否有易于使用的 c++ 库来对矩形矩阵进行“瘦”QR 分解？
Eigen 似乎只支持全 Q 矩阵。我可以获取完整的 Q 并丢弃一些列，但是不计算它们是否更有效？

Answer 1

Newmat完全按照您的意愿行事。

要将 A 分解为 QR，你可以这样做:

Matrix Q = A;
UpperTriangularMatrix R;
QRZ(Q, R)

如果 A 是 3x5 矩阵，则 R 将为 3x3，Q 也将为 3x5。