c++ - 问题与 eigs_sym 获取具有最小幅度的特征值

Question

我正在尝试获取有限数量的平方对称矩阵的最小量值的特征值。 为此，我首先使用 Armadillo 文档中的示例 (http://arma.sourceforge.net/docs.html#eigs_sym):

sp_mat A = sprandu<sp_mat>(1000, 1000, 0.1);
sp_mat B = A.t()*A;
arma::vec eigval;
mat eigvec;
eigs_sym(eigval, eigvec, B, 10, "sm");//i add "sm" to get the eigen-  
                                      //values of smallest magnitude
cout<<eigval<<endl; 

这里我得到一个错误，说 ddcomposition 失败 [无法收敛]。

但是，当我这样调用 eigs_sym 时:

eigs_sym(eigval, eigvec, B, 10); //obtain the eigen-values with 
                                 //LARGEST magnitude (default call)

这个工作很好，我得到了预期的结果:

1.1596e+02
1.1680e+02
1.1785e+02
1.1815e+02
1.1927e+02
1.2017e+02
1.2108e+02
1.2256e+02
1.2323e+02
2.5413e+03

我在 Ubuntu 操作系统上，这是我的 .pro 文件 (Qt):

 LIBS += -lgsl -lgslcblas -lX11 -lpthread -llapack  -lm -fopenmp   
         -larmadillo

有解决这个问题的想法吗？

谢谢

Answer 1

我通过选择更多的特征值来提取来解决这个问题。 显然，要提取的特征值数量较少会使特征求解器注意收敛。如果你更换

eigs_sym(eigval, eigvec, B, 10,"sm")

由

eigs_sym(eigval, eigvec, B, 100,"sm")

这会起作用。