konrad.hinsen at laposte.net wrote: > How about > > A[:, :, NewAxis]*A[:, NewAxis, :] > > That works for your example at least. I am not quite sure why you use > matrixmultiply in your definition as there doesn't seem to be any > summation. You're right, it was my mistake to use matrixmultiply. Thanks a lot.