📂ジュリア

Juliaで配列を列ごとに計算する方法

コード

配列を列ごと、または行ごとに計算する方法には、map、ブロードキャスト、コンプリヘンションがあります。

julia> using Statistics

julia> X = stack([i*ones(8) for i ∈ 1:9], dims=2)
8×9 Matrix{Float64}:
 1.0  2.0  3.0  4.0  5.0  6.0  7.0  8.0  9.0
 1.0  2.0  3.0  4.0  5.0  6.0  7.0  8.0  9.0
 1.0  2.0  3.0  4.0  5.0  6.0  7.0  8.0  9.0
 1.0  2.0  3.0  4.0  5.0  6.0  7.0  8.0  9.0
 1.0  2.0  3.0  4.0  5.0  6.0  7.0  8.0  9.0
 1.0  2.0  3.0  4.0  5.0  6.0  7.0  8.0  9.0
 1.0  2.0  3.0  4.0  5.0  6.0  7.0  8.0  9.0
 1.0  2.0  3.0  4.0  5.0  6.0  7.0  8.0  9.0

# map
julia> map(mean, eachcol(X))
9-element Vector{Float64}:
 1.0
 2.0
 3.0
 4.0
 5.0
 6.0
 7.0
 8.0
 9.0

julia> map(eachcol(X)) do col
           mean(col)
       end
9-element Vector{Float64}:
 1.0
 2.0
 3.0
 4.0
 5.0
 6.0
 7.0
 8.0
 9.0

# 브로드캐스트
julia> mean.(eachcol(X))
9-element Vector{Float64}:
 1.0
 2.0
 3.0
 4.0
 5.0
 6.0
 7.0
 8.0
 9.0

# 컴프리헨션
julia> [mean(col) for col ∈ eachcol(X)]
9-element Vector{Float64}:
 1.0
 2.0
 3.0
 4.0
 5.0
 6.0
 7.0
 8.0
 9.0

違い¹

通常、マップとコンプリヘンションの間にはパフォーマンスの違いはありません。二つの違いは返り値にあります。マップはイテレータをタプルで渡す場合にはタプルを返します。

julia> X = (1, 2, 3, 4, 5)
(1, 2, 3, 4, 5)

julia> [x^2 for x ∈ X]
5-element Vector{Int64}:
  1
  4
  9
 16
 25

julia> map(x->x^2, X)
(1, 4, 9, 16, 25)

コンプリヘンションでタプルを得るには、

julia> (x^2 for x ∈ X).iter
(1, 2, 3, 4, 5)

環境

• OS: Windows11
• Version: Julia 1.10.0