351
|
1 |
(*
|
|
2 |
* Hedgewars, a Worms style game
|
|
3 |
* Copyright (c) 2007 Andrey Korotaev <unC0Rr@gmail.com>
|
|
4 |
*
|
|
5 |
* This program is free software; you can redistribute it and/or modify
|
|
6 |
* it under the terms of the GNU General Public License as published by
|
|
7 |
* the Free Software Foundation; version 2 of the License
|
|
8 |
*
|
|
9 |
* This program is distributed in the hope that it will be useful,
|
|
10 |
* but WITHOUT ANY WARRANTY; without even the implied warranty of
|
|
11 |
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
|
|
12 |
* GNU General Public License for more details.
|
|
13 |
*
|
|
14 |
* You should have received a copy of the GNU General Public License
|
|
15 |
* along with this program; if not, write to the Free Software
|
|
16 |
* Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA
|
|
17 |
*)
|
|
18 |
|
|
19 |
unit uFloat;
|
|
20 |
interface
|
|
21 |
|
|
22 |
{$IFDEF FPC}
|
|
23 |
{$ifdef FPC_LITTLE_ENDIAN}
|
|
24 |
type hwFloat = record
|
|
25 |
isNegative: boolean;
|
|
26 |
case byte of
|
|
27 |
0: (Frac, Round: Longword);
|
|
28 |
1: (QWordValue : QWord);
|
|
29 |
end;
|
|
30 |
{$else FPC_LITTLE_ENDIAN}
|
|
31 |
type hwFloat = record
|
|
32 |
isNegative: boolean;
|
|
33 |
case byte of
|
|
34 |
0: (Round, Frac: Longword);
|
|
35 |
1: (QWordValue : QWord);
|
|
36 |
end;
|
|
37 |
{$endif FPC_LITTLE_ENDIAN}
|
|
38 |
|
498
|
39 |
function int2hwFloat (i: LongInt) : hwFloat;
|
351
|
40 |
|
|
41 |
operator + (z1, z2: hwFloat) z : hwFloat;
|
|
42 |
operator - (z1, z2: hwFloat) z : hwFloat;
|
|
43 |
operator - (z1: hwFloat) z : hwFloat;
|
|
44 |
|
|
45 |
operator * (z1, z2: hwFloat) z : hwFloat;
|
|
46 |
operator * (z1: hwFloat; z2: LongInt) z : hwFloat;
|
|
47 |
operator / (z1, z2: hwFloat) z : hwFloat;
|
498
|
48 |
operator / (z1: hwFloat; z2: LongInt) z : hwFloat;
|
351
|
49 |
|
|
50 |
operator < (z1, z2: hwFloat) b : boolean;
|
|
51 |
operator > (z1, z2: hwFloat) b : boolean;
|
|
52 |
|
|
53 |
function cstr(z: hwFloat): string;
|
371
|
54 |
function hwRound(t: hwFloat): LongInt;
|
351
|
55 |
function hwAbs(t: hwFloat): hwFloat;
|
|
56 |
function hwSqr(t: hwFloat): hwFloat;
|
357
|
57 |
function hwSqrt(t: hwFloat): hwFloat;
|
351
|
58 |
function Distance(dx, dy: hwFloat): hwFloat;
|
498
|
59 |
function DistanceI(dx, dy: LongInt): hwFloat;
|
357
|
60 |
function AngleSin(Angle: Longword): hwFloat;
|
|
61 |
function AngleCos(Angle: Longword): hwFloat;
|
498
|
62 |
function SignAs(num, signum: hwFloat): hwFloat;
|
351
|
63 |
|
|
64 |
const _1div1024: hwFloat = (isNegative: false; QWordValue: 4194304);
|
|
65 |
_1div10000: hwFloat = (isNegative: false; QWordValue: 429496);
|
|
66 |
_1div50000: hwFloat = (isNegative: false; QWordValue: 85899);
|
|
67 |
_1div100000: hwFloat = (isNegative: false; QWordValue: 42950);
|
|
68 |
_1div3: hwFloat = (isNegative: false; QWordValue: 1431655766);
|
|
69 |
hwPi: hwFloat = (isNegative: false; QWordValue: 13493037704);
|
|
70 |
_0_000004: hwFloat = (isNegative: false; QWordValue: 17179);
|
|
71 |
_0_0002: hwFloat = (isNegative: false; QWordValue: 858993);
|
|
72 |
_0_001: hwFloat = (isNegative: false; QWordValue: 4294967);
|
|
73 |
_0_003: hwFloat = (isNegative: false; QWordValue: 12884902);
|
|
74 |
_0_004: hwFloat = (isNegative: false; QWordValue: 17179869);
|
|
75 |
_0_005: hwFloat = (isNegative: false; QWordValue: 21474836);
|
|
76 |
_0_01: hwFloat = (isNegative: false; QWordValue: 42949673);
|
|
77 |
_0_02: hwFloat = (isNegative: false; QWordValue: 85899345);
|
|
78 |
_0_03: hwFloat = (isNegative: false; QWordValue: 128849018);
|
|
79 |
_0_08: hwFloat = (isNegative: false; QWordValue: 343597383);
|
365
|
80 |
_0_1: hwFloat = (isNegative: false; QWordValue: 429496730);
|
351
|
81 |
_0_15: hwFloat = (isNegative: false; QWordValue: 644245094);
|
|
82 |
_0_2: hwFloat = (isNegative: false; QWordValue: 858993459);
|
|
83 |
_0_25: hwFloat = (isNegative: false; QWordValue: 1073741824);
|
|
84 |
_0_3: hwFloat = (isNegative: false; QWordValue: 1288490189);
|
|
85 |
_0_35: hwFloat = (isNegative: false; QWordValue: 1503238553);
|
|
86 |
_0_4: hwFloat = (isNegative: false; QWordValue: 1717986918);
|
|
87 |
_0_45: hwFloat = (isNegative: false; QWordValue: 1932735283);
|
|
88 |
_0_5: hwFloat = (isNegative: false; QWordValue: 2147483648);
|
|
89 |
_0_55: hwFloat = (isNegative: false; QWordValue: 2362232012);
|
|
90 |
_0_6: hwFloat = (isNegative: false; QWordValue: 2576980377);
|
358
|
91 |
_0_7: hwFloat = (isNegative: false; QWordValue: 3006477107);
|
351
|
92 |
_0_8: hwFloat = (isNegative: false; QWordValue: 3435973837);
|
|
93 |
_0_84: hwFloat = (isNegative: false; QWordValue: 3607772528);
|
|
94 |
_0_87: hwFloat = (isNegative: false; QWordValue: 3736621547);
|
|
95 |
_0_9: hwFloat = (isNegative: false; QWordValue: 3865470566);
|
|
96 |
_0_93: hwFloat = (isNegative: false; QWordValue: 3994319585);
|
|
97 |
_0_96: hwFloat = (isNegative: false; QWordValue: 4123168604);
|
|
98 |
_0_995: hwFloat = (isNegative: false; QWordValue: 4273492459);
|
|
99 |
_0_999: hwFloat = (isNegative: false; QWordValue: 4290672328);
|
|
100 |
_1_9: hwFloat = (isNegative: false; QWordValue: 8160437862);
|
498
|
101 |
_0: hwFloat = (isNegative: false; QWordValue: 0);
|
|
102 |
_1: hwFloat = (isNegative: false; QWordValue: 4294967296);
|
|
103 |
_2: hwFloat = (isNegative: false; QWordValue: 4294967296 * 2);
|
|
104 |
_3: hwFloat = (isNegative: false; QWordValue: 4294967296 * 3);
|
|
105 |
_4: hwFloat = (isNegative: false; QWordValue: 4294967296 * 4);
|
|
106 |
_5: hwFloat = (isNegative: false; QWordValue: 4294967296 * 5);
|
|
107 |
_6: hwFloat = (isNegative: false; QWordValue: 4294967296 * 6);
|
|
108 |
_10: hwFloat = (isNegative: false; QWordValue: 4294967296 * 10);
|
|
109 |
_16: hwFloat = (isNegative: false; QWordValue: 4294967296 * 16);
|
|
110 |
_19: hwFloat = (isNegative: false; QWordValue: 4294967296 * 19);
|
|
111 |
_20: hwFloat = (isNegative: false; QWordValue: 4294967296 * 20);
|
|
112 |
_25: hwFloat = (isNegative: false; QWordValue: 4294967296 * 25);
|
|
113 |
_30: hwFloat = (isNegative: false; QWordValue: 4294967296 * 30);
|
|
114 |
_128: hwFloat = (isNegative: false; QWordValue: 4294967296 * 128);
|
|
115 |
_450: hwFloat = (isNegative: false; QWordValue: 4294967296 * 450);
|
|
116 |
_1024: hwFloat = (isNegative: false; QWordValue: 4398046511104);
|
|
117 |
_2048: hwFloat = (isNegative: false; QWordValue: 8796093022208);
|
|
118 |
_10000: hwFloat = (isNegative: false; QWordValue: 4294967296 * 10000);
|
351
|
119 |
|
|
120 |
cLittle: hwFloat = (isNegative: false; QWordValue: 1);
|
|
121 |
cHHKick: hwFloat = (isNegative: false; QWordValue: 128849018);
|
|
122 |
{$ENDIF}
|
|
123 |
|
|
124 |
{$IFNDEF FPC}
|
|
125 |
type hwFloat = Extended;
|
|
126 |
{$ENDIF}
|
|
127 |
|
|
128 |
implementation
|
|
129 |
uses uConsts;
|
|
130 |
|
|
131 |
{$IFDEF FPC}
|
|
132 |
|
498
|
133 |
function int2hwFloat (i: LongInt) : hwFloat;
|
351
|
134 |
begin
|
498
|
135 |
int2hwFloat.isNegative:= i < 0;
|
|
136 |
int2hwFloat.Round:= abs(i);
|
|
137 |
int2hwFloat.Frac:= 0
|
351
|
138 |
end;
|
|
139 |
|
|
140 |
operator + (z1, z2: hwFloat) z : hwFloat;
|
|
141 |
begin
|
|
142 |
if z1.isNegative = z2.isNegative then
|
|
143 |
begin
|
|
144 |
z.isNegative:= z1.isNegative;
|
|
145 |
z.QWordValue:= z1.QWordValue + z2.QWordValue
|
|
146 |
end
|
|
147 |
else
|
|
148 |
if z1.QWordValue > z2.QWordValue then
|
|
149 |
begin
|
|
150 |
z.isNegative:= z1.isNegative;
|
|
151 |
z.QWordValue:= z1.QWordValue - z2.QWordValue
|
|
152 |
end else
|
|
153 |
begin
|
|
154 |
z.isNegative:= z2.isNegative;
|
|
155 |
z.QWordValue:= z2.QWordValue - z1.QWordValue
|
|
156 |
end
|
|
157 |
end;
|
|
158 |
|
|
159 |
operator - (z1, z2: hwFloat) z : hwFloat;
|
|
160 |
begin
|
|
161 |
if z1.isNegative = z2.isNegative then
|
|
162 |
if z1.QWordValue > z2.QWordValue then
|
|
163 |
begin
|
|
164 |
z.isNegative:= z1.isNegative;
|
|
165 |
z.QWordValue:= z1.QWordValue - z2.QWordValue
|
|
166 |
end else
|
|
167 |
begin
|
|
168 |
z.isNegative:= not z2.isNegative;
|
|
169 |
z.QWordValue:= z2.QWordValue - z1.QWordValue
|
|
170 |
end
|
|
171 |
else begin
|
|
172 |
z.isNegative:= z1.isNegative;
|
|
173 |
z.QWordValue:= z1.QWordValue + z2.QWordValue
|
|
174 |
end
|
|
175 |
end;
|
|
176 |
|
|
177 |
operator - (z1: hwFloat) z : hwFloat;
|
|
178 |
begin
|
|
179 |
z:= z1;
|
|
180 |
z.isNegative:= not z.isNegative
|
|
181 |
end;
|
|
182 |
|
|
183 |
|
|
184 |
operator * (z1, z2: hwFloat) z : hwFloat;
|
|
185 |
begin
|
|
186 |
z.isNegative:= z1.isNegative xor z2.isNegative;
|
|
187 |
z.QWordValue:= QWord(z1.Round) * z2.Frac +
|
|
188 |
QWord(z1.Frac) * z2.Round +
|
|
189 |
((QWord(z1.Frac) * z2.Frac) shr 32);
|
|
190 |
z.Round:= z.Round + QWord(z1.Round) * z2.Round;
|
|
191 |
end;
|
|
192 |
|
|
193 |
operator * (z1: hwFloat; z2: LongInt) z : hwFloat;
|
|
194 |
begin
|
|
195 |
z.isNegative:= z1.isNegative xor (z2 < 0);
|
|
196 |
z2:= abs(z2);
|
355
|
197 |
z.QWordValue:= z1.QWordValue * z2
|
351
|
198 |
end;
|
|
199 |
|
|
200 |
operator / (z1, z2: hwFloat) z : hwFloat;
|
|
201 |
var t: hwFloat;
|
|
202 |
begin
|
|
203 |
z.isNegative:= z1.isNegative xor z2.isNegative;
|
|
204 |
z.Round:= z1.QWordValue div z2.QWordValue;
|
|
205 |
t:= z1 - z2 * z.Round;
|
|
206 |
if t.QWordValue = 0 then
|
|
207 |
z.Frac:= 0
|
|
208 |
else
|
|
209 |
begin
|
|
210 |
while ((t.QWordValue and $8000000000000000) = 0) and
|
|
211 |
((z2.QWordValue and $8000000000000000) = 0) do
|
|
212 |
begin
|
|
213 |
t.QWordValue:= t.QWordValue shl 1;
|
|
214 |
z2.QWordValue:= z2.QWordValue shl 1
|
|
215 |
end;
|
|
216 |
z.Frac:= (t.QWordValue) div (z2.Round)
|
|
217 |
end
|
|
218 |
end;
|
|
219 |
|
498
|
220 |
operator / (z1: hwFloat; z2: LongInt) z : hwFloat;
|
|
221 |
begin
|
|
222 |
z.isNegative:= z1.isNegative xor (z2 < 0);
|
|
223 |
z2:= abs(z2);
|
|
224 |
z.QWordValue:= z1.QWordValue div z2
|
|
225 |
end;
|
|
226 |
|
351
|
227 |
operator < (z1, z2: hwFloat) b : boolean;
|
|
228 |
begin
|
|
229 |
if z1.isNegative <> z2.isNegative then
|
|
230 |
b:= z1.isNegative
|
|
231 |
else
|
|
232 |
if z1.QWordValue = z2.QWordValue then
|
|
233 |
b:= false
|
|
234 |
else
|
|
235 |
b:= (z1.QWordValue < z2.QWordValue) xor z1.isNegative
|
|
236 |
end;
|
|
237 |
|
|
238 |
operator > (z1, z2: hwFloat) b : boolean;
|
|
239 |
begin
|
|
240 |
if z1.isNegative <> z2.isNegative then
|
|
241 |
b:= z2.isNegative
|
|
242 |
else
|
|
243 |
if z1.QWordValue = z2.QWordValue then
|
|
244 |
b:= false
|
|
245 |
else
|
|
246 |
b:= (z1.QWordValue > z2.QWordValue) xor z2.isNegative
|
|
247 |
end;
|
|
248 |
|
|
249 |
function cstr(z: hwFloat): string;
|
|
250 |
var tmpstr: string;
|
|
251 |
begin
|
|
252 |
str(z.Round, cstr);
|
|
253 |
if z.Frac <> 0 then
|
|
254 |
begin
|
|
255 |
str(z.Frac / $100000000:1:15, tmpstr);
|
|
256 |
delete(tmpstr, 1, 2);
|
|
257 |
cstr:= cstr + '.' + tmpstr
|
|
258 |
end;
|
|
259 |
if z.isNegative then cstr:= '-' + cstr
|
|
260 |
end;
|
|
261 |
|
371
|
262 |
function hwRound(t: hwFloat): LongInt;
|
351
|
263 |
begin
|
|
264 |
if t.isNegative then hwRound:= -t.Round
|
|
265 |
else hwRound:= t.Round
|
|
266 |
end;
|
|
267 |
|
|
268 |
function hwAbs(t: hwFloat): hwFloat;
|
|
269 |
begin
|
|
270 |
hwAbs:= t;
|
|
271 |
hwAbs.isNegative:= false
|
|
272 |
end;
|
|
273 |
|
|
274 |
function hwSqr(t: hwFloat): hwFloat;
|
|
275 |
begin
|
|
276 |
hwSqr:= t * t
|
|
277 |
end;
|
|
278 |
|
357
|
279 |
function hwSqrt(t: hwFloat): hwFloat;
|
|
280 |
begin
|
|
281 |
hwSqrt.isNegative:= false;
|
|
282 |
hwSqrt.QWordValue:= Round(sqrt(1.0 / $100000000 * (t.QWordValue)) * $100000000)
|
|
283 |
end;
|
|
284 |
|
351
|
285 |
function Distance(dx, dy: hwFloat): hwFloat;
|
|
286 |
var x, y: hwFloat;
|
|
287 |
Result: hwFloat;
|
|
288 |
begin
|
|
289 |
x:= dx * dx;
|
|
290 |
y:= dy * dy;
|
|
291 |
Result:= x + y;
|
357
|
292 |
Result.QWordValue:= Round(sqrt(1.0 / $100000000 * (Result.QWordValue)) * $100000000);
|
351
|
293 |
Distance:= Result
|
|
294 |
end;
|
|
295 |
|
498
|
296 |
function DistanceI(dx, dy: LongInt): hwFloat;
|
|
297 |
begin
|
|
298 |
DistanceI:= Distance(int2hwFloat(dx), int2hwFloat(dy))
|
|
299 |
end;
|
|
300 |
|
|
301 |
function SignAs(num, signum: hwFloat): hwFloat;
|
|
302 |
begin
|
|
303 |
SignAs:= num;
|
|
304 |
SignAs.isNegative:= signum.isNegative
|
|
305 |
end;
|
|
306 |
|
357
|
307 |
{$INCLUDE SinTable.inc}
|
|
308 |
|
|
309 |
function AngleSin(Angle: Longword): hwFloat;
|
351
|
310 |
begin
|
|
311 |
AngleSin.isNegative:= false;
|
357
|
312 |
if Angle < 1024 then AngleSin.QWordValue:= SinTable[Angle]
|
|
313 |
else AngleSin.QWordValue:= SinTable[2048 - Angle]
|
351
|
314 |
end;
|
|
315 |
|
357
|
316 |
function AngleCos(Angle: Longword): hwFloat;
|
351
|
317 |
begin
|
357
|
318 |
AngleCos.isNegative:= Angle > 1024;
|
|
319 |
if Angle < 1024 then AngleCos.QWordValue:= SinTable[1024 - Angle]
|
|
320 |
else AngleCos.QWordValue:= SinTable[Angle - 1024]
|
351
|
321 |
end;
|
|
322 |
|
|
323 |
{$ENDIF}
|
|
324 |
|
|
325 |
end.
|