访问矩阵时SPLAT操作员的性能

Question

假设我们有一个矩阵和要访问该矩阵的指数的矩阵和向量：

matrix = Matrix{Float64}(undef, 5000, 4000)
point = [1244, 3353]

我们可以使用Splat Operator（...）扩展点向量以访问矩阵，或者我们可以明确使用点[1 ]和Point [2]作为指数：

matrix[point...]
vs
matrix[point[1], point[2]]

第一个更优雅，但显然也更慢：

@btime $matrix[$point...]
70.083 ns (4 allocations: 64 bytes)

@btime $matrix[$point[1], $point[2]]
1.800 ns (0 allocations: 0 bytes)

这种差异来自何处？有没有办法使SPLAT操作员更具性能，或者如果我关心性能，我应该简单地使用其他解决方案？

Answer 1

元组似乎与明确的破坏一样快：（

julia> @btime $matrix[point...] setup=(point = [1244, 3353])
  91.024 ns (4 allocations: 64 bytes)
0.0

julia> @btime $matrix[point[1], point[2]] setup=(point = [1244, 3353])
  2.530 ns (0 allocations: 0 bytes)
0.0

julia> @btime $matrix[point...] setup=(point = (1244, 3353))
  2.624 ns (0 allocations: 0 bytes)
0.0

julia> @btime $matrix[point[1], point[2]] setup=(point = (1244, 3353))
  2.255 ns (0 allocations: 0 bytes)
0.0

我认为后三个之间的区别大多是微不足道的 - 在较早的试验中，第二个“赢了”。

搭配 如果您直接插入元组，就不足为奇了。