1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
|
program Plasmapropagation;
{$mode objfpc}{$H+}
{$INCLUDE Compilerswitches.inc}
uses
{$IFDEF UNIX}{$IFDEF UseCThreads}
cthreads,
{$ENDIF}{$ENDIF}
Classes, SysUtils, CustApp,
{ you can add units after this }
math, Physikunit, protokollunit, lowlevelunit, fftw,
compilerswitchupdateunit;
type
{ TPlasmapropagation }
TPlasmapropagation = class(TCustomApplication)
protected
procedure DoRun; override;
public
constructor Create(TheOwner: TComponent); override;
destructor Destroy; override;
end;
{ TPlasmapropagation }
procedure TPlasmapropagation.DoRun;
var
simulation: tSimulation;
start,zeitPhysik,zeitDatei: double;
prot: tProtokollant;
s,t,u: string;
begin
prot:=tProtokollant.create('error');
if paramcount<>1 then begin
prot.schreibe('Bitte genau einen Parameter übergeben, nämlich die Parameterdatei!',true);
prot.destroyAll;
halt(1);
end;
start:=now;
s:=extractfilepath(paramstr(0))+'fftWisdom.dat';
if fileexists(s) then
if fftw_import_wisdom_from_filename(pChar(s))=0 then begin
prot.schreibe('Kann vorhandene Weisheitsdatei '''+s+''' nicht lesen!');
prot.destroyAll;
halt(1);
end;
zeitDatei:=0;
zeitPhysik:=0;
simulation:=tSimulation.create(paramstr(1),prot,'simulation');
if fftw_export_wisdom_to_filename(pChar(s))=0 then begin
prot.schreibe('Kann Weisheitsdatei '''+s+''' nicht schreiben!');
prot.destroyAll;
halt(1);
end;
while simulation.iteriereSchritt(start,zeitPhysik,zeitDatei) do ;
case errorCode of
2: prot.schreibe('Simulation wurde auf halbem Wege abgebrochen wegen Überlaufs.',true);
3: prot.schreibe('Simulation wurde auf halbem Wege abgebrochen wegen Benutzereingriffs.',true);
end{of case};
simulation.free;
prot.schreibe('fertig!',true);
s:=timetostr(now-start);
t:=timetostr(zeitDatei);
u:=timetostr(zeitPhysik);
while length(s)<max(length(t),length(u)) do s:=' '+s;
while length(t)<max(length(s),length(u)) do t:=' '+t;
while length(u)<max(length(s),length(t)) do u:=' '+u;
prot.schreibe('Das hat '+s+' gedauert,',true);
prot.schreibe(' davon '+t+' für Dateizugriffe',true);
prot.schreibe('und nur '+u+' für die eigentliche Physik!',true);
prot.free;
// stop program loop
Terminate;
if errorCode=1 then errorCode:=0;
end;
constructor TPlasmapropagation.Create(TheOwner: TComponent);
begin
inherited Create(TheOwner);
errorCode:=1;
StopOnException:=True;
end;
destructor TPlasmapropagation.Destroy;
begin
inherited Destroy;
end;
var
Application: TPlasmapropagation;
begin
Application:=TPlasmapropagation.Create(nil);
Application.Run;
Application.Free;
Halt(errorCode);
end.
|
