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
244
245
246
247
248
249
250
251
252
253
254
|
unit matheunit;
{$mode objfpc}{$H+}
interface
uses
Classes, SysUtils, Gmp, Math;
type
tExtPoint = record
x,y: extended;
end;
tExtPointArray = array of tExtPoint;
pTExtPointArray = ^tExtPointArray;
tLongintArray = array of longint;
pTLongintArray = ^tLongintArray;
tExtendedArray = array of extended;
function plus(a,b: tExtPoint): tExtPoint;
function durch(a: tExtPoint; b: extended): tExtPoint;
function myFrac(x: extended): extended;
function mpfToStr(f: mpf_t): string;
function signSqr(x: extended): extended; inline;
function mpfMyRoot(rad: mpf_t; wzlExp: int64): extended;
function myTimeToStr(t: extended): string;
function cmpStr(s1,s2: string): longint;
function mitte(s1,s2: string): string;
function myFloatToStr(x: extended): string;
function myStrToFloat(s: string): extended;
procedure copyArray(i: tExtPointArray; out o: tExtPointArray); overload;
procedure copyArray(i: tLongintArray; out o: tLongintArray); overload;
procedure copyArray(i: tExtendedArray; out o: tExtendedArray); overload;
implementation
function plus(a,b: tExtPoint): tExtPoint;
begin
result.x:=a.x+b.x;
result.y:=a.y+b.y;
end;
function durch(a: tExtPoint; b: extended): tExtPoint;
begin
result.x:=a.x/b;
result.y:=a.y/b;
end;
function myFrac(x: extended): extended;
begin
result:=frac(x);
while result<0 do
result:=result+1;
end;
function mpfToStr(f: mpf_t): string;
var
ex: int64;
off: byte;
begin
result:=mpf_get_str(nil,ex,10,0,f);
off:=1+byte(pos('-',result)=1);
if result='' then
result:='0'
else if ex=1 then
result:=copy(result,1,off)+','+copy(result,off+1,length(result)-off)
else
result:=copy(result,1,off)+','+copy(result,off+1,length(result)-off)+' * 10^'+inttostr(ex-1);
end;
function signSqr(x: extended): extended;
begin
result:=sign(x)*sqr(x);
end;
function mpfMyRoot(rad: mpf_t; wzlExp: int64): extended;
var
ex: int64;
begin
result:=power(mpf_get_d_2exp(ex,rad),1/wzlExp);
result:=result*power(2,ex/wzlExp);
end;
function myTimeToStr(t: extended): string;
var
tim: int64;
begin
tim:=floor(t*24*60*60);
result:=inttostr(tim mod 10)+'s';
tim:=tim div 10;
if tim=0 then exit;
result:=inttostr(tim mod 6)+result;
tim:=tim div 6;
if tim=0 then exit;
result:=inttostr(tim mod 10)+'min '+result;
tim:=tim div 10;
if tim=0 then exit;
result:=inttostr(tim mod 6)+result;
tim:=tim div 6;
if tim=0 then exit;
result:=inttostr(tim mod 24)+'h '+result;
tim:=tim div 24;
if tim=0 then exit;
result:=' '+result;
if (tim mod 7)<>1 then
result:='e'+result;
result:=inttostr(tim mod 7)+'Tag'+result;
tim:=tim div 7;
if tim=0 then exit;
result:=' '+result;
if tim<>1 then
result:='n'+result;
result:=inttostr(tim)+'Woche'+result;
end;
function cmpStr(s1,s2: string): longint;
var
i: longint;
begin
for i:=1 to min(length(s1),length(s2)) do
if s1[i]<>s2[i] then begin
result:=2*byte(s1[i]>s2[i])-1;
exit;
end;
if length(s1)<>length(s2) then begin
result:=2*byte(length(s1)>length(s2))-1;
exit;
end;
result:=0;
end;
function mitte(s1,s2: string): string;
var
i: longint;
w,nw: word;
begin
setlength(result,max(length(s1),length(s2)));
w:=0;
for i:=length(result) downto 1 do begin // result:= "s1+s2";
if i<=length(s1) then
w:=w+byte(s1[i]);
if i<=length(s2) then
w:=w+byte(s2[i]);
result[i]:=char(w and $ff);
w:=w shr 8;
end;
result:=char(w)+result;
w:=0;
for i:=1 to length(result) do begin
nw:=byte(odd(byte(result[i])+w));
result[i]:=char((byte(result[i])+w) div 2);
w:=nw shl 8;
end;
if w<>0 then
result:=result+char(w div 2);
if result[1]<>#0 then begin
writeln('Fehler bei der Mittenfindeung!');
halt;
end;
delete(result,1,1);
end;
function myFloatToStr(x: extended): string;
var
i,e: longint;
begin
e:=0;
if x<0 then begin
result:='-';
x:=-x;
end
else
result:='';
if x=0 then begin
result:='0';
exit;
end;
while x<1 do begin
dec(e);
x:=x*10;
end;
while x>=10 do begin
inc(e);
x:=x/10;
end;
result:=result+char(ord('0')+floor(x))+'.';
for i:=0 to 20 do begin
x:=(x-floor(x))*10;
result:=result+char(ord('0')+floor(x));
end;
if e<>0 then
result:=result+'E'+inttostr(e);
end;
function myStrToFloat(s: string): extended;
var
i,e: longint;
neg: boolean;
begin
if pos('E',s)>0 then begin
e:=strtoint(rightStr(s,length(s)-pos('E',s)));
delete(s,pos('E',s),length(s));
end
else e:=0;
if pos('.',s)=0 then begin
result:=strtoint(s)*power(10,e);
exit;
end;
neg:=leftStr(s,1)='-';
if neg then
delete(s,1,1);
if pos('.',s)=2 then begin
result:=0;
for i:=length(s) downto 3 do
result:=result/10 + ord(s[i])-ord('0');
result:=result/10 + ord(s[1])-ord('0');
end
else result:=strtofloat(s);
result:=result*power(10,e);
if neg then
result:=-result;
end;
procedure copyArray(i: tExtPointArray; out o: tExtPointArray);
var
j: longint;
begin
setlength(o,length(i));
for j:=0 to length(o)-1 do begin
o[j].x:=i[j].x;
o[j].y:=i[j].y;
end;
end;
procedure copyArray(i: tLongintArray; out o: tLongintArray);
var
j: longint;
begin
setlength(o,length(i));
for j:=0 to length(o)-1 do
o[j]:=i[j];
end;
procedure copyArray(i: tExtendedArray; out o: tExtendedArray);
var
j: longint;
begin
setlength(o,length(i));
for j:=0 to length(o)-1 do
o[j]:=i[j];
end;
end.
|
