Armadillo eigs_sym (A, k) для комплексной разреженной матрицы A

Question

Чтобы найти 10 наименьших собственных значений разреженной матрицы 'A', приведенный ниже минимальный код работает хорошо:

g++ -std=c++17  -o test_sparse.o -c test_sparse.cpp
g++ -std=c++17  -o myapp test_sparse.o -larmadillo -larpack

#include <armadillo>
#include <iostream>
int main(){
    arma::SpMat<double> A = arma::sprandu(100,100,0.1) ;
    A = A.t()*A ;
    arma::dvec e = arma::eigs_sym(A,10,"sm") ;
    std::cout << e ; 
    return 0 ;
}

Но когда я меняю A на сложную разреженную матрицу, например:

#include <armadillo>
#include <iostream>
#include <complex>
int main(){
    arma::SpMat<arma::cx_double> A = arma::sprandu<arma::SpMat<arma::cx_double>>(100,100,0.1) ;
    A = A.t()*A ;
    arma::dvec e = arma::eigs_sym(A,1,"sm") ;
    std::cout << e ; 
    return 0 ;
}

с теми же флагами компиляции я получаю следующую ошибку функции отсутствия соответствия:

g++ -std=c++17  -o test_sparse.o -c test_sparse.cpp
test_sparse.cpp:8:43: error: no matching function for call to ‘eigs_sym(arma::SpMat<std::complex<double> >&, int, const char [3])’
    8 |     arma::dvec e = arma::eigs_sym(A,1,"sm") ;                    ^
make: *** [Makefile:47: test_sparse.o] Error 1

Я знаю из http://arma.sourceforge.net/docs.html#config_hpp что

ARMA_USE_ARPACK Разрешить использование ARPACK или высокоскоростной замены ARPACK. Armadillo требует ARPACK для собственного разложения сложных разреженных матриц, т.е. eigs_gen (), eigs_sym () и svds ()

поэтому я меняю файл config.hpp, и вот соответствующая строка в моем config.hppфайле:

#if !defined(ARMA_USE_NEWARP)
#define ARMA_USE_NEWARP
#endif
#if !defined(ARMA_USE_ARPACK)
#define ARMA_USE_ARPACK
#endif
#if !defined(ARMA_USE_SUPERLU)
#define ARMA_USE_SUPERLU
#endif

подробнее: я могу запустить arpack из gfortran без проблем.

Есть идеи, как сделать трюк? Заранее благодарю за помощь.

Answer 1

Это внутреннее ограничение библиотеки. Согласно документации (выделено мной):

eigs_sym ограниченное количество собственных значений и собственных векторов разреженной симметричной вещественной матрицы

eigs_gen ограниченное количество собственных значений и собственных векторов разреженной квадратной матрицы общего вида

Лучше использовать то, eigs_genчто позволяет создавать сложные матрицы. Или лучше преобразовать матрицу в плотную и использовать eig_sym.