diff options
| author | Erich Eckner <git@eckner.net> | 2015-10-05 13:58:58 +0200 |
|---|---|---|
| committer | Erich Eckner <git@eckner.net> | 2015-10-05 14:03:14 +0200 |
| commit | 88699eb92b33f51062b19105aded78ac7c3e2bf8 (patch) | |
| tree | 245bbeb2d78e17d7d5add5582c24196fe64fc11b | |
| parent | 8eb4bb7cd4c67415a6fa482b21918ea3dec97ee8 (diff) | |
| download | Plasmapropagation-88699eb92b33f51062b19105aded78ac7c3e2bf8.tar.xz | |
negativeDichtenüberwachung geändert und aktiviert
| -rw-r--r-- | Physikunit.pas | 51 |
1 files changed, 10 insertions, 41 deletions
diff --git a/Physikunit.pas b/Physikunit.pas index 4b73d8b..0e8075b 100644 --- a/Physikunit.pas +++ b/Physikunit.pas @@ -8,7 +8,7 @@ unit Physikunit; {$DEFINE Zeitschrittueberwachung} {$DEFINE Dichteueberwachung} -{ $DEFINE negativeDichteueberwachung} +{$DEFINE negativeDichteueberwachung} interface @@ -1194,46 +1194,15 @@ begin aktuelleFelder:=1-aktuelleFelder; {$IFDEF negativeDichteueberwachung} - for i:=0 to length(felders[aktuelleFelder].inhalt)-1 do - for j:=0 to felders[aktuelleFelder].matAnz-1 do - if felders[aktuelleFelder].inhalt[i].matWerte[j,mfN,false]<0 then begin - pro:=tProtokollant.create(prot,'iteriereSchritt'); - pro.schreibe('n<0 bei:',true); - pro.schreibe(' t = ' +floattostr(t)+',',true); - pro.schreibe(' i = ' +inttostr(i)+' (x = '+floattostr(xl+i*dX)+'),',true); - pro.schreibe(' alte Werte:',true); - with felders[1-aktuelleFelder] do begin - pro.schreibe(' v_l = ' +floattostr(inhalt[i-1].matWerte[j,mfDPsiDX,false]* - inhalt[i-1].matWerte[j,mfIGamma,false])+',',true); - pro.schreibe(' n_l = ' +floattostr(inhalt[i-1].matWerte[j,mfN,false])+',',true); - pro.schreibe(' n_l'' = '+floattostr(inhalt[i-1].matWerte[j,mfN,true])+',',true); - pro.schreibe(' v = ' +floattostr(inhalt[i].matWerte[j,mfDPsiDX,false]* - inhalt[i].matWerte[j,mfIGamma,false])+',',true); - pro.schreibe(' n = ' +floattostr(inhalt[i].matWerte[j,mfN,false])+',',true); - pro.schreibe(' n'' = '+floattostr(inhalt[i].matWerte[j,mfN,true])+',',true); - pro.schreibe(' v_r = ' +floattostr(inhalt[i+1].matWerte[j,mfDPsiDX,false]* - inhalt[i+1].matWerte[j,mfIGamma,false])+',',true); - pro.schreibe(' n_r = ' +floattostr(inhalt[i+1].matWerte[j,mfN,false])+',',true); - pro.schreibe(' n_r'' = '+floattostr(inhalt[i+1].matWerte[j,mfN,true]),true); - end; - pro.schreibe(' neue Werte:',true); - with felders[aktuelleFelder] do begin - pro.schreibe(' v_l = ' +floattostr(inhalt[i-1].matWerte[j,mfDPsiDX,false]* - inhalt[i-1].matWerte[j,mfIGamma,false])+',',true); - pro.schreibe(' n_l = ' +floattostr(inhalt[i-1].matWerte[j,mfN,false])+',',true); - pro.schreibe(' n_l'' = '+floattostr(inhalt[i-1].matWerte[j,mfN,true])+',',true); - pro.schreibe(' v = ' +floattostr(inhalt[i].matWerte[j,mfDPsiDX,false]* - inhalt[i].matWerte[j,mfIGamma,false])+',',true); - pro.schreibe(' n = ' +floattostr(inhalt[i].matWerte[j,mfN,false])+',',true); - pro.schreibe(' n'' = '+floattostr(inhalt[i].matWerte[j,mfN,true])+',',true); - pro.schreibe(' v_r = ' +floattostr(inhalt[i+1].matWerte[j,mfDPsiDX,false]* - inhalt[i+1].matWerte[j,mfIGamma,false])+',',true); - pro.schreibe(' n_r = ' +floattostr(inhalt[i+1].matWerte[j,mfN,false])+',',true); - pro.schreibe(' n_r'' = '+floattostr(inhalt[i+1].matWerte[j,mfN,true]),true); - end; - pro.free; - abbrechen; - end; + for i:=0 to length(felders[aktuelleFelder].impulsraum)-1 do + for j:=0 to felders[aktuelleFelder].aX*felders[aktuelleFelder].aP-1 do + 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(' i = ' +inttostr(i)+' (x = '+floattostr(xl+i*dX)+'),',true); + abbrechen; + end; {$ENDIF} t:=t+dT; end; |