Dieter Kaiser wrote:
> Hello Ray,
>
> no problem. I am always a bit in fear that the code will not work on other
> systems. Perhaps in the feature I will install additionally a Linux to do even
> more tests with other compilers. But now my computer is too old and my harddisks
> are full.
I think you only need to have what you're comfortable with. :-)
Anyway, I've updated the code once again (but didn't get your
rtest_expintegral.mac change today). I removed the binary directories
and rebuilt from scratch. Now the errors are different. Clisp passes
all the tests, as does gcl. CMUCL and Ecl fail quite a few.
I won't list them all. I think if I can understand one of the failures,
the rest will be obvious. So, one cmucl failure is rtest_expintegral
problem 43:
********************** Problem 43 ***************
Input:
ev(test_value(%, expintegral_e(3, 0.5), 15), numer)
Result:
.3926990816987241
With expintexpand:erf, we have ratsimp(expintegral_e(3,1/2)) ->
(sqrt(%e)*(2*log(2)+log(-1/2)-2*expintegral_ei(-1/2)-log(-2))+4)/(16*sqrt(%e));
Numerically evaluating, we get .2216043642751784 - .3926990816987241 %i.
But expintegral_e(3,0.5) -> .2216043642751785.
More investigation needed....
Ray