A (quite) accurate ln1+x function, or "how close can you get" part II
02-01-2019, 04:26 PM
Post: #5
 Albert Chan Senior Member Posts: 684 Joined: Jul 2018
RE: A (quite) accurate ln1+x function, or "how close can you get" part II
(01-31-2019 07:04 PM)Albert Chan Wrote:  Excess ULP error is due to correction *lowering* decimal exponent.

To avoid excess ULP error, we like correction same sign as X

Y = 1+X, rounded-toward 1.0

log1p(X) ~ LN(Y) - (Y-1-X)/Y

Previous example, log1p(X = -0.00099950016) :
Y = round-toward-1 of 1+X = 0.9990004999 (10 digits)

log1p(X) ~ LN(Y) - (Y-1-X)/Y
= -9.999999333e-4 - 6.006003001e-11
= -9.999999934e-4 (all digits correct)
 « Next Oldest | Next Newest »

 Messages In This Thread A (quite) accurate ln1+x function, or "how close can you get" part II - Dieter - 04-09-2014, 06:44 PM RE: A (quite) accurate ln1+x function, or "how close can you get" part II - htom trites - 04-10-2014, 04:47 AM RE: A (quite) accurate ln1+x function, or "how close can you get" part II - Dieter - 04-11-2014, 07:01 PM RE: A (quite) accurate ln1+x function, or "how close can you get" part II - Albert Chan - 01-31-2019, 07:04 PM RE: A (quite) accurate ln1+x function, or "how close can you get" part II - Albert Chan - 02-01-2019 04:26 PM

User(s) browsing this thread: 1 Guest(s)