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

oneMKL potrf on GPU reports the failing column relative to its 32-column panel

Đang mở
  • Có gì không ổn?
#755 0 bình luận 0 reaction 0 người được giao Xem trên GitHub

Maintainer thường phản hồi trong vòng 1 ngày

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

Đánh giá

Độ khó
4/5
Thời gian dự kiến
1-2 ngày
Mức phù hợp với người mới
55/100
Loại issue
Lỗi
Độ rõ ràng
Khá rõ ràng
Mức độ hoạt động
Sôi nổi
Công nghệ
julia
Lĩnh vực
backend, testing

Hướng nghiên cứu

Start with the onemkl<T>potrf shim and the LinearAlgebra.LAPACK.potrf! changes referenced in #754; compare the raw oneMKL info with the failing diagonal entry in the reproducer. Check the existing @test_broken coverage and add a focused regression test for the affected GPU types. Done means the reported failure column is correct, with a raw-shim test retaining coverage of oneMKL's behavior if a workaround is added.

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

Mô tả

On the GPU, oneMKL's potrf returns the wrong info when the matrix stops being positive definite after column 32. This is oneMKL 2026.1, through oneAPI_Support_jll 0.11.0, on an Intel Data Center GPU Max 1100.

  • Real types get the column's position inside its 32-column panel, i.e. (j - 1) % 32 + 1.
  • Complex types get panel * 1024 + position.

Columns 1 to 32 are reported correctly. The shim onemkl<T>potrf returns computation_error::info() unchanged, so the wrong value comes from oneMKL itself, not from oneAPI.jl.

Reproducer

using oneAPI, LinearAlgebra
using oneAPI: Support, sycl_queue, global_queue

n, j = 100, 50
G = randn(n, n); S = Matrix(Hermitian(G * G' / n + n * I)); S[j, j] = -1
A = oneMatrix(S)
q = sycl_queue(global_queue(oneAPI.context(A), oneAPI.device(A)))
lw = Support.onemklDpotrf_scratchpad_size(q, 'L', n, n)
Support.onemklDpotrf(q, 'L', n, A, n, oneVector{Float64}(undef, lw), lw)   # returns 18, expected 50

Observed values for an HPD matrix with A[j, j] = -1. Results are the same for 'L' and 'U'.

type n j info
Float32 / Float64 50 17 17
Float32 / Float64 50 33 1
Float32 / Float64 100 50 18
Float32 / Float64 100 100 4
Float32 / Float64 257 257 1
ComplexF32 / ComplexF64 100 50 1042
ComplexF32 / ComplexF64 100 97 3073
ComplexF32 / ComplexF64 257 254 7198

In every case the factor itself still follows LAPACK potf2. The first diagonal entry that is not a finite positive number is exactly column j, so a caller can recover the right column from the diagonal. getrf's info is correct. The environment variables OverrideDefaultFP64Settings and IGC_EnableDPEmulation make no difference.

Impact. cholesky, issuccess and isposdef only check whether info is nonzero, so they are unaffected. Only the exact column number is wrong, for example in a PosDefException. #754 makes LinearAlgebra.LAPACK.potrf! return info and marks these cases @test_broken.

Possible handling

  1. Report the bug to Intel, which is the real fix.
  2. Optionally work around it in Julia: when info > 0, take the first non-positive diagonal entry as the column. This costs one small diagonal copy on the failure path only, since the shim already synchronizes. If you prefer to carry a workaround, a separate test on the raw shim would keep tracking the oneMKL bug.
Ngôn ngữ chính
Julia
Star
218
Fork
36
Merge trung bình
1 ngày 6 giờ
Pull request đã merge (30 ngày)
65

Chuẩn bị môi trường

Dự án này không cung cấp dev container, Dockerfile hay hướng dẫn đóng góp, nên bạn cần tự thiết lập môi trường: hãy bắt đầu từ README và xem hướng dẫn đóng góp lần đầu của chúng tôi để biết các bước chung.

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 JuliaGPU/oneAPI.jl

Tất cả issue của JuliaGPU/oneAPI.jl

Issue tương tự

Thêm issue về Julia

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.