wrong eigh eigenvectors for complex valued matrix

Aperta
#307 0 commenti 0 reazioni 0 assegnatari Vedi su GitHub

Nessuno ha ancora preso questa issue.

Valutazione

Difficoltà
4/5
Tempo stimato
3-5 giorni
Idoneità per principianti
42/100
Tipo di issue
Bug
Chiarezza
Abbastanza chiara
Stato di attività
Ferma
Stack tecnologico
rust
Ambito
hpc

Direzione di ricerca

Start by reproducing the provided complex-valued matrix example and trace the implementation behind eigh with UPLO::Lower. Compare the eigenvector result for the matrix and its Hermitian conjugate, then add a regression test based on the shown assertions; done means the original matrix's eigenvectors satisfy the eigenvalue equation.

Scritto dal modello di indicizzazione a partire dal testo della issue.

Descrizione

The resulting eigenvectors for eigh of complex valued matrix is wrong.

If eigh is called with the hermitian conjugate of the matrix,
the resulting eigenvectors are the eigenvectors of the original matrix.

   let values = vec![
        c64::new(1.0, 0.0), c64::new(0.0, -2.0),
        c64::new(0.0, 2.0), c64::new(5.0, 0.0),
    ];
    let matrix = Array2::<c64>::from_shape_vec([2; 2], values).unwrap();
    // assertion passed
    let (eigvalues, eigvecs) = matrix.t().mapv(|v| v.conj()).eigh(UPLO::Lower).unwrap();
    for i in 0 .. 2 {
        let dot = matrix.dot(&eigvecs.column(i));
        let mul = &eigvecs.column(i).mapv(|c| c*eigvalues[i]);
        assert_close_l2!(&dot, &mul, 1.0e-10);
    }
    // assertion failed.
    let (eigvalues, eigvecs) = matrix.eigh(UPLO::Lower).unwrap();
    for i in 0 .. 2 {
        let dot = matrix.dot(&eigvecs.column(i));
        let mul = &eigvecs.column(i).mapv(|c| c*eigvalues[i]);
        assert_close_l2!(&dot, &mul, 1.0e-10);
    }
Lingua principale
Rust
Stelle
452
Fork
95
Nessuna PR unita negli ultimi 30g

Guida per i contributori

Nessuna guida per i contributori indicizzata per questo repository

Come iniziare

  1. Leggi tutta la issue e poi la guida ai contributi del progetto.
  2. Commenta sulla issue per dire che te ne occupi tu — evita che due persone facciano lo stesso lavoro.
  3. Fai un fork del repository e lavora su un branch.
  4. Apri una pull request che faccia riferimento al numero della issue.

Altre issue di rust-ndarray/ndarray-linalg

Tutte le issue di rust-ndarray/ndarray-linalg

Issue simili

Altre issue su Rust

Ricevi le nuove issue nella tua casella

Un breve riepilogo di issue GitHub adatte ai principianti.