def _sq_dist(self, a, b=None):
        '''Compute a matrix of all pairwise squared distances
        between two sets of vectors, stored in the row of the two matrices:
        a (of size n by D) and b (of size m by D).'''
        n = a.shape[0]
        D = a.shape[1]
        m = n
        if b is None:
            b = a.transpose()
        else:
            m = b.shape[0]
            b = b.transpose()
        C = np.zeros((n,m))
        for d in range(0,D):
            tt  = a[:,d]
            tt  = tt.reshape(n,1)
            tem = np.kron(np.ones((1,m)), tt)
            tem = np.kron(np.ones((n,1)), b[d,:]) - tem
            C   = C + tem * tem
        return C
dataset/ETHPy150Open marionmari/pyGPs/pyGPs/Core/cov.py/Kernel._sq_dist/swapped operands of `-`