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

perf(forest): thread the reference BlockForest path — leaf sweeps and exchange descriptor loops

Đang mở
#83 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
48/100
Loại issue
Tái cấu trúc
Độ 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, performance

Hướng nghiên cứu

Start with _forest_sweep_leaves! and _forest_adjoint_sweep_leaves! in src/operators/forest.jl, then inspect the three descriptor loops in _exchange_storage! in src/transfer.jl. Check whether descriptors in each exchange phase can alias destinations and compare the executor decision with #51. Run examples/niederer_benchmark.jl and verify diffusion scales with -t while the adjoint sweep preserves the stated inner-product equality bit-for-bit against the serial path.

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

Mô tả

performance

Carved out of #49 (candidate 1) so it can be tracked as one unit of work.

_forest_sweep_leaves! (src/operators/forest.jl), _forest_adjoint_sweep_leaves!, and the three descriptor loops in _exchange_storage! (src/transfer.jl) are plain serial for loops; on CPU the packed kernels fall back to the same broadcasts. Leaves are disjoint and the exchange phases are ordered to be independent within a phase, so this looks close to embarrassingly parallel. In the Niederer benchmark the reaction term is threaded in example code while diffusion is not, which is a large part of why the two phases scale so differently (#49: diffusion 0.44× on 5.3× fewer cells).

Check first

Whether descriptors within one exchange phase can alias the same destination block. If they can, the threaded sweep needs a partition by destination or atomics; if they cannot, a Threads.@threads over descriptors is enough.

Constraints

  • Same executor decision as #51 (KA CPU backend vs a threaded executor at the apply! seam); pick once for both single-grid and forest paths.
  • Must not disturb the _exchange_storage! / _bc_storage! rule seams: no allocation inside them, no Field-typed rule arguments.
  • Acceptance: the examples/niederer_benchmark.jl phase table shows the diffusion phase scaling with -t, and the adjoint sweep keeps ⟨Lx,y⟩ = ⟨x,Lᵀy⟩ bit-for-bit against the serial path.

🤖 Beep boop — Claude split this off #49; Kyle is not forking himself.

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

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 RallypointOne/MatrixFreeOperators.jl

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