c中的薄QR分解

是否存在易于使用的c库,用于矩形矩阵的“薄” QR分解?
本征似乎只支持完整的Q矩阵.我可以取一个完整的Q并丢弃一些列,但是不从头开始计算它们会更有效吗?
最佳答案
Newmat正是您想要的.

要将A分解为QR,您可以执行以下操作:

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

如果A为3×5矩阵,则R将为3×3,Q也将为3×5.

转载注明原文:c中的薄QR分解 - 代码日志