Armadillo eigs_sym (A, k) untuk matriks jarang kompleks A

Question

Untuk mencari 10 nilai eigen terkecil dari matriks renggang 'A', kode minimal di bawah ini berfungsi dengan baik:

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

Tetapi ketika saya mengubah A menjadi matriks renggang yang kompleks, seperti:

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

dengan bendera kompilasi yang sama, saya mendapatkan kesalahan fungsi tidak cocok berikut ini:

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

Saya tahu dari http://arma.sourceforge.net/docs.html#config_hpp bahwa

ARMA_USE_ARPACK Mengaktifkan penggunaan ARPACK, atau pengganti ARPACK berkecepatan tinggi. Armadillo membutuhkan ARPACK untuk dekomposisi eigen dari matriks jarang kompleks, yaitu. eigs_gen (), eigs_sym () dan svds ()

jadi saya mengubah file config.hpp dan berikut adalah baris yang sesuai di config.hppfile saya :

#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

info lebih lanjut: Saya dapat menjalankan arpack dari gfortran tanpa masalah.

Ada ide bagaimana melakukan triknya? Terima kasih atas bantuan Anda sebelumnya.

Answer 1

Ini adalah batasan inheren dari perpustakaan. Menurut dokumentasi (penekanan saya):

eigs_sym jumlah terbatas nilai eigen & vektor eigen dari matriks riil simetris renggang

eigs_gen jumlah terbatas nilai eigen & vektor eigen dari matriks kuadrat umum renggang

Anda sebaiknya menggunakan eigs_genyang memungkinkan untuk matriks kompleks. Atau Anda harus mengubah matriks menjadi matriks padat dan menggunakan eig_sym.