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