On 4/21/09, Mathew Yeates <[email protected]> wrote: > Hi > I posted something about this earlier > > Say I have 2 arrays X and Y with shapes (N,3) where N is large > I am doing the following > > for row in range(N): > result=polyfit(X[row,:],Y[row,:],1,full=True) # fit 3 points with a line > > This takes forever and I was hoping to find a way to speed things up. > But now I'm starting to wonder if this pointless. If the routine "poly > fit takes a long time, when compared with the time for a Python > function call, then things can't be sped up. > > Any comments?
If you remove the mean from x and y (along axis = 1) then can't you just do something like (x*y).sum(1) / (x*x).sum(axis=1) _______________________________________________ Numpy-discussion mailing list [email protected] http://mail.scipy.org/mailman/listinfo/numpy-discussion
