Extremely slow BLAS muladd when α,β are Int
メンテナーはふだん 1 日以内に返信
まだ誰も着手していません。
評価
- 難易度
- 4/5
- 見積もり時間
- 3〜5日
- 初心者へのやさしさ
- 38/100
- issue の種類
- バグ
- 明瞭さ
- 明確に書かれている
- 活発さ
- 停滞
- 技術スタック
- julia
- 領域
- performance
調査の方向性
まず、issue にある再現可能なベンチマークを実行し、default_blasmul!、tiled_blasmul!、および MulAdd wrapper を調べます。整数の α と β に対する5引数の LinearAlgebra.mul! と、これらの挙動および実行時間を比較します。報告された低速化が説明され、影響を受ける操作が同等のパフォーマンスを持つか、文書化された制限があれば完了です。
索引モデルが issue の本文から書いたものです。
説明
Hi, perhaps I do not understand something about how ArrayLayouts are supposed to work and be used, but it seems that BLAS matrix multiplication (both default_blasmul! and tiled_blasmul!) is very slow in comparison to 5-argument LinearAlgebra.mul!.
Below is a reproducible example.
using LinearAlgebra
import ArrayLayouts: default_blasmul!, tiled_blasmul!, tile_size
A, B, C = (rand(1000, 1000) for _ in 1:3)
# mul! is fast
mul!(copy(C), A, B, 1, 1) ≈ A*B + C # true
@btime mul!(copy(C), A, B, 1, 1);
# 6.448 ms (4 allocations: 7.63 MiB)
# default_blasmul! is super slow!!
default_blasmul!(1, A, B, 1, copy(C)) ≈ A*B + C # true
@btime default_blasmul!(1, A, B, 1, copy(C));
# 852.139 ms (3 allocations: 7.63 MiB)
# tiled_blasmul! is also very slow!
ts = tile_size(eltype(A), eltype(B), eltype(C))
tiled_blasmul!(ts, 1, A, B, 1, copy(C)) ≈ A*B + C # true
@btime tiled_blasmul!(ts, 1, A, B, 1, copy(C));
# 402.555 ms (12 allocations: 7.66 MiB)
And this of course affects operations using MulAdd wrapper. What am I doing wrong? Or is it really a bug?
Version info
Julia Version 1.11.3
Commit d63adeda50d (2025-01-21 19:42 UTC)
Build Info:
Official https://julialang.org/ release
Platform Info:
OS: Linux (x86_64-linux-gnu)
CPU: 16 × AMD Ryzen 9 5900HX with Radeon Graphics
WORD_SIZE: 64
LLVM: libLLVM-16.0.6 (ORCJIT, znver3)
Threads: 1 default, 0 interactive, 1 GC (on 16 virtual cores)
Environment:
JULIA_PROJECT = @.
JULIA_PKG_PRESERVE_TIERED_INSTALLED = true
JULIA_REVISE = manual
PKG status
[dce04be8] ArgCheck v2.4.0
[7d9fca2a] Arpack v0.5.4
[4c555306] ArrayLayouts v1.11.0
[86223c79] Graphs v1.12.0
⌃ [0b1a1467] KrylovKit v0.9.3
⌃ [5078a376] LazyArrays v2.4.0
[2ab3a3ac] LogExpFunctions v0.3.29
[5a0628fe] Supposition v0.3.5
[37e2e46d] LinearAlgebra v1.11.0
[2f01184e] SparseArrays v1.11.0
[f489334b] StyledStrings v1.11.0
[8dfed614] Test v1.11.0
- 主要言語
- Julia
- スター
- 57
- フォーク
- 10
- 平均マージ
- 1日 9時間
- マージ済み PR（30日）
- 11
環境構築
このプロジェクトには開発コンテナ、Dockerfile、コントリビューションガイドがありません。まず README を読み、一般的な手順ははじめてのコントリビューションガイドを参照してください。
はじめの一歩
- issue を最後まで読み、次にプロジェクトのコントリビューションガイドを読みます。
- 着手することを issue にコメントします — 二人が同じ作業をするのを防げます。
- リポジトリをフォークし、ブランチを切って変更します。
- issue 番号を参照したプルリクエストを送ります。
JuliaLinearAlgebra/ArrayLayouts.jl のほかの issue
-
難易度 4/5 3〜5日 初心者へのやさしさ 45/100
JuliaLinearAlgebra/ArrayLayouts.jl#277 ·
メンテナーはふだん 1 日以内に返信
-
難易度 3/5 1〜2日 初心者へのやさしさ 48/100
JuliaLinearAlgebra/ArrayLayouts.jl#243 ·
メンテナーはふだん 1 日以内に返信
-
難易度 5/5 1週間以上 初心者へのやさしさ 15/100
JuliaLinearAlgebra/ArrayLayouts.jl#185 ·
メンテナーはふだん 1 日以内に返信
-
難易度 2/5 1〜3時間 初心者へのやさしさ 35/100
JuliaLinearAlgebra/ArrayLayouts.jl#157 · コメント 2 件 ·
メンテナーはふだん 1 日以内に返信
-
難易度 5/5 1週間以上 初心者へのやさしさ 28/100
JuliaLinearAlgebra/ArrayLayouts.jl#155 · コメント 1 件 ·
メンテナーはふだん 1 日以内に返信
JuliaLinearAlgebra/ArrayLayouts.jl の issue をすべて見る
似ている issue
-
bug
難易度 2/5 1〜3時間 初心者へのやさしさ 61/100
-
難易度 2/5 1〜3時間 初心者へのやさしさ 83/100
-
bug
難易度 2/5 1〜3時間 初心者へのやさしさ 80/100
Deltares/Wflow.jl#1145 · コメント 1 件 ·
メンテナーはふだん 1 日以内に返信
-
難易度 2/5 1〜3時間 初心者へのやさしさ 72/100
JuliaDiff/SparseMatrixColorings.jl#356 ·
メンテナーはふだん 1 日以内に返信
-
難易度 2/5 1〜3時間 初心者へのやさしさ 72/100
JuliaDecisionFocusedLearning/ImplicitDifferentiation.jl#231 ·