Hanni Ali wrote: > I'm afraid the matrix is not available anywhere and I would not be > able to make it available. >
Forget it, Sebastian is right. I was confused by the range of the error, but the error between two floating point numbers is indeed 1e-7 for float on most runtimes (FLT_EPS in C, which is the minimal value such as 1. + FLT_EPS != 1.) I am not sure whether we should define the accumulator to a double in the float case by default; generally, you use float for saving memory and computing speed, and you lose quite a bit speed-wise by using a double accumulator. FWIW, matlab behaves as well (not that it is a justification by itself, but at least it should not surprise matlab users). cheers, David _______________________________________________ Numpy-discussion mailing list Numpy-discussion@scipy.org http://projects.scipy.org/mailman/listinfo/numpy-discussion
