test suite failures with cvs Maxima



Hello Robert,

I have today updated my CVS. On my system (Windows XP with GCL 2.6.8) I get test
failures in rtest_hypgeo.lisp and rtesthyp.lisp.

First: The problems with rtest_hypgeo.test can almost be solved when indroducing
a dummy test at the place where the loading of orthopoly is killed. All tests
are well numbered. This scheme and the expected failures get out of order
killing the line. Should I commit the dummy test, to reorder the numbering
scheme?

Second: There remains one open failure. Pochhammer(1,i) no longer simplifies to
i!. After loading orthopoly by hand, it works again, the Pochhammer symbol now
simplifies.

Third:
I have the 3 failures (13,14,66) in rtesthyp.mac too. The reason is that e.g. in
the following example

   ratsimp(hgfred([-4],[1/2],z^2/2)-hermite(8,z/sqrt(2))/1680);

hgfred returns an expression in terms of the Hermite function. But now the
Hermite function we use to compare will automatically expanded in a series. 

So, I think all problems in GCL arise because of the changed orthopoly loading
scheme.

Dieter Kaiser

-----Urspr?ngliche Nachricht-----
Von: maxima-bounces at math.utexas.edu [mailto:maxima-bounces at math.utexas.edu] Im
Auftrag von Robert Dodier
Gesendet: Sonntag, 27. Juli 2008 03:50
An: Maxima mailing list
Betreff: [Maxima] test suite failures with cvs Maxima

Hello,

There are some failures in the test suite when running with
Clisp, CMUCL, SBCL, and GCL on Linux.
All 4 report failures in rtesthyp and CMUCL and Clisp report
some other errors as well. I've attached an abbreviated output
log from run_testsuite. Maybe someone can take a look at
these errors.

best

Robert Dodier

PS.
-----------------------------------------------------------------------
Using Lisp CMU Common Lisp 18e

Running tests in rtest_trig:
********************** Problem 60 ***************
Input:
block([buggy:[]],
      for f in alltrig do
          for p in pts do
              block([e,fop:op(f)],
                    e:buildq([p,fop],
                             lambda([],
                                    cabs(
                                     float(rectform(fop(p)))
                                      -fop(float(p))))),
                    if e() > 1.e-13 then buggy:cons(e,buggy)),
      buggy)


Result:
The number 0.5*%i+0.5 isn't in the domain of acot
?error\-catch

This differed from the expected result:
[]

Running tests in rtesthyp:
********************** Problem 13 ***************
Input:
ratsimp(hgfred([-4],[1/2],z^2/2)-hermite(8,z/sqrt(2))/1680)


Result:
-(16*z^8-448*z^6+3360*z^4-6720*z^2-hermite(8,z/sqrt(2))+1680)
 /1680

This differed from the expected result:
0

********************** Problem 14 ***************
Input:
ratsimp(hgfred([-4],[3/2],z^2/2)-sqrt(2)*hermite(9,z/sqrt(2))
                                 /30240
                                 /z)


Result:
-(16*sqrt(2)*z^9-576*sqrt(2)*z^7+6048*sqrt(2)*z^5
                -20160*sqrt(2)*z^3+15120*sqrt(2)*z
                -hermite(9,z/sqrt(2)))
 /(15120*sqrt(2)*z)

This differed from the expected result:
0

********************** Problem 66 ***************
Input:
ratsimp(hgfred([-2],[3/2],z)-hermite(5,sqrt(z))/(120*sqrt(z)))


Result:
-(sqrt(z)*(32*z^2-160*z+120)-hermite(5,sqrt(z)))/(120*sqrt(z))

This differed from the expected result:
0

Running tests in rtest8:
********************** Problem 98 ***************
Input:
ev(e4,foo(u) = u)


Result:
[2.0,0.0,25,0]

This differed from the expected result:
[1.999999999999999,2.220446049250313e-16,25,0]

********************** Problem 107 ***************
Input:
ev(e7,alfa = 2,vita = 1)


Result:
[.05296102778655729,4.85722573273506e-17,50,0]

This differed from the expected result:
[.05296102778655729,5.551115123125782e-17,50,0]

105/107 tests passed.

The following 2 problems failed: (98 107)
Running tests in rtest4:
********************** Problem 86 ***************
Input:
is(errcatch(rat(x^2^128)) = [])


Result:
true

This differed from the expected result:
false


-----------------------------------------------------------------------
Using Lisp SBCL 1.0

Running tests in rtesthyp:
********************** Problem 13 ***************
Input:
ratsimp(hgfred([-4],[1/2],z^2/2)-hermite(8,z/sqrt(2))/1680)


Result:
-(16*z^8-448*z^6+3360*z^4-6720*z^2-hermite(8,z/sqrt(2))+1680)
 /1680

This differed from the expected result:
0

********************** Problem 14 ***************
Input:
ratsimp(hgfred([-4],[3/2],z^2/2)-sqrt(2)*hermite(9,z/sqrt(2))
                                 /30240
                                 /z)


Result:
-(16*sqrt(2)*z^9-576*sqrt(2)*z^7+6048*sqrt(2)*z^5
                -20160*sqrt(2)*z^3+15120*sqrt(2)*z
                -hermite(9,z/sqrt(2)))
 /(15120*sqrt(2)*z)

This differed from the expected result:
0

********************** Problem 66 ***************
Input:
ratsimp(hgfred([-2],[3/2],z)-hermite(5,sqrt(z))/(120*sqrt(z)))


Result:
-(sqrt(z)*(32*z^2-160*z+120)-hermite(5,sqrt(z)))/(120*sqrt(z))

This differed from the expected result:
0


-----------------------------------------------------------------------
Using Lisp GNU Common Lisp (GCL) GCL 2.6.7 (aka GCL)

Running tests in rtesthyp:
********************** Problem 13 ***************
Input:
ratsimp(hgfred([-4],[1/2],z^2/2)-hermite(8,z/sqrt(2))/1680)


Result:
-(16*z^8-448*z^6+3360*z^4-6720*z^2-hermite(8,z/sqrt(2))+1680)
 /1680

This differed from the expected result:
0

********************** Problem 14 ***************
Input:
ratsimp(hgfred([-4],[3/2],z^2/2)-sqrt(2)*hermite(9,z/sqrt(2))
                                 /30240
                                 /z)


Result:
-(16*sqrt(2)*z^9-576*sqrt(2)*z^7+6048*sqrt(2)*z^5
                -20160*sqrt(2)*z^3+15120*sqrt(2)*z
                -hermite(9,z/sqrt(2)))
 /(15120*sqrt(2)*z)

This differed from the expected result:
0

********************** Problem 66 ***************
Input:
ratsimp(hgfred([-2],[3/2],z)-hermite(5,sqrt(z))/(120*sqrt(z)))


Result:
-(sqrt(z)*(32*z^2-160*z+120)-hermite(5,sqrt(z)))/(120*sqrt(z))

This differed from the expected result:
0


-----------------------------------------------------------------------
Using Lisp CLISP 2.38 (2006-01-24)

Running tests in rtest14:
********************** Problem 231 ***************
Input:
test_bessel(bessel_i(-2.5,-1.5),-.8015666610717216*%i,14)


Result:
0.234773439454848

This differed from the expected result:
true

********************** Problem 232 ***************
Input:
test_bessel(bessel_i(-2.5,1.5),.8015666610717216,14)


Result:
.2347734394548481

This differed from the expected result:
true

********************** Problem 242 ***************
Input:
test_bessel(bessel_i(2.5,-1.0),.05709890920304825*%i,15)


Result:
.01672388330721784

This differed from the expected result:
true

287/290 tests passed (not counting 2 expected errors).

The following 3 problems failed: (231 232 242)
Running tests in rtesthyp:
********************** Problem 13 ***************
Input:
ratsimp(hgfred([-4],[1/2],z^2/2)-hermite(8,z/sqrt(2))/1680)


Result:
-(16*z^8-448*z^6+3360*z^4-6720*z^2-hermite(8,z/sqrt(2))+1680)
 /1680

This differed from the expected result:
0

********************** Problem 14 ***************
Input:
ratsimp(hgfred([-4],[3/2],z^2/2)-sqrt(2)*hermite(9,z/sqrt(2))
                                 /30240
                                 /z)


Result:
-(16*sqrt(2)*z^9-576*sqrt(2)*z^7+6048*sqrt(2)*z^5
                -20160*sqrt(2)*z^3+15120*sqrt(2)*z
                -hermite(9,z/sqrt(2)))
 /(15120*sqrt(2)*z)

This differed from the expected result:
0

********************** Problem 66 ***************
Input:
ratsimp(hgfred([-2],[3/2],z)-hermite(5,sqrt(z))/(120*sqrt(z)))


Result:
-(sqrt(z)*(32*z^2-160*z+120)-hermite(5,sqrt(z)))/(120*sqrt(z))

This differed from the expected result:
0
_______________________________________________
Maxima mailing list
Maxima at math.utexas.edu
http://www.math.utexas.edu/mailman/listinfo/maxima