armadillo c ++ lu分解：我做错了吗？

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;
}

当我运行此代码时，它返回：warning：lu（）：分解失败。但是，matlab确实可以分解这个矩阵。

犰狳是否经过广泛使用和测试，或者我应该小心使用它吗？

Answer 1

我检查了源代码，它似乎依赖于部分旋转，已知只有方形可逆矩阵才能稳定。

我的矩阵不是正方形，所以我出去了。

来源：Eigen library documentation