From ab9808b04e22fae1fb867703be71a537ecf604b9 Mon Sep 17 00:00:00 2001 From: Erich Eckner Date: Mon, 24 Aug 2015 16:26:55 +0200 Subject: Tagesendstand, immernoch langsam, aber vmtl. ist die Methode einfach langsam --- Physikunit.pas | 14 ++---- Plasmapropagation.lps | 116 +++++++++++++++++++++++++------------------------- input.plap | 2 +- 3 files changed, 62 insertions(+), 70 deletions(-) diff --git a/Physikunit.pas b/Physikunit.pas index ef77ee4..556334d 100644 --- a/Physikunit.pas +++ b/Physikunit.pas @@ -1142,10 +1142,7 @@ begin // Ableitungen nach PX - write('.'); - len:=length(inhalt[0].phasenraum); - write(len); setlength(res,len); setlength(ims,len); for i:=0 to length(inhalt)-1 do @@ -1168,16 +1165,14 @@ begin end; fft(res,ims,true); for j:=0 to length(inhalt[i].phasenraum)-1 do - inhalt[i].phasenraum[j].grad[tlc,0]:=res[j]*iDP*0; + inhalt[i].phasenraum[j].grad[tlc,0]:=res[j]*iDP; end; // Ableitungen nach X - write('.'); len:=length(inhalt); - write(len); - setlength(res,2*len); - setlength(ims,2*len); + setlength(res,len); + setlength(ims,len); for i:=0 to length(inhalt[0].phasenraum)-1 do for tlc:=0 to length(teilchen)-1 do begin for j:=0 to length(inhalt)-1 do begin @@ -1196,9 +1191,8 @@ begin end; fft(res,ims,true); for j:=0 to length(inhalt)-1 do - inhalt[j].phasenraum[i].grad[tlc,1]:=res[j]*iDX*0; + inhalt[j].phasenraum[i].grad[tlc,1]:=res[j]*iDX; end; - write('>'); end; procedure tFelder.berechneAbleitungen(dX,iDX: extended); diff --git a/Plasmapropagation.lps b/Plasmapropagation.lps index 557ace6..c266d7b 100644 --- a/Plasmapropagation.lps +++ b/Plasmapropagation.lps @@ -16,10 +16,10 @@ - - - - + + + + @@ -27,15 +27,14 @@ - + - - - + + @@ -45,14 +44,14 @@ - + - + @@ -66,10 +65,9 @@ - - - - + + + @@ -155,123 +153,123 @@ - + - + - - + + - - + + - - + + - - + + - - + + - - + + - - + + - - + + - - + + - + - + - + - - + + - - + + - - + + - - + + - - + + - - + + - + - + - + - + - + - + - + - + - + - + diff --git a/input.plap b/input.plap index f83acba..3196fee 100644 --- a/input.plap +++ b/input.plap @@ -10,7 +10,7 @@ allgemein ortsschritt 10^-2 * λ zeitschritt 10^-2 * T maximalimpuls 5 - impulsschritt 10^-2 + impulsschritt 2*10^-2 zeit 2 * T !setze $breite: (5 * λ) breite $breite -- cgit v1.2.3-70-g09d2