The doubly compensated sum of a matrix.
s = accsum_dblcompsum ( x ) [s,e] = accsum_dblcompsum ( x )
a n-by-m matrix of doubles
a 1-by-1 matrix of doubles, the highest significant digits of the sum
a 1-by-1, matrix of doubles, the lowest significant digits of the sum
Returns the sum of x. Algorithm 4.2 in SNAA (Doubly Compensated Summation) Due to Priest (and Kahan) Assumes that Scilab uses round-to-nearest.
"Stability and numerical accuracy of algorithms", Nicolas Higham
"Handbook of Floating Point Computations", Muller et al