+1
Thanks

Accumulation of double precision vector

dattorro 5 years ago updated by Pavel Holoborodko 5 years ago 1

Quadruple precision floating-point accumulation of a double precision vector  p  is emulated numerically in Matlab:

function t = orosumvec(p, recurs)
   x = cumsum(p);
   z = x - p;
   u = (z - x) + p;
   v = [0
        x(1:end-1)] - z;
   if ~recurs
      t = sum([u
               v
               x(end)]);
   else
      t = orosumvec([u
                     v
                     x(end)], 0);
   end
end

This is verified by Advanpix Multiprecision Toolbox.


Matlab's built-in variable precision arithmetic, vpa() from Mathworks Symbolic Math Toolbox, 

produces erroneous results in 2018One recursive call is necessary and sufficient.

Thanks

Thank you very much for the report!

Could you please provide full example, so that other people can easily copy and reproduce the situation?

If you want, we can add this example to the VPA vs. MCT comparison here:

Symbolic Math Toolbox (VPA) vs. Multiprecision Computing Toolbox

We will also include backlink to your page, of course.