# - - parameters - - - - - - - - - - - 
# Physical setting 
system["L"] = (4, 4, 4, 4)
system["β"] = 5.7
system["NC"] = 3
system["Nthermalization"] = 0
system["Nsteps"] = 5
system["initial"] = "cold"
system["initialtrj"] = 1
system["update_method"] = "HMC"
system["Nwing"] = 1
	
# Physical setting(fermions)
system["quench"] = true
system["Dirac_operator"] = nothing
wilson["Clover_coefficient"] = 0
wilson["r"] = 1
wilson["hop"] = 0
staggered["Nf"] = 0
staggered["mass"] = 0
system["BoundaryCondition"] = [1, 1, 1, -1]
	
# System Control
system["log_dir"] = "./logs"
system["logfile"] = "HMC_L04040404_beta5.7_quenched_su3.txt"
system["saveU_dir"] = ""
system["saveU_format"] = nothing
system["saveU_every"] = 1
system["verboselevel"] = 2
system["randomseed"] = 111
measurement["measurement_basedir"] = "./measurements"
measurement["measurement_dir"] = "HMC_L04040404_beta5.7_quenched_su3"
	
# HMC related
md["Δτ"] = 0.06666666666666667
md["SextonWeingargten"] = false
md["N_SextonWeingargten"] = 2
md["MDsteps"] = 15
cg["eps"] = 1.0e-19
cg["MaxCGstep"] = 3000
	
# Action parameter for SLMC
actions["use_autogeneratedstaples"] = false
actions["couplingcoeff"] = Any[]
actions["couplinglist"] = Any[]
actions["coupling_loops"] = nothing
	
# Measurement set
measurement["measurement_methods"] = Dict[ 
  Dict{Any,Any}("methodname" => "Chiral_condensate",
    "eps" => 1.0e-19,
    "fermiontype" => "Staggered",
    "Nf" => 4,
    "mass" => 0.5,
    "measure_every" => 5,
    "MaxCGstep" => 3000
  ),
  Dict{Any,Any}("methodname" => "Polyakov_loop",
    "measure_every" => 1
  ),
  Dict{Any,Any}("methodname" => "Topological_charge",
    "Nflowsteps" => 4,
    "numflow" => 10,
    "measure_every" => 10,
    "eps_flow" => 0.01
  ),
  Dict{Any,Any}("methodname" => "Pion_correlator",
    "eps" => 1.0e-19,
    "fermiontype" => "Wilson",
    "r" => 1,
    "measure_every" => 5,
    "hop" => 0.141139,
    "MaxCGstep" => 3000
  ),
  Dict{Any,Any}("methodname" => "Plaquette",
    "measure_every" => 1
  )
]
	
# - - - - - - - - - - - - - - - - - - -
# /Users/yuki/gitback/LatticeQCD.jl/test
# 2021-10-26T09:57:54.950
Julia Version 1.6.2
Commit 1b93d53fc4 (2021-07-14 15:36 UTC)
Platform Info:
  OS: macOS (x86_64-apple-darwin18.7.0)
  CPU: Apple M1
  WORD_SIZE: 64
  LIBM: libopenlibm
  LLVM: libLLVM-11.0.1 (ORCJIT, westmere)
Environment:
  JULIA_SSL_NO_VERIFY_HOSTS = github.com
  JULIA_LOAD_PATH = @:/var/folders/kq/sxg320kn3hvdxkxky0c3dxwh0000gn/T/jl_x3hqgd
-----------------
Chiral condensate for Nf = 4, Dirac_operator = Staggered, factor=1.0 is multiplied.
# 0 1 1.3369615388151035 # itrj irand chiralcond
# 0 2 1.3161174483831688 # itrj irand chiralcond
# 0 3 1.3210328731178065 # itrj irand chiralcond
# 0 4 1.2845825229987118 # itrj irand chiralcond
# 0 5 1.33804701577134 # itrj irand chiralcond
# 0 6 1.340461155983293 # itrj irand chiralcond
# 0 7 1.312524152263332 # itrj irand chiralcond
# 0 8 1.2953855253473705 # itrj irand chiralcond
# 0 9 1.305101094752209 # itrj irand chiralcond
# 0 10 1.2996309880364953 # itrj irand chiralcond
0 1.314984431546883 # pbp Nr=10
-----------------
-----------------
0 3.0 0.0 # poly
-----------------
-----------------
# epsilon for the Wilson flow is 0.01
0 0.0 1.0 1.0 -0.0 -0.0 0.0 #flow itrj flowtime plaq E Qplaq Qclov Qimproved
0 0.04 NaN NaN NaN NaN NaN #flow itrj flowtime plaq E Qplaq Qclov Qimproved
0 0.08 NaN NaN NaN NaN NaN #flow itrj flowtime plaq E Qplaq Qclov Qimproved
0 0.12 NaN NaN NaN NaN NaN #flow itrj flowtime plaq E Qplaq Qclov Qimproved
0 0.16 NaN NaN NaN NaN NaN #flow itrj flowtime plaq E Qplaq Qclov Qimproved
0 0.2 NaN NaN NaN NaN NaN #flow itrj flowtime plaq E Qplaq Qclov Qimproved
0 0.24 NaN NaN NaN NaN NaN #flow itrj flowtime plaq E Qplaq Qclov Qimproved
0 0.28 NaN NaN NaN NaN NaN #flow itrj flowtime plaq E Qplaq Qclov Qimproved
0 0.32 NaN NaN NaN NaN NaN #flow itrj flowtime plaq E Qplaq Qclov Qimproved
0 0.36 NaN NaN NaN NaN NaN #flow itrj flowtime plaq E Qplaq Qclov Qimproved
0 0.4 NaN NaN NaN NaN NaN #flow itrj flowtime plaq E Qplaq Qclov Qimproved
-----------------
-----------------
Hadron spectrum started
1 1
Hadron spectrum: Inversion 1/12 is done
1 2
Hadron spectrum: Inversion 2/12 is done
1 3
Hadron spectrum: Inversion 3/12 is done
1 4
Hadron spectrum: Inversion 4/12 is done
2 1
Hadron spectrum: Inversion 5/12 is done
2 2
Hadron spectrum: Inversion 6/12 is done
2 3
Hadron spectrum: Inversion 7/12 is done
2 4
Hadron spectrum: Inversion 8/12 is done
3 1
Hadron spectrum: Inversion 9/12 is done
3 2
Hadron spectrum: Inversion 10/12 is done
3 3
Hadron spectrum: Inversion 11/12 is done
3 4
Hadron spectrum: Inversion 12/12 is done
Hadron spectrum: Reconstruction
Hadron spectrum end
0 118.73224788364023 4.042622131015124 3.974503959071931e-25 4.042622131015124 #pioncorrelator
-----------------
-----------------
0 1.0 # plaq
-----------------
Acceptance 1/1 : 100.0 %
-----------------
1 1.6844293149677714 -0.014131801865616809 # poly
-----------------
-----------------
1 0.8352308729144909 # plaq
-----------------
-------------------------------------
Acceptance 2/2 : 100.0 %
-----------------
2 1.1058082016305715 0.003884846423294677 # poly
-----------------
-----------------
2 0.7364714493418552 # plaq
-----------------
-------------------------------------
Acceptance 3/3 : 100.0 %
-----------------
3 1.0861809827852005 0.061650613266930665 # poly
-----------------
-----------------
3 0.6620724930664146 # plaq
-----------------
-------------------------------------
Acceptance 4/4 : 100.0 %
-----------------
4 0.9659827553766934 0.0011714698229286285 # poly
-----------------
-----------------
4 0.6394146161193309 # plaq
-----------------
-------------------------------------
Acceptance 5/5 : 100.0 %
-----------------
Chiral condensate for Nf = 4, Dirac_operator = Staggered, factor=1.0 is multiplied.
# 5 1 1.2894975040422194 # itrj irand chiralcond
# 5 2 1.2521277422730637 # itrj irand chiralcond
# 5 3 1.2946588365023775 # itrj irand chiralcond
# 5 4 1.325515451522386 # itrj irand chiralcond
# 5 5 1.3171820065631015 # itrj irand chiralcond
# 5 6 1.2852731474608794 # itrj irand chiralcond
# 5 7 1.3331232944930973 # itrj irand chiralcond
# 5 8 1.3330965223362194 # itrj irand chiralcond
# 5 9 1.2596141495142785 # itrj irand chiralcond
# 5 10 1.313567471228398 # itrj irand chiralcond
5 1.300365612593602 # pbp Nr=10
-----------------
-----------------
5 0.8328693656663123 0.0885139312603251 # poly
-----------------
-----------------
Hadron spectrum started
1 1
Hadron spectrum: Inversion 1/12 is done
1 2
Hadron spectrum: Inversion 2/12 is done
1 3
Hadron spectrum: Inversion 3/12 is done
1 4
Hadron spectrum: Inversion 4/12 is done
2 1
Hadron spectrum: Inversion 5/12 is done
2 2
Hadron spectrum: Inversion 6/12 is done
2 3
Hadron spectrum: Inversion 7/12 is done
2 4
Hadron spectrum: Inversion 8/12 is done
3 1
Hadron spectrum: Inversion 9/12 is done
3 2
Hadron spectrum: Inversion 10/12 is done
3 3
Hadron spectrum: Inversion 11/12 is done
3 4
Hadron spectrum: Inversion 12/12 is done
Hadron spectrum: Reconstruction
Hadron spectrum end
5 137.25290161749606 0.7938690280905328 0.28022062911910994 1.341196091976362 #pioncorrelator
-----------------
-----------------
5 0.6173953730210353 # plaq
-----------------
-------------------------------------
# - - parameters - - - - - - - - - - - 
# Physical setting 
system["L"] = (4, 4, 4, 4)
system["β"] = 5.7
system["NC"] = 3
system["Nthermalization"] = 0
system["Nsteps"] = 5
system["initial"] = "cold"
system["initialtrj"] = 1
system["update_method"] = "HMC"
system["Nwing"] = 1
	
# Physical setting(fermions)
system["quench"] = true
system["Dirac_operator"] = nothing
wilson["Clover_coefficient"] = 0
wilson["r"] = 1
wilson["hop"] = 0
staggered["Nf"] = 0
staggered["mass"] = 0
system["BoundaryCondition"] = [1, 1, 1, -1]
	
# System Control
system["log_dir"] = "./logs"
system["logfile"] = "HMC_L04040404_beta5.7_quenched_su3.txt"
system["saveU_dir"] = ""
system["saveU_format"] = nothing
system["saveU_every"] = 1
system["verboselevel"] = 2
system["randomseed"] = 111
measurement["measurement_basedir"] = "./measurements"
measurement["measurement_dir"] = "HMC_L04040404_beta5.7_quenched_su3"
	
# HMC related
md["Δτ"] = 0.06666666666666667
md["SextonWeingargten"] = false
md["N_SextonWeingargten"] = 2
md["MDsteps"] = 15
cg["eps"] = 1.0e-19
cg["MaxCGstep"] = 3000
	
# Action parameter for SLMC
actions["use_autogeneratedstaples"] = false
actions["couplingcoeff"] = Any[]
actions["couplinglist"] = Any[]
actions["coupling_loops"] = nothing
	
# Measurement set
measurement["measurement_methods"] = Dict[ 
  Dict{Any,Any}("methodname" => "Chiral_condensate",
    "eps" => 1.0e-19,
    "fermiontype" => "Staggered",
    "Nf" => 4,
    "mass" => 0.5,
    "measure_every" => 5,
    "MaxCGstep" => 3000
  ),
  Dict{Any,Any}("methodname" => "Polyakov_loop",
    "measure_every" => 1
  ),
  Dict{Any,Any}("methodname" => "Topological_charge",
    "Nflowsteps" => 4,
    "numflow" => 10,
    "measure_every" => 10,
    "eps_flow" => 0.01
  ),
  Dict{Any,Any}("methodname" => "Pion_correlator",
    "eps" => 1.0e-19,
    "fermiontype" => "Wilson",
    "r" => 1,
    "measure_every" => 5,
    "hop" => 0.141139,
    "MaxCGstep" => 3000
  ),
  Dict{Any,Any}("methodname" => "Plaquette",
    "measure_every" => 1
  )
]
	
# - - - - - - - - - - - - - - - - - - -
# /Users/yuki/gitback/LatticeQCD.jl/test
# 2021-10-26T10:44:45.777
Julia Version 1.6.2
Commit 1b93d53fc4 (2021-07-14 15:36 UTC)
Platform Info:
  OS: macOS (x86_64-apple-darwin18.7.0)
  CPU: Apple M1
  WORD_SIZE: 64
  LIBM: libopenlibm
  LLVM: libLLVM-11.0.1 (ORCJIT, westmere)
Environment:
  JULIA_SSL_NO_VERIFY_HOSTS = github.com
  JULIA_LOAD_PATH = @:/var/folders/kq/sxg320kn3hvdxkxky0c3dxwh0000gn/T/jl_akKdPj
-----------------
Chiral condensate for Nf = 4, Dirac_operator = Staggered, factor=1.0 is multiplied.
# 0 1 1.3369615388151035 # itrj irand chiralcond
# 0 2 1.3161174483831688 # itrj irand chiralcond
# 0 3 1.3210328731178065 # itrj irand chiralcond
# 0 4 1.2845825229987118 # itrj irand chiralcond
# 0 5 1.33804701577134 # itrj irand chiralcond
# 0 6 1.340461155983293 # itrj irand chiralcond
# 0 7 1.312524152263332 # itrj irand chiralcond
# 0 8 1.2953855253473705 # itrj irand chiralcond
# 0 9 1.305101094752209 # itrj irand chiralcond
# 0 10 1.2996309880364953 # itrj irand chiralcond
0 1.314984431546883 # pbp Nr=10
-----------------
-----------------
0 3.0 0.0 # poly
-----------------
-----------------
# epsilon for the Wilson flow is 0.01
0 0.0 1.0 1.0 -0.0 -0.0 0.0 #flow itrj flowtime plaq E Qplaq Qclov Qimproved
0 0.04 NaN NaN NaN NaN NaN #flow itrj flowtime plaq E Qplaq Qclov Qimproved
0 0.08 NaN NaN NaN NaN NaN #flow itrj flowtime plaq E Qplaq Qclov Qimproved
0 0.12 NaN NaN NaN NaN NaN #flow itrj flowtime plaq E Qplaq Qclov Qimproved
0 0.16 NaN NaN NaN NaN NaN #flow itrj flowtime plaq E Qplaq Qclov Qimproved
0 0.2 NaN NaN NaN NaN NaN #flow itrj flowtime plaq E Qplaq Qclov Qimproved
0 0.24 NaN NaN NaN NaN NaN #flow itrj flowtime plaq E Qplaq Qclov Qimproved
0 0.28 NaN NaN NaN NaN NaN #flow itrj flowtime plaq E Qplaq Qclov Qimproved
0 0.32 NaN NaN NaN NaN NaN #flow itrj flowtime plaq E Qplaq Qclov Qimproved
0 0.36 NaN NaN NaN NaN NaN #flow itrj flowtime plaq E Qplaq Qclov Qimproved
0 0.4 NaN NaN NaN NaN NaN #flow itrj flowtime plaq E Qplaq Qclov Qimproved
-----------------
-----------------
Hadron spectrum started
1 1
Hadron spectrum: Inversion 1/12 is done
1 2
Hadron spectrum: Inversion 2/12 is done
1 3
Hadron spectrum: Inversion 3/12 is done
1 4
Hadron spectrum: Inversion 4/12 is done
2 1
Hadron spectrum: Inversion 5/12 is done
2 2
Hadron spectrum: Inversion 6/12 is done
2 3
Hadron spectrum: Inversion 7/12 is done
2 4
Hadron spectrum: Inversion 8/12 is done
3 1
Hadron spectrum: Inversion 9/12 is done
3 2
Hadron spectrum: Inversion 10/12 is done
3 3
Hadron spectrum: Inversion 11/12 is done
3 4
Hadron spectrum: Inversion 12/12 is done
Hadron spectrum: Reconstruction
Hadron spectrum end
0 118.73224788364023 4.042622131015124 3.974503959071931e-25 4.042622131015124 #pioncorrelator
-----------------
-----------------
0 1.0 # plaq
-----------------
Acceptance 1/1 : 100.0 %
-----------------
1 1.4488798388323665 -0.026943853252128072 # poly
-----------------
-----------------
1 0.8260478740475447 # plaq
-----------------
-------------------------------------
Acceptance 2/2 : 100.0 %
-----------------
2 1.174497471420054 -0.0479070151822465 # poly
-----------------
-----------------
2 0.7230547717069059 # plaq
-----------------
-------------------------------------
Acceptance 3/3 : 100.0 %
-----------------
3 0.8719477256759247 0.01335178016553522 # poly
-----------------
-----------------
3 0.6636818524769045 # plaq
-----------------
-------------------------------------
Acceptance 4/4 : 100.0 %
-----------------
4 0.7750264876644488 -0.1430700062252773 # poly
-----------------
-----------------
4 0.6449549071431334 # plaq
-----------------
-------------------------------------
Acceptance 5/5 : 100.0 %
-----------------
Chiral condensate for Nf = 4, Dirac_operator = Staggered, factor=1.0 is multiplied.
# 5 1 1.2990728707564398 # itrj irand chiralcond
# 5 2 1.2557880582089342 # itrj irand chiralcond
# 5 3 1.3809365114450536 # itrj irand chiralcond
# 5 4 1.3714152395010701 # itrj irand chiralcond
# 5 5 1.3308636454136042 # itrj irand chiralcond
# 5 6 1.348665299105648 # itrj irand chiralcond
# 5 7 1.3318453929582827 # itrj irand chiralcond
# 5 8 1.3214503489721934 # itrj irand chiralcond
# 5 9 1.287673492891771 # itrj irand chiralcond
# 5 10 1.3384003921839163 # itrj irand chiralcond
5 1.3266111251436914 # pbp Nr=10
-----------------
-----------------
5 0.7933562443217119 -0.0282676008655901 # poly
-----------------
-----------------
Hadron spectrum started
1 1
Hadron spectrum: Inversion 1/12 is done
1 2
Hadron spectrum: Inversion 2/12 is done
1 3
Hadron spectrum: Inversion 3/12 is done
1 4
Hadron spectrum: Inversion 4/12 is done
2 1
Hadron spectrum: Inversion 5/12 is done
2 2
Hadron spectrum: Inversion 6/12 is done
2 3
Hadron spectrum: Inversion 7/12 is done
2 4
Hadron spectrum: Inversion 8/12 is done
3 1
Hadron spectrum: Inversion 9/12 is done
3 2
Hadron spectrum: Inversion 10/12 is done
3 3
Hadron spectrum: Inversion 11/12 is done
3 4
Hadron spectrum: Inversion 12/12 is done
Hadron spectrum: Reconstruction
Hadron spectrum end
5 135.0351034213489 0.7904500504738583 0.42462041171896375 0.9817811700058765 #pioncorrelator
-----------------
-----------------
5 0.6198464934205753 # plaq
-----------------
-------------------------------------
