Steven D'Aprano <st...@pearwood.info> writes: >>> protection against overflow of intermediate values. >> Might be simplest to just add the logarithms. > Simplest, but least accurate, even with Kahan summation or equivalent.
Well the idea was to avoid overflow first, then hold on to whatever precision you have after that. Is there a way to use extended (e.g. 80-bit) floats in numpy/scipy/etc.? -- https://mail.python.org/mailman/listinfo/python-list
