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
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
|
{
This file is part of the Free Pascal run time library.
Copyright (c) 1999-2001 by the Free Pascal development team
Implementation of mathematical routines (for extended type)
See the file COPYING.FPC, included in this distribution,
for details about the copyright.
This program is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
**********************************************************************}
{****************************************************************************
FPU Control word
****************************************************************************}
procedure Set8087CW(cw:word);
begin
{ pic-safe ; cw will not be a regvar because it's accessed from }
{ assembler }
default8087cw:=cw;
asm
fnclex
fldcw cw
end;
end;
function Get8087CW:word;assembler;
asm
pushl $0
fnstcw (%esp)
popl %eax
end;
procedure SetSSECSR(w : dword);
begin
mxcsr:=w;
asm
ldmxcsr w
end;
end;
function GetSSECSR : dword;
var
_w : dword;
begin
asm
stmxcsr _w
end;
result:=_w;
end;
{****************************************************************************
EXTENDED data type routines
****************************************************************************}
{$define FPC_SYSTEM_HAS_PI}
function fpc_pi_real : ValReal;compilerproc;
begin
{ Function is handled internal in the compiler }
runerror(207);
result:=0;
end;
{$define FPC_SYSTEM_HAS_ABS}
function fpc_abs_real(d : ValReal) : ValReal;compilerproc;
begin
{ Function is handled internal in the compiler }
runerror(207);
result:=0;
end;
{$define FPC_SYSTEM_HAS_SQR}
function fpc_sqr_real(d : ValReal) : ValReal;compilerproc;
begin
{ Function is handled internal in the compiler }
runerror(207);
result:=0;
end;
{$define FPC_SYSTEM_HAS_SQRT}
function fpc_sqrt_real(d : ValReal) : ValReal;compilerproc;
begin
{ Function is handled internal in the compiler }
runerror(207);
result:=0;
end;
{$define FPC_SYSTEM_HAS_ARCTAN}
function fpc_arctan_real(d : ValReal) : ValReal;compilerproc;
begin
{ Function is handled internal in the compiler }
runerror(207);
result:=0;
end;
{$define FPC_SYSTEM_HAS_LN}
function fpc_ln_real(d : ValReal) : ValReal;compilerproc;
begin
{ Function is handled internal in the compiler }
runerror(207);
result:=0;
end;
{$define FPC_SYSTEM_HAS_SIN}
function fpc_sin_real(d : ValReal) : ValReal;compilerproc;
begin
{ Function is handled internal in the compiler }
runerror(207);
result:=0;
end;
{$define FPC_SYSTEM_HAS_COS}
function fpc_cos_real(d : ValReal) : ValReal;compilerproc;
begin
{ Function is handled internal in the compiler }
runerror(207);
result:=0;
end;
{$define FPC_SYSTEM_HAS_EXP}
function fpc_exp_real(d : ValReal) : ValReal;assembler;compilerproc;
var
cw1,cw2: word;
asm
// comes from DJ GPP
fldt d
fldl2e
fmulp %st,%st(1)
fstcw CW1
fstcw CW2
fwait
andw $0xf3ff,CW2
orw $0x0400,CW2
fldcw CW2
fld %st(0)
frndint
fldcw CW1
fxch %st(1)
fsub %st(1),%st
f2xm1
fld1
faddp %st,%st(1)
fscale
fstp %st(1)
end;
{$define FPC_SYSTEM_HAS_FRAC}
function fpc_frac_real(d : ValReal) : ValReal;assembler;compilerproc;
asm
subl $4,%esp
fnstcw (%esp)
fwait
movw (%esp),%cx
orw $0x0f00,(%esp)
fldcw (%esp)
fldt d
frndint
fldt d
fsub %st(1),%st
fstp %st(1)
movw %cx,(%esp)
fldcw (%esp)
end;
{$define FPC_SYSTEM_HAS_INT}
function fpc_int_real(d : ValReal) : ValReal;assembler;compilerproc;
asm
subl $4,%esp
fnstcw (%esp)
fwait
movw (%esp),%cx
orw $0x0f00,(%esp)
fldcw (%esp)
fwait
fldt d
frndint
fwait
movw %cx,(%esp)
fldcw (%esp)
end;
{$define FPC_SYSTEM_HAS_TRUNC}
function fpc_trunc_real(d : ValReal) : int64;assembler;compilerproc;
asm
subl $12,%esp
fldt d
fnstcw (%esp)
movw (%esp),%cx
orw $0x0f00,(%esp)
fldcw (%esp)
movw %cx,(%esp)
fistpq 4(%esp)
fldcw (%esp)
fwait
movl 4(%esp),%eax
movl 8(%esp),%edx
end;
{$define FPC_SYSTEM_HAS_ROUND}
{ keep for bootstrapping with 2.0.x }
function fpc_round_real(d : ValReal) : int64;compilerproc;assembler;
var
res : int64;
asm
fldt d
fistpq res
fwait
movl res,%eax
movl res+4,%edx
end;
{$define FPC_SYSTEM_HAS_POWER}
function power(bas,expo : ValReal) : ValReal;
begin
if bas=0 then
begin
if expo<>0 then
power:=0.0
else
HandleError(207);
end
else if expo=0 then
power:=1
else
{ bas < 0 is not allowed when doing roots }
if (bas<0) and (frac(expo) <> 0) then
handleerror(207)
else
begin
power:=exp(ln(abs(bas))*expo);
if (bas < 0) and
odd(trunc(expo)) then
begin
power := -power;
end;
end;
end;
|
