paweld/fpc
1
0
Fork 0
mirror of https://gitlab.com/freepascal.org/fpc/source.git synced 2025-04-21 21:51:42 +02:00
Code Issues Packages Projects Releases Wiki Activity

*** empty log message ***

Browse Source
This commit is contained in:
florian 2005-01-12 20:16:59 +00:00
parent abe11d29c2
commit a433f722be
1 changed files with 132 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

132
tests/test/units/math/ttrig1.pp Normal file
View File

@ -0,0 +1,132 @@
program test03;
uses
{$IFDEF CPU86}
CMem,
{$ENDIF}
SysUtils,Math;
const
dim = 36;
MaxFloat = 1.1000000000000000E+4932;
// from dcc
sinus : array[1..dim] of double = (
1.7364817766693035E-0001, 3.4202014332566873E-0001, 5.0000000000000000E-0001, 6.4278760968653933E-0001,
7.6604444311897804E-0001, 8.6602540378443865E-0001, 9.3969262078590838E-0001, 9.8480775301220806E-0001,
1.0000000000000000E+0000, 9.8480775301220806E-0001, 9.3969262078590838E-0001, 8.6602540378443865E-0001,
7.6604444311897804E-0001, 6.4278760968653933E-0001, 5.0000000000000000E-0001, 3.4202014332566873E-0001,
1.7364817766693035E-0001,-5.4210108624275222E-0020,-1.7364817766693035E-0001,-3.4202014332566873E-0001,
-5.0000000000000000E-0001,-6.4278760968653933E-0001,-7.6604444311897804E-0001,-8.6602540378443865E-0001,
-9.3969262078590838E-0001,-9.8480775301220806E-0001,-1.0000000000000000E+0000,-9.8480775301220806E-0001,
-9.3969262078590838E-0001,-8.6602540378443865E-0001,-7.6604444311897804E-0001,-6.4278760968653933E-0001,
-5.0000000000000000E-0001,-3.4202014332566873E-0001,-1.7364817766693035E-0001, 1.0842021724855044E-0019
);
var
i : integer;
Delta,Ref,Value : Double;
begin
for i:=1 to dim do
begin
// Generate constant array above output
write(sin(i*10/180.0*pi),',');
if i mod 4 = 0 then writeln;
end;
writeln('Testing SIN');
for i:=1 to dim do
begin
Ref:=sinus[i];
Value:=sin(i*10/180.0*pi);
Delta := Value - Ref;
if Abs(Delta) > 1E-15 then
begin
writeln(' Error for Sin(',i*10,') was:',Value,' should be:',Ref) ;
halt(1);
end;
end;
writeln('Testing COS');
for i:=1 to dim do
begin
Ref := sin(pi/2-i*10/180*pi);
Value := cos(i*10/180*pi);
Delta := Value - Ref;
if Abs(Delta) > 1E-15 then
begin
writeln(' Error for Cos(',i*10,') was:',Value,' should be:',Ref) ;
halt(1);
end;
end;
writeln('Testing TAN');
for i:=1 to dim do
begin
if i=9 then Ref := MaxFloat
else if i=27 then Ref := -Maxfloat
else Ref:=sin(i*10/180*pi)/cos(i*10/180*pi);
Value := tan(i*10/180*pi);
Delta := Value - Ref;
if Abs(Delta) > 1E-15 then
begin
writeln(' Error for Tan(',i*10,') was:',Value,' should be:',Ref) ;
halt(1);
end;
end;
writeln('Testing ARCTAN...');
for i:=1 to 8 do
begin
Ref := i*10;
Value := arctan(tan(i*10/180*pi))/pi*180;
Delta := Value - Ref;
if Abs(Delta) > 1E-14 then
begin
writeln(' Error for ArcTan(',i*10,') was:',Value,' should be:',Ref);
halt(1);
end;
end;
for i:=-1 downto -8 do
begin
Ref := i*10;
Value := arctan(tan(i*10/180*pi))/pi*180;
Delta := Value - Ref;
if Abs(Delta) > 1E-14 then
begin
writeln(' Error for ArcTan(',i*10,') was:',Value,' should be:',Ref);
halt(1);
end;
end;
writeln('Tan +/- 90 deg test:');
writeln('tan(89.999):',tan(89.999/180*pi));
writeln('tan(90.000):',tan(90.000/180*pi));
writeln('tan(-89.999):',tan(-89.999/180*pi));
writeln('tan(-90.000):',tan(-90.000/180*pi));
writeln('ArcTan2 kwadrants:');
writeln('Kwadrant 1 ( 1, 1):',arctan2( 1, 1)/pi*180:5:1,' deg');
writeln('Kwadrant 2 (-1, 1):',arctan2( 1,-1)/pi*180:5:1,' deg');
writeln('Kwadrant 3 (-1,-1):',arctan2(-1,-1)/pi*180:5:1,' deg');
writeln('Kwadrant 4 ( 1,-1):',arctan2(-1, 1)/pi*180:5:1,' deg');
writeln('ArcTan2 special cases:');
writeln('Kwadrant X ( 0, 0):',arctan2( 0, 0)/pi*180:5:1,' deg');
writeln('Kwadrant 1 ( 1, 0):',arctan2( 0, 1)/pi*180:5:1,' deg');
writeln('Kwadrant 2 ( 0, 1):',arctan2( 1, 0)/pi*180:5:1,' deg');
writeln('Kwadrant 3 (-1, 0):',arctan2( 0,-1)/pi*180:5:1,' deg');
writeln('Kwadrant 4 ( 0,-1):',arctan2(-1, 0)/pi*180:5:1,' deg');
end.