Hacktoberfest 2026:メンテナが10月に向けて印を付けた、オープンで初心者向けの issue。 Hacktoberfest の issue を見る

Wrong result from a reduction with a branching operator on Julia 1.12 with --check-bounds=yes

クローズ
#678 コメント 0 件 リアクション 0 件 担当者 0 名 GitHub で見る

メンテナーはふだん 1 日以内に返信

まだ誰も着手していません。

評価

難易度
5/5
見積もり時間
1週間以上
初心者へのやさしさ
15/100
issue の種類
バグ
明瞭さ
おおむね明確
活発さ
活発
領域
compilers

調査の方向性

Start with the reproducer in the issue, running it on Julia 1.12 with --check-bounds=yes, then compare against --check-bounds=auto. The failing test is the strided reductions testset in test/array.jl:223. The bug sits in the SPIR-V/LLVM path for a branching reducer over a named tuple on 2-D inputs, not in the test itself. Done means the reproducer returns true on 1.12 with bounds checking on, and the testset passes there.

索引モデルが issue の本文から書いたものです。

説明

With the GPUArrays 12 port (#665), == between two equal 64×64 oneArrays returns false. This only happens on Julia 1.12 and with bounds checking enabled, which Pkg.test does by default. It makes the strided reductions testset in test/array.jl:223 fail on 1.12. CI only tests 1.10 and 1.13, so it doesn't show up there.

GPUArrays' == reduces over a named tuple of two Bools (one tracking missing, one the result), and its reducer branches on the missing flag. Reduced:

using oneAPI
A = rand(Int32(1):Int32(4), 64, 64); dA = oneArray(A); dB = oneArray(A)
init = (; is_missing=false, is_equal=true)
mapper(a, b) = (; is_missing=false, is_equal=a == b)
branching(a, b) = (a.is_missing || b.is_missing) ? (; is_missing=true, is_equal=false) :
                                                  (; is_missing=false, is_equal=a.is_equal & b.is_equal)
branchless(a, b) = (; is_missing=a.is_missing | b.is_missing, is_equal=a.is_equal & b.is_equal)
@show mapreduce(mapper, branching, dA, dB; init)    # is_equal = false with --check-bounds=yes on 1.12
@show mapreduce(mapper, branchless, dA, dB; init)   # is_equal = true
@show mapreduce(mapper, branching, oneArray(ones(Int32, 4096)), oneArray(ones(Int32, 4096)); init)  # 1-D: true
--check-bounds=yes --check-bounds=auto
Julia 1.12.7 wrong (deterministic) correct
Julia 1.13.1 correct correct
Julia 1.10.12 correct correct
  • Only 2-D inputs fail. 1-D vectors of 1 to 4096 elements are correct.
  • The missing-tracking mapper doesn't matter. Replacing the branching reducer with a branchless one fixes it.
  • It doesn't depend on AcceleratedKernels (0.5.2, 0.5.3 and 0.5.4 all fail) or GPUCompiler (2.13.0, 2.13.1 and 2.13.2 all fail).
  • The same code on CUDA.jl with Julia 1.12 and --check-bounds=yes is correct.
  • oneAPI's main branch (GPUArrays 11) isn't affected, but it uses different reduction kernels.

The SPIR-V path for this kernel seems the likely culprit, on Julia 1.12's LLVM 18 or in IGC.

Environment: Intel Iris Xe Graphics (driver 1.3.38308), NEO 26.18.38308, libigc 2.34.4, SPIRV_LLVM_Backend 23.1.1, SPIRV_LLVM_Translator 23.1.1, Julia 1.12.7 (LLVM 18.1.7), oneAPI.jl from #665 (2.9.2 + GPUArrays 12.1.0, AcceleratedKernels 0.5.4, KernelAbstractions 0.9.44, GPUCompiler 2.13.2).

主要言語
Julia
スター
217
フォーク
37
平均マージ
1日 7時間
マージ済み PR(30日)
33

環境構築

このプロジェクトには開発コンテナ、Dockerfile、コントリビューションガイドがありません。まず README を読み、一般的な手順ははじめてのコントリビューションガイドを参照してください。

はじめの一歩

  1. issue を最後まで読み、次にプロジェクトのコントリビューションガイドを読みます。
  2. 着手することを issue にコメントします — 二人が同じ作業をするのを防げます。
  3. リポジトリをフォークし、ブランチを切って変更します。
  4. issue 番号を参照したプルリクエストを送ります。

JuliaGPU/oneAPI.jl のほかの issue

JuliaGPU/oneAPI.jl の issue をすべて見る

似ている issue

Julia の issue をもっと見る

新しい issue をメールで受け取る

初心者向けの GitHub issue を短くまとめたダイジェスト。