1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
|
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;
Ref2,Value2,Dummy : Float;
begin
for i:=1 to dim do
begin
// Generate constant array above output
write(sin(i*10/float(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/float(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/float(180)*pi);
Value := cos(i*10/float(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/float(180)*pi)/cos(i*10/float(180)*pi);
Value := tan(i*10/float(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/float(180)*pi))/pi*float(180);
Delta := Value - Ref;
if Abs(Delta) > {$ifdef cpuarm} 1E-13 {$else} 1E-14 {$endif} 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/float(180)*pi))/pi*float(180);
Delta := Value - Ref;
if Abs(Delta) > {$ifdef cpuarm} 1E-13 {$else} 1E-14 {$endif} then
begin
writeln(' Error for ArcTan(',i*10,') was:',Value,' should be:',Ref);
halt(1);
end;
end;
writeln('Testing SINCOS');
for i:=1 to dim do
begin
sincos(pi/2-i*10/float(180)*pi,Ref2,Dummy);
sincos(i*10/float(180)*pi,Dummy,Value2);
Delta := Value2 - Ref2;
if Abs(Delta) > 1E-15 then
begin
writeln(' Error for Cos(',i*10,') was:',Value2,' should be:',Ref2) ;
halt(1);
end;
end;
writeln('Tan +/- 90 deg test:');
writeln('tan(89.999):',tan(89.999/float(180)*pi));
writeln('tan(90.000):',tan(90.000/float(180)*pi));
writeln('tan(-89.999):',tan(-89.999/float(180)*pi));
writeln('tan(-90.000):',tan(-90.000/float(180)*pi));
writeln('ArcTan2 kwadrants:');
writeln('Kwadrant 1 ( 1, 1):',arctan2( 1, 1)/pi*float(180):5:1,' deg');
writeln('Kwadrant 2 (-1, 1):',arctan2( 1,-1)/pi*float(180):5:1,' deg');
writeln('Kwadrant 3 (-1,-1):',arctan2(-1,-1)/pi*float(180):5:1,' deg');
writeln('Kwadrant 4 ( 1,-1):',arctan2(-1, 1)/pi*float(180):5:1,' deg');
writeln('ArcTan2 special cases:');
writeln('Kwadrant X ( 0, 0):',arctan2( 0, 0)/pi*float(180):5:1,' deg');
writeln('Kwadrant 1 ( 1, 0):',arctan2( 0, 1)/pi*float(180):5:1,' deg');
writeln('Kwadrant 2 ( 0, 1):',arctan2( 1, 0)/pi*float(180):5:1,' deg');
writeln('Kwadrant 3 (-1, 0):',arctan2( 0,-1)/pi*float(180):5:1,' deg');
writeln('Kwadrant 4 ( 0,-1):',arctan2(-1, 0)/pi*float(180):5:1,' deg');
end.
|
