2010/8/3 matt_in_nyc <mkran...@gmail.com>: > X is an m by n matrix, and I want to store outer products of the form Y[i] = > numpy.outer(X[i,:], X[i,:]), leading to the relation Y[i,j,k] = > X[i,j]*X[i,k] for i = 0,...,m-1; j,k = 0,...,n-1. I am trying to think of > how to do this using tensordot, but so far I am finding no inspiration.
What about: X[:, :, numpy.newaxis] * X[:, numpy.newaxis, :] ? untested. Friedrich _______________________________________________ NumPy-Discussion mailing list NumPy-Discussion@scipy.org http://mail.scipy.org/mailman/listinfo/numpy-discussion