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