独立滚动矩阵行

Question

我有一个矩阵（准确地说是2d numpy ndarray）：

A = np.array([[4, 0, 0],
              [1, 2, 3],
              [0, 0, 5]])

我想A根据另一个数组中的roll值独立滚动每一行：

r = np.array([2, 0, -1])

也就是说，我要这样做：

print np.array([np.roll(row, x) for row,x in zip(A, r)])

[[0 0 4]
 [1 2 3]
 [0 5 0]]

有办法有效地做到这一点吗？也许使用花哨的索引技巧？

Answer 1

确保您可以使用高级索引来做到这一点，这是否是最快的方法可能取决于您的数组大小（如果行很大，则可能不是）：

rows, column_indices = np.ogrid[:A.shape[0], :A.shape[1]]

# Use always a negative shift, so that column_indices are valid.
# (could also use module operation)
r[r < 0] += A.shape[1]
column_indices = column_indices - r[:, np.newaxis]

result = A[rows, column_indices]