Skip to content
Snippets Groups Projects
Commit f699c5a3 authored by Lukas Bentkamp's avatar Lukas Bentkamp Committed by Cristian Lalescu
Browse files

WIP: Gaussian field convergence test

parent f93dc4d7
No related branches found
No related tags found
No related merge requests found
......@@ -58,9 +58,10 @@ def main():
def main_convergence_test():
n_arr = [64, 128, 256, 512]
diss =
energy =
etaK =
dk_arr = [1., 0.5, 0.25, 0.125]
diss = 0.007036634455033392
energy = 0.0417604575121591
etaK = 0.2157937040813811
Lint = energy**(3./2.)/diss
c_L, c_eta = calculate_constants(epsilon=diss, Lint=Lint, etaK=etaK)
......@@ -68,7 +69,7 @@ def main_convergence_test():
bin_no = 100
if not os.path.exists(simname+'.h5'):
for n in n_arr:
for n, dk in zip(n_arr, dk_arr):
for run in range(3):
rseed = int(time.time())
simname = f'conv_test_n{n}_run{run}'
......@@ -79,6 +80,9 @@ def main_convergence_test():
'--nx', str(n),
'--ny', str(n),
'--nz', str(n),
'--dkx', str(dk),
'--dky', str(dk),
'--dkz', str(dk),
'--simname', simname,
'--np', '4',
'--environment', 'short',
......@@ -96,7 +100,7 @@ def main_convergence_test():
sys.argv[1:])
else:
plt.figure()
# TODO
for n in n_arr:
for run in range(3):
df = h5py.File(simname + '.h5', 'r')
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment