Hacktoberfest 2026: những issue maintainer đã đánh dấu cho tháng Mười, đang mở và phù hợp người mới. Xem issue Hacktoberfest

wrong eigh eigenvectors for complex valued matrix

Đang mở
#307 0 bình luận 0 reaction 0 người được giao Xem trên GitHub

Chưa có ai nhận issue này.

Đánh giá

Độ khó
4/5
Thời gian dự kiến
3-5 ngày
Mức phù hợp với người mới
42/100
Loại issue
Lỗi
Độ rõ ràng
Khá rõ ràng
Mức độ hoạt động
Đình trệ
Công nghệ
rust
Lĩnh vực
hpc

Hướng nghiên cứu

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.

Do mô hình lập chỉ mục viết ra từ nội dung của issue.

Mô tả

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);
    }
Ngôn ngữ chính
Rust
Star
452
Fork
95
Không có pull request nào được merge trong 30 ngày

Hướng dẫn đóng góp

Chưa lập chỉ mục được hướng dẫn đóng góp cho kho mã nguồn này

Bắt đầu từ đâu

  1. Đọc hết issue, rồi đọc hướng dẫn đóng góp của dự án.
  2. Bình luận trên issue rằng bạn sẽ nhận — tránh hai người làm cùng một việc.
  3. Fork repository và làm thay đổi trên một nhánh.
  4. Mở pull request có tham chiếu số hiệu của issue.

Issue khác của rust-ndarray/ndarray-linalg

Tất cả issue của rust-ndarray/ndarray-linalg

Issue tương tự

Thêm issue về Rust

Nhận issue mới trong hộp thư của bạn

Bản tóm tắt ngắn những issue GitHub phù hợp với người mới.