# - - parameters - - - - - - - - - - - 
# Physical setting 
system["L"] = (4, 4, 4, 4)
system["β"] = 9.0
system["NC"] = 4
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_su4.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_su4"
	
# 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:58:08.177
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.7452094905710938 # itrj irand chiralcond
# 0 2 1.744817153119431 # itrj irand chiralcond
# 0 3 1.7132852134634409 # itrj irand chiralcond
# 0 4 1.7628154356130994 # itrj irand chiralcond
# 0 5 1.74883561832629 # itrj irand chiralcond
# 0 6 1.7548897726389254 # itrj irand chiralcond
# 0 7 1.733162945797312 # itrj irand chiralcond
# 0 8 1.7587966170550384 # itrj irand chiralcond
# 0 9 1.8079376241328518 # itrj irand chiralcond
# 0 10 1.8094920687428437 # itrj irand chiralcond
0 1.7579241939460328 # pbp Nr=10
-----------------
-----------------
0 4.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 1.0 1.0 -0.0 -0.0 0.0 #flow itrj flowtime plaq E Qplaq Qclov Qimproved
0 0.08 1.0 1.0 -0.0 -0.0 0.0 #flow itrj flowtime plaq E Qplaq Qclov Qimproved
0 0.12 1.0 1.0 -0.0 -0.0 0.0 #flow itrj flowtime plaq E Qplaq Qclov Qimproved
0 0.16 1.0 1.0 -0.0 -0.0 0.0 #flow itrj flowtime plaq E Qplaq Qclov Qimproved
0 0.2 1.0 1.0 -0.0 -0.0 0.0 #flow itrj flowtime plaq E Qplaq Qclov Qimproved
0 0.24 1.0 1.0 -0.0 -0.0 0.0 #flow itrj flowtime plaq E Qplaq Qclov Qimproved
0 0.28 1.0 1.0 -0.0 -0.0 0.0 #flow itrj flowtime plaq E Qplaq Qclov Qimproved
0 0.32 1.0 1.0 -0.0 -0.0 0.0 #flow itrj flowtime plaq E Qplaq Qclov Qimproved
0 0.36 1.0 1.0 -0.0 -0.0 0.0 #flow itrj flowtime plaq E Qplaq Qclov Qimproved
0 0.4 1.0 1.0 -0.0 -0.0 0.0 #flow itrj flowtime plaq E Qplaq Qclov Qimproved
-----------------
-----------------
Hadron spectrum started
1 1
Hadron spectrum: Inversion 1/16 is done
1 2
Hadron spectrum: Inversion 2/16 is done
1 3
Hadron spectrum: Inversion 3/16 is done
1 4
Hadron spectrum: Inversion 4/16 is done
2 1
Hadron spectrum: Inversion 5/16 is done
2 2
Hadron spectrum: Inversion 6/16 is done
2 3
Hadron spectrum: Inversion 7/16 is done
2 4
Hadron spectrum: Inversion 8/16 is done
3 1
Hadron spectrum: Inversion 9/16 is done
3 2
Hadron spectrum: Inversion 10/16 is done
3 3
Hadron spectrum: Inversion 11/16 is done
3 4
Hadron spectrum: Inversion 12/16 is done
4 1
Hadron spectrum: Inversion 13/16 is done
4 2
Hadron spectrum: Inversion 14/16 is done
4 3
Hadron spectrum: Inversion 15/16 is done
4 4
Hadron spectrum: Inversion 16/16 is done
Hadron spectrum: Reconstruction
Hadron spectrum end
0 211.079551793159 7.186883788470843 7.065784816127768e-25 7.186883788470843 #pioncorrelator
-----------------
-----------------
0 1.0 # plaq
-----------------
Acceptance 1/1 : 100.0 %
-----------------
1 1.9698434498220048 0.005788167429629749 # poly
-----------------
-----------------
1 0.8146770820561052 # plaq
-----------------
-------------------------------------
Acceptance 2/2 : 100.0 %
-----------------
2 1.5953507028356364 -0.031071735653988033 # poly
-----------------
-----------------
2 0.6952060303710503 # plaq
-----------------
-------------------------------------
Acceptance 3/3 : 100.0 %
-----------------
3 1.242885437887568 0.04130781481091159 # poly
-----------------
-----------------
3 0.60942385332367 # plaq
-----------------
-------------------------------------
Acceptance 4/4 : 100.0 %
-----------------
4 0.784464418883137 -0.16539557396038684 # poly
-----------------
-----------------
4 0.5543264357400718 # plaq
-----------------
-------------------------------------
Acceptance 5/5 : 100.0 %
-----------------
Chiral condensate for Nf = 4, Dirac_operator = Staggered, factor=1.0 is multiplied.
# 5 1 1.8498738638599748 # itrj irand chiralcond
# 5 2 1.8319947153407354 # itrj irand chiralcond
# 5 3 1.7924796441385225 # itrj irand chiralcond
# 5 4 1.8888935761028773 # itrj irand chiralcond
# 5 5 1.784558795396484 # itrj irand chiralcond
# 5 6 1.8161734804173306 # itrj irand chiralcond
# 5 7 2.000073559434345 # itrj irand chiralcond
# 5 8 1.7151776391092306 # itrj irand chiralcond
# 5 9 1.7563591244296022 # itrj irand chiralcond
# 5 10 1.8362984590249218 # itrj irand chiralcond
5 1.827188285725402 # pbp Nr=10
-----------------
-----------------
5 0.7742708986781698 -0.09824472386053917 # poly
-----------------
-----------------
Hadron spectrum started
1 1
Hadron spectrum: Inversion 1/16 is done
1 2
Hadron spectrum: Inversion 2/16 is done
1 3
Hadron spectrum: Inversion 3/16 is done
1 4
Hadron spectrum: Inversion 4/16 is done
2 1
Hadron spectrum: Inversion 5/16 is done
2 2
Hadron spectrum: Inversion 6/16 is done
2 3
Hadron spectrum: Inversion 7/16 is done
2 4
Hadron spectrum: Inversion 8/16 is done
3 1
Hadron spectrum: Inversion 9/16 is done
3 2
Hadron spectrum: Inversion 10/16 is done
3 3
Hadron spectrum: Inversion 11/16 is done
3 4
Hadron spectrum: Inversion 12/16 is done
4 1
Hadron spectrum: Inversion 13/16 is done
4 2
Hadron spectrum: Inversion 14/16 is done
4 3
Hadron spectrum: Inversion 15/16 is done
4 4
Hadron spectrum: Inversion 16/16 is done
Hadron spectrum: Reconstruction
Hadron spectrum end
5 250.9688754478505 1.7116768952953496 0.5752874912465761 2.2481126575747408 #pioncorrelator
-----------------
-----------------
5 0.4987986141745169 # plaq
-----------------
-------------------------------------
# - - parameters - - - - - - - - - - - 
# Physical setting 
system["L"] = (4, 4, 4, 4)
system["β"] = 9.0
system["NC"] = 4
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_su4.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_su4"
	
# 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:57.066
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.7452094905710938 # itrj irand chiralcond
# 0 2 1.744817153119431 # itrj irand chiralcond
# 0 3 1.7132852134634409 # itrj irand chiralcond
# 0 4 1.7628154356130994 # itrj irand chiralcond
# 0 5 1.74883561832629 # itrj irand chiralcond
# 0 6 1.7548897726389254 # itrj irand chiralcond
# 0 7 1.733162945797312 # itrj irand chiralcond
# 0 8 1.7587966170550384 # itrj irand chiralcond
# 0 9 1.8079376241328518 # itrj irand chiralcond
# 0 10 1.8094920687428437 # itrj irand chiralcond
0 1.7579241939460328 # pbp Nr=10
-----------------
-----------------
0 4.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 1.0 1.0 -0.0 -0.0 0.0 #flow itrj flowtime plaq E Qplaq Qclov Qimproved
0 0.08 1.0 1.0 -0.0 -0.0 0.0 #flow itrj flowtime plaq E Qplaq Qclov Qimproved
0 0.12 1.0 1.0 -0.0 -0.0 0.0 #flow itrj flowtime plaq E Qplaq Qclov Qimproved
0 0.16 1.0 1.0 -0.0 -0.0 0.0 #flow itrj flowtime plaq E Qplaq Qclov Qimproved
0 0.2 1.0 1.0 -0.0 -0.0 0.0 #flow itrj flowtime plaq E Qplaq Qclov Qimproved
0 0.24 1.0 1.0 -0.0 -0.0 0.0 #flow itrj flowtime plaq E Qplaq Qclov Qimproved
0 0.28 1.0 1.0 -0.0 -0.0 0.0 #flow itrj flowtime plaq E Qplaq Qclov Qimproved
0 0.32 1.0 1.0 -0.0 -0.0 0.0 #flow itrj flowtime plaq E Qplaq Qclov Qimproved
0 0.36 1.0 1.0 -0.0 -0.0 0.0 #flow itrj flowtime plaq E Qplaq Qclov Qimproved
0 0.4 1.0 1.0 -0.0 -0.0 0.0 #flow itrj flowtime plaq E Qplaq Qclov Qimproved
-----------------
-----------------
Hadron spectrum started
1 1
Hadron spectrum: Inversion 1/16 is done
1 2
Hadron spectrum: Inversion 2/16 is done
1 3
Hadron spectrum: Inversion 3/16 is done
1 4
Hadron spectrum: Inversion 4/16 is done
2 1
Hadron spectrum: Inversion 5/16 is done
2 2
Hadron spectrum: Inversion 6/16 is done
2 3
Hadron spectrum: Inversion 7/16 is done
2 4
Hadron spectrum: Inversion 8/16 is done
3 1
Hadron spectrum: Inversion 9/16 is done
3 2
Hadron spectrum: Inversion 10/16 is done
3 3
Hadron spectrum: Inversion 11/16 is done
3 4
Hadron spectrum: Inversion 12/16 is done
4 1
Hadron spectrum: Inversion 13/16 is done
4 2
Hadron spectrum: Inversion 14/16 is done
4 3
Hadron spectrum: Inversion 15/16 is done
4 4
Hadron spectrum: Inversion 16/16 is done
Hadron spectrum: Reconstruction
Hadron spectrum end
0 211.079551793159 7.186883788470843 7.065784816127768e-25 7.186883788470843 #pioncorrelator
-----------------
-----------------
0 1.0 # plaq
-----------------
Acceptance 1/1 : 100.0 %
-----------------
1 1.9125686807318039 0.055277432999890484 # poly
-----------------
-----------------
1 0.8070721354546924 # plaq
-----------------
-------------------------------------
Acceptance 2/2 : 100.0 %
-----------------
2 1.3480973198729955 0.1276238889951535 # poly
-----------------
-----------------
2 0.6880533658299977 # plaq
-----------------
-------------------------------------
Acceptance 3/3 : 100.0 %
-----------------
3 0.9177916772923059 -0.024637446631856557 # poly
-----------------
-----------------
3 0.6052695266631649 # plaq
-----------------
-------------------------------------
Acceptance 4/4 : 100.0 %
-----------------
4 0.8240368261297066 -0.1132414579201632 # poly
-----------------
-----------------
4 0.5464523971449331 # plaq
-----------------
-------------------------------------
Acceptance 5/5 : 100.0 %
-----------------
Chiral condensate for Nf = 4, Dirac_operator = Staggered, factor=1.0 is multiplied.
# 5 1 1.8127839461066846 # itrj irand chiralcond
# 5 2 1.8477245363474801 # itrj irand chiralcond
# 5 3 1.8636336246089307 # itrj irand chiralcond
# 5 4 1.844268949001634 # itrj irand chiralcond
# 5 5 1.8480520382743 # itrj irand chiralcond
# 5 6 1.8116093121553194 # itrj irand chiralcond
# 5 7 1.8162187426380214 # itrj irand chiralcond
# 5 8 1.8265317862670016 # itrj irand chiralcond
# 5 9 1.8502615703730956 # itrj irand chiralcond
# 5 10 1.9164061289604213 # itrj irand chiralcond
5 1.843749063473289 # pbp Nr=10
-----------------
-----------------
5 0.660149485244372 0.12329463752684058 # poly
-----------------
-----------------
Hadron spectrum started
1 1
Hadron spectrum: Inversion 1/16 is done
1 2
Hadron spectrum: Inversion 2/16 is done
1 3
Hadron spectrum: Inversion 3/16 is done
1 4
Hadron spectrum: Inversion 4/16 is done
2 1
Hadron spectrum: Inversion 5/16 is done
2 2
Hadron spectrum: Inversion 6/16 is done
2 3
Hadron spectrum: Inversion 7/16 is done
2 4
Hadron spectrum: Inversion 8/16 is done
3 1
Hadron spectrum: Inversion 9/16 is done
3 2
Hadron spectrum: Inversion 10/16 is done
3 3
Hadron spectrum: Inversion 11/16 is done
3 4
Hadron spectrum: Inversion 12/16 is done
4 1
Hadron spectrum: Inversion 13/16 is done
4 2
Hadron spectrum: Inversion 14/16 is done
4 3
Hadron spectrum: Inversion 15/16 is done
4 4
Hadron spectrum: Inversion 16/16 is done
Hadron spectrum: Reconstruction
Hadron spectrum end
5 250.79624911126703 1.7852633865054777 0.5135247054474883 2.4344163113053914 #pioncorrelator
-----------------
-----------------
5 0.4968305279850038 # plaq
-----------------
-------------------------------------
