巨大な疎な対称行列のためのSpectraライブラリの中で最も速い固有値ソルバはどれですか？

Question

大規模な疎な対称行列（通常20,000〜60,000のオーダー）の固有値を解きたい。私はSymEigsShiftSolver/SymEigsSolverを活用しようとしましたが、私の固有値を生成するにはかなりの時間がかかります。スピードアップの方法はありますか？私は何を目指しているのかのデモコードを添付しました。私は何かを逃している？

メインINT（）{

Eigen::MatrixXd A = Eigen::MatrixXd::Random(10000, 10000); 
Eigen::MatrixXd M = A.transpose() * A; 

// Matrix operation objects 
DenseGenMatProd<double> op_largest(M); 
DenseSymShiftSolve<double> op_smallest(M); 

// Construct solver object, requesting the smallest 10 eigenvalues 
SymEigsShiftSolver< double, LARGEST_MAGN, DenseSymShiftSolve<double> > 
    eigs_smallest(&op_smallest, 10, 30, 0.0); 

eigs_smallest.init(); 
eigs_smallest.compute(); 
std::cout << "Smallest 10 Eigenvalues :\n" << 
    eigs_smallest.eigenvalues() << std::endl; 

return 0; 

}

Answer 1

次のコード例は、密行列を融合されます。まず、Eigen::SparseMatrix<double>を詳細にthereとして適切に組み立てることです。次に、Spectra側でそれぞれSparseSymShiftSolveヘルパーを使用します。

また、コンパイラの最適化をオンにしてコンパイルすることを忘れないでください。このようなヘッダー専用のlibでは重要です。