From: Dominique O. <dom...@gm...> - 2009-11-17 21:27:38
|
On Tue, Nov 17, 2009 at 1:24 PM, Fergus Gallagher <fe...@go...> wrote: > On Tue, Nov 17, 2009 at 01:09:24PM -0500, Dominique Orban wrote: >> >> If you really need A*A^T explicitly, I suggest building A^T and then using >> dot(). >> > > I think we do. Our matrices are ~ 10 Million square, so A^T not > cheap doing it either explictly as A[i,j]=A[j,i] (ouch) or as > A^T = spmatrix.dot(A,I) OK. What I meant was to construct B=A^T to start with, instead of A. Then you can say dot(B,B) which will give you A*A^T. -- Dominique |