朱莉娅:阵列中的作业

时间:2015-04-20 19:48:53

标签: multidimensional-array variable-assignment julia

为数组索引多个级别时,它可以正常工作。但是当我用它来分配值时,却没有。有谁知道为什么A下面没有变化?

In  [4]: A = rand(6)

Out [4]: 6-element Array{Float64,1}:
 0.111552
 0.155126
 0.78485 
 0.147477
 0.362078
 0.959022

In  [5]: A[3:5][[true,false,true]]

Out [5]: 2-element Array{Float64,1}:
 0.78485 
 0.362078

In  [6]: A[3:5][[true,false,true]] = [99, 999]

Out [6]: 2-element Array{Int64,1}:
  99
 999

In  [7]: A

Out [7]: 6-element Array{Float64,1}:
 0.111552
 0.155126
 0.78485 
 0.147477
 0.362078
 0.959022

1 个答案:

答案 0 :(得分:5)

这是因为按范围和向量索引数组会返回带有输出的新数组(而不是原始数组中的视图)。您的陈述等同于以下内容:

julia> A = rand(6)
6-element Array{Float64,1}:
 0.806919
 0.445286
 0.882625
 0.556251
 0.719156
 0.276755

julia> B = A[3:5]
3-element Array{Float64,1}:
 0.882625
 0.556251
 0.719156

julia> B[[true,false,true]] = [99, 999]
2-element Array{Int64,1}:
  99
 999

julia> A'
1x6 Array{Float64,2}:
 0.806919  0.445286  0.882625  0.556251  0.719156  0.276755

julia> B'
1x3 Array{Float64,2}:
 99.0  0.556251  999.0

你实际上可以看到这就是Julia通过它的一些表达式实用程序所做的事情。注意明确的括号 - 它调用setindex!关于索引的结果,已经制作了副本。 (GenSym()是指定临时变量的内部方式):

julia> :(A[3:5][[true,false,true]] = [99, 999])
:((A[3:5])[[true,false,true]] = [99,999])

julia> expand(:(A[3:5][[true,false,true]] = [99, 999]))
:(begin
        GenSym(0) = (top(vect))(99,999)
        setindex!(getindex(A,colon(3,5)),GenSym(0),(top(vect))(true,false,true))
        return GenSym(0)
    end)

目标是最终让所有数组索引返回视图而不是副本,但这仍然是一项正在进行的工作。