Я пытаюсь использовать функционал cumprod numpy, но для матриц. Имея матрицу, заполненную однородными преобразованиями (здесь случайными для простоты), я хочу накапливать преобразования, используя матричное умножение numpy.
n = 1000
data = np.random([n, 4,4])
dataAcc = cumProd(data, np.matmul)
dataA cc должен содержать следующее:
dataAcc[0,:,:] = data[0,:,:]
dataAcc[1,:,:] = data[0,:,:] @ data[1,:,:] = dataAcc[0,:,:] @ data[1,:,:]
dataAcc[2,:,:] = dataAcc[1,:,:] @ data[2,:,:]
...
Есть ли способ сделать это, используя быстрые numpy функции или аналогичные? Не хочу использовать петли ...
Ура