ENH: Twice as fast by using np.einsum