Julia 中矩阵 `按列` 存储, 以及 for 循环指标变化次序

Julia 的矩阵,例如 A[i,j,k] 是按照 column 存储的,最内层的指标 i 变化最快,所以如果要 for 循环矩阵,一定要使得 i 指标变化最快,才是最优的,例如 (这里补充下,对于 for a = A, b = B, ... end 这种循环,b 是变化最快的,即 b 是循环的最内层):

A = zeros(5,6,7);

I = 1:5; J = 1:6; K = 1:7;

for k = K, j = J, i = I
    A[i,j,k] = 0.1
end

更多的事例可以参考:Access arrays in memory order, along columns

你可能感兴趣的:(Julia)