double-precision logs with error bounds
multiply by inverse