kees-study/cds-monte-carlo-methods
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

09: Do things after file things

Browse Source
This commit is contained in:
Kees van Kempen
2022-11-17 17:19:31 +01:00
parent 124d2d077e
commit 641bdf75a4
1 changed files with 13 additions and 12 deletions
Download Patch File Download Diff File Expand all files Collapse all files

25
Exercise sheet 9/latticescalar.py
View File

@ -91,18 +91,7 @@ def main():
else:
output_filename = args.f
# Measure
# TODO: Does mean_magn need to be a list?
mean_magn = []
phi_state = np.zeros((width,width,width,width))
run_scalar_MH(phi_state,lamb,kappa,delta,equil_sweeps * num_sites)
magnetizations = np.empty(measurements)
for i in range(measurements):
run_scalar_MH(phi_state,lamb,kappa,delta,measure_sweeps * num_sites)
magnetizations[i] = np.mean(phi_state)
mean, err = batch_estimate(np.abs(magnetizations),lambda x:np.mean(x),10)
mean_magn.append([mean,err])
last_output_time = time.time()
with h5py.File(output_filename,'a') as f:
if not "mean-magn" in f:
dataset = f.create_dataset("mean-magn", chunks=True, data=mean_magn)
@ -118,5 +107,17 @@ def main():
dataset.attrs["start time"] = starttime
dataset.attrs["stop time"] = time.asctime()
# Measure
# TODO: Does mean_magn need to be a list?
mean_magn = []
phi_state = np.zeros((width,width,width,width))
run_scalar_MH(phi_state,lamb,kappa,delta,equil_sweeps * num_sites)
magnetizations = np.empty(measurements)
for i in range(measurements):
run_scalar_MH(phi_state,lamb,kappa,delta,measure_sweeps * num_sites)
magnetizations[i] = np.mean(phi_state)
mean, err = batch_estimate(np.abs(magnetizations),lambda x:np.mean(x),10)
mean_magn.append([mean,err])
if __name__ == "__main__":
main()