
Re: An independent integration test suite
Posted:
Jul 21, 2013 1:30 AM


On Friday, July 19, 2013 11:34:11 AM UTC4, clicl...@freenet.de wrote: > daly@axiomdeveloper.org schrieb: > > > > > > Albert, > > > > > > I'm getting different answers than you for some of the problems. > > > In Axiom, if I do > > > > > > t0:=1/((1x)^(9/2)*(1+x)^(5/2)) > > > t1:=integrate(t0,x) > > > t2:=D(t1,x) > > > t3:=t0t2 > > > > > > I get 0 as a result. However, when I take the derivative of your > > > result and difference it from your initial equation I get a > > > nonconstant result dependent on x. > > > > > > Perhaps you could check your answer in other systems and let me > > > know if you agree. > > > > > > I will be posting the first set of results shortly. > > > > > > > This integrand appears on page 1102 of rich1b.input.pdf. Its > > antiderivative > > > >  (8*x^5  16*x^4  4*x^3 + 24*x^2  9*x  6) > > /(21*(1  x)^(7/2)*(x + 1)^(3/2)) > > > > is simply an algebraic function, so should be behave the same on all > > systems. Indeed your expression d000 at the bottom of page 1103 > > simplifies to zero for arbitrary complex x. > > > > You are having more serious problems with the integrands on pages 994, > > 996, 998, 1000, 1004, 1008, 1017, 1019, 1020, 1021, 1022, 1107, 1109, > > and maybe elsewhere. Can FriCAS handle these (correctly) already? > > > > Martin.
I have completed the Rich's test set (a+b x)^m (A+B x) (d+e x)^p for Axiom. The results are at http://axiomdeveloper.org/axiomwebsite/CATS
Just a side note, the "source" format is a pamphlet file which is pure latex. Axiom uses these literate programming documents. The executable code lives in a "chunk" environment (similar to a "verbatim" environment). This file format is used so the expected results are cached in readable form.
To run these files (e.g for rich2.input.pamphlet), start Axiom and type )tangle rich2 )read rich2 The first command reads rich2.input.pamphlet and creates rich2.input The second command reads and executes rich2.input
You can also check the results for regression failures with the command )regress rich2 This command reads a file "rich2.output" (which is the spool file output requested by the )spool command in the rich2.input file).
Failing stanzas (not marked "ok to fail") will generate a "FAILED" line. A summary of FAILED lines is given at the end.
Some of the tests do not compare with the stored results (the R comments) because Axiom uses a random number generator and may get different constants. Only explicitly "FAILED" lines are interesting.
I'll post again as I generate more results.
Tim

