在Julia中交换行的最简单方法

在Julia中,就像在Numpy中一样,使用方括号对数组进行索引,而在Matlab中,索引使用正则括号。

所以当你写作的时候 A([k r],:) ,它被解释为函数调用。将其更改为

A[[k r], :] = A[[r k], :]

作品但是用矢量索引比用矩阵索引更习惯,所以添加逗号,如下所示:

A[[k, r], :] = A[[r, k], :]

你说“最简单”,但如果你关心性能,那么下面的速度会快得多,而且不使用内存:

function swaprows!(A, k, r)
    for i in axes(A, 2)
        (A[k, i], A[r, i]) = (A[r, i], A[k, i])
    end
end