Блог

Чтобы найти 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 ;
}
  

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

 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> >amp;, 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 без проблем.

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

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

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

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

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