Jeroen Demeyer on Tue, 12 Nov 2013 16:43:14 +0100


[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]

Re: Large errors for eint()


On 2013-11-09 23:36, Karim Belabas wrote:
1) For eint1, it was an oversight (1-word cancellation just before
returning the result, plus a few other places where we computed with one
word less than intended). This one is now fixed in commit f0ed0cd2
Shouldn't the result of eint1() be truncated to the original precision?

gp> x=15.0;
gp> \x
[&=0000000001ae3d90] REAL(lg=4,CLONE):0500000000000004 (+,expo=3):6000000000000003 f000000000000000 0000000000000000

gp> eint1(x);
gp> \x
[&=00000000019eefb0] REAL(lg=5,CLONE):0500000000000005 (+,expo=-26):5fffffffffffffe6 a4cf12d33ea678a5 77d1130620dd3589 3573133270000000