Question

设A是一个形状（a，b，c）的numpy数组，B是一个numpy形状的数组（a＆＃39;，b，c）。设f（A_，B_）是将形状（b，c）的numpy数组A_和形状（b，c）的numpy数组B_映射到实数的函数。我想构造一个形状（a，a＆＃39;）的numpy数组C，其中通过将f应用于第一个索引上的切片来给出条目。天真的解决方案是

    A=np.reshape(range(2*3*4), (2,3,4))
    B=np.reshape(range(3*3*4), (3,3,4))
    C=np.empty((2,3))

    def f(A_,B_):
      return np.prod(A_)+np.prod(B_)

    for i in range(A.shape[0]):
        for j in range(B.shape[0]):
           C[i,j]=f(A[i],B[j])

将C返回为

 [[  0.00000000e+00,   6.47647525e+14,   3.99703747e+17],
  [  6.47647525e+14,   1.29529505e+15,   4.00351395e+17]]

我要将它应用于更大的阵列A，B，其中f的计算成本很高（f以上只是一个玩具示例）。我通常尽量避免以元素方式访问numpy数组，但在上述情况下，我不知道如何实现这一点。

Answer 1

对于示例中的尺寸：

A2 = np.prod(A, axis=2).prod(axis=1)
B2 = np.prod(B, axis=2).prod(axis=1)
Bv, Av = np.meshgrid(B2, A2)
C2 = Av + Bv

array([[                 0,    647647525324800, 399703747322880000],
       [   647647525324800,   1295295050649600, 400351394848204800]])

元素函数在两个不同形状的numpy数组中的条目

1 个答案: