用数组索引分配矩阵值

Question

我想访问一个数组索引为 D 的numpy数组 M ，并使用另一个矩阵将其存储到同一数组 M 尺寸为 D x D 的strong> N 。

因此，我尝试了以下操作。我真的不想循环这样做。我来自Matlab，它与下面代码的第6行一样简单。

M= numpy.random.rand(6,6)
M= numpy.random.rand(6,6)  #M: 6x6
indices = [0,2,5]          #D=3
N = numpy.random.rand(3,3) #N: 3x3
#Don't work:
M[indices,indices]=M[indices,indices]+N 
M[tuple([indices,indices])]=M[tuple([indices,indices])]+N

通常 M 的格式如下：

  

array([r r ... r]) #in most attempts

因此，由于尺寸不同，口译员无法向 M 广播 N 。