Julia，@view 以防止生成新数组

Question

我的老朋友叫 Pop 的 3d 数组，我想删除列 (d2)，当 d1 末尾的值 == 1 时跨越所有 d3，所以我有这个代码：

Pop[end, :, 1] .!=1

我认为在此之前添加 @view 会将更改写回 Pop，而不是在内存中生成额外的副本。此代码工作正常

@view(Pop[ :, Pop[end, :, 1] .!=1, :])

但它不会改变名为 Pop 的原始 3d 数组。我可以

Pop = @view(Pop[ :, Pop[end, :, 1] .!=1, :])

但我相信这会在内存中创建另一个我试图避免的副本。我错过了什么简单的语法？谢谢。 J

Answer 1

@view 不会修改原始数组的大小，它提供了一个“视图”（例如，在您的情况下省略了一些列）。我觉得

没什么问题

Pop = view(Pop, :, Pop[end, :, 1] .≠ 1, :)

因为现在 Pop 是对旧的完整 Pop 的视图，但它的行为类似于一个数组，因此您可以修改它的条目，例如，您可以执行类似

julia> using Random # using a fixed seed for reproducibility

julia> Random.seed!(0) ;

julia> Pop = rand(1:5, (2,4,2)) # original Pop
2×4×2 Array{Int64,3}:
[:, :, 1] =
 4  3  5  5
 1  1  3  5

[:, :, 2] =
 2  2  3  1
 2  5  1  1

julia> Pop[end,:,1] .≠ 1 # columns to keep
4-element BitArray{1}:
 0
 0
 1
 1

julia> Pop = view(Pop, :, Pop[end, :, 1] .≠ 1, :) # make it a view
2×2×2 view(::Array{Int64,3}, :, [3, 4], :) with eltype Int64:
[:, :, 1] =
 5  5
 3  5

[:, :, 2] =
 3  1
 1  1

julia> Pop[end,:,1] .= 1 ; # use your new view to manipulate data

julia> Pop # view of the modified Pop
2×2×2 view(::Array{Int64,3}, :, [3, 4], :) with eltype Int64:
[:, :, 1] =
 5  5
 1  1

[:, :, 2] =
 3  1
 1  1

julia> Pop.parent # original full Pop (now Pop.parent) has been modified
2×4×2 Array{Int64,3}:
[:, :, 1] =
 4  3  5  5
 1  1  1  1

[:, :, 2] =
 2  2  3  1
 2  5  1  1

Julia，@view 以防止生成新数组

Julia, @view to prevent production of a new array

arrays

julia

row-removal