[Maxima] Test suite failures with cvs version?
raymond.toy at ericsson.com
Fri Oct 24 13:13:21 CDT 2008
Dieter Kaiser wrote:
> Hello Ray,
> thank you very much for your help. I would have some problems to find so fast
> the reason for the problems.
I've been working on this off and on since yesterday. I'm not very fast
> Here is a version which work for GCL and CLISP and no longer needs to calculate
> the log function:
> (defun expintegral-ei (z)
> (- (expintegral-e 1 (- z)))
> ; (- (* 0.5 (- (log z) (log (/ 1 z)))) (log (- z)))
> (if (and (= (imagpart z) 0) (> (realpart z) 0))
> ;; Positive real value. Add a phase factor -%pi*%i
> (complex 0 (- (float pi)))
> ;; For all other values. No phase factor.
> With this change I have no problems with the testsuite.
> Could you verify, that this would work with other compilers?
This works well with cmucl and ecl, and the original issue is gone. But
there are issues. In rtest_expintegral, I get a failure in problem 149.
The result is
[628.3185307179587, [0.5, 1.140302841043172 - 12.56637061435917 %i,
1.140302841, 12.56637061435917], [0.49,
1.137167432518459 - 12.82282715750936 %i, 1.137167432, 12.82282715750936],
[0.48, 1.134046737613986 - 13.08996938995747 %i, 1.134046738,
13.08996938995747], [0.47, 1.130940670735224 - 13.36847937697784 %i,
1.130940671, 13.36847937697784], [0.46,
instead of . I haven't looked at this any further.
I'm still getting failures in problem 143, 147, and others where the
result is error-catch instead of . I haven't looked at these yet, but
these existed before.
Ecl works ok, but I also get error-catch issues on problems 153, 159,
More information about the Maxima