BTW, i have implemented sumOfLogs [1], it is more precise then everything else.
Thanks, I was going to use 32.64 fixed point for my program but now I think I will get better precision modifying that to use a 32 bit exponent field with a 1.64 fixed point field (to represent [1,2) normalized number.)
The [0.5,1) normalization due to frexp was confusing at first glance.
