def euclidean_dist_matrix(data_1, data_2):
    """
    Returns matrix of pairwise, squared Euclidean distances
    """
    norms_1 = (data_1 ** 2).sum(axis=1)
    norms_2 = (data_2 ** 2).sum(axis=1)
    return np.abs(norms_1.reshape(-1, 1) + 2 * np.dot(data_1, data_2.T) - norms_2)
dataset/ETHPy150Open gmum/pykernels/pykernels/utils.py/euclidean_dist_matrix/swapped operands of `-`