在Julia中评估具有矢量化函数的组合

Question

在Julia中，具有点.的矢量化函数用于按元素操作。

运行f.(x)意味着f(x[1])，f(x[2])，......依次执行

但是，假设我有一个函数，它接受两个参数，说g(x,y)

我想g(x[1],y[1])，g(x[2],y[1])，g(x[3],y[1])，... ，g(x[1],y[2]),g(x[2],y[2]),g(x[3],y[2])，...

有没有什么办法来评估x和的所有组合？

Answer 1

马特的答案是好的，但我想提供使用数组理解的选择：

julia> x = 1:5 
     y = 10:10:50 
     [i + j for i in x, j in y] 
5×5 Array{Int64,2}: 
11 21 31 41 51 
12 22 32 42 52 
13 23 33 43 53 
14 24 34 44 54 
15 25 35 45 55 

在我看来，阵列的理解往往更容易阅读，比broadcast和reshape更加灵活。

Answer 2

是的，重塑y，使它正交于x。 .矢量化使用broadcast来完成其工作。我把它想象成在所有其他维度上“挤压”单体尺寸。

这意味着，对于矢量x和y，你可以简单地重塑他们的一个评价x和y所有组合的产品：

julia> x = 1:5 
     y = 10:10:50 
     (+).(x, reshape(y, 1, length(y))) 
5×5 Array{Int64,2}: 
11 21 31 41 51 
12 22 32 42 52 
13 23 33 43 53 
14 24 34 44 54 
15 25 35 45 55 

请注意，数组的形状的方向相匹配争论; x跨越行，并且y跨越列，因为它被转置为单行矩阵。