fuegeSortiertHinzu neu in lowlevelunit.pas

1 files changed, 24 insertions, 0 deletions

diff --git a/lowlevelunit.pas b/lowlevelunit.pas
index 52786db..940e4bc 100644
--- a/lowlevelunit.pas
+++ b/lowlevelunit.pas
@@ -101,6 +101,7 @@ function myUtf8Encode(s: string): string;
 function permutation(len: longint): tLongintArray; overload;
 function permutation(len,rest,maxBasis: longint; basen: tLongintArray): tLongintArray; overload;
 procedure llPermutation(len,offset: longint; var ar: tLongintArray); inline;
+procedure fuegeSortiertHinzu(x: extended; var xa: tExtendedArray);
 
 var
   base64Chars: array[0..63] of char;
@@ -923,6 +924,29 @@ begin
       fehler('permutation: '+inttostr(i)+' wurde nicht verteilt!');
 end;
 
+procedure fuegeSortiertHinzu(x: extended; var xa: tExtendedArray);
+var
+  mi,ma,i: longint;
+begin
+  mi:=0;
+  ma:=length(xa)-1;
+  while mi<=ma do begin
+    i:=(mi+ma) div 2;
+    if x<xa[i] then
+      ma:=i-1
+    else if x>xa[i] then
+      mi:=i+1
+    else
+      exit;
+  end;
+  if mi<>ma+1 then
+    fehler('Bisektion fehlgeschlagen! ('+inttostr(mi)+' = mi <> ma+1 = '+inttostr(ma+1)+')');
+  setlength(xa,length(xa)+1);
+  for i:=length(xa)-1 downto mi+1 do
+    xa[i]:=xa[i-1];
+  xa[mi]:=x;
+end;
+
 var
   b: byte;