summaryrefslogtreecommitdiff
diff options
context:
space:
mode:
-rw-r--r--Physikunit.pas602
-rw-r--r--Plasmapropagation.lpi1
-rw-r--r--Plasmapropagation.lps193
-rw-r--r--input.epost10
-rw-r--r--input.plap24
-rw-r--r--linearkombination.inc2
-rw-r--r--linearkombinationen.inc84
7 files changed, 523 insertions, 393 deletions
diff --git a/Physikunit.pas b/Physikunit.pas
index 727f0dc..4861a14 100644
--- a/Physikunit.pas
+++ b/Physikunit.pas
@@ -26,6 +26,8 @@ type
tEMQuellGroesze = (
eqRho,eqJX,eqJY
);
+ tMaterieFeldGroesze = (mfPY);
+ tTeilchenSpeziesGroeszen = (tsgMasse,tsgIQdrMasse,tsgLadung);
tEMQuellen = array[tEMQuellGroesze] of extended;
const
@@ -68,12 +70,14 @@ type
tTeilchenSpezies = class
private
- ortsGrenzen: array of extended;
- dichteStuecke: array of string;
+ ortsGrenzen: array of extended;
+ dichteStuecke: array of string;
public
- nMin,nMax, // minimale und maximale Massendichte in nc
- spezifischeLadung: extended; // q/m in e/me
- constructor create;
+ nMin,nMax,breite: extended; // minimale und maximale Massendichte in nc; Breite der Impulsverteilung
+ eigenschaften: array[tTeilchenSpeziesGroeszen] of extended; // q in e; m in me
+ diffusion: array[0..1] of extended; // Diffusionsterm in x- und p-Richtung (= maximales d(ln n)/d[xp] * iD[XP])
+ constructor create; overload;
+ constructor create(original: tTeilchenSpezies); overload;
destructor destroy; override;
procedure dump(prot: tProtokollant; prefix: string);
function gibDichte(x: extended; kvs: tKnownValues): extended; // Massendichte
@@ -85,19 +89,20 @@ type
tImpulsPunkt = class(tObject) // Besetzungsdichten verschiedener Teilchenspezies für ein Paar (r;p)
private
raumpunkt: tRaumPunkt;
- nachbarn: array[0..2,boolean] of tImpulsPunkt; // x-, px-, py-Nachbarn (jeweil in Richtung - und +)
+ nachbarn: array[0..1,boolean] of tImpulsPunkt; // x- und px-Nachbarn (jeweil in Richtung - und +)
public
- werte: array of array[boolean] of extended; // Besetzungsdichten und deren Ableitungen für verschiedene Teilchenspezies
- constructor create(rp: tRaumPunkt; linksRaum,linksImpuls,untenImpuls: tImpulsPunkt; anzTeilchen: longint);
+ werte: array of array[boolean] of extended; // Besetzungsdichten und deren Ableitungen für verschiedene Teilchenspezies
+ iMGamma: array of extended; // 1/m/gamma (= v/p)
+ constructor create(rp: tRaumPunkt; linksRaum,linksImpuls: tImpulsPunkt; anzTeilchen: longint);
destructor destroy; override;
{$DEFINE LiKotImpulsPunktHeader}
{$INCLUDE linearkombinationen.inc}
{$UNDEF LiKotImpulsPunktHeader}
- procedure nichtnegativieren;
- procedure akkumuliereEMQuellen(var emQuellen: tEMQuellen; vx,vy,dVP: extended);
- function gibMirWas(var defizit: extended; entfernung, teilchen,richtung: longint; positiv: boolean): boolean; // ich habe jemanden erreicht
- procedure setzeNull;
- procedure berechneAbleitungen(pX,pY,iDX,iDPX,iDPY: extended);
+ function nichtnegativieren: extended; inline;
+ procedure akkumuliereEMQuellen(var emQuellen: tEMQuellen; pX,dVP: extended); inline;
+ function gibMirWas(var defizit: extended; entfernung, teilchen,richtung: longint; positiv: boolean): boolean; inline; // ich habe jemanden erreicht
+ procedure setzeNull; inline;
+ procedure berechneAbleitungen(pX,iDX,iDPX: extended); inline;
end;
{ tRaumPunkt }
@@ -106,56 +111,55 @@ type
private
felder: tFelder;
public
- phasenraum: array of tImpulsPunkt;
- // Phasenraumbesetzungsdichten
- emFelder: array[tEMFeldGroesze,boolean] of extended;
- // EM-Felder und deren Zeitableitungen
- // dPhiDX und B[xyz] haben keine sinnvolle Ableitung hier!
- emQuellen: array[tEMQuellGroesze] of extended;
- // Quellen für die EM-Felder
- lN,rN: tRaumPunkt;
- // Nachbarpunkte
- dPX,dPY: extended;
- // Impulsdiskretisierung
- aPX,aPY: longint;
- // Anzahl der Impulsbins
-
- constructor create(linkerNachbar: tRaumPunkt; derBesitzer: tFelder; teilchen,_aPX,_aPY: longint; _dPX,_dPY: extended);
+ phasenraum: array of tImpulsPunkt; // Phasenraumbesetzungsdichten
+ emFelder: array[tEMFeldGroesze,boolean] of extended; // EM-Felder und deren Zeitableitungen
+ // dPhiDX und B[xyz] haben keine sinnvolle Ableitung hier!
+ emQuellen: array[tEMQuellGroesze] of extended; // Quellen für die EM-Felder
+ matFelder: array of array[tMaterieFeldGroesze,boolean] of extended; // impulsunabhängige Materiefelder getrennt nach Teilchenspezies
+ lN,rN: tRaumPunkt; // Nachbarpunkte
+ dP: extended; // Impulsdiskretisierung
+ aP: longint; // Anzahl der Impulsbins
+
+ constructor create(linkerNachbar: tRaumPunkt; derBesitzer: tFelder; teilchen,_aP: longint; _dP: extended);
destructor destroy; override;
- procedure berechneEMFelder(dX: extended; rechts: boolean); overload;
- procedure berechneEMFelder(dX,iDX: extended); overload;
- procedure berechneEMAbleitungen(iDX: extended);
- procedure berechnePhasenraumAbleitungen(iDX: extended);
+ procedure berechneEMFelder(dX: extended; rechts: boolean); overload; inline;
+ procedure berechneEMFelder(dX,iDX: extended); overload; inline;
+ procedure berechneEMQuellen;
+ procedure berechneEMAbleitungen(iDX: extended); inline;
+ procedure berechnePhasenraumAbleitungen(iDX: extended); inline;
{$DEFINE LiKotRaumPunktHeader}
{$INCLUDE linearkombinationen.inc}
{$UNDEF LiKotRaumPunktHeader}
- procedure nichtnegativieren;
- function impulsIntegral(teilchen: longint; emQ: tEMQuellGroesze): extended;
- procedure initialisiereDichte(teilchen: longint; n: extended);
- procedure reflektiereTeilchen(rechts: boolean);
- procedure setzeNull;
+ function nichtnegativieren: extended; inline;
+ function impulsIntegral(teilchen: longint; emQ: tEMQuellGroesze): extended; overload; inline;
+ function impulsIntegral(teilchen: longint): extended; overload; inline;
+ procedure initialisiereDichte(teilchen: longint; breite,n: extended); inline;
+ procedure reflektiereTeilchen(rechts: boolean); inline;
+ procedure setzeNull; inline;
+ procedure berechnePY; inline;
end;
{ TFelder }
tFelder = class(tObject) // repräsentiert eine Simulationsbox von Feldern inklusive deren Ableitungen
private
- matAnz: longint;
- spezLadungen,
- massen: tExtendedArray;
- lichters: array[boolean] of tMyStringlist;
+ teilchen: array of tTeilchenSpezies;
+ lichters: array[boolean] of tMyStringlist;
+ massen: array of extended; // Gesamtmassen je Teilchenspezies zur Überwachung der Erhaltungsgrößen
+ gesamtDefizit: extended;
procedure setzeRaender(iDX: extended);
public
inhalt: array of tRaumPunkt;
gitter: tGitter;
- constructor create(groesse: longint; teilchen: array of tTeilchenSpezies; lichter: tMyStringList; parent: tGitter; aPX,aPY: longint; dPX,dPY: extended);
+ constructor create(groesse: longint; _teilchen: array of tTeilchenSpezies; lichter: tMyStringList; parent: tGitter; aP: longint; dP: extended);
destructor destroy; override;
- procedure berechneAbleitungen(dX,iDX: extended);
+ procedure berechneAbleitungen(dX,iDX: extended); inline;
{$DEFINE LiKotFelderHeader}
{$INCLUDE linearkombinationen.inc}
{$UNDEF LiKotFelderHeader}
- procedure setzeNull;
+ procedure setzeNull; inline;
+ procedure dumpErhaltungsgroeszen(prot: tProtokollant);
end;
{ tGitter }
@@ -174,9 +178,10 @@ type
aktuelleFelder: longint;
felders: array of tFelder; // mehrere komplette Simulationsboxen von Feldern, benötigt um Zwischenschritte für die Zeitentwicklung zu speichern
- constructor create(derBesitzer: tSimulation; size: longint; deltaX: extended; bekannteWerte: tKnownValues; teilchen: array of tTeilchenSpezies; lichter: tMyStringlist; zv: tZeitverfahren; protokollant: tProtokollant; name: string; aPX,aPY: longint; dPX,dPY: extended);
+ constructor create(derBesitzer: tSimulation; size: longint; deltaX: extended; bekannteWerte: tKnownValues; teilchen: array of tTeilchenSpezies; lichter: tMyStringlist; zv: tZeitverfahren; protokollant: tProtokollant; name: string; aP: longint; dP: extended);
destructor destroy; override;
procedure iteriereSchritt(dT: extended);
+ procedure dumpErhaltungsgroeszen;
end;
{ tSimulation }
@@ -349,7 +354,7 @@ procedure tAusgabeDatei.speichereWerte(gitter: tGitter; sT, sDX: extended);
var i,cnt: longint;
sX,cX,val: extended;
begin
- if (teilchen>=gitter.felders[gitter.aktuelleFelder].matAnz) then begin
+ if (teilchen>=length(gitter.felders[gitter.aktuelleFelder].teilchen)) then begin
pro.schreibe('Teilchen '+inttostr(teilchen+1)+' gibt es nicht, da kann ich auch nichts speichern!',true);
pro.destroyall;
halt(1);
@@ -424,6 +429,8 @@ end;
// tTeilchenSpezies ************************************************************
constructor tTeilchenSpezies.create;
+var
+ egr: tTeilchenSpeziesGroeszen;
begin
inherited create;
fillchar(ortsGrenzen,sizeof(ortsGrenzen),#0);
@@ -432,7 +439,26 @@ begin
setlength(dichteStuecke,0);
nMin:=1E-3;
nMax:=nMin;
- spezifischeLadung:=0;
+ for egr:=low(tTeilchenSpeziesGroeszen) to high(tTeilchenSpeziesGroeszen) do
+ eigenschaften[egr]:=byte(egr in [tsgMasse,tsgIQdrMasse]);
+ diffusion[0]:=0;
+ diffusion[1]:=0;
+ breite:=1;
+end;
+
+constructor tTeilchenSpezies.create(original: tTeilchenSpezies);
+var
+ egr: tTeilchenSpeziesGroeszen;
+ i: longint;
+begin
+ create;
+ nMin:=original.nMin;
+ nMax:=original.nMax;
+ for egr:=low(tTeilchenSpeziesGroeszen) to high(tTeilchenSpeziesGroeszen) do
+ eigenschaften[egr]:=original.eigenschaften[egr];
+ for i:=0 to 1 do
+ diffusion[i]:=original.diffusion[i];
+ breite:=original.breite;
end;
destructor tTeilchenSpezies.destroy;
@@ -446,8 +472,11 @@ procedure tTeilchenSpezies.dump(prot: tProtokollant; prefix: string);
var
i: longint;
begin
- prot.schreibe(prefix+'nMax = '+floattostr(nMax)+' * nc');
- prot.schreibe(prefix+'q/m = '+floattostr(-spezifischeLadung)+' e/me');
+ prot.schreibe(prefix+'nMax = '+floattostr(nMax)+' * nc');
+ prot.schreibe(prefix+'m = '+floattostr(eigenschaften[tsgMasse])+' me');
+ prot.schreibe(prefix+'q = '+floattostr(eigenschaften[tsgLadung])+' e');
+ prot.schreibe(prefix+'dLnN/dX < '+floattostr(diffusion[0])+' / dX');
+ prot.schreibe(prefix+'dLnN/dP < '+floattostr(diffusion[1])+' / dP');
prot.schreibe(prefix+'n(x):');
for i:=0 to length(dichteStuecke)-1 do begin
prot.schreibe(prefix+' '+dichteStuecke[i]);
@@ -520,7 +549,7 @@ end;
// tImpulsPunkt ****************************************************************
-constructor tImpulsPunkt.create(rp: tRaumPunkt; linksRaum,linksImpuls,untenImpuls: tImpulsPunkt; anzTeilchen: longint);
+constructor tImpulsPunkt.create(rp: tRaumPunkt; linksRaum,linksImpuls: tImpulsPunkt; anzTeilchen: longint);
var
i: longint;
b: boolean;
@@ -532,13 +561,15 @@ begin
for i:=0 to length(werte)-1 do
for b:=false to true do
werte[i,b]:=0;
+ fillchar(iMGamma,sizeof(iMGamma),#0);
+ setlength(iMGamma,length(werte));
+ for i:=0 to length(iMGamma)-1 do
+ iMGamma[i]:=0;
nachbarn[0,false]:=linksRaum;
nachbarn[0,true]:=nil;
nachbarn[1,false]:=linksImpuls;
nachbarn[1,true]:=nil;
- nachbarn[2,false]:=untenImpuls;
- nachbarn[2,true]:=nil;
- for i:=0 to 2 do
+ for i:=0 to 1 do
if assigned(nachbarn[i,false]) then
nachbarn[i,false].nachbarn[i,true]:=self;
end;
@@ -548,33 +579,36 @@ var
i: longint;
b: boolean;
begin
- for i:=0 to 3 do
+ for i:=0 to 1 do
for b:=false to true do
if assigned(nachbarn[i,b]) then
nachbarn[i,b].nachbarn[i,not b]:=nil;
setlength(werte,0);
+ setlength(iMGamma,0);
inherited destroy;
end;
-procedure tImpulsPunkt.nichtnegativieren; // Dichten nicht negativ machen
+function tImpulsPunkt.nichtnegativieren: extended; // Dichten nicht negativ machen
var
- i,j,entfernung: longint;
- defizit: extended;
- pro: tProtokollant;
- jemandErreicht,b: boolean;
+ i,richtung,entfernung: longint;
+ defizit: extended;
+ pro: tProtokollant;
+ jemandErreicht,vorZurueck: boolean;
begin
+ result:=0;
for i:=0 to length(werte)-1 do
if werte[i,false]<0 then begin
defizit:=-werte[i,false];
+ result:=result+defizit;
werte[i,false]:=0;
entfernung:=1;
jemandErreicht:=true;
while (defizit>0) and jemandErreicht do begin
jemandErreicht:=false;
- for j:=0 to 3 do
- for b:=false to true do
- if assigned(nachbarn[j,b]) then
- jemandErreicht:=nachbarn[j,b].gibMirWas(defizit,i,entfernung,j,b) or jemandErreicht;
+ for richtung:=0 to 1 do
+ for vorZurueck:=false to true do
+ if assigned(nachbarn[richtung,vorZurueck]) then
+ jemandErreicht:=nachbarn[richtung,vorZurueck].gibMirWas(defizit,entfernung,i,richtung,vorZurueck) or jemandErreicht;
inc(entfernung);
end;
@@ -586,22 +620,24 @@ begin
end;
end;
-procedure tImpulsPunkt.akkumuliereEMQuellen(var emQuellen: tEMQuellen; vx,vy,dVP: extended);
+procedure tImpulsPunkt.akkumuliereEMQuellen(var emQuellen: tEMQuellen; pX,dVP: extended);
var
i: longint;
tmp: extended;
begin
for i:=0 to length(werte)-1 do begin
- tmp:=werte[i,false] * raumpunkt.felder.spezLadungen[i] * dVP;
- emQuellen[eqRho]:=emQuellen[eqRho] + tmp;
- emQuellen[eqJX]:= emQuellen[eqJX] + tmp * vx;
- emQuellen[eqJY]:= emQuellen[eqJY] + tmp * vy;
+ iMGamma[i]:=1/(sqrt(1 + (sqr(pX) + sqr(raumpunkt.matFelder[i,mfPY,false]))*raumpunkt.felder.teilchen[i].eigenschaften[tsgIQdrMasse])*raumpunkt.felder.teilchen[i].eigenschaften[tsgMasse]);
+ tmp:=werte[i,false] * raumpunkt.felder.teilchen[i].eigenschaften[tsgLadung] * dVP;
+ emQuellen[eqRho]:=emQuellen[eqRho] + tmp; // * iGamma[i] ??? TODO ... falls, dann aber ohne M!
+ emQuellen[eqJX]:= emQuellen[eqJX] + tmp * pX * iMGamma[i];
+ emQuellen[eqJY]:= emQuellen[eqJY] + tmp * raumpunkt.matFelder[i,mfPY,false] * iMGamma[i];
end;
end;
-function tImpulsPunkt.gibMirWas(var defizit: extended; entfernung, teilchen,richtung: longint; positiv: boolean): boolean; // ich habe jemanden erreicht
+function tImpulsPunkt.gibMirWas(var defizit: extended; entfernung,teilchen,richtung: longint; positiv: boolean): boolean; // ich habe jemanden erreicht
var
i: longint;
+ b: boolean;
begin
result:=entfernung<=0;
if result then begin
@@ -612,17 +648,16 @@ begin
end
else begin
werte[teilchen,false]:=werte[teilchen,false]-defizit;
- defizit:=0;;
+ defizit:=0;
end;
end;
end
else
- if defizit>0 then begin
- result:=false;
- for i:=richtung to 4 do
- if assigned(nachbarn[i,positiv]) then
- result:=nachbarn[i,positiv].gibMirWas(defizit,entfernung-1,teilchen,i,positiv) or result;
- end;
+ if defizit>0 then
+ for i:=richtung to 1 do
+ for b:=false to i>richtung do
+ if assigned(nachbarn[i,positiv xor b]) then
+ result:=nachbarn[i,positiv xor b].gibMirWas(defizit,entfernung-1,teilchen,i,positiv xor b) or result;
end;
procedure tImpulsPunkt.setzeNull;
@@ -630,41 +665,61 @@ var
i: longint;
abl: boolean;
begin
- for i:=0 to length(werte)-1 do
+ for i:=0 to length(werte)-1 do begin
for abl:=false to true do
werte[i,abl]:=0;
+ iMGamma[i]:=0;
+ end;
end;
-procedure tImpulsPunkt.berechneAbleitungen(pX,pY,iDX,iDPX,iDPY: extended);
+procedure tImpulsPunkt.berechneAbleitungen(pX,iDX,iDPX: extended);
var
- igam: extended;
- i: longint;
+ i: longint;
+ tx,tp: extended;
begin
- igam:=1/sqrt(1+sqr(pX)+sqr(pY));
- pX:=pX*igam;
- pY:=pY*igam;
-
// df/dt = - v Nabla f - q(E + v/c x B) Nabla_p f
- for i:=0 to length(werte)-1 do
+
+ for i:=0 to length(werte)-1 do begin
+ tX:= // der Teil vor "Nabla f"
+ - pX * iMGamma[i];
+ tP:= // der Teil vor "Nabla_p f"
+ - raumpunkt.felder.teilchen[i].eigenschaften[tsgLadung] *
+ ( raumpunkt.emFelder[efEX,false]
+ + raumpunkt.matFelder[i,mfPY,false] * iMGamma[i] * raumpunkt.emFelder[efBZ,false]);
werte[i,true]:=
- - pX * (nachbarn[0,true].werte[i,false] - nachbarn[0,false].werte[i,false]) * iDX/2
- - raumpunkt.felder.spezLadungen[i] * (
- (raumpunkt.emFelder[efEX,false] + pY*raumpunkt.emFelder[efBZ,false])
- * (nachbarn[1,true].werte[i,false] - nachbarn[1,false].werte[i,false]) * iDPX/2
- -pX*raumpunkt.emFelder[efBZ,false]
- * (nachbarn[2,true].werte[i,false] - nachbarn[2,false].werte[i,false]) * iDPY/2
- );
+ tX * (
+ ( nachbarn[0,true].werte[i,false] - nachbarn[0,false].werte[i,false])
+ - sign(tX)*raumpunkt.felder.teilchen[i].diffusion[0] * (
+ nachbarn[0,true].werte[i,false] +
+ nachbarn[0,false].werte[i,false] -
+ 2*werte[i,false]
+ )
+ ) * iDX/2
+ + tP * (
+ ( nachbarn[1,true].werte[i,false] - nachbarn[1,false].werte[i,false])
+ - sign(tP)*raumpunkt.felder.teilchen[i].diffusion[1] * (
+ nachbarn[1,true].werte[i,false] +
+ nachbarn[1,false].werte[i,false] -
+ 2*werte[i,false]
+ )
+ )* iDPX/2;
+ end;
+
+ // Der jeweils zweite Summand entspricht (in etwa) einer thermischen
+ // Verschmierung und soll die numerische Stabilität bis zu
+ // d(ln n)/d[xp] * d[XP] = "diffusion" gewährleisten.
end;
// tRaumPunkt ******************************************************************
-constructor tRaumPunkt.create(linkerNachbar: tRaumPunkt; derBesitzer: tFelder; teilchen,_aPX,_aPY: longint; _dPX,_dPY: extended);
+constructor tRaumPunkt.create(linkerNachbar: tRaumPunkt; derBesitzer: tFelder; teilchen,_aP: longint; _dP: extended);
var
- emF: tEMFeldGroesze;
- emQ: tEMQuellGroesze;
- abl: boolean;
- i,j: longint;
- xN,pxN,pyN: tImpulsPunkt;
+ emF: tEMFeldGroesze;
+ emQ: tEMQuellGroesze;
+ maF: tMaterieFeldGroesze;
+ abl: boolean;
+ i: longint;
+ xN,pN: tImpulsPunkt;
begin
inherited create;
felder:=derBesitzer;
@@ -673,32 +728,31 @@ begin
for emF:=low(tEMFeldGroesze) to high(tEMFeldGroesze) do
for abl:=false to true do
emFelder[emF,abl]:=0;
+ fillchar(matFelder,sizeof(matFelder),#0);
+ setlength(matFelder,teilchen);
+ for i:=0 to length(matFelder)-1 do
+ for maF:=low(tMaterieFeldGroesze) to high(tMaterieFeldGroesze) do
+ for abl:=false to true do
+ matFelder[i,maF,abl]:=0;
lN:=linkerNachbar;
rN:=nil;
if assigned(lN) then
lN.rN:=self;
- aPX:=_aPX+byte(not odd(_aPX));
- aPY:=_aPY+byte(not odd(_aPY));
- dPX:=_dPX;
- dPY:=_dPY;
+ aP:=_aP+byte(not odd(_aP));
+ dP:=_dP;
fillchar(phasenraum,sizeof(phasenraum),#0);
- setlength(phasenraum,aPX*aPY);
- for i:=0 to aPX-1 do
- for j:=0 to aPY-1 do begin
- if assigned(linkerNachbar) then
- xN:=linkerNachbar.phasenraum[i+aPX*j]
- else
- xN:=nil;
- if i>0 then
- pxN:=phasenraum[i-1+aPX*j]
- else
- pxN:=nil;
- if j>0 then
- pyN:=phasenraum[i+aPX*(j-1)]
- else
- pyN:=nil;
- phasenraum[i+aPX*j]:=tImpulsPunkt.create(self,xN,pxN,pyN,teilchen);
- end;
+ setlength(phasenraum,aP);
+ for i:=0 to aP-1 do begin
+ if assigned(linkerNachbar) then
+ xN:=linkerNachbar.phasenraum[i]
+ else
+ xN:=nil;
+ if i>0 then
+ pN:=phasenraum[i-1]
+ else
+ pN:=nil;
+ phasenraum[i]:=tImpulsPunkt.create(self,xN,pN,teilchen);
+ end;
end;
destructor tRaumPunkt.destroy;
@@ -747,17 +801,19 @@ begin
(rN.emFelder[efAY,false]-lN.emFelder[efAY,false])*iDX/2;
end;
-procedure tRaumPunkt.berechneEMAbleitungen(iDX: extended); // Zeitableitungen der EM-Potentiale berechnen
+procedure tRaumPunkt.berechneEMQuellen;
var
- i,j: longint;
+ i: longint;
emQ: tEMQuellGroesze;
begin
for emQ:=low(tEMQuellGroesze) to high(tEMQuellGroesze) do
emQuellen[emQ]:=0;
- for i:=0 to aPX-1 do
- for j:=0 to aPY-1 do
- phasenraum[i+aPX*j].akkumuliereEMQuellen(emQuellen,(i-aPX/2)*dPX,(j-aPY/2)*dPY,dPX*dPY);
+ for i:=0 to aP-1 do
+ phasenraum[i].akkumuliereEMQuellen(emQuellen,(i-aP/2)*dP,dP);
+end;
+procedure tRaumPunkt.berechneEMAbleitungen(iDX: extended); // Zeitableitungen der EM-Potentiale berechnen
+begin
// d2A/dt2 = Laplace(A) - j ( - dE/dt wird auf dA/dt direkt aufgeschlagen !!! )
emFelder[efDAXDT,true]:=
( rN.emFelder[efAX,false] - 2*emFelder[efAX,false] + lN.emFelder[efAX,false] )*sqr(iDX)
@@ -773,25 +829,25 @@ end;
procedure tRaumPunkt.berechnePhasenraumAbleitungen(iDX: extended);
var
- i,j: longint;
+ i: longint;
begin
// df/dt = - v Nabla f - q(E + v/c x B) Nabla_p f
- for i:=1 to aPX-2 do
- for j:=1 to aPY-2 do
- phasenraum[i+aPX*j].berechneAbleitungen((i-(aPX div 2))*dPX,(j-(aPY div 2))*dPY,iDX,1/dPX,1/dPY);
+ for i:=1 to aP-2 do // Rand wird ausgelassen!
+ phasenraum[i].berechneAbleitungen((i-(aP div 2))*dP,iDX,1/dP);
end;
-procedure tRaumPunkt.nichtnegativieren; // Dichten nicht negativ machen
+function tRaumPunkt.nichtnegativieren: extended; // Dichten nicht negativ machen
var
- i: longint;
+ i: longint;
begin
+ result:=0;
for i:=0 to length(phasenraum)-1 do
- phasenraum[i].nichtnegativieren;
+ result:=result+phasenraum[i].nichtnegativieren;
end;
function tRaumPunkt.impulsIntegral(teilchen: longint; emQ: tEMQuellGroesze): extended;
var
- i,j: longint;
+ i: longint;
begin
if teilchen<0 then begin
result:=emQuellen[emQ]; // das ist leicht :-)
@@ -802,48 +858,71 @@ begin
case emQ of
eqRho:
- for i:=0 to aPX-1 do
- for j:=0 to aPY-1 do
- result:=
- result +
- phasenraum[i+aPX*j].werte[teilchen,false]*dPX*dPY;
- eqJX:
- for i:=0 to aPX-1 do
- for j:=0 to aPY-1 do
- result:=
- result +
- phasenraum[i+aPX*j].werte[teilchen,false]*dPX*dPY * (i-aPX/2)*dPX;
- eqJY:
- for i:=0 to aPX-1 do
- for j:=0 to aPY-1 do
- result:=
- result +
- phasenraum[i+aPX*j].werte[teilchen,false]*dPX*dPY * (j-aPY/2)*dPY;
+ for i:=0 to aP-1 do
+ result:=
+ result +
+ phasenraum[i].werte[teilchen,false]*dP;
+ eqJX: begin
+ for i:=0 to aP-1 do
+ result:=
+ result +
+ phasenraum[i].werte[teilchen,false]*dP * (i-aP/2)*dP
+ / sqrt(1 + (sqr((i-aP/2)*dP) + sqr(matFelder[teilchen,mfPY,false]))*felder.teilchen[teilchen].eigenschaften[tsgIQdrMasse]);
+ result:=result / felder.teilchen[teilchen].eigenschaften[tsgMasse];
+ end;
+ eqJY: begin
+ for i:=0 to aP-1 do
+ result:=
+ result +
+ phasenraum[i].werte[teilchen,false]*dP
+ / sqrt(1 + (sqr((i-aP/2)*dP) + sqr(matFelder[teilchen,mfPY,false]))*felder.teilchen[teilchen].eigenschaften[tsgIQdrMasse]);
+ result:=result * matFelder[teilchen,mfPY,false] / felder.teilchen[teilchen].eigenschaften[tsgMasse];
+ end;
end{of case};
+ result:=result * felder.teilchen[teilchen].eigenschaften[tsgLadung];
end;
-procedure tRaumPunkt.initialisiereDichte(teilchen: longint; n: extended);
+function tRaumPunkt.impulsIntegral(teilchen: longint): extended;
var
- i: longint;
+ i,j: longint;
begin
+ result:=0;
+
+ if teilchen<0 then begin
+ for j:=0 to length(matFelder)-1 do
+ for i:=0 to aP-1 do
+ result:=result + phasenraum[i].werte[j,false]*dP;
+ end
+ else
+ for i:=0 to aP-1 do
+ result:=result + phasenraum[i].werte[teilchen,false]*dP;
+end;
+
+procedure tRaumPunkt.initialisiereDichte(teilchen: longint; breite,n: extended);
+var
+ i: longint;
+ ges: extended;
+begin
+ ges:=0;
+ for i:=0 to length(phasenraum)-1 do begin
+ phasenraum[i].werte[teilchen,false]:=power(1/2,sqr((i-aP/2)*dP*2/breite));
+ ges:=ges + phasenraum[i].werte[teilchen,false];
+ end;
+ ges:=n/ges/dP;
for i:=0 to length(phasenraum)-1 do
- if (aPX div 2) + aPX*(aPY div 2) = i then
- phasenraum[i].werte[teilchen,false]:=n
- else
- phasenraum[i].werte[teilchen,false]:=0;
+ phasenraum[i].werte[teilchen,false]:=phasenraum[i].werte[teilchen,false] * ges;
end;
procedure tRaumPunkt.reflektiereTeilchen(rechts: boolean);
var
- i,j,k: longint;
+ i,j: longint;
begin
- for i:=(aPX div 2 + 1)*byte(rechts) to (aPX div 2)-1 + (apx-(aPX div 2))*byte(rechts) do
- for j:=0 to aPY-1 do
- for k:=0 to felder.matAnz-1 do
- if phasenraum[i+aPX*j].werte[k,false]>0 then begin
- phasenraum[aPX-1-i+aPX*j].werte[k,false]:=phasenraum[i+aPX*j].werte[k,false];
- phasenraum[i+aPX*j].werte[k,false]:=0;
- end;
+ for i:=(aP div 2 + 1)*byte(rechts) to (aP div 2)-1 + (aP-(aP div 2))*byte(rechts) do
+ for j:=0 to length(felder.teilchen)-1 do
+ if phasenraum[i].werte[j,false]>0 then begin
+ phasenraum[aP-1-i].werte[j,false]:=phasenraum[aP-1-i].werte[j,false]+phasenraum[i].werte[j,false];
+ phasenraum[i].werte[j,false]:=0;
+ end;
end;
procedure tRaumPunkt.setzeNull;
@@ -851,6 +930,7 @@ var
i: longint;
emF: tEMFeldGroesze;
emQ: tEMQuellGroesze;
+ maF: tMaterieFeldGroesze;
abl: boolean;
begin
for i:=0 to length(phasenraum)-1 do
@@ -860,6 +940,18 @@ begin
for emF:=low(tEMFeldGroesze) to high(tEMFeldGroesze) do
for abl:=false to true do
emFelder[emF,abl]:=0;
+ for i:=0 to length(matFelder)-1 do
+ for maF:=low(tMaterieFeldGroesze) to high(tMaterieFeldGroesze) do
+ for abl:=false to true do
+ matFelder[i,maF,abl]:=0;
+end;
+
+procedure tRaumPunkt.berechnePY;
+var
+ i: longint;
+begin
+ for i:=0 to length(matFelder)-1 do
+ matFelder[i,mfPY,false]:=-felder.teilchen[i].eigenschaften[tsgLadung] * emFelder[efAY,false];
end;
// linearkombinationsspezifische Methoden von tImpulsPunkt, tRaumPunkt und tFelder
@@ -870,7 +962,7 @@ end;
// tFelder *********************************************************************
-constructor tFelder.create(groesse: longint; teilchen: array of tTeilchenSpezies; lichter: tMyStringList; parent: tGitter; aPX,aPY: longint; dPX,dPY: extended);
+constructor tFelder.create(groesse: longint; _teilchen: array of tTeilchenSpezies; lichter: tMyStringList; parent: tGitter; aP: longint; dP: extended);
var
i,j: longint;
rechts: boolean;
@@ -878,25 +970,25 @@ var
begin
inherited create;
gitter:=parent;
- matAnz:=length(teilchen);
- fillchar(spezLadungen,sizeof(spezLadungen),#0);
- setlength(spezLadungen,matAnz);
- for i:=0 to length(spezLadungen)-1 do
- spezLadungen[i]:=teilchen[i].spezifischeLadung;
+ gesamtDefizit:=0;
+ fillchar(teilchen,sizeof(teilchen),#0);
+ setlength(teilchen,length(_teilchen));
+ for i:=0 to length(teilchen)-1 do
+ teilchen[i]:=tTeilchenSpezies.create(_teilchen[i]);
fillchar(inhalt,sizeof(inhalt),#0);
setlength(inhalt,groesse+4); // zwei Felder links und rechts extra für Randbedingungen
- inhalt[0]:=tRaumPunkt.create(nil,self,length(teilchen),aPX,aPY,dPX,dPY);
+ inhalt[0]:=tRaumPunkt.create(nil,self,length(teilchen),aP,dP);
for i:=1 to length(inhalt)-1 do
- inhalt[i]:=tRaumPunkt.create(inhalt[i-1],self,length(teilchen),aPX,aPY,dPX,dPY);
+ inhalt[i]:=tRaumPunkt.create(inhalt[i-1],self,length(teilchen),aP,dP);
fillchar(massen,sizeof(massen),#0);
setlength(massen,length(teilchen));
for i:=0 to length(massen)-1 do
massen[i]:=0;
for i:=0 to length(inhalt)-1 do
for j:=0 to length(teilchen)-1 do begin
- dens:=teilchen[j].gibDichte(gitter.xl+(i-2)*gitter.dX,parent.kvs);
- inhalt[i].initialisiereDichte(j,dens);
- massen[j]:=massen[j]+dens*gitter.dX;
+ dens:=_teilchen[j].gibDichte(gitter.xl+(i-2)*gitter.dX,parent.kvs);
+ inhalt[i].initialisiereDichte(j,_teilchen[j].breite,dens);
+ massen[j]:=massen[j]+dens*gitter.dX*teilchen[j].eigenschaften[tsgMasse];
end;
for rechts:=false to true do begin
lichters[rechts]:=tMyStringlist.create(nil,'');
@@ -912,7 +1004,10 @@ destructor tFelder.destroy;
var
i: longint;
begin
- setlength(spezLadungen,0);
+ setlength(massen,0);
+ for i:=0 to length(teilchen)-1 do
+ teilchen[i].free;
+ setlength(teilchen,0);
for i:=0 to length(inhalt)-1 do
inhalt[i].free;
setlength(inhalt,0);
@@ -948,19 +1043,26 @@ var
i: longint;
b: boolean;
begin
- for b:=false to true do
+ for b:=false to true do // PX ggf. korrigieren
inhalt[(length(inhalt)-1)*byte(b)].reflektiereTeilchen(b);
- inhalt[0].berechneEMFelder(dX,false);
+ for i:=0 to length(inhalt)-1 do // PY berechnen
+ inhalt[i].berechnePY;
+
+ for i:=0 to length(inhalt)-1 do // rho, j und nebenbei iGamma berechnen
+ inhalt[i].berechneEMQuellen;
+
+ inhalt[0].berechneEMFelder(dX,false); // E und B aus Rho und A berechnen
for i:=1 to length(inhalt)-2 do
inhalt[i].berechneEMFelder(dX,iDX);
inhalt[length(inhalt)-1].berechneEMFelder(dX,false);
- setzeRaender(iDX);
+ setzeRaender(iDX); // Randbedingungen für die Ableitungen des A-Felds setzen
- for i:=1 to length(inhalt)-2 do
+ for i:=1 to length(inhalt)-2 do // sonstige Ableitungen des A-Feldes berechnen
inhalt[i].berechneEMAbleitungen(iDX);
- for i:=1 to length(inhalt)-2 do
+
+ for i:=1 to length(inhalt)-2 do // Ableitungen der Phasenraumbesetzungsdichten berechnen (Rand wird ausgelassen!)
inhalt[i].berechnePhasenraumAbleitungen(iDX);
end;
@@ -972,9 +1074,27 @@ begin
inhalt[i].setzeNull;
end;
+procedure tFelder.dumpErhaltungsgroeszen(prot: tProtokollant);
+var
+ i,j: longint;
+ dens: extended;
+begin
+ dens:=0;
+ for i:=0 to length(massen)-1 do
+ dens:=dens+massen[i]/teilchen[i].eigenschaften[tsgMasse];
+ prot.schreibe('Gesamtdefizit: '+floattostr(gesamtDefizit)+' (Anteil '+floattostr(gesamtDefizit/dens)+')',true);
+ for i:=0 to length(massen)-1 do begin
+ dens:=0;
+ for j:=0 to length(inhalt)-1 do
+ dens:=dens+inhalt[j].impulsIntegral(i);
+ dens:=dens*teilchen[i].eigenschaften[tsgMasse]*gitter.dX;
+ prot.schreibe('n['+inttostr(i+1)+'] = '+floattostr(dens)+' (relative Abweichung: '+floattostr(dens/massen[i]-1)+')',true);
+ end;
+end;
+
// tGitter *********************************************************************
-constructor tGitter.create(derBesitzer: tSimulation; size: longint; deltaX: extended; bekannteWerte: tKnownValues; teilchen: array of tTeilchenSpezies; lichter: tMyStringlist; zv: tZeitverfahren; protokollant: tProtokollant; name: string; aPX,aPY: longint; dPX,dPY: extended);
+constructor tGitter.create(derBesitzer: tSimulation; size: longint; deltaX: extended; bekannteWerte: tKnownValues; teilchen: array of tTeilchenSpezies; lichter: tMyStringlist; zv: tZeitverfahren; protokollant: tProtokollant; name: string; aP: longint; dP: extended);
var
i: longint;
begin
@@ -1001,7 +1121,7 @@ begin
xl:=dX/2;
for i:=0 to length(felders)-1 do
- felders[i]:=tFelder.create(size,teilchen,lichter,self,aPX,aPY,dPX,dPY);
+ felders[i]:=tFelder.create(size,teilchen,lichter,self,aP,dP);
aktuelleFelder:=0;
t:=0;
@@ -1111,20 +1231,25 @@ begin
t:=t+dT;
end;
+procedure tGitter.dumpErhaltungsgroeszen;
+begin
+ felders[aktuelleFelder].dumpErhaltungsgroeszen(prot);
+end;
+
// tSimulation *****************************************************************
constructor tSimulation.create(inName: string; protokollant: tProtokollant; name: string);
var
- ifile: tMyStringlist;
- zeitverfahren: tZeitverfahren;
- s,t,aSuffix,aPrefix: string;
- deltaX,breite,dPX,dPY,pXMax,pYMax: extended;
- i: longint;
- kvs: tKnownValues;
- teilchen: array of tTeilchenSpezies;
- lichter: tMyStringlist;
- pro: tProtokollant;
- na: pSigActionRec;
+ ifile: tMyStringlist;
+ zeitverfahren: tZeitverfahren;
+ s,t,aSuffix,aPrefix: string;
+ dX,breite,dP,pMax: extended;
+ i: longint;
+ kvs: tKnownValues;
+ teilchen: array of tTeilchenSpezies;
+ lichter: tMyStringlist;
+ pro: tProtokollant;
+ na: pSigActionRec;
begin
inherited create;
prot:=tProtokollant.create(protokollant,name);
@@ -1140,12 +1265,12 @@ begin
// Standardeinstellungen Bereich 'allgemein'
zeitverfahren:=zfRungeKuttaVier;
- deltaX:=1e-2;
+ dX:=1e-2;
kvs.add('λ',1);
kvs.add('T',1);
kvs.add('ω',2*pi);
- kvs.add('dX',deltaX);
- kvs.add('q',1);
+ kvs.add('dX',dX);
+ kvs.add('qe',1);
kvs.add('me',1);
dT:=-1;
sDT:=-1;
@@ -1154,10 +1279,8 @@ begin
breite:=10.0;
fortschrittsAnzeige:=false;
gotSighup:=false;
- dPX:=-1;
- dPY:=-1;
- pXMax:=3;
- pYMax:=3;
+ dP:=-1;
+ pMax:=3;
// Standardeinstellungen Bereich 'ausgaben'
aPrefix:=extractfilepath(paramstr(0));
@@ -1213,8 +1336,8 @@ begin
continue;
end;
if startetMit('ortsschritt ',s) then begin
- deltaX:=exprtofloat(false,s,kvs,nil);
- kvs.add('dX',deltaX);
+ dX:=exprtofloat(false,s,kvs,nil);
+ kvs.add('dX',dX);
continue;
end;
if startetMit('zeitschritt ',s) then begin
@@ -1223,29 +1346,11 @@ begin
continue;
end;
if startetMit('impulsschritt ',s) then begin
- if startetMit('x ',s) then begin
- dPX:=exprtofloat(false,s,kvs,nil);
- continue;
- end;
- if startetMit('y ',s) then begin
- dPY:=exprtofloat(false,s,kvs,nil);
- continue;
- end;
- dPX:=exprtofloat(false,s,kvs,nil);
- dPY:=dPX;
+ dP:=exprtofloat(false,s,kvs,nil);
continue;
end;
if startetMit('maximalimpuls ',s) then begin
- if startetMit('x ',s) then begin
- pXMax:=exprtofloat(false,s,kvs,nil);
- continue;
- end;
- if startetMit('y ',s) then begin
- pYMax:=exprtofloat(false,s,kvs,nil);
- continue;
- end;
- pXMax:=exprtofloat(false,s,kvs,nil);
- pYMax:=pXMax;
+ pMax:=exprtofloat(false,s,kvs,nil);
continue;
end;
if startetMit('zeit ',s) then begin
@@ -1340,9 +1445,15 @@ begin
halt(1);
end;
if s='teilchen'+inttostr(length(teilchen))+'Ende' then break;
- if startetMit('spezifische Ladung ',s) then begin
- teilchen[length(teilchen)-1].spezifischeLadung:=-exprToFloat(false,s,kvs,nil);
- kvs.add('qm'+inttostr(length(teilchen)),teilchen[length(teilchen)-1].spezifischeLadung);
+ if startetMit('ladung ',s) then begin
+ teilchen[length(teilchen)-1].eigenschaften[tsgLadung]:=-exprToFloat(false,s,kvs,nil);
+ kvs.add('q'+inttostr(length(teilchen)),teilchen[length(teilchen)-1].eigenschaften[tsgLadung]);
+ continue;
+ end;
+ if startetMit('masse ',s) then begin
+ teilchen[length(teilchen)-1].eigenschaften[tsgMasse]:=exprToFloat(false,s,kvs,nil);
+ teilchen[length(teilchen)-1].eigenschaften[tsgIQdrMasse]:=1/sqr(teilchen[length(teilchen)-1].eigenschaften[tsgMasse]);
+ kvs.add('m'+inttostr(length(teilchen)),teilchen[length(teilchen)-1].eigenschaften[tsgMasse]);
continue;
end;
if startetMit('maximaldichte ',s) then begin
@@ -1359,6 +1470,18 @@ begin
teilchen[length(teilchen)-1].liesDichteFunktion(s,ifile,kvs,teilchen,prot);
continue;
end;
+ if startetMit('maximales dLnN/dX',s) then begin
+ teilchen[length(teilchen)-1].diffusion[0]:=exprToFloat(false,s,kvs,nil);
+ continue;
+ end;
+ if startetMit('maximales dLnN/dP',s) then begin
+ teilchen[length(teilchen)-1].diffusion[1]:=exprToFloat(false,s,kvs,nil);
+ continue;
+ end;
+ if startetMit('impulsbreite ',s) then begin
+ teilchen[length(teilchen)-1].breite:=exprToFloat(false,s,kvs,nil);
+ continue;
+ end;
prot.schreibe('Unbekannter Befehl '''+s+''' in Parameterdatei '''+inName+''' im Bereich teilchen!',true);
prot.destroyall;
@@ -1381,18 +1504,22 @@ begin
end;
if dT<0 then
- dT:=deltaX/10;
+ dT:=dX/10;
if sDT<0 then
sDT:=dT;
sDT:=1/round(1/sDT);
sT:=sDT/2;
if sDX<0 then
- sDX:=deltaX;
+ sDX:=dX;
+
+ if dP<0 then
+ dP:=dX;
+
+ for i:=0 to length(teilchen)-1 do begin
+ teilchen[i].diffusion[0]:=teilchen[i].diffusion[0] * dX;
+ teilchen[i].diffusion[1]:=teilchen[i].diffusion[1] * dP;
+ end;
- if dPX<0 then
- dPX:=deltaX;
- if dPY<0 then
- dPY:=sqrt(deltaX);
pro:=tProtokollant.create(prot,'create');
case zeitverfahren of
@@ -1411,13 +1538,11 @@ begin
else
pro.schreibe('Iteration mittels unbekanntem Verfahren');
end{of case};
- pro.schreibe('dX = '+floattostr(deltaX));
- pro.schreibe('dT = '+floattostr(dT));
- pro.schreibe('breite = '+floattostr(breite)+' (= '+inttostr(round(breite/deltaX)+1)+' Punkte)');
- pro.schreibe('dPX = '+floattostr(dPX));
- pro.schreibe('pXMax = '+floattostr(pXMax));
- pro.schreibe('dPY = '+floattostr(dPY));
- pro.schreibe('pYMax = '+floattostr(pYMax));
+ pro.schreibe('dX = '+floattostr(dX));
+ pro.schreibe('dT = '+floattostr(dT));
+ pro.schreibe('breite = '+floattostr(breite)+' (= '+inttostr(round(breite/dX)+1)+' Punkte)');
+ pro.schreibe('dP = '+floattostr(dP));
+ pro.schreibe('pMax = '+floattostr(pMax)+' (= '+inttostr(2*round(pMax/dP)+1)+' Punkte)');
pro.schreibe('bekannte Werte:');
kvs.dump(pro,' ');
if length(teilchen)>0 then begin
@@ -1456,7 +1581,7 @@ begin
setlength(sighupSimulationen,length(sighupSimulationen)+1);
sighupSimulationen[length(sighupSimulationen)-1]:=self;
- gitter:=tGitter.create(self,round(breite/deltaX)+1,deltaX,kvs,teilchen,lichter,zeitverfahren,prot,'gitter',2*round(pXMax/dPX)+1,2*round(pYMax/dPY)+1,dPX,dPY);
+ gitter:=tGitter.create(self,round(breite/dX)+1,dX,kvs,teilchen,lichter,zeitverfahren,prot,'gitter',2*round(pMax/dP)+1,dP);
for i:=0 to length(teilchen)-1 do
teilchen[i].free;
setlength(teilchen,0);
@@ -1509,6 +1634,7 @@ begin
prot.schreibe(timetostr(now-start)+' ('+floattostr(zeitPhysik/max(1e-11,zeitPhysik+zeitDatei))+')',true);
prot.schreibe('ETA: '+timetostr((now-start)*(tEnde-Gitter.t)/max(Gitter.t,dT)),true);
prot.schreibe('aktueller Zeitschritt: '+floattostr(dT)+'T',true);
+ gitter.dumpErhaltungsgroeszen;
end;
result:=gitter.t<tEnde;
diff --git a/Plasmapropagation.lpi b/Plasmapropagation.lpi
index fcb19b7..d9909e2 100644
--- a/Plasmapropagation.lpi
+++ b/Plasmapropagation.lpi
@@ -29,6 +29,7 @@
<RunParams>
<local>
<FormatVersion Value="1"/>
+ <CommandLineParams Value="input.plap"/>
<LaunchingApplication PathPlusParams="/usr/bin/screen -T 'Lazarus Run Output' -e $(LazarusDir)/tools/runwait.sh $(TargetCmdLine)"/>
</local>
</RunParams>
diff --git a/Plasmapropagation.lps b/Plasmapropagation.lps
index adb0f8b..d47f765 100644
--- a/Plasmapropagation.lps
+++ b/Plasmapropagation.lps
@@ -7,8 +7,7 @@
<Unit0>
<Filename Value="Plasmapropagation.lpr"/>
<IsPartOfProject Value="True"/>
- <TopLine Value="54"/>
- <CursorPos X="37" Y="34"/>
+ <CursorPos Y="5"/>
<UsageCount Value="200"/>
<Loaded Value="True"/>
</Unit0>
@@ -17,12 +16,12 @@
<IsPartOfProject Value="True"/>
<IsVisibleTab Value="True"/>
<EditorIndex Value="1"/>
- <TopLine Value="1376"/>
- <CursorPos Y="1396"/>
- <FoldState Value=" T3m00/5 pjdlc0A41212133131112 p0tP0601221[C4a07 T0\QPc07]dL"/>
- <UsageCount Value="155"/>
+ <TopLine Value="180"/>
+ <CursorPos Y="901"/>
+ <FoldState Value=" T3m50/5 pjdlc0G[945jZ0N4 piZjP0A2]9Yj90C7]RDjf0E5]IfjW0E[j4Ejg082]93jH031 pp4pD0G2 T0$2Pc071F311K"/>
+ <UsageCount Value="167"/>
<Bookmarks Count="1">
- <Item0 Y="1065"/>
+ <Item0 Y="1185"/>
</Bookmarks>
<Loaded Value="True"/>
</Unit1>
@@ -31,246 +30,244 @@
<IsPartOfProject Value="True"/>
<EditorIndex Value="-1"/>
<CursorPos X="15" Y="46"/>
- <UsageCount Value="117"/>
+ <UsageCount Value="129"/>
</Unit2>
<Unit3>
<Filename Value="input.plap"/>
<IsPartOfProject Value="True"/>
- <EditorIndex Value="4"/>
- <CursorPos X="44" Y="23"/>
- <UsageCount Value="116"/>
+ <EditorIndex Value="2"/>
+ <CursorPos X="20" Y="11"/>
+ <UsageCount Value="128"/>
<Loaded Value="True"/>
<DefaultSyntaxHighlighter Value="None"/>
</Unit3>
<Unit4>
<Filename Value="linearkombination.inc"/>
<IsPartOfProject Value="True"/>
- <EditorIndex Value="3"/>
- <TopLine Value="51"/>
- <CursorPos X="62" Y="77"/>
- <UsageCount Value="64"/>
- <Loaded Value="True"/>
+ <EditorIndex Value="-1"/>
+ <TopLine Value="286"/>
+ <CursorPos Y="324"/>
+ <UsageCount Value="76"/>
</Unit4>
<Unit5>
+ <Filename Value="linearkombinationen.inc"/>
+ <IsPartOfProject Value="True"/>
+ <EditorIndex Value="-1"/>
+ <TopLine Value="40"/>
+ <CursorPos X="234" Y="65"/>
+ <UsageCount Value="39"/>
+ </Unit5>
+ <Unit6>
<Filename Value="input.epost"/>
- <EditorIndex Value="5"/>
+ <EditorIndex Value="3"/>
<TopLine Value="22"/>
- <CursorPos X="36" Y="29"/>
- <UsageCount Value="92"/>
+ <CursorPos X="38" Y="29"/>
+ <UsageCount Value="98"/>
<Loaded Value="True"/>
<DefaultSyntaxHighlighter Value="None"/>
- </Unit5>
- <Unit6>
+ </Unit6>
+ <Unit7>
<Filename Value="../units/matheunit.pas"/>
<EditorIndex Value="-1"/>
<TopLine Value="366"/>
<CursorPos X="27" Y="390"/>
<FoldState Value=" T3i905B pj0jV034 piaj60U511+"/>
- <UsageCount Value="15"/>
- </Unit6>
- <Unit7>
- <Filename Value="../units/lowlevelunit.pas"/>
- <EditorIndex Value="-1"/>
- <TopLine Value="46"/>
- <CursorPos X="3" Y="12"/>
- <UsageCount Value="8"/>
+ <UsageCount Value="14"/>
</Unit7>
<Unit8>
+ <Filename Value="../units/lowlevelunit.pas"/>
+ <EditorIndex Value="-1"/>
+ <TopLine Value="549"/>
+ <CursorPos X="41" Y="567"/>
+ <UsageCount Value="9"/>
+ </Unit8>
+ <Unit9>
<Filename Value="../units/mystringlistunit.pas"/>
<EditorIndex Value="-1"/>
<TopLine Value="367"/>
<CursorPos X="17" Y="390"/>
<FoldState Value=" T3i3075 piZjD0WQ"/>
- <UsageCount Value="6"/>
- </Unit8>
- <Unit9>
+ <UsageCount Value="5"/>
+ </Unit9>
+ <Unit10>
<Filename Value="../epost/werteunit.pas"/>
<EditorIndex Value="-1"/>
<TopLine Value="950"/>
<CursorPos X="30" Y="1054"/>
- <UsageCount Value="3"/>
- </Unit9>
- <Unit10>
+ <UsageCount Value="2"/>
+ </Unit10>
+ <Unit11>
<Filename Value="../epost/typenunit.pas"/>
<EditorIndex Value="-1"/>
<TopLine Value="347"/>
<CursorPos X="62" Y="358"/>
- <UsageCount Value="3"/>
- </Unit10>
- <Unit11>
+ <UsageCount Value="2"/>
+ </Unit11>
+ <Unit12>
<Filename Value="../units/systemunit.pas"/>
<EditorIndex Value="-1"/>
<CursorPos X="3" Y="79"/>
- <UsageCount Value="3"/>
- </Unit11>
- <Unit12>
+ <UsageCount Value="2"/>
+ </Unit12>
+ <Unit13>
<Filename Value="/usr/lib/fpc/src/rtl/inc/objpash.inc"/>
<EditorIndex Value="-1"/>
<TopLine Value="232"/>
<CursorPos X="23" Y="192"/>
- <UsageCount Value="3"/>
- </Unit12>
- <Unit13>
- <Filename Value="rk14.inc"/>
- <EditorIndex Value="-1"/>
- <UsageCount Value="6"/>
+ <UsageCount Value="2"/>
</Unit13>
<Unit14>
- <Filename Value="rk1210.inc"/>
+ <Filename Value="rk14.inc"/>
<EditorIndex Value="-1"/>
- <TopLine Value="492"/>
- <CursorPos X="28" Y="565"/>
- <UsageCount Value="7"/>
+ <UsageCount Value="5"/>
</Unit14>
<Unit15>
- <Filename Value="rk1412.inc"/>
+ <Filename Value="rk1210.inc"/>
<EditorIndex Value="-1"/>
- <TopLine Value="945"/>
- <CursorPos X="5" Y="985"/>
- <UsageCount Value="7"/>
+ <TopLine Value="492"/>
+ <CursorPos X="28" Y="565"/>
+ <UsageCount Value="6"/>
</Unit15>
<Unit16>
- <Filename Value="rk108.inc"/>
+ <Filename Value="rk1412.inc"/>
<EditorIndex Value="-1"/>
- <CursorPos X="28" Y="8"/>
- <UsageCount Value="7"/>
+ <TopLine Value="755"/>
+ <CursorPos X="3" Y="775"/>
+ <UsageCount Value="9"/>
</Unit16>
<Unit17>
- <Filename Value="rk3_8.inc"/>
+ <Filename Value="rk108.inc"/>
<EditorIndex Value="-1"/>
- <CursorPos X="41" Y="6"/>
- <UsageCount Value="10"/>
+ <CursorPos X="28" Y="8"/>
+ <UsageCount Value="6"/>
</Unit17>
<Unit18>
- <Filename Value="linearkombinationen.inc"/>
- <IsPartOfProject Value="True"/>
- <EditorIndex Value="2"/>
- <TopLine Value="89"/>
- <CursorPos X="77" Y="102"/>
- <UsageCount Value="27"/>
- <Loaded Value="True"/>
+ <Filename Value="rk3_8.inc"/>
+ <EditorIndex Value="-1"/>
+ <CursorPos X="41" Y="6"/>
+ <UsageCount Value="9"/>
</Unit18>
</Units>
<JumpHistory Count="30" HistoryIndex="29">
<Position1>
<Filename Value="Physikunit.pas"/>
- <Caret Line="33" Column="34"/>
+ <Caret Line="680" TopLine="595"/>
</Position1>
<Position2>
<Filename Value="Physikunit.pas"/>
- <Caret Line="4" Column="37"/>
+ <Caret Line="78" Column="140" TopLine="55"/>
</Position2>
<Position3>
<Filename Value="Physikunit.pas"/>
- <Caret Line="207" Column="20" TopLine="184"/>
+ <Caret Line="1467" Column="88" TopLine="1438"/>
</Position3>
<Position4>
<Filename Value="Physikunit.pas"/>
- <Caret Line="601" TopLine="561"/>
+ <Caret Line="1463" Column="37" TopLine="1443"/>
</Position4>
<Position5>
<Filename Value="Physikunit.pas"/>
- <Caret Line="100" Column="41" TopLine="75"/>
+ <Caret Line="1469" Column="22" TopLine="1444"/>
</Position5>
<Position6>
<Filename Value="Physikunit.pas"/>
- <Caret Line="523" Column="82" TopLine="521"/>
+ <Caret Line="554" Column="31" TopLine="220"/>
</Position6>
<Position7>
<Filename Value="Physikunit.pas"/>
- <Caret Line="91" Column="73" TopLine="78"/>
+ <Caret Line="1479" Column="57" TopLine="1457"/>
</Position7>
<Position8>
<Filename Value="Physikunit.pas"/>
- <Caret Line="597" Column="78" TopLine="559"/>
+ <Caret Line="682" TopLine="336"/>
</Position8>
<Position9>
<Filename Value="Physikunit.pas"/>
- <Caret Line="646" TopLine="602"/>
+ <Caret Line="920" Column="83" TopLine="810"/>
</Position9>
<Position10>
<Filename Value="Physikunit.pas"/>
- <Caret Line="656" TopLine="628"/>
+ <Caret Line="912" Column="37" TopLine="855"/>
</Position10>
<Position11>
<Filename Value="Physikunit.pas"/>
- <Caret Line="707" Column="74" TopLine="687"/>
+ <Caret Line="137" Column="61" TopLine="120"/>
</Position11>
<Position12>
<Filename Value="Physikunit.pas"/>
- <Caret Line="80" TopLine="80"/>
+ <Caret Line="915" Column="54" TopLine="778"/>
</Position12>
<Position13>
<Filename Value="Physikunit.pas"/>
- <Caret Line="685" Column="31" TopLine="654"/>
+ <Caret Line="908" Column="68" TopLine="821"/>
</Position13>
<Position14>
<Filename Value="Physikunit.pas"/>
- <Caret Line="123" Column="111" TopLine="97"/>
+ <Caret Line="997" Column="52" TopLine="965"/>
</Position14>
<Position15>
<Filename Value="Physikunit.pas"/>
- <Caret Line="700" Column="60" TopLine="669"/>
+ <Caret Line="1484" Column="55" TopLine="1458"/>
</Position15>
<Position16>
<Filename Value="Physikunit.pas"/>
- <Caret Line="746" TopLine="658"/>
+ <Caret Line="75" TopLine="49"/>
</Position16>
<Position17>
<Filename Value="Physikunit.pas"/>
- <Caret Line="759" Column="93" TopLine="739"/>
+ <Caret Line="434" TopLine="353"/>
</Position17>
<Position18>
<Filename Value="Physikunit.pas"/>
- <Caret Line="772" TopLine="748"/>
+ <Caret Line="464" TopLine="428"/>
</Position18>
<Position19>
<Filename Value="Physikunit.pas"/>
- <Caret Line="781" Column="102" TopLine="761"/>
+ <Caret Line="1521" Column="60" TopLine="1379"/>
</Position19>
<Position20>
<Filename Value="Physikunit.pas"/>
- <Caret Line="794" Column="8" TopLine="784"/>
+ <Caret Line="175" TopLine="175"/>
</Position20>
<Position21>
<Filename Value="Physikunit.pas"/>
- <Caret Line="830" Column="37" TopLine="810"/>
+ <Caret Line="1247" Column="5" TopLine="1234"/>
</Position21>
<Position22>
<Filename Value="Physikunit.pas"/>
- <Caret Line="841" Column="20" TopLine="821"/>
+ <Caret Line="1269" Column="5" TopLine="1250"/>
</Position22>
<Position23>
<Filename Value="Physikunit.pas"/>
- <Caret Line="873" Column="149" TopLine="868"/>
+ <Caret Line="1273" Column="17" TopLine="1253"/>
</Position23>
<Position24>
<Filename Value="Physikunit.pas"/>
- <Caret Line="152" Column="127" TopLine="132"/>
+ <Caret Line="1341" Column="25" TopLine="1320"/>
</Position24>
<Position25>
<Filename Value="Physikunit.pas"/>
- <Caret Line="890" Column="83" TopLine="870"/>
+ <Caret Line="1517" Column="10" TopLine="1488"/>
</Position25>
<Position26>
<Filename Value="Physikunit.pas"/>
- <Caret Line="930" Column="24" TopLine="910"/>
+ <Caret Line="1520" Column="59" TopLine="1500"/>
</Position26>
<Position27>
<Filename Value="Physikunit.pas"/>
- <Caret Line="977" TopLine="966"/>
+ <Caret Line="1544" Column="77" TopLine="1522"/>
</Position27>
<Position28>
<Filename Value="Physikunit.pas"/>
- <Caret Line="177" Column="243" TopLine="157"/>
+ <Caret Line="1546" Column="63" TopLine="1539"/>
</Position28>
<Position29>
<Filename Value="Physikunit.pas"/>
- <Caret Line="1127" Column="38" TopLine="1100"/>
+ <Caret Line="1585" Column="51" TopLine="1539"/>
</Position29>
<Position30>
<Filename Value="Physikunit.pas"/>
- <Caret Line="1161" TopLine="1139"/>
+ <Caret Line="908" Column="67" TopLine="815"/>
</Position30>
</JumpHistory>
</ProjectSession>
diff --git a/input.epost b/input.epost
index c6f8d9a..96ec796 100644
--- a/input.epost
+++ b/input.epost
@@ -26,8 +26,8 @@ Threadanzahl: 11
parallel lesen
-!setze $symFelder: DPSIDX1 DPHIDX AY DAYDT
-!setze $asyFelder: N1
+!setze $symFelder: JX1 AY DAYDT # EX
+!setze $asyFelder: RHO1
!Schleife: $Feld: $asyFelder $symFelder
@@ -41,10 +41,10 @@ Ende
!Schleifenende
-!setze $symBerFelder: DPSIDX1_N1
+!setze $symBerFelder: VX1
-Multipliziere DPSIDX1 mal N1
- Name: DPSIDX1_N1
+Teile JX1 durch RHO1
+ Name: VX1
Ende
!Schleife: $symmetrie: symmetrisch asymmetrisch
diff --git a/input.plap b/input.plap
index 8e6fe6c..a6905dd 100644
--- a/input.plap
+++ b/input.plap
@@ -7,11 +7,11 @@ allgemein
runge-Kutta-3/8
# runge-Kutta-4
# euler-Vorwärts
- ortsschritt 10^-2 * λ
- zeitschritt 10^-2 * T
-# minZeitschritt 10^-2 * T
+ ortsschritt 5*10^-3 * λ
+ zeitschritt 5*10^-3 * T
+ maximalimpuls 10
+ impulsschritt 5 * 10^-3
zeit 40 * T
-# diffusionsterm 2 # max. p/Lp
!setze $breite: (5 * λ)
breite $breite
mit Fortschrittsanzeige
@@ -20,7 +20,7 @@ allgemeinEnde
ausgaben
prefix /home_raid/erich/Dokumente/Prograemmchen/Plasmapropagation/Daten/
suffix _test.dat
- felder AX,AY,dAYDT,EX,BY,BZ,Rho1,Rho2
+ felder AX,AY,dAYDT,EX,BY,BZ,Rho1,Rho2,JX1
ausgabenEnde
!setze $tFwhm: (2.5 * T)
@@ -31,9 +31,12 @@ licht von links 2 * 2^(-2*((t-$tMitte)/$tFwhm)^2) * ω * cos(ω*t) # Zei
!setze $IonenMassenFaktor: (1836.15267245 + 1838.68366158)
teilchen1
- spezifische Ladung -q/me
+ ladung -qe
+ masse me
maximaldichte 10
minimaldichte 0 # 10^-2
+ maximales dLnN/dX 10
+ maximales dLnN/dP 10
!setze $profilbreite: (4 * λ)
!setze $randbreite: (0.1 * λ)
verteilung stückweise
@@ -50,9 +53,12 @@ teilchen1
teilchen1Ende
teilchen2
- spezifische Ladung q/$IonenMassenFaktor/me
- maximaldichte nmax1*$IonenMassenFaktor
- minimaldichte nmin1*$IonenMassenFaktor
+ ladung qe
+ masse $IonenMassenFaktor * me
+ maximaldichte nmax1
+ minimaldichte nmin1
+ maximales dLnN/dX 10
+ maximales dLnN/dP 10
# verteilung wie teilchen1
verteilung stückweise
0
diff --git a/linearkombination.inc b/linearkombination.inc
index 0f74497..85ed695 100644
--- a/linearkombination.inc
+++ b/linearkombination.inc
@@ -320,6 +320,6 @@ begin
{$IFDEF lkA32},fak32
{$ENDIF}{$ENDIF}{$ENDIF}{$ENDIF}{$ENDIF}{$ENDIF}{$ENDIF}{$ENDIF}{$ENDIF}{$ENDIF}{$ENDIF}{$ENDIF}{$ENDIF}{$ENDIF}{$ENDIF}{$ENDIF}{$ENDIF}{$ENDIF}{$ENDIF}{$ENDIF});
for i:=0 to length(inhalt)-1 do
- inhalt[i].nichtnegativieren;
+ gesamtDefizit:=gesamtDefizit+inhalt[i].nichtnegativieren;
end;
diff --git a/linearkombinationen.inc b/linearkombinationen.inc
index 724c695..6717141 100644
--- a/linearkombinationen.inc
+++ b/linearkombinationen.inc
@@ -2,51 +2,51 @@
// interface *******************************************************************
{$IFDEF LiKotImpulsPunktHeader}
-procedure liKo(in1,in2: tImpulsPunkt; fak2: extended); overload; // Werte werden auf (in1 + \sum_i faki * ini') gesetzt
-procedure liKo(in1,in2,in3: tImpulsPunkt; fak2,fak3: extended); overload;
-procedure liKo(in1,in2,in3,in4: tImpulsPunkt; fak2,fak3,fak4: extended); overload;
-procedure liKo(in1,in2,in3,in4,in5: tImpulsPunkt; fak2,fak3,fak4,fak5: extended); overload;
-procedure liKo(in1,in2,in3,in4,in5,in6: tImpulsPunkt; fak2,fak3,fak4,fak5,fak6: extended); overload;
-procedure liKo(in1,in2,in3,in4,in5,in6,in7: tImpulsPunkt; fak2,fak3,fak4,fak5,fak6,fak7: extended); overload;
-procedure liKo(in1,in2,in3,in4,in5,in6,in7,in8: tImpulsPunkt; fak2,fak3,fak4,fak5,fak6,fak7,fak8: extended); overload;
-procedure liKo(in1,in2,in3,in4,in5,in6,in7,in8,in9: tImpulsPunkt; fak2,fak3,fak4,fak5,fak6,fak7,fak8,fak9: extended); overload;
-procedure liKo(in1,in2,in3,in4,in5,in6,in7,in8,in9,in10: tImpulsPunkt; fak2,fak3,fak4,fak5,fak6,fak7,fak8,fak9,fak10: extended); overload;
-procedure liKo(in1,in2,in3,in4,in5,in6,in7,in8,in9,in10,in11: tImpulsPunkt; fak2,fak3,fak4,fak5,fak6,fak7,fak8,fak9,fak10,fak11: extended); overload;
-procedure liKo(in1,in2,in3,in4,in5,in6,in7,in8,in9,in10,in11,in12: tImpulsPunkt; fak2,fak3,fak4,fak5,fak6,fak7,fak8,fak9,fak10,fak11,fak12: extended); overload;
-procedure liKo(in1,in2,in3,in4,in5,in6,in7,in8,in9,in10,in11,in12,in13,in14: tImpulsPunkt; fak2,fak3,fak4,fak5,fak6,fak7,fak8,fak9,fak10,fak11,fak12,fak13,fak14: extended); overload;
-procedure liKo(in1,in2,in3,in4,in5,in6,in7,in8,in9,in10,in11,in12,in13,in14,in15: tImpulsPunkt; fak2,fak3,fak4,fak5,fak6,fak7,fak8,fak9,fak10,fak11,fak12,fak13,fak14,fak15: extended); overload;
-procedure liKo(in1,in2,in3,in4,in5,in6,in7,in8,in9,in10,in11,in12,in13,in14,in15,in16: tImpulsPunkt; fak2,fak3,fak4,fak5,fak6,fak7,fak8,fak9,fak10,fak11,fak12,fak13,fak14,fak15,fak16: extended); overload;
-procedure liKo(in1,in2,in3,in4,in5,in6,in7,in8,in9,in10,in11,in12,in13,in14,in15,in16,in17: tImpulsPunkt; fak2,fak3,fak4,fak5,fak6,fak7,fak8,fak9,fak10,fak11,fak12,fak13,fak14,fak15,fak16,fak17: extended); overload;
-procedure liKo(in1,in2,in3,in4,in5,in6,in7,in8,in9,in10,in11,in12,in13,in14,in15,in16,in17,in18: tImpulsPunkt; fak2,fak3,fak4,fak5,fak6,fak7,fak8,fak9,fak10,fak11,fak12,fak13,fak14,fak15,fak16,fak17,fak18: extended); overload;
-procedure liKo(in1,in2,in3,in4,in5,in6,in7,in8,in9,in10,in11,in12,in13,in14,in15,in16,in17,in18,in19: tImpulsPunkt; fak2,fak3,fak4,fak5,fak6,fak7,fak8,fak9,fak10,fak11,fak12,fak13,fak14,fak15,fak16,fak17,fak18,fak19: extended); overload;
-procedure liKo(in1,in2,in3,in4,in5,in6,in7,in8,in9,in10,in11,in12,in13,in14,in15,in16,in17,in18,in19,in20,in21,in22: tImpulsPunkt; fak2,fak3,fak4,fak5,fak6,fak7,fak8,fak9,fak10,fak11,fak12,fak13,fak14,fak15,fak16,fak17,fak18,fak19,fak20,fak21,fak22: extended); overload;
-procedure liKo(in1,in2,in3,in4,in5,in6,in7,in8,in9,in10,in11,in12,in13,in14,in15,in16,in17,in18,in19,in20,in21,in22,in23: tImpulsPunkt; fak2,fak3,fak4,fak5,fak6,fak7,fak8,fak9,fak10,fak11,fak12,fak13,fak14,fak15,fak16,fak17,fak18,fak19,fak20,fak21,fak22,fak23: extended); overload;
-procedure liKo(in1,in2,in3,in4,in5,in6,in7,in8,in9,in10,in11,in12,in13,in14,in15,in16,in17,in18,in19,in20,in21,in22,in23,in24,in25,in26,in27,in28,in29,in30,in31: tImpulsPunkt; fak2,fak3,fak4,fak5,fak6,fak7,fak8,fak9,fak10,fak11,fak12,fak13,fak14,fak15,fak16,fak17,fak18,fak19,fak20,fak21,fak22,fak23,fak24,fak25,fak26,fak27,fak28,fak29,fak30,fak31: extended); overload;
-procedure liKo(in1,in2,in3,in4,in5,in6,in7,in8,in9,in10,in11,in12,in13,in14,in15,in16,in17,in18,in19,in20,in21,in22,in23,in24,in25,in26,in27,in28,in29,in30,in31,in32: tImpulsPunkt; fak2,fak3,fak4,fak5,fak6,fak7,fak8,fak9,fak10,fak11,fak12,fak13,fak14,fak15,fak16,fak17,fak18,fak19,fak20,fak21,fak22,fak23,fak24,fak25,fak26,fak27,fak28,fak29,fak30,fak31,fak32: extended); overload;
+procedure liKo(in1,in2: tImpulsPunkt; fak2: extended); overload; inline; // Werte werden auf (in1 + \sum_i faki * ini') gesetzt
+procedure liKo(in1,in2,in3: tImpulsPunkt; fak2,fak3: extended); overload; inline;
+procedure liKo(in1,in2,in3,in4: tImpulsPunkt; fak2,fak3,fak4: extended); overload; inline;
+procedure liKo(in1,in2,in3,in4,in5: tImpulsPunkt; fak2,fak3,fak4,fak5: extended); overload; inline;
+procedure liKo(in1,in2,in3,in4,in5,in6: tImpulsPunkt; fak2,fak3,fak4,fak5,fak6: extended); overload; inline;
+procedure liKo(in1,in2,in3,in4,in5,in6,in7: tImpulsPunkt; fak2,fak3,fak4,fak5,fak6,fak7: extended); overload; inline;
+procedure liKo(in1,in2,in3,in4,in5,in6,in7,in8: tImpulsPunkt; fak2,fak3,fak4,fak5,fak6,fak7,fak8: extended); overload; inline;
+procedure liKo(in1,in2,in3,in4,in5,in6,in7,in8,in9: tImpulsPunkt; fak2,fak3,fak4,fak5,fak6,fak7,fak8,fak9: extended); overload; inline;
+procedure liKo(in1,in2,in3,in4,in5,in6,in7,in8,in9,in10: tImpulsPunkt; fak2,fak3,fak4,fak5,fak6,fak7,fak8,fak9,fak10: extended); overload; inline;
+procedure liKo(in1,in2,in3,in4,in5,in6,in7,in8,in9,in10,in11: tImpulsPunkt; fak2,fak3,fak4,fak5,fak6,fak7,fak8,fak9,fak10,fak11: extended); overload; inline;
+procedure liKo(in1,in2,in3,in4,in5,in6,in7,in8,in9,in10,in11,in12: tImpulsPunkt; fak2,fak3,fak4,fak5,fak6,fak7,fak8,fak9,fak10,fak11,fak12: extended); overload; inline;
+procedure liKo(in1,in2,in3,in4,in5,in6,in7,in8,in9,in10,in11,in12,in13,in14: tImpulsPunkt; fak2,fak3,fak4,fak5,fak6,fak7,fak8,fak9,fak10,fak11,fak12,fak13,fak14: extended); overload; inline;
+procedure liKo(in1,in2,in3,in4,in5,in6,in7,in8,in9,in10,in11,in12,in13,in14,in15: tImpulsPunkt; fak2,fak3,fak4,fak5,fak6,fak7,fak8,fak9,fak10,fak11,fak12,fak13,fak14,fak15: extended); overload; inline;
+procedure liKo(in1,in2,in3,in4,in5,in6,in7,in8,in9,in10,in11,in12,in13,in14,in15,in16: tImpulsPunkt; fak2,fak3,fak4,fak5,fak6,fak7,fak8,fak9,fak10,fak11,fak12,fak13,fak14,fak15,fak16: extended); overload; inline;
+procedure liKo(in1,in2,in3,in4,in5,in6,in7,in8,in9,in10,in11,in12,in13,in14,in15,in16,in17: tImpulsPunkt; fak2,fak3,fak4,fak5,fak6,fak7,fak8,fak9,fak10,fak11,fak12,fak13,fak14,fak15,fak16,fak17: extended); overload; inline;
+procedure liKo(in1,in2,in3,in4,in5,in6,in7,in8,in9,in10,in11,in12,in13,in14,in15,in16,in17,in18: tImpulsPunkt; fak2,fak3,fak4,fak5,fak6,fak7,fak8,fak9,fak10,fak11,fak12,fak13,fak14,fak15,fak16,fak17,fak18: extended); overload; inline;
+procedure liKo(in1,in2,in3,in4,in5,in6,in7,in8,in9,in10,in11,in12,in13,in14,in15,in16,in17,in18,in19: tImpulsPunkt; fak2,fak3,fak4,fak5,fak6,fak7,fak8,fak9,fak10,fak11,fak12,fak13,fak14,fak15,fak16,fak17,fak18,fak19: extended); overload; inline;
+procedure liKo(in1,in2,in3,in4,in5,in6,in7,in8,in9,in10,in11,in12,in13,in14,in15,in16,in17,in18,in19,in20,in21,in22: tImpulsPunkt; fak2,fak3,fak4,fak5,fak6,fak7,fak8,fak9,fak10,fak11,fak12,fak13,fak14,fak15,fak16,fak17,fak18,fak19,fak20,fak21,fak22: extended); overload; inline;
+procedure liKo(in1,in2,in3,in4,in5,in6,in7,in8,in9,in10,in11,in12,in13,in14,in15,in16,in17,in18,in19,in20,in21,in22,in23: tImpulsPunkt; fak2,fak3,fak4,fak5,fak6,fak7,fak8,fak9,fak10,fak11,fak12,fak13,fak14,fak15,fak16,fak17,fak18,fak19,fak20,fak21,fak22,fak23: extended); overload; inline;
+procedure liKo(in1,in2,in3,in4,in5,in6,in7,in8,in9,in10,in11,in12,in13,in14,in15,in16,in17,in18,in19,in20,in21,in22,in23,in24,in25,in26,in27,in28,in29,in30,in31: tImpulsPunkt; fak2,fak3,fak4,fak5,fak6,fak7,fak8,fak9,fak10,fak11,fak12,fak13,fak14,fak15,fak16,fak17,fak18,fak19,fak20,fak21,fak22,fak23,fak24,fak25,fak26,fak27,fak28,fak29,fak30,fak31: extended); overload; inline;
+procedure liKo(in1,in2,in3,in4,in5,in6,in7,in8,in9,in10,in11,in12,in13,in14,in15,in16,in17,in18,in19,in20,in21,in22,in23,in24,in25,in26,in27,in28,in29,in30,in31,in32: tImpulsPunkt; fak2,fak3,fak4,fak5,fak6,fak7,fak8,fak9,fak10,fak11,fak12,fak13,fak14,fak15,fak16,fak17,fak18,fak19,fak20,fak21,fak22,fak23,fak24,fak25,fak26,fak27,fak28,fak29,fak30,fak31,fak32: extended); overload; inline;
{$ENDIF}
{$IFDEF LiKotRaumPunktHeader}
-procedure liKo(in1,in2: tRaumPunkt; fak2: extended); overload; // Werte werden auf (in1 + \sum_i faki * ini') gesetzt
-procedure liKo(in1,in2,in3: tRaumPunkt; fak2,fak3: extended); overload;
-procedure liKo(in1,in2,in3,in4: tRaumPunkt; fak2,fak3,fak4: extended); overload;
-procedure liKo(in1,in2,in3,in4,in5: tRaumPunkt; fak2,fak3,fak4,fak5: extended); overload;
-procedure liKo(in1,in2,in3,in4,in5,in6: tRaumPunkt; fak2,fak3,fak4,fak5,fak6: extended); overload;
-procedure liKo(in1,in2,in3,in4,in5,in6,in7: tRaumPunkt; fak2,fak3,fak4,fak5,fak6,fak7: extended); overload;
-procedure liKo(in1,in2,in3,in4,in5,in6,in7,in8: tRaumPunkt; fak2,fak3,fak4,fak5,fak6,fak7,fak8: extended); overload;
-procedure liKo(in1,in2,in3,in4,in5,in6,in7,in8,in9: tRaumPunkt; fak2,fak3,fak4,fak5,fak6,fak7,fak8,fak9: extended); overload;
-procedure liKo(in1,in2,in3,in4,in5,in6,in7,in8,in9,in10: tRaumPunkt; fak2,fak3,fak4,fak5,fak6,fak7,fak8,fak9,fak10: extended); overload;
-procedure liKo(in1,in2,in3,in4,in5,in6,in7,in8,in9,in10,in11: tRaumPunkt; fak2,fak3,fak4,fak5,fak6,fak7,fak8,fak9,fak10,fak11: extended); overload;
-procedure liKo(in1,in2,in3,in4,in5,in6,in7,in8,in9,in10,in11,in12: tRaumPunkt; fak2,fak3,fak4,fak5,fak6,fak7,fak8,fak9,fak10,fak11,fak12: extended); overload;
-procedure liKo(in1,in2,in3,in4,in5,in6,in7,in8,in9,in10,in11,in12,in13,in14: tRaumPunkt; fak2,fak3,fak4,fak5,fak6,fak7,fak8,fak9,fak10,fak11,fak12,fak13,fak14: extended); overload;
-procedure liKo(in1,in2,in3,in4,in5,in6,in7,in8,in9,in10,in11,in12,in13,in14,in15: tRaumPunkt; fak2,fak3,fak4,fak5,fak6,fak7,fak8,fak9,fak10,fak11,fak12,fak13,fak14,fak15: extended); overload;
-procedure liKo(in1,in2,in3,in4,in5,in6,in7,in8,in9,in10,in11,in12,in13,in14,in15,in16: tRaumPunkt; fak2,fak3,fak4,fak5,fak6,fak7,fak8,fak9,fak10,fak11,fak12,fak13,fak14,fak15,fak16: extended); overload;
-procedure liKo(in1,in2,in3,in4,in5,in6,in7,in8,in9,in10,in11,in12,in13,in14,in15,in16,in17: tRaumPunkt; fak2,fak3,fak4,fak5,fak6,fak7,fak8,fak9,fak10,fak11,fak12,fak13,fak14,fak15,fak16,fak17: extended); overload;
-procedure liKo(in1,in2,in3,in4,in5,in6,in7,in8,in9,in10,in11,in12,in13,in14,in15,in16,in17,in18: tRaumPunkt; fak2,fak3,fak4,fak5,fak6,fak7,fak8,fak9,fak10,fak11,fak12,fak13,fak14,fak15,fak16,fak17,fak18: extended); overload;
-procedure liKo(in1,in2,in3,in4,in5,in6,in7,in8,in9,in10,in11,in12,in13,in14,in15,in16,in17,in18,in19: tRaumPunkt; fak2,fak3,fak4,fak5,fak6,fak7,fak8,fak9,fak10,fak11,fak12,fak13,fak14,fak15,fak16,fak17,fak18,fak19: extended); overload;
-procedure liKo(in1,in2,in3,in4,in5,in6,in7,in8,in9,in10,in11,in12,in13,in14,in15,in16,in17,in18,in19,in20,in21,in22: tRaumPunkt; fak2,fak3,fak4,fak5,fak6,fak7,fak8,fak9,fak10,fak11,fak12,fak13,fak14,fak15,fak16,fak17,fak18,fak19,fak20,fak21,fak22: extended); overload;
-procedure liKo(in1,in2,in3,in4,in5,in6,in7,in8,in9,in10,in11,in12,in13,in14,in15,in16,in17,in18,in19,in20,in21,in22,in23: tRaumPunkt; fak2,fak3,fak4,fak5,fak6,fak7,fak8,fak9,fak10,fak11,fak12,fak13,fak14,fak15,fak16,fak17,fak18,fak19,fak20,fak21,fak22,fak23: extended); overload;
-procedure liKo(in1,in2,in3,in4,in5,in6,in7,in8,in9,in10,in11,in12,in13,in14,in15,in16,in17,in18,in19,in20,in21,in22,in23,in24,in25,in26,in27,in28,in29,in30,in31: tRaumPunkt; fak2,fak3,fak4,fak5,fak6,fak7,fak8,fak9,fak10,fak11,fak12,fak13,fak14,fak15,fak16,fak17,fak18,fak19,fak20,fak21,fak22,fak23,fak24,fak25,fak26,fak27,fak28,fak29,fak30,fak31: extended); overload;
-procedure liKo(in1,in2,in3,in4,in5,in6,in7,in8,in9,in10,in11,in12,in13,in14,in15,in16,in17,in18,in19,in20,in21,in22,in23,in24,in25,in26,in27,in28,in29,in30,in31,in32: tRaumPunkt; fak2,fak3,fak4,fak5,fak6,fak7,fak8,fak9,fak10,fak11,fak12,fak13,fak14,fak15,fak16,fak17,fak18,fak19,fak20,fak21,fak22,fak23,fak24,fak25,fak26,fak27,fak28,fak29,fak30,fak31,fak32: extended); overload;
+procedure liKo(in1,in2: tRaumPunkt; fak2: extended); overload; inline; // Werte werden auf (in1 + \sum_i faki * ini') gesetzt
+procedure liKo(in1,in2,in3: tRaumPunkt; fak2,fak3: extended); overload; inline;
+procedure liKo(in1,in2,in3,in4: tRaumPunkt; fak2,fak3,fak4: extended); overload; inline;
+procedure liKo(in1,in2,in3,in4,in5: tRaumPunkt; fak2,fak3,fak4,fak5: extended); overload; inline;
+procedure liKo(in1,in2,in3,in4,in5,in6: tRaumPunkt; fak2,fak3,fak4,fak5,fak6: extended); overload; inline;
+procedure liKo(in1,in2,in3,in4,in5,in6,in7: tRaumPunkt; fak2,fak3,fak4,fak5,fak6,fak7: extended); overload; inline;
+procedure liKo(in1,in2,in3,in4,in5,in6,in7,in8: tRaumPunkt; fak2,fak3,fak4,fak5,fak6,fak7,fak8: extended); overload; inline;
+procedure liKo(in1,in2,in3,in4,in5,in6,in7,in8,in9: tRaumPunkt; fak2,fak3,fak4,fak5,fak6,fak7,fak8,fak9: extended); overload; inline;
+procedure liKo(in1,in2,in3,in4,in5,in6,in7,in8,in9,in10: tRaumPunkt; fak2,fak3,fak4,fak5,fak6,fak7,fak8,fak9,fak10: extended); overload; inline;
+procedure liKo(in1,in2,in3,in4,in5,in6,in7,in8,in9,in10,in11: tRaumPunkt; fak2,fak3,fak4,fak5,fak6,fak7,fak8,fak9,fak10,fak11: extended); overload; inline;
+procedure liKo(in1,in2,in3,in4,in5,in6,in7,in8,in9,in10,in11,in12: tRaumPunkt; fak2,fak3,fak4,fak5,fak6,fak7,fak8,fak9,fak10,fak11,fak12: extended); overload; inline;
+procedure liKo(in1,in2,in3,in4,in5,in6,in7,in8,in9,in10,in11,in12,in13,in14: tRaumPunkt; fak2,fak3,fak4,fak5,fak6,fak7,fak8,fak9,fak10,fak11,fak12,fak13,fak14: extended); overload; inline;
+procedure liKo(in1,in2,in3,in4,in5,in6,in7,in8,in9,in10,in11,in12,in13,in14,in15: tRaumPunkt; fak2,fak3,fak4,fak5,fak6,fak7,fak8,fak9,fak10,fak11,fak12,fak13,fak14,fak15: extended); overload; inline;
+procedure liKo(in1,in2,in3,in4,in5,in6,in7,in8,in9,in10,in11,in12,in13,in14,in15,in16: tRaumPunkt; fak2,fak3,fak4,fak5,fak6,fak7,fak8,fak9,fak10,fak11,fak12,fak13,fak14,fak15,fak16: extended); overload; inline;
+procedure liKo(in1,in2,in3,in4,in5,in6,in7,in8,in9,in10,in11,in12,in13,in14,in15,in16,in17: tRaumPunkt; fak2,fak3,fak4,fak5,fak6,fak7,fak8,fak9,fak10,fak11,fak12,fak13,fak14,fak15,fak16,fak17: extended); overload; inline;
+procedure liKo(in1,in2,in3,in4,in5,in6,in7,in8,in9,in10,in11,in12,in13,in14,in15,in16,in17,in18: tRaumPunkt; fak2,fak3,fak4,fak5,fak6,fak7,fak8,fak9,fak10,fak11,fak12,fak13,fak14,fak15,fak16,fak17,fak18: extended); overload; inline;
+procedure liKo(in1,in2,in3,in4,in5,in6,in7,in8,in9,in10,in11,in12,in13,in14,in15,in16,in17,in18,in19: tRaumPunkt; fak2,fak3,fak4,fak5,fak6,fak7,fak8,fak9,fak10,fak11,fak12,fak13,fak14,fak15,fak16,fak17,fak18,fak19: extended); overload; inline;
+procedure liKo(in1,in2,in3,in4,in5,in6,in7,in8,in9,in10,in11,in12,in13,in14,in15,in16,in17,in18,in19,in20,in21,in22: tRaumPunkt; fak2,fak3,fak4,fak5,fak6,fak7,fak8,fak9,fak10,fak11,fak12,fak13,fak14,fak15,fak16,fak17,fak18,fak19,fak20,fak21,fak22: extended); overload; inline;
+procedure liKo(in1,in2,in3,in4,in5,in6,in7,in8,in9,in10,in11,in12,in13,in14,in15,in16,in17,in18,in19,in20,in21,in22,in23: tRaumPunkt; fak2,fak3,fak4,fak5,fak6,fak7,fak8,fak9,fak10,fak11,fak12,fak13,fak14,fak15,fak16,fak17,fak18,fak19,fak20,fak21,fak22,fak23: extended); overload; inline;
+procedure liKo(in1,in2,in3,in4,in5,in6,in7,in8,in9,in10,in11,in12,in13,in14,in15,in16,in17,in18,in19,in20,in21,in22,in23,in24,in25,in26,in27,in28,in29,in30,in31: tRaumPunkt; fak2,fak3,fak4,fak5,fak6,fak7,fak8,fak9,fak10,fak11,fak12,fak13,fak14,fak15,fak16,fak17,fak18,fak19,fak20,fak21,fak22,fak23,fak24,fak25,fak26,fak27,fak28,fak29,fak30,fak31: extended); overload; inline;
+procedure liKo(in1,in2,in3,in4,in5,in6,in7,in8,in9,in10,in11,in12,in13,in14,in15,in16,in17,in18,in19,in20,in21,in22,in23,in24,in25,in26,in27,in28,in29,in30,in31,in32: tRaumPunkt; fak2,fak3,fak4,fak5,fak6,fak7,fak8,fak9,fak10,fak11,fak12,fak13,fak14,fak15,fak16,fak17,fak18,fak19,fak20,fak21,fak22,fak23,fak24,fak25,fak26,fak27,fak28,fak29,fak30,fak31,fak32: extended); overload; inline;
{$ENDIF}
{$IFDEF LiKotFelderHeader}