# - - parameters - - - - - - - - - - - 
# Physical setting 
system["L"] = (4, 4, 4, 4)
system["β"] = 1.9
system["NC"] = 2
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_su2.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_su2"
	
# 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:42.601
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 0.8997670230155203 # itrj irand chiralcond
# 0 2 0.888283077955815 # itrj irand chiralcond
# 0 3 0.8675279843280307 # itrj irand chiralcond
# 0 4 0.8679008104614023 # itrj irand chiralcond
# 0 5 0.8647324335715463 # itrj irand chiralcond
# 0 6 0.8854424656132709 # itrj irand chiralcond
# 0 7 0.8902772446764718 # itrj irand chiralcond
# 0 8 0.8807265955754637 # itrj irand chiralcond
# 0 9 0.8806530145717533 # itrj irand chiralcond
# 0 10 0.8699537884622627 # itrj irand chiralcond
0 0.8795264438231538 # pbp Nr=10
-----------------
-----------------
0 2.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/8 is done
1 2
Hadron spectrum: Inversion 2/8 is done
1 3
Hadron spectrum: Inversion 3/8 is done
1 4
Hadron spectrum: Inversion 4/8 is done
2 1
Hadron spectrum: Inversion 5/8 is done
2 2
Hadron spectrum: Inversion 6/8 is done
2 3
Hadron spectrum: Inversion 7/8 is done
2 4
Hadron spectrum: Inversion 8/8 is done
Hadron spectrum: Reconstruction
Hadron spectrum end
0 52.76988794828185 1.7967209471177712 1.766446204031999e-25 1.7967209471177712 #pioncorrelator
-----------------
-----------------
0 1.0 # plaq
-----------------
Acceptance 1/1 : 100.0 %
-----------------
1 1.021073029385355 0.0 # poly
-----------------
-----------------
1 0.6941894930053661 # plaq
-----------------
-------------------------------------
Acceptance 2/2 : 100.0 %
-----------------
2 0.5303186307234687 0.0 # poly
-----------------
-----------------
2 0.5968992770023934 # plaq
-----------------
-------------------------------------
Acceptance 3/3 : 100.0 %
-----------------
3 0.7323656700197926 0.0 # poly
-----------------
-----------------
3 0.5680765565783865 # plaq
-----------------
-------------------------------------
Acceptance 4/4 : 100.0 %
-----------------
4 0.6669580613484427 0.0 # poly
-----------------
-----------------
4 0.5734360882646219 # plaq
-----------------
-------------------------------------
Acceptance 5/5 : 100.0 %
-----------------
Chiral condensate for Nf = 4, Dirac_operator = Staggered, factor=1.0 is multiplied.
# 5 1 0.9485657558971221 # itrj irand chiralcond
# 5 2 0.9034796642285211 # itrj irand chiralcond
# 5 3 0.8963442679586214 # itrj irand chiralcond
# 5 4 0.9506246557014337 # itrj irand chiralcond
# 5 5 0.9106702987578454 # itrj irand chiralcond
# 5 6 0.9519605663368266 # itrj irand chiralcond
# 5 7 0.9339468898888418 # itrj irand chiralcond
# 5 8 0.9437917333797714 # itrj irand chiralcond
# 5 9 0.9421981884826576 # itrj irand chiralcond
# 5 10 1.002230786435428 # itrj irand chiralcond
5 0.938381280706707 # pbp Nr=10
-----------------
-----------------
5 0.3021143350816606 0.0 # poly
-----------------
-----------------
Hadron spectrum started
1 1
Hadron spectrum: Inversion 1/8 is done
1 2
Hadron spectrum: Inversion 2/8 is done
1 3
Hadron spectrum: Inversion 3/8 is done
1 4
Hadron spectrum: Inversion 4/8 is done
2 1
Hadron spectrum: Inversion 5/8 is done
2 2
Hadron spectrum: Inversion 6/8 is done
2 3
Hadron spectrum: Inversion 7/8 is done
2 4
Hadron spectrum: Inversion 8/8 is done
Hadron spectrum: Reconstruction
Hadron spectrum end
5 67.02498906054976 1.4468431351769575 0.4994023656748559 1.2654306029346762 #pioncorrelator
-----------------
-----------------
5 0.5430145443592769 # plaq
-----------------
-------------------------------------
# - - parameters - - - - - - - - - - - 
# Physical setting 
system["L"] = (4, 4, 4, 4)
system["β"] = 1.9
system["NC"] = 2
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_su2.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_su2"
	
# 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:34.022
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 0.8997670230155203 # itrj irand chiralcond
# 0 2 0.888283077955815 # itrj irand chiralcond
# 0 3 0.8675279843280307 # itrj irand chiralcond
# 0 4 0.8679008104614023 # itrj irand chiralcond
# 0 5 0.8647324335715463 # itrj irand chiralcond
# 0 6 0.8854424656132709 # itrj irand chiralcond
# 0 7 0.8902772446764718 # itrj irand chiralcond
# 0 8 0.8807265955754637 # itrj irand chiralcond
# 0 9 0.8806530145717533 # itrj irand chiralcond
# 0 10 0.8699537884622627 # itrj irand chiralcond
0 0.8795264438231538 # pbp Nr=10
-----------------
-----------------
0 2.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/8 is done
1 2
Hadron spectrum: Inversion 2/8 is done
1 3
Hadron spectrum: Inversion 3/8 is done
1 4
Hadron spectrum: Inversion 4/8 is done
2 1
Hadron spectrum: Inversion 5/8 is done
2 2
Hadron spectrum: Inversion 6/8 is done
2 3
Hadron spectrum: Inversion 7/8 is done
2 4
Hadron spectrum: Inversion 8/8 is done
Hadron spectrum: Reconstruction
Hadron spectrum end
0 52.76988794828185 1.7967209471177712 1.766446204031999e-25 1.7967209471177712 #pioncorrelator
-----------------
-----------------
0 1.0 # plaq
-----------------
Acceptance 1/1 : 100.0 %
-----------------
1 0.940402830164371 0.0 # poly
-----------------
-----------------
1 0.6727258187413361 # plaq
-----------------
-------------------------------------
Acceptance 2/2 : 100.0 %
-----------------
2 0.6514590146470397 0.0 # poly
-----------------
-----------------
2 0.5723482574960691 # plaq
-----------------
-------------------------------------
Acceptance 3/3 : 100.0 %
-----------------
3 0.3823777915954255 0.0 # poly
-----------------
-----------------
3 0.5317955892555435 # plaq
-----------------
-------------------------------------
Acceptance 4/4 : 100.0 %
-----------------
4 0.4215896063043498 0.0 # poly
-----------------
-----------------
4 0.543134859891627 # plaq
-----------------
-------------------------------------
Acceptance 5/5 : 100.0 %
-----------------
Chiral condensate for Nf = 4, Dirac_operator = Staggered, factor=1.0 is multiplied.
# 5 1 0.9887203451078492 # itrj irand chiralcond
# 5 2 0.9504835441719351 # itrj irand chiralcond
# 5 3 0.9744318235652712 # itrj irand chiralcond
# 5 4 0.9039171130056395 # itrj irand chiralcond
# 5 5 0.9714905211645191 # itrj irand chiralcond
# 5 6 0.883771997361961 # itrj irand chiralcond
# 5 7 0.9911297224031665 # itrj irand chiralcond
# 5 8 0.9794198438662888 # itrj irand chiralcond
# 5 9 0.9277094879481179 # itrj irand chiralcond
# 5 10 1.0411313374238589 # itrj irand chiralcond
5 0.9612205736018609 # pbp Nr=10
-----------------
-----------------
5 0.22816957969053114 0.0 # poly
-----------------
-----------------
Hadron spectrum started
1 1
Hadron spectrum: Inversion 1/8 is done
1 2
Hadron spectrum: Inversion 2/8 is done
1 3
Hadron spectrum: Inversion 3/8 is done
1 4
Hadron spectrum: Inversion 4/8 is done
2 1
Hadron spectrum: Inversion 5/8 is done
2 2
Hadron spectrum: Inversion 6/8 is done
2 3
Hadron spectrum: Inversion 7/8 is done
2 4
Hadron spectrum: Inversion 8/8 is done
Hadron spectrum: Reconstruction
Hadron spectrum end
5 64.77476376029514 1.353616174646661 0.6665467762825079 1.7088502548862055 #pioncorrelator
-----------------
-----------------
5 0.511296176942659 # plaq
-----------------
-------------------------------------
