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

lu_instance/qr_instance run a full factorization for FixedSizeArrays

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

@devmotion đang làm issue này rồi.

Từ ngày 7/10/2026.

  • #511 của @devmotion — đang mở

Đánh giá

Độ khó
2/5
Thời gian dự kiến
1-3 giờ
Mức phù hợp với người mới
72/100
Loại issue
Tính năng
Độ rõ ràng
Đặc tả rõ ràng
Mức độ hoạt động
Sôi nổi
Công nghệ
julia
Lĩnh vực
performance

Hướng nghiên cứu

Đọc các tệp có sẵn ext/StaticArraysCoreExt.jl và ext/FillArraysExt.jl (hoặc các tệp tương đương trong ext/) làm mẫu, sau đó thêm một FixedSizeArraysExt định nghĩa các phương pháp rẻ lu_instance / qr_instance (và các *_instance liên quan) cho FixedSizeMatrix. Xác nhận rằng kiểu của mỗi instance khớp với những gì lu(a; check=false) / qr(a) thực sự trả về cho một FixedSizeMatrix — ví dụ pivot là một FixedSizeVector — bằng cách chạy các lệnh gọi đó trong REPL. Xong nghĩa là ArrayInterface.lu_instance(F) ở thang độ nanogiây và benchmark trong issue không còn cho thấy chênh lệch 50×; kiểm tra test/ để xem các phương pháp mở rộng được kiểm thử ở đâu.

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

Mô tả

lu_instance and qr_instance only have cheap methods for Matrix. A FixedSizeMatrix from FixedSizeArrays.jl falls back to lu(a, check = false) / qr(a) and runs a full factorization. LinearSolve's init calls these, so building a cache costs about 50× more than for a Matrix of the same size:

using ArrayInterface, FixedSizeArrays, LinearSolve, Chairmarks
A = randn(64, 64); F = FixedSizeMatrix(A)
@b ArrayInterface.lu_instance(A)                           # 24 ns
@b ArrayInterface.lu_instance(F)                           # 19 μs
@b LinearSolve.init(LinearProblem(similar(A), similar(A))) # 3.1 μs
@b LinearSolve.init(LinearProblem(similar(F), similar(F))) # 168 μs

Suggested fix: add an FixedSizeArraysExt (like the existing StaticArraysCore and FillArrays extensions) with cheap *_instance methods for FixedSizeMatrix. The instances have to match the types the real factorizations return; for example, lu on a FixedSizeMatrix returns a FixedSizeVector pivot vector.

ArrayInterface 7.30.2, FixedSizeArrays 1.3.0, LinearSolve 5.18.2, Julia 1.13.1. Found via ExponentialUtilities.jl, where the default exponential! is about 4× slower on small FixedSizeMatrix inputs than on Matrix.

Ngôn ngữ chính
Julia
Star
140
Fork
45
Merge trung bình
1 phút
Pull request đã merge (30 ngày)
1

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 JuliaArrays/ArrayInterface.jl

Tất cả issue của JuliaArrays/ArrayInterface.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.