log10(x) = log10e * log(x). The formula would work if the RHS were evaluated in extended precision with an extended precision log(). This actually happened with the i387 log() because it returns excess precision. Found by: ucbtest
Obtained from: netlib.att.com