在Julia中获得6维矩阵中子集的均值的最快方法

Question

我有一个Float64数组x，其中6个维度为[10,10,100,4,4,1000]。

我想取这个数组子集的均值，例如：

x = rand(10,10,100,4,4,1000)   #Dummy array
r = squeeze(mean(x[:,:,1:80,:,:,56:800],(1,2,3,4,5)),(1,2,3,4,5))

这需要很长时间才能完成，我假设由于矩阵索引中的复制过程。

有更快的方法吗？

编辑：我很高兴修改数组本身，并想知道是否有多维等效的pop!？