A Self Compensated Sum algorithm
s = accsum_scs ( x ) [s,e] = accsum_scs ( x )
a m-by-n, 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
A Self Compensated Sum algorithm. Uses accsum_fasttwosum.
"Stability and numerical accuracy of algorithms", Nicolas Higham
"Handbook of Floating Point Computations", Muller et al
https://hpcrd.lbl.gov/SCG/ocean/NRS/ECMWF/img14.htm
https://hpcrd.lbl.gov/SCG/ocean/NRS/SCSsum.F