Artifact 304ce03966ae5a8fff363c77d63c803e424f135ce975e7ab0b2bfa3ec0fa0df6:
- Executable file
r36/XMPL/MATH.TST
— part of check-in
[f2fda60abd]
at
2011-09-02 18:13:33
on branch master
— Some historical releases purely for archival purposes
git-svn-id: https://svn.code.sf.net/p/reduce-algebra/code/trunk/historical@1375 2bfe0521-f11c-4a00-b80e-6202646ff360 (user: arthurcnorman@users.sourceforge.net, size: 10922) [annotate] [blame] [check-ins using] [more...]
% Test of Math Package. % Author: Stanley L. Kameny <stan%valley.uucp@rand.org>. % Copyright (c) 1987, 1988, 1989, 1991 Stanley L. Kameny. % All Rights Reserved. %********************************************************************* %** %** This math package will compute the floating point values of ** %** the usual elementary functions, namely: ** %** sin asin sind asind sinh asinh ** %** cos acos cosd acosd cosh acosh ** %** tan atan tand atand tanh atanh ** %** cot acot cotd acotd coth acoth ** %** sec asec secd asecd sech asech ** %** csc acsc cscd acscd csch acsch ** %** atan2 atan2d ** %** exp ln sqrt ** %** expt log cbrt ** %** logb hypot ** %** log10 floor ** %** ceiling ** %** round ** %** ** %** All functions are computed to the accuracy of the floating- ** %** point precision of the system set up at the time. ** %** ** %********************************************************************* %** File #1===Trig Function Tests=== %** Trig functions are tested in both degrees and radians modes. %********************************************************************* symbolic; math!!label; symbolic procedure terpr(i,j); if remainder(i,j)=0 then terpri()$ on rounded; % We need !!plumin, etc. % #1: sind**2+cosd**2 test: ideal answers 1.0 for all args. for i:=0:45 do <<write " ",i,"->",sind float i**2+cosd float i**2; terpr(i,4)>>; % #2: quadrant test of sind, cosd: proper answers + +,+ -,- -,- +. begin scalar a; a:= sind 45.0; for i:= 0.0:3.0 do <<write " ",sind(i*90+45)/a," ", cosd (i*90+45)/a;terpri()>> end$ % #3: scaling test: all values should be 1 exactly. begin scalar a; a:= cosd 60.0; % for i:= -10.0:10.0 do write fix(cosd(60+i*360)/a)," " for i:= -10.0:10.0 do write round(cosd(60+i*360)/a)," " end$ % #4: test of radians -> degrees evaluation: ideal values 1.0. array a(6)$ begin for i:=1:6 do a(i):=sind(15.0*i); for i:=1:6 do <<write sin(!!pii2*i/6.0)/a(i)," "; terpr(i,3)>> end$ % #5: test of tand*cotd: ideal values 1.0. begin for i:=5 step 5 until 85 do <<write tand float i*cotd float i," "; terpr(i,25)>>; terpri() end$ % #6: test of secd**2-tand**2: ideal values 1.0. begin for i:=5 step 5 until 85 do <<write secd float i**2-tand float i**2," "; terpr(i,25)>> end$ % #7: test of cscd**2-cotd**2: ideal values 1.0. begin for i:=5 step 5 until 85 do <<write cscd float i**2-cotd float i**2," "; terpr(i,25)>> end$ % #8: test of asind+acosd: ideal values 1.0. begin write "sind and cosd"; terpri(); for i:=-10:10 do <<write (asind(0.1*i)+acosd(0.1*i))/90," "; terpr(i,5)>>; write "sin and cos";terpri(); for i:=-10:10 do <<write (acos(0.1*i)+asin(0.1*i))/!!pii2," "; terpr(i,5)>> end$ % #9: test of atand+acotd: ideal values 1.0. begin scalar x; write "tand, atand and acotd"; terpri(); for i:=-80 step 10 until 80 do <<x:=tand float i; write (atand x+acotd x)/90," "; terpr(i,50)>>; terpri(); write "tan, atan and acot";terpri(); for i:=-80 step 10 until 80 do <<x:= tan (!!pii2*i/90.0); write (atan x+acot x)/!!pii2," "; terpr(i,50)>> end$ % #10: test of atand tand: ideal values i for i:=-9:89. begin for i:=-9:89 do <<write " ",i,"->",if i=0 then 1.0 else atand tand float i; terpr(i,4)>> end$ % #11: test of acot cotd: ideal values 10*i for i:=1:17. begin for i:=10 step 10 until 170 do <<write " ",i,"->",acotd cotd i; terpr(i,40)>>; terpri();terpri() end$ % #12: test of asind sind: ideal values 10*i for i:=-9:9. begin for i:=-90 step 10 until 90 do <<write " ",i,"->",asind sind float i; terpr(i,40)>> end$ % #13: test of acosd cosd: ideal values 10*i for i:=1:18. begin for i:=10 step 10 until 180 do <<write " ",i,"->",acosd cosd float i; terpr(i,40)>> end$ % #14: test of acscd cscd: ideal values 10*i for i:=-9:9, except % error for i=0. begin for i:=-90 step 10 until 90 do <<write " ",i,"->",if i=0 then "error" else acscd cscd float i; terpr(i,40)>> end$ % #15: test of asecd secd: ideal values 10*i for i :=0:18. except % error for i=9. begin for i:=0 step 10 until 180 do <<write" ",i,"->",if i=90 then "error" else asecd secd float i; terpr(i,40)>> end$ %********************************************************************* %** ===Exp,Log,Sqrt,Cbrt, and Expt Function tests=== %********************************************************************* % #16: test of properties of exp function: ideal results 1.0. array b(5)$ begin scalar x; x:=0; write "multiplicative property";terpri(); for i:=0:5 do b(i):=1+i/6.0; for i:=0:5 do for j:=i:5 do <<write " ",exp (b(i)+b(j))/exp(b(i))/exp(b(j)); terpr(x:=x+1,5)>> end$ % #17: various properties of exp: ideal results 1.0. begin write "inverse property"$ terpri()$ for i:=1:5 do write " ",exp(b(i))*exp(-b(i));terpri(); write "squares"; terpri(); for i:=-10:10 do <<write " ",sqrt(exp(0.2*i))/exp(0.1*i); terpr(i,5)>>; write "cubes"; terpri(); for i:=-10:10 do <<write " ",cbrt(exp(0.3*i))/exp(0.1*i); terpr(i,5)>> end$ % #18: test of log exp: ideal results 1.0. begin for i:=-5:5 do <<write if i=0 then "0/0" else (log exp float i)/i," "; terpr(i,5)>> end$ % #19: test of log10 expt(10.0,i): ideal results 1.0. begin scalar i; write "small values i:=-5:5"; terpri(); for j:=-5:5 do <<write if j neq 0 then log10 float expt(10.0,j)/j else "zero"," "; terpr(j,5)>>; write "large i=2**j where j:=0:6"; terpri(); for j:=0:5 do <<write (log10 float expt(10.0,2**j))/2**j," "; terpr(j,5)>>; terpri(); write "noninteger values of i=j/10.0 where j:=1:20";terpri(); for j:=1:20 do <<i:=j/10.0; write (log10 float expt(10,i))/i," "; terpr(j,5)>> end$ % #20: test of properties of expt(x,i)*(expt(x,-i). ideal result 1.0. begin integer j; for x:=2:6 do for i:=2:6 do <<write expt(float x,i)*expt(float x,-i)," "; terpr(j:=j+1,5)>> end$ % #21: test of expt(-x,i)/expt(x,i) for fractional i. begin integer j,k; write "odd numerator. ideal result -1.0"; terpri(); for i:=1:10 do <<k:=(2*i-1.0)/(8*i+1); write rexpt(-8,k)/rexpt(8,k)," "; terpr(j:=j+1,5)>>; write "even numerator. ideal result 1.0"; terpri(); for i:=1:10 do <<k:=(2.0*i)/(8*i+1); write rexpt(-8,k)/rexpt(8,k)," "; terpr(j:=j+1,5)>> end$ % #22: test of properties of ln or log or logb: % inverse argument: ideal result -1.0. begin integer x; for i:=2:5 do for j:= 2:10 do <<x:=x+1; write logb(float i,float j)/logb(float i,1.0/j)," "; terpr(x,5)>> end$ % #23: test of log(a*b) = log a+log b: ideal result 1.0. begin integer x; for i:=1:5 do for j:=i:5 do <<write log (i*j*0.01)/(log (i*0.1)+log(j*0.1))," "; terpr(x:=x+1,5)>> end$ % #24:test of sqrt x*sqrt x/x for x:=5i*(5i/3)**i where i:=1:20 % (test values strictly arbitrary): ideal results 1.0. begin scalar x,s; for i:=1:20 do <<x:= 5.0*i;s:=sqrt(x:=x*(expt(x/3,i))); write s*s/x," "; terpr(i,5)>> end$ % #25: test of cbrt x**3/x for x:=5i*(5i/3)**i where i:=-9:10 % (test values strictly arbitrary):ideal results 1.0. begin scalar x,s; for i:=-9:10 do <<x:= 5.0*i; if i neq 0 then s:= cbrt(x:=x*(expt(x/3,i))); write if i=0 then 1 else s*s*s/x," "; terpr(i,5)>> end$ %********************************************************************* %** ===Hyperbolic Function Tests=== %********************************************************************* % #26: test of sinh x+ cosh x= exp x: ideal results 1.0. begin scalar x; for i:=1:10 do <<x:=ln float i$ write (sinh x+cosh x)/exp x," "$ terpr(i,5)>> end$ % #27: test of cosh x-sinh x= exp(-x): ideal results 1.0. begin scalar x; for i:=1:10 do <<x:=ln float i$ write(cosh x-sinh x)*exp x," "$ terpr(i,5)>> end$ % #28: test of (cosh x)**2-(sinh x)**2: ideal results 1.0. begin scalar x$ for i:=1:10 do <<x:=ln float i$write(cosh x)**2-(sinh x)**2," "; terpr(i,5)>> end$ % #29: test of tanh*cosh/sinh: ideal results 1.0. begin scalar x; for i:=1:20 do <<x:=ln(i*0.1); write if i=10 then 1 else tanh x*cosh x/sinh x," "; terpr(i,5)>> end$ % #30: test of tanh*coth: ideal results 1.0. begin scalar x; for i:=1:20 do <<x:=ln(i*0.1); write if i=10 then 1 else tanh x*coth x," "; terpr(i,5)>> end$ % #31: test of sech*cosh: ideal results 1.0. begin scalar x; for i:=1:20 do <<x:=ln(i*0.1); write sech x*cosh x," "; terpr(i,5)>> end$ % #32: test of csch*sinh: ideal results 1.0. begin scalar x; for i:=1:20 do <<x:=ln(i*0.1); write if i=10 then 1 else csch x*sinh x," "; terpr(i,5)>> end$ % #33: test of asinh sinh: ideal results 1.0. begin scalar x; for i:=1:20 do <<x:=ln(i*0.1); write if i=10 then 1 else asinh sinh x/x," "; terpr(i,5)>> end$ % #34: test of acosh cosh: ideal results 1.0. However, acosh x % loses accuracy as x -> 1 since d/dx cosh x -> 0. begin scalar x; for i:=1:20 do <<x:=ln(1+i*0.05); write acosh cosh x/x," "; terpr(i,5)>> end$ % #35: test of cosh acosh:ideal results 1.0. begin scalar x; for i:=1:50 do <<x:=1+i/25.0; write (cosh acosh x)/x," "; terpr(i,5)>> end$ % #36: test of atanh tanh: ideal results 1.0. begin scalar x; for i:=1:20 do <<x:=ln(i*0.1); write if i=10 then 1 else (atanh tanh x)/x," "; terpr(i,5)>> end$ % #37: test of acoth coth: ideal results 1.0. begin scalar x; for i:=1:20 do <<x:=ln(i*0.1); write if i=10 then 1 else (acoth coth x)/x," "; terpr(i,5)>> end$ % #38: test of asech sech: ideal results 1.0. However, asech x % loses accuracy as x -> 1 since d/dx sech x -> 0. begin scalar x; for i:=1:20 do <<x:=ln(1+i*0.05); write (asech sech x)/x," "; terpr(i,5)>> end$ % #39: test of acsch csch: ideal results 1.0. begin scalar x; for i:=1:20 do <<x:=ln(i*0.1); write if i=10 then 1 else (acsch csch x)/x," "; terpr(i,5)>> end$ % End of Test. end;