-
Notifications
You must be signed in to change notification settings - Fork 1
/
LQCDAnalyser.py
115 lines (96 loc) · 4.49 KB
/
LQCDAnalyser.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
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
113
114
115
#!/usr/bin/env python2
from pre_analysis.pre_analyser import pre_analysis
from post_analysis.post_analyser import post_analysis
from tools.folderreadingtools import get_num_observables
import copy
import os
from analysis_batches.main_analysis import main_analysis
from analysis_batches.distribution_analysis import distribution_analysis
from analysis_batches.topc_modes_analysis import topc_modes_analysis
from analysis_batches.thermalization_comparison_analysis \
import thermalization_analysis
from analysis_batches.b645_L32x32_analysis import beta645_L32_analysis
from analysis_batches.scaling_analysis import scaling_analysis
from analysis_batches.ABC_analysis import ABC_analysis
from analysis_batches.BCD_analysis import BCD_analysis
def main():
# Printing settings
section_seperator = "="*160
# Default params
run_pre_analysis = True
run_post_analysis = True
only_generate_data = False
observables = None
post_analysis_data_type = ["bootstrap", "unanalyzed", "blocked_bootstrap",
"jackknife", "blocked", "bootstrap_time_series"]
# Overriding params for what to run
run_pre_analysis = False
# run_post_analysis = False
only_generate_data = True
# Observables selection, full
observables = ["plaq", "topc", "topc2", "topc4", "topcr", "topsus",
"topsusqtq0", "qtq0e", "qtq0eff", "topcMC"]
# The essential observables we are inspecting
observables = ["topc", "topcr", "topsus", "topsusqtq0", "qtq0e", "qtq0eff"]
# observables = ["topsusqtq0"]
observables = ["qtq0e"]
# observables = ["topcr", "qtq0e", "qtq0eff", "energy"]
# observables = ["topcMC"]
# observables = ["energy"]
# observables = ["qtq0eff"]
# observables = ["topc"]
# observables = ["topcMC"]
# observables = ["topsus"]
# observables = ["topcr"]
observables = ["topsus", "topsusqtq0"]
# observables = ["topsus"]
# observables = ["topsusqtq0"]
# Sets the post analysis type to use
post_analysis_data_type = ["bootstrap_time_series", "bootstrap"]
# post_analysis_data_type = ["bootstrap_time_series"]
# post_analysis_data_type = ["bootstrap"]
# Full analysis of *all* elements available
beta645_L32_analysis(run_pre_analysis=run_pre_analysis,
run_post_analysis=run_post_analysis,
only_generate_data=only_generate_data,
observables=observables,
post_analysis_data_type=post_analysis_data_type,
include_b645x48xx3x96=True)
main_analysis(run_pre_analysis=False,
run_post_analysis=run_post_analysis,
only_generate_data=only_generate_data,
observables=observables,
post_analysis_data_type=post_analysis_data_type)
beta645_L32_analysis(run_pre_analysis=False,
run_post_analysis=run_post_analysis,
only_generate_data=only_generate_data,
observables=observables,
post_analysis_data_type=post_analysis_data_type,
include_b645x48xx3x96=False)
# # Only run these in presence of topological susceptibility
# if ("topsus" in observables) or ("energy" in observables):
# observables = list(set(observables) & set(["topsus", "energy"]))
# ABC_analysis(run_pre_analysis=False,
# run_post_analysis=run_post_analysis,
# only_generate_data=only_generate_data,
# observables=observables,
# post_analysis_data_type=post_analysis_data_type)
# BCD_analysis(run_pre_analysis=False,
# run_post_analysis=run_post_analysis,
# only_generate_data=only_generate_data,
# observables=observables,
# post_analysis_data_type=post_analysis_data_type,
# include_b645x48xx3x96=False)
# BCD_analysis(run_pre_analysis=False,
# run_post_analysis=run_post_analysis,
# only_generate_data=only_generate_data,
# observables=observables,
# post_analysis_data_type=post_analysis_data_type,
# include_b645x48xx3x96=True)
# distribution_analysis()
# topc_modes_analysis()
# lattice_updates_analysis()
# thermalization_analysis()
# scaling_analysis()
if __name__ == '__main__':
main()