c++ - C++ でのシン QR 分解

Question

長方形行列の「薄い」QR 分解のための使いやすい c++ ライブラリはありますか?
Eigen は完全な Q 行列のみをサポートしているようです。完全な Q を取得して一部の列を破棄することはできますが、最初から計算しない方が効率的でしょうか?

Answer 1

Newmatはまさにあなたが望むことを行います。

A を QR に分解するには、次のようにします。

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

A が 3x5 行列の場合、R は 3x3 になり、Q も 3x5 になります。