[Numpy-discussion] Efficient orthogonalisation with scipy/numpy

josef.pktd at gmail.com josef.pktd at gmail.com
Tue Jan 19 16:19:19 EST 2010

On Tue, Jan 19, 2010 at 4:11 PM, Gael Varoquaux
<gael.varoquaux at normalesup.org> wrote:
>> you can  combine x, and y for one call to leastsq, if it makes a difference
>> linalg.lstsq(confounds.T, [x,y])  #format? columnstack?
> Indeed! Thank you Joseph. That's a gain of 10% in the total computation
> time of my algorithm (and 20% on the partial correlation calculation).

if you have z=[x,y] stacked, just one call to dot might also help for
the correlation

zz = dot(z.T, z)

You might be able to do everything in a stacked version/

there is no ph in Josef (unless you talk about my french father-in-law)


> Gaël
> _______________________________________________
> NumPy-Discussion mailing list
> NumPy-Discussion at scipy.org
> http://mail.scipy.org/mailman/listinfo/numpy-discussion

More information about the NumPy-Discussion mailing list