Physikunit.pas aufgefrischt

1 files changed, 15 insertions, 3 deletions

diff --git a/Physikunit.pas b/Physikunit.pas
index 680419b..935dcad 100644
--- a/Physikunit.pas
+++ b/Physikunit.pas
@@ -6,7 +6,6 @@ unit Physikunit;
   {$ERROR This program can be compiled only on/for Unix/Linux based systems.}
 {$ENDIF}
 
-{$DEFINE Zeitschrittueberwachung}
 {$DEFINE Dichteueberwachung}
 {$DEFINE negativeDichteueberwachung}
 { $DEFINE exzessiveArrayBereichsTests}
@@ -133,7 +132,7 @@ type
   {$IFDEF exzessiveArrayBereichsTests}
     procedure pruefeArrayEnden(fehler: string);
   {$ENDIF}
-    procedure nichtnegativieren;
+    procedure nichtnegativieren; inline;
     procedure berechnePhasenraumAbleitungen; inline;
   public
     emFelder:  array[tEMFeldGroesze,boolean] of pDouble;         // EM-Felder und deren Zeitableitungen
@@ -1070,6 +1069,18 @@ begin
         (impulsraum[i,false]+j+(aX-1)*aP)^:=0;
       end;
 
+  for i:=0 to length(teilchen)-1 do
+    for j:=0 to aX-1 do begin
+      if (impulsraum[i,false]+j*aP)^>0 then begin // unterer Rand
+        (impulsraum[i,false]+1+j*aP)^:=(impulsraum[i,false]+1+j*aP)^+(impulsraum[i,false]+j*aP)^;
+        (impulsraum[i,false]+j*aP)^:=0;
+      end;
+      if (impulsraum[i,false]+aP-1+j*aP)^>0 then begin // oberer Rand
+        (impulsraum[i,false]+aP-2+j*aP)^:=(impulsraum[i,false]+aP-2+j*aP)^+(impulsraum[i,false]+aP-1+j*aP)^;
+        (impulsraum[i,false]+aP-1+j*aP)^:=0;
+      end;
+    end;
+
   // PY und iMGamma berechnen
   for i:=0 to length(matFelder)-1 do
     for j:=0 to aX-1 do begin
@@ -1613,8 +1624,9 @@ begin
       if (felders[aktuelleFelder].impulsraum[i,false]+j)^<0 then begin
         pro:=tProtokollant.create(prot,'iteriereSchritt');
         pro.schreibe('n<0 bei:',true);
-        pro.schreibe(' t    = ' +floattostr(t)+',',true);
+        pro.schreibe(' t    = ' +floattostr(t)+' (aktuelleFelder = '+inttostr(aktuelleFelder)+'),',true);
         pro.schreibe(' i    = ' +inttostr(i)+' (x = '+floattostr(xl+i*dX)+'),',true);
+        pro.schreibe(' n    = ' +floattostr((felders[aktuelleFelder].impulsraum[i,false]+j)^),true);
         abbrechen;
       end;
   {$ENDIF}