On Sun, Jul 10, 2011 at 8:26 PM, Roy Lowrance <[email protected]>wrote:
> I have a 1D float64 array ts. I want to square each element, so I compute > x = ts * ts > > I get a floating point overflow error. > > However, when I access each element separately and multiple, I get no > error: > for i in ts.shape[0]: > > Data please. The element by element squaring is handled by python, not numpy, so I expect numpy and python handle the errors differently. Catching floating point errors is a bit unreliable in any case. What OS/compiler are you using? Are you running 32 bit or 64 bit? Chuck
_______________________________________________ NumPy-Discussion mailing list [email protected] http://mail.scipy.org/mailman/listinfo/numpy-discussion
