From: Barton W. <wil...@us...> - 2007-05-25 19:23:08
|
Update of /cvsroot/maxima/maxima/tests In directory sc8-pr-cvs16.sourceforge.net:/tmp/cvs-serv31741/tests Modified Files: rtestint.mac rtestode.mac rtestsum.mac Log Message: o make tests less dependent on ratsimp and evaluation o prepend kill(all) to rtestint o correct some float results in rtestsum (last few digits missing) Index: rtestint.mac =================================================================== RCS file: /cvsroot/maxima/maxima/tests/rtestint.mac,v retrieving revision 1.21 retrieving revision 1.22 diff -u -d -r1.21 -r1.22 --- rtestint.mac 20 May 2007 22:31:39 -0000 1.21 +++ rtestint.mac 25 May 2007 19:23:02 -0000 1.22 @@ -1,3 +1,6 @@ +(kill(all),0); +0$ + /* * This is a set of test integrals for testing various parts of the * integration routines. Index: rtestode.mac =================================================================== RCS file: /cvsroot/maxima/maxima/tests/rtestode.mac,v retrieving revision 1.14 retrieving revision 1.15 diff -u -d -r1.14 -r1.15 --- rtestode.mac 21 Nov 2006 05:16:56 -0000 1.14 +++ rtestode.mac 25 May 2007 19:23:02 -0000 1.15 @@ -12,7 +12,7 @@ ode2(x^2*'diff(y,x)+3*x*y=sin(x)/x, y, x); y = (%c-cos(x))/x^3; ic1(%, x=1, y=1); -y = -(cos(x)-cos(1)-1)/x^3; +y = -((cos(x)-cos(1)-1)/x^3); method; linear; @@ -23,7 +23,7 @@ ic2(soln, x=0, y=1, diff(y,x)=3); y = -sin(x)+cos(x)+4*x; bc2(soln, x=0, y=3, x=2, y=1); -y = -(3*cos(2)+7)*sin(x)/sin(2) + 3*cos(x) + 4*x; +y = -((3*cos(2)+7)*sin(x)/sin(2)) + 3*cos(x) + 4*x; ode2((3*x^2+4*x+2)=(2*y-1)*'diff(y,x), y, x); y^2-y = x^3+2*x^2+2*x+%c; @@ -43,7 +43,7 @@ exp(2*y)/y; ode2( 'diff(y,x)=(y/x)^2+2*(y/x), y, x); --(x*y+x^2)/y = %c; +-((x*y+x^2)/y) = %c; method; exact; @@ -114,7 +114,7 @@ eq: 'diff(y,x,2)+x*'diff(y,x)+exp(-x^2)*y=0; 'diff(y,x,2)+x*'diff(y,x,1)+%e^-x^2*y = 0; ans:ode2(eq,y,x); -y = %k1*sin(sqrt(2)*sqrt(%pi)*erf(x/sqrt(2))/2)+%k2*cos(sqrt(2)*sqrt(%pi)*erf(x/sqrt(2))/2); +y = %k1*sin((1/2) * sqrt(2)*sqrt(%pi)*erf(x/sqrt(2)))+%k2*cos((1/2) * sqrt(2)*sqrt(%pi)*erf(x/sqrt(2))); is(ratsimp(ev(eq,ans,diff))); true; method; @@ -156,7 +156,7 @@ eqn3: 'diff(f(x),x,2)+f(x)=2*x; 'diff(f(x),x,2)+f(x)=2*x; desolve(eqn3,f(x)); -f(x) = sin(x)*(at('diff(f(x),x,1),x = 0)-2)+f(0)*cos(x)+2*x; +''(f(x) = sin(x)*(at('diff(f(x),x,1),x = 0)-2)+f(0)*cos(x)+2*x); /* Examples mentioned in bug report [ 1063454 ] bug in ode2 * First one was reported to fail in CMUCL with "run out of heap" message. Index: rtestsum.mac =================================================================== RCS file: /cvsroot/maxima/maxima/tests/rtestsum.mac,v retrieving revision 1.11 retrieving revision 1.12 diff -u -d -r1.11 -r1.12 --- rtestsum.mac 8 May 2007 03:37:14 -0000 1.11 +++ rtestsum.mac 25 May 2007 19:23:02 -0000 1.12 @@ -49,7 +49,7 @@ /* further examples in the same vein: part of summand depends on index, part doesn't. */ sum (sin(g(j)) - %pi^j + cos(h(k)) - %e^l, j, 1, m), simpsum; -m*cos(h(k)) - m*%e^l + 'sum (sin(g(j)), j, 1, m) - (%pi^(m + 1) - %pi)/(%pi - 1); +m*cos(h(k)) - m*%e^l + 'sum (sin(g(j)), j, 1, m) - ((%pi^(m + 1) - %pi)/(%pi - 1)); sum (sum (f(j) + g(k), j, 1, m), k, 1, n), simpsum; m * 'sum (g(k), k, 1, n) + n * 'sum (f(j), j, 1, m); @@ -66,7 +66,7 @@ /* see SF bug report # 649428 */ (factsum3 (mt, ej) := sum ((-1)^(k + 1)/(k*mt^k)*sum((1 - l)^k - l^k, l, 1, ej),k,1,inf), t% : taylor (factsum3 (mt, ej), [mt, 0, 3, asymp])); -('sum(-2*l+1,l,1,ej))/mt-('sum(-2*l+1,l,1,ej))/(2*mt^2) +('sum(-2*l^3+3*l^2-3*l+1,l,1,ej)) /(3*mt^3); +sum(-2*l+1,l,1,ej)/mt-(sum(-2*l+1,l,1,ej)/(2*mt^2)) + sum(-2*l^3+3*l^2-3*l+1,l,1,ej) /(3*mt^3); ''t%, simpsum; -(6*ej^2*mt^2-3*ej^2*mt+ej^4+ej^2)/(6*mt^3); @@ -142,10 +142,10 @@ 1618091/196035840; %, numer; -0.00825405701325; +0.0082540570132481894; u, a=2, simpsum, numer; -0.00825405701325; +0.0082540570132481911; u, a=2, simpsum, bfloat; 8.25405701324819B-3; |