The twosum sum of a and b.
s = accsum_twosum ( a , b ) [s,t] = accsum_twosum ( a , b )
a 1-by-1 matrix of doubles
a 1-by-1 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 a and b so that a + b = s + t exactly. where s contains the highest significant digits and t the lowest significant digits. Algorithm 4.4 in HFCC Due to Knuth No assumption on a , b Assumes that Scilab uses round-to-nearest.
"Stability and numerical accuracy of algorithms", Nicolas Higham
"Handbook of Floating Point Computations", Muller et al