Re: Same .m file: different results with different versions of Octave

[Judd Storrs]

I'm not sure that's really a good fix for the problem. GSL's approach
is fundamentally superior because it avoids cancellation entirely by
factoring the results into a product. All of the additions in GSL's
implementation are between probably positive values that cannot have
cancellation problems. glibc's algorithim is comparatively naive and
should just be replaced but I haven't had a chance to verify the math
behind the GSL terms yet.


--judd