-
Notifications
You must be signed in to change notification settings - Fork 1
/
RunSingle.py
82 lines (70 loc) · 2.45 KB
/
RunSingle.py
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
'''
RunSingle
Start a simulation of the agent-based model by launching cpp-code 'swarmdyn'
with the here defined parameters.
Also visualizes the simulation by calling the script "AnimateRun.py".
v0.1, 13.5.2020
(C) 2020 Pascal Klamser
This program is free software: you can redistribute it and/or modify
it under the terms of the GNU General Public License as published by
the Free Software Foundation, either version 3 of the License, or
(at your option) any later version.
This program is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
GNU General Public License for more details.
You should have received a copy of the GNU General Public License
along with this program. If not, see <http://www.gnu.org/licenses/>.
'''
import os
import time as pytime
import matplotlib
if __name__ == '__main__':
matplotlib.use('TkAgg')
import SwarmDynByPy as swarmPy
import AnimateRun
import json
from pathlib import Path
def main():
# Input Parameters
#########################################
trans_time = 20
record_time = 120
outpath = ''
dic = dict()
dic = swarmPy.get_base_params(record_time, trans_time=trans_time)
# changes from base-parameters
dicC = dict()
dicC['N'] = 400
dicC['beta'] = 2 # relaxation coefficient
f_name = ''
keys = sorted(list(dicC.keys()))
for k in keys:
f_name += '{}{}'.format(k, dicC[k])
dic[k] = dicC[k]
print(f_name)
# Generate and Run Command
#########################################
swarmPy.solve_parameter_dependencies(dic)
f_para = Path.cwd() / 'parameters.json'
json.dump(dic, f_para.open('w', encoding='utf-8'), indent=4, ensure_ascii=False)
command = swarmPy.dic2swarmdyn_command(dic)
print(command)
t0 = pytime.time()
os.system('make cl;')
os.system(command)
t1 = pytime.time()
print(t1-t0)
size = None
if 'BC' in dicC.keys():
size = dicC['size']
AnimateRun.main(None, size=size)
# reproducible movies:
if False:
outpath = Path('/home/klamser/Seafile/PoSI_EmergentStructure/Data/Videos_simu/xso_BoundaryN400')
outpath = Path('/home/klamser/Seafile/PoSI_EmergentStructure/Data/SimuPPK/')
outpath /= f_name
swarmPy.reproducible(outpath)
return
if __name__ == '__main__':
main()