On Wed, Jun 9, 2010 at 12:57 PM, "V. Armando Solé" <s...@esrf.fr> wrote: > greg whittier wrote: >> a = np.ones((400, 500000), dtype=np.float32) >> c = np.dot(a, a.T) >> >> > In such cases I create a matrix of zeros with the final size and I fill > it with a loop of dot products of smaller chunks of the original a matrix. > > The MDP package also does something similar. > > Armando >
Thanks. I've done that as a workaround, but I was hoping there was a more elegant, "native" solution. Thanks for the pointer to MDP by the way! Looks very interesting. _______________________________________________ NumPy-Discussion mailing list NumPy-Discussion@scipy.org http://mail.scipy.org/mailman/listinfo/numpy-discussion