Armadilloを使用して固有ペアを計算するにはどうすればよいですか？

Question

私はアルマジロを持って、私は「テスト」私がコンパイルしようとすると、それ

int main(){ 
    double myMatrix[6][6]; 
    for (int i = 0; i < 6; i++){ 
     for (int j = 0; j < 6; j++){ 
      myMatrix[i][j] = i+2*j; 
     } 
    } 

    mat ARMA_L; 
    vec ARMA_eigenval; 
    mat ARMA_eigenvec; 

    for (int i = 0; i < 6; i++){ 
     for (int j = 0; j < 6; j++){ 
      ARMA_L(i,j) = myMatrix[i][j]; 
     } 
    } 
    eig_gen(ARMA_eigenval, ARMA_eigenvec, ARMA_L); 
    return 0; 
} 

で実行をやっている自分自身を慣れたかった、それは（eig_genへの呼び出しのための機能がありませんと言う） が、ドキュメントしかし、私がやっているべきことをやっているようだ。

Answer 1

eig_gen()機能に間違ったタイプを使用しています。 ARMA_eigenvalとARMA_eigenvecはそれぞれcx_vecとcx_mat次のようになります。

cx_vec ARMA_eigenval; 
cx_mat ARMA_eigenvec; 

documentationの例を参照してください：

mat A = randu<mat>(10,10); 

cx_vec eigval; 
cx_mat eigvec; 

eig_gen(eigval, eigvec, A);