Date: 2007-02-15 14:19:40
> Did you try to multiply a column major A' with a row major A? Then you
> can compute the columns of A'A one by one - each as a linear
> combination of a few columns of A'. (call axpy_prod(AT,x,y,false) with
> AT=a column major copy of A', x = row(A, i) )
I'm confused! A is mxn and A' is nxm. Multiplying A' by a row in A has a
mismatch in dimension: (nxm) * (nx1).