I have adjusted the accuracy of some tests. Please report when there are more
problems.
Dieter Kaiser
-----Urspr?ngliche Nachricht-----
Von: maxima-bounces at math.utexas.edu [mailto:maxima-bounces at math.utexas.edu] Im
Auftrag von Dieter Kaiser
Gesendet: Montag, 5. Januar 2009 22:12
An: 'Raymond Toy'
Cc: maxima at math.utexas.edu
Betreff: Re: [Maxima] Testsuite failures
Hello Ray,
thank you for your report. I will have a look at the numerical accuracy of the
functions.
Dieter Kaiser
-----Urspr?ngliche Nachricht-----
Von: maxima-bounces at math.utexas.edu [mailto:maxima-bounces at math.utexas.edu] Im
Auftrag von Raymond Toy
Gesendet: Montag, 5. Januar 2009 17:40
An: Maxima List
Betreff: [Maxima] Testsuite failures
I'm getting the following failures in the testsuite:
********************** Problem 303 ***************
Input:
integrate(bessel_y(11, x), x)
Result:
error-catch
I think this due to a typo in the integral formula for bessel_y.
There's a call to simplify with two args instead of one. I think the
second arg should be removed, or simplify should by simplifya.
Clisp (2.47) fails rtest_gamma:
********************** Problem 354 ***************
Input:
relerror(gamma_incomplete(0.5, 0.1 %i - 2.5),
1.00288947695445 - 9.842709236689227 %i,
1.0E-15)
Result:
1.105194365996218E-15
Cmucl fails rtest_gamma:
********************** Problem 371 ***************
Input:
relerror(gamma_incomplete(0.5, - 5 %i),
0.364419841063559 - .2436855906381129 %i,
1.0e-15)
Result:
2.788636698038084e-15
********************** Problem 381 ***************
Input:
relerror(gamma_incomplete(0.5, 5 %i), .2436855906381129 %i
+ 0.364419841063559,
1.0e-15)
Result:
2.788636698038084e-15
********************** Problem 388 ***************
Input:
relerror(gamma_incomplete(0.5, 2.0 %i - 2),
- .6696404542328767 %i - 4.636214977691262,
2.0e-15)
Result:
2.051455149042451e-15
********************** Problem 389 ***************
Input:
relerror(gamma_incomplete(0.5, 3.0 %i - 3),
8.244293668813985 %i - 6.274836016533672,
2.5e-15)
Result:
3.682253398201579e-15
********************** Problem 391 ***************
Input:
relerror(gamma_incomplete(0.5, 5.0 %i - 5),
9.813229347687738 %i + 57.54071128202554,
2.0e-15)
Result:
2.449911020455208e-15
********************** Problem 396 ***************
Input:
relerror(gamma_incomplete(0.5, 2.0 %i - 2),
- .6696404542328767 %i - 4.636214977691262,
2.0e-15)
Result:
2.051455149042451e-15
********************** Problem 397 ***************
Input:
relerror(gamma_incomplete(0.5, 3.0 %i - 3),
8.244293668813985 %i - 6.274836016533672,
3.0e-15)
Result:
3.682253398201579e-15
********************** Problem 399 ***************
Input:
relerror(gamma_incomplete(0.5, 5.0 %i - 5),
9.813229347687738 %i + 57.54071128202554,
2.0e-15)
Result:
2.449911020455208e-15
Gcl fails rtest_gamma:
********************** Problem 592 ***************
Input:
relerror(fresnel_s(%i + 5.0),
3.84389093777759513198736581956757227475265101347572827364567521b4 %i
- 2.044525505063210590493330126425293361797143144299005524393297869b5,
2.e-14)
Result:
2.542490092925679718842752779534059543174617462409733517737362792b-14
********************** Problem 614 ***************
Input:
relerror(fresnel_c(%i + 5.0),
2.044530505063119004499308897858846121088557663560705550579996093b5 %i
+ 3.84394093777740143202918713550472184235160647415045418329908291b4,
2.e-14)
I guess the tolerance needs to be increased a bit. The CMUCL failures
may indicate an issue with some CMUCL routines since CMUCL fails more
tests than the other lisps.
Ray
_______________________________________________
Maxima mailing list
Maxima at math.utexas.edu
http://www.math.utexas.edu/mailman/listinfo/maxima
_______________________________________________
Maxima mailing list
Maxima at math.utexas.edu
http://www.math.utexas.edu/mailman/listinfo/maxima