アルマジロC++ル分解：間違っていますか？

Question

私はいくつかの本当に基本的なマトリックス操作のためにarmadilloライブラリを使用しています。だから私は3x5の行列を持っており、私はLUの分解を実行したい。

#include "myIncludes.hpp" 
#include <armadillo> 

using namespace std; 
using namespace arma; 


int main() 
{ 
mat L, U; 
mat A(3, 5); A = "1 1 1 0 0; 1 1 1 0 0; 2 1 2 0 1"; 

arma::lu(L,U,A); 

return 0; 
} 

このコードを実行すると、次のコードが返されます。警告：lu（）：分解に失敗しました。しかし、matlabは実際にこの行列を分解することができます。

アルマジロは広く使用され、テストされていますか？それとも注意して使用する必要がありますか？

Answer 1

ソースコードをチェックしたところ、正方形の可逆行列に対してのみ安定であることが知られている部分ピボットに依存しているようです。

私の行列は正方形ではないので、私は外です。

出典：Eigen library documentation