Wenn es nur um Summation geht, kannst Du meine aus meiner
AMath-Unit die Funktionen
Delphi-Quellcode:
function sum2x(const a: array of extended; n: integer): extended;
{-Compute accurate sum(a[i], i=0..n-1) of extended vector}
function sum2x(const a: array of extended; n: integer): extended;
{-Compute accurate sum(a[i], i=0..n-1) of extended vector}
benutzen, die Genaugkeit wird dadurch ca verdoppelt. Intern verwenden sie den Algorithmus von "T. Ogita, S.M. Rump, and S. Oishi, Accurate sum and dot product".