[Numpy-discussion] performance issue (again)

[Keith Goodman]

On Wed, Apr 22, 2009 at 8:48 AM, Mathew Yeates <myeates at jpl.nasa.gov> wrote:
> well, this isn't a perfect solution. polyfit  is better because it
> determines rank based on condition values. Finds the eigenvalues ...
> etc. But, unless it can vectorized without Python looping, it's too slow
> for me to use

I liked your "sheer genius" comment better.

Yeah, maybe use polyfit only for those cases where abs((x*y).sum(1) /
(x*x).sum(1)) is "large"? And ignore the slope calculation where
(x*x).sum(1) is "small".