Compensated summation in sum?

Does Matlab use some algorithm such as compensated summation to reduce the roundoff error in sum?

Answer 1

[Jan]

There is no compensation in Matlab's SUM. For small arrays (there was a limit at 88999 elements, but this might change with the Matlab release) the sum is computed directly. For bigger arrays the sum is divided in parts and distributed over different threads.

If you need a compensated sum, try:

http://www.mathworks.com/matlabcentral/fileexchange/26800-xsum

This offers a Kahan-compensation, Knuth's method with intermediate 128 and 192 bit precision, an 80bit accumulator and Knuth with 160 bits (the last two are not supported by all compilers and platforms).

[Andrew Newell]

Thanks!