summaryrefslogtreecommitdiff
path: root/Plasmapropagation.lpr
blob: 478901f3069966075920a26a88ce09f34386e34f (plain)
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
program Plasmapropagation;

{$mode objfpc}{$H+}

uses
  {$IFDEF UNIX}{$IFDEF UseCThreads}
  cthreads,
  {$ENDIF}{$ENDIF}
  Classes, SysUtils, CustApp,
  { you can add units after this }
  math, Physikunit, protokollunit, lowlevelunit, fftw;

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.