python - numpy 点积的步骤

Question

我正在尝试将我的点积分解为步骤。就我而言，只需 2 个步骤。

>>> a=np.array([[1,2,3],[4,5,6],[7,8,9],[10,11,12],[13,14,15],[16,17,18]])
>>> a.dot(a.T)
array([[ 14,  32,  50,  68,  86, 104],
   [ 32,  77, 122, 167, 212, 257],
   [ 50, 122, 194, 266, 338, 410],
   [ 68, 167, 266, 365, 464, 563],
   [ 86, 212, 338, 464, 590, 716],
   [104, 257, 410, 563, 716, 869]])

我能够获得第一和第四象限，但不确定如何获得第二和第三象限

>>> a[0:3].dot(a[0:3].T)
array([[ 14,  32,  50],
       [ 32,  77, 122],
       [ 50, 122, 194]])


>>> a[3:].dot(a[3:].T)
array([[365, 464, 563],
       [464, 590, 716],
       [563, 716, 869]])

Answer 1

您想要做的是一个循环的 GEMM。以下是执行此操作的一些快速代码:

def loop_gemm(a, b, c=None, chunksize=100):

    size_i = a.shape[0]
    size_zip = a.shape[1]

    size_j = b.shape[1]
    size_alt_zip = b.shape[0]

    if size_zip != size_alt_zip:
        ValueError("Loop GEMM zip index is not of the same size for both tensors")

    if c is None:
        c = np.zeros((size_i, size_j))

    istart = 0
    for i in range(int(np.ceil(size_i / float(chunksize)))):

        left_slice = slice(istart, istart+chunksize)
        left_view = a[left_slice]  

        jstart = 0
        for j in range(int(np.ceil(size_j / float(chunksize)))):

            right_slice = slice(jstart, jstart+chunksize)
            right_view = b[:, right_slice]

            c[left_slice, right_slice] = np.dot(left_view, right_view)
            jstart += chunksize

        istart += chunksize

    return c

我们看到它在这里工作:

a = np.array([[1,2,3],[4,5,6],[7,8,9],[10,11,12],[13,14,15],[16,17,18]])
loop_gemm(a, a.T, chunksize=2)
array([[  14.,   32.,   50.,   68.,   86.,  104.],
       [  32.,   77.,  122.,  167.,  212.,  257.],
       [  50.,  122.,  194.,  266.,  338.,  410.],
       [  68.,  167.,  266.,  365.,  464.,  563.],
       [  86.,  212.,  338.,  464.,  590.,  716.],
       [ 104.,  257.,  410.,  563.,  716.,  869.]])

您还可以遍历邮政编码并缩放不同部分以保留 GEMM 的完整功能。

但是，除非您正在做一些特定的事情，否则这通常不是一个好方法:

a = np.random.rand(1E3, 1E4)
b = np.random.rand(1E4, 1E3)

%timeit np.dot(a,b)
10 loops, best of 3: 137 ms per loop

%timeit loop_gemm(a,b)
1 loops, best of 3: 311 ms per loop