mirror of
https://igit.ific.uv.es/fernando.p.csic.es/sfcf.jl.git
synced 2025-11-30 01:46:52 +01:00
Minor changes.
This commit is contained in:
Fernando P. Panadero
parent
7d14ab241a
commit
5790c8f863
3 changed files with 87 additions and 80 deletions
|
|
@ -28,35 +28,38 @@ tolerance = 1.0e-13
|
||||||
maxiter = 1000
|
maxiter = 1000
|
||||||
|
|
||||||
Reading gauge field with LGPU native convention from: ./cnfg/randfield...
|
Reading gauge field with LGPU native convention from: ./cnfg/randfield...
|
||||||
|
MD5 checksum of gauge config: 89c5bd25d1cd2c74e59f0ecfcf95ebef
|
||||||
|
|
||||||
|
|
||||||
Gauge meassurements:
|
Gauge meassurements:
|
||||||
|
|
||||||
Plaquette: 0.3306405864818687
|
Plaquette: 0.3306405864818687
|
||||||
Topological charge: -0.057705150807682376
|
Topological charge: -0.057705150807682376
|
||||||
SF coupling (dS / deta): 6.346811365600632
|
SF coupling (dS / deta): 6.346811365600632
|
||||||
|
SF coupling (d^2S / deta dnu): 26.166706988703776
|
||||||
|
|
||||||
|
|
||||||
Computing propagators from t=0 to the bulk
|
Computing propagators from t=0 to the bulk
|
||||||
CG converged for c=1 and s=1 after 50 iterations with absolute residue 0.0009654241012493543
|
CG converged for c=1 and s=1 after 50 iterations with absolute residue 0.0009654241012493528
|
||||||
CG converged for c=2 and s=1 after 50 iterations with absolute residue 0.0009792505060407602
|
CG converged for c=2 and s=1 after 50 iterations with absolute residue 0.0009792505060407621
|
||||||
CG converged for c=3 and s=1 after 50 iterations with absolute residue 0.0009410537206344275
|
CG converged for c=3 and s=1 after 50 iterations with absolute residue 0.000941053720634425
|
||||||
CG converged for c=1 and s=2 after 50 iterations with absolute residue 0.0009706380660725444
|
CG converged for c=1 and s=2 after 50 iterations with absolute residue 0.000970638066072543
|
||||||
CG converged for c=2 and s=2 after 50 iterations with absolute residue 0.0009348723257035598
|
CG converged for c=2 and s=2 after 50 iterations with absolute residue 0.0009348723257035582
|
||||||
CG converged for c=3 and s=2 after 50 iterations with absolute residue 0.0009646157878217441
|
CG converged for c=3 and s=2 after 50 iterations with absolute residue 0.000964615787821745
|
||||||
Computing propagators from t=T to the bulk
|
Computing propagators from t=T to the bulk
|
||||||
CG converged for c=1 and s=1 after 50 iterations with absolute residue 0.0009409155699679223
|
CG converged for c=1 and s=1 after 50 iterations with absolute residue 0.00094091556996792
|
||||||
CG converged for c=2 and s=1 after 50 iterations with absolute residue 0.0009591325218333536
|
CG converged for c=2 and s=1 after 50 iterations with absolute residue 0.0009591325218333531
|
||||||
CG converged for c=3 and s=1 after 50 iterations with absolute residue 0.000993180324656404
|
CG converged for c=3 and s=1 after 50 iterations with absolute residue 0.000993180324656406
|
||||||
CG converged for c=1 and s=2 after 50 iterations with absolute residue 0.0009696740933237723
|
CG converged for c=1 and s=2 after 50 iterations with absolute residue 0.0009696740933237713
|
||||||
CG converged for c=2 and s=2 after 50 iterations with absolute residue 0.0009753894671580818
|
CG converged for c=2 and s=2 after 50 iterations with absolute residue 0.0009753894671580808
|
||||||
CG converged for c=3 and s=2 after 50 iterations with absolute residue 0.0009677391305951339
|
CG converged for c=3 and s=2 after 50 iterations with absolute residue 0.0009677391305951349
|
||||||
|
|
||||||
Computing correlators
|
Computing correlators
|
||||||
|
|
||||||
Computing f1...
|
Computing f1...
|
||||||
|
|
||||||
f1 =
|
f1 =
|
||||||
7.968153189691683e-10
|
7.968153189691686e-10 + 0.0im
|
||||||
|
|
||||||
Computing k1...
|
Computing k1...
|
||||||
|
|
||||||
|
|
@ -70,115 +73,117 @@ fP =
|
||||||
0.39164297606334497 + 0.0im
|
0.39164297606334497 + 0.0im
|
||||||
0.03838253084495811 + 0.0im
|
0.03838253084495811 + 0.0im
|
||||||
0.0036829603591832434 + 0.0im
|
0.0036829603591832434 + 0.0im
|
||||||
0.00034990928932360324 + 0.0im
|
0.0003499092893236032 + 0.0im
|
||||||
3.1848709670960246e-5 + 0.0im
|
3.1848709670960246e-5 + 0.0im
|
||||||
2.83087866678673e-6 + 0.0im
|
2.8308786667867295e-6 + 0.0im
|
||||||
2.552038846719931e-7 + 0.0im
|
2.552038846719931e-7 + 0.0im
|
||||||
|
|
||||||
Computing fA...
|
Computing fA...
|
||||||
|
|
||||||
fA =
|
fA =
|
||||||
0.0 + 0.0im
|
0.0 + 0.0im
|
||||||
-0.28022427481948287 + 3.9466441386130056e-20im
|
-0.28022427481948287 + 2.968585782330306e-20im
|
||||||
-0.028432536079057023 - 9.446402595355576e-22im
|
-0.028432536079057023 + 1.435985543392056e-21im
|
||||||
-0.002717462099554097 + 4.301339185275744e-23im
|
-0.002717462099554098 + 1.8611563782443123e-23im
|
||||||
-0.0002553997594898297 - 2.8822074468644558e-24im
|
-0.0002553997594898297 + 2.1067256225959924e-23im
|
||||||
-2.3190475107607993e-5 - 1.34416849539867e-24im
|
-2.3190475107607993e-5 + 3.723928343622517e-24im
|
||||||
-2.078551613006522e-6 - 3.0342745858421e-25im
|
-2.0785516130065226e-6 - 1.0905213153775267e-25im
|
||||||
-1.8801834874904572e-7 - 1.7443292336247013e-26im
|
-1.880183487490457e-7 + 3.6085653248429964e-26im
|
||||||
|
|
||||||
Computing gP...
|
Computing gP...
|
||||||
|
|
||||||
gP =
|
gP =
|
||||||
0.0 + 0.0im
|
0.0 + 0.0im
|
||||||
2.538597063821986e-7 + 0.0im
|
2.5385970638219865e-7 + 0.0im
|
||||||
2.8229409631152233e-6 + 0.0im
|
2.8229409631152233e-6 + 0.0im
|
||||||
3.096526353494133e-5 + 0.0im
|
3.0965263534941325e-5 + 0.0im
|
||||||
0.0003331843094846162 + 0.0im
|
0.0003331843094846162 + 0.0im
|
||||||
0.0036254847521787863 + 0.0im
|
0.0036254847521787863 + 0.0im
|
||||||
0.03837127818992878 + 0.0im
|
0.03837127818992878 + 0.0im
|
||||||
0.39410791485294716 + 0.0im
|
0.39410791485294727 + 0.0im
|
||||||
|
|
||||||
Computing gA...
|
Computing gA...
|
||||||
|
|
||||||
gA =
|
gA =
|
||||||
0.0 + 0.0im
|
0.0 + 0.0im
|
||||||
-1.8696507552685477e-7 - 6.096317075547979e-27im
|
-1.8696507552685472e-7 - 6.203602158658037e-27im
|
||||||
-2.0658936442177055e-6 - 3.842068152788416e-26im
|
-2.0658936442177055e-6 - 1.7405427012796404e-25im
|
||||||
-2.2674085372631826e-5 - 8.56261180963095e-25im
|
-2.2674085372631822e-5 + 2.1422685395416303e-24im
|
||||||
-0.00024433351353064797 + 9.060212646869881e-24im
|
-0.000244333513530648 + 5.8458404852771e-23im
|
||||||
-0.0026482638745495414 + 2.161009350294785e-22im
|
-0.0026482638745495414 + 3.7698756417104236e-22im
|
||||||
-0.028194731180271545 + 5.375019620369574e-21im
|
-0.028194731180271545 + 6.309733712554495e-21im
|
||||||
-0.28340818434504744 + 2.511320339710992e-20im
|
-0.28340818434504744 + 1.9051623868321333e-20im
|
||||||
|
|
||||||
Computing kV...
|
Computing kV...
|
||||||
|
|
||||||
kV =
|
kV =
|
||||||
0.0 + 0.0im
|
0.0 + 0.0im
|
||||||
0.33632172225753315 - 1.0340860564396276e-20im
|
0.33632172225753315 + 1.0940842235368113e-21im
|
||||||
0.03186802210515694 - 9.926167350636303e-23im
|
0.03186802210515694 + 1.0742318888355314e-21im
|
||||||
0.0029308316253718663 - 1.1084220208210557e-22im
|
0.0029308316253718663 - 9.760731228125728e-23im
|
||||||
0.0002663612299706971 + 3.171720661258017e-23im
|
0.00026636122997069706 + 1.672886624241357e-23im
|
||||||
2.3221979105382962e-5 - 1.3495537858449775e-24im
|
2.3221979105382955e-5 + 1.2580038482577253e-24im
|
||||||
1.9759435694009697e-6 + 3.989827059409007e-25im
|
1.9759435694009697e-6 - 2.610519543848177e-25im
|
||||||
1.7038427603748213e-7 + 2.5655860266867883e-26im
|
1.7038427603748213e-7 + 1.5226067285305813e-26im
|
||||||
|
|
||||||
Computing lV...
|
Computing lV...
|
||||||
|
|
||||||
lV =
|
lV =
|
||||||
0.0 + 0.0im
|
0.0 + 0.0im
|
||||||
1.692538762882869e-7 - 1.2470313189733742e-26im
|
1.692538762882869e-7 + 3.2328571710497395e-26im
|
||||||
1.972297868357158e-6 + 1.9703431420432253e-25im
|
1.9722978683571574e-6 + 2.1049754032009423e-25im
|
||||||
2.2655589268307515e-5 + 4.094974855373194e-24im
|
2.2655589268307515e-5 - 3.08684848382419e-24im
|
||||||
0.0002539424459369091 - 2.540133797714921e-23im
|
0.0002539424459369091 + 3.9704669402545317e-23im
|
||||||
0.0028788471359388916 - 4.1193594505140788e-22im
|
0.002878847135938892 + 5.79026428787119e-23im
|
||||||
0.031803280548268274 + 2.5322755818956662e-21im
|
0.03180328054826828 + 2.4749243927586577e-21im
|
||||||
0.3390357554359771 - 7.90564084104014e-21im
|
0.339035755435977 + 2.883441324611509e-20im
|
||||||
|
|
||||||
Computing kT...
|
Computing kT...
|
||||||
|
|
||||||
kT =
|
kT =
|
||||||
0.0 + 0.0im
|
0.0 + 0.0im
|
||||||
-0.30173882074287517 + 5.293955920339405e-22im
|
-0.30173882074287517 + 9.03501810404587e-21im
|
||||||
-0.028918463611070797 - 1.548482106699267e-21im
|
-0.0289184636110708 + 1.6234798155707435e-21im
|
||||||
-0.002655774263798736 + 6.81045371001993e-23im
|
-0.002655774263798736 + 2.451211881865472e-22im
|
||||||
-0.00024049035597286907 + 3.794691060087014e-23im
|
-0.00024049035597286904 + 5.2388105461691765e-24im
|
||||||
-2.0961481048454104e-5 - 1.4734154661100807e-24im
|
-2.0961481048454104e-5 - 1.861156378244313e-24im
|
||||||
-1.7864824922599887e-6 + 1.9037001727701538e-25im
|
-1.786482492259989e-6 + 1.9925574651342477e-25im
|
||||||
-1.543158528170733e-7 + 2.862618452866008e-26im
|
-1.543158528170733e-7 - 4.2577452591128694e-26im
|
||||||
|
|
||||||
Computing lT...
|
Computing lT...
|
||||||
|
|
||||||
lT =
|
lT =
|
||||||
0.0 + 0.0im
|
0.0 + 0.0im
|
||||||
-1.5314862576847136e-7 + 3.022494262990794e-26im
|
-1.531486257684714e-7 + 5.2170001198616245e-27im
|
||||||
-1.781979951825098e-6 + 3.204247815553729e-26im
|
-1.781979951825098e-6 - 8.427979548473216e-26im
|
||||||
-2.0486213262577787e-5 - 4.26515003347655e-25im
|
-2.0486213262577787e-5 + 1.921471631242971e-24im
|
||||||
-0.00022976131918157333 - 4.718376077437893e-23im
|
-0.00022976131918157338 - 2.1437764208665974e-23im
|
||||||
-0.002600005011614236 + 1.9521462456251453e-22im
|
-0.0026000050116142363 - 5.018229049488372e-23im
|
||||||
-0.028781733302634345 - 3.211666591672554e-21im
|
-0.02878173330263435 - 1.5793635162345812e-21im
|
||||||
-0.3047801284340767 + 6.988021814847987e-21im
|
-0.3047801284340767 + 1.4470146182260945e-20im
|
||||||
|
Creating new BDIO output file ./output/Template_run.bdio
|
||||||
────────────────────────────────────────────────────────────────────────────────
|
────────────────────────────────────────────────────────────────────────────────
|
||||||
Time Allocations
|
Time Allocations
|
||||||
─────────────────────── ────────────────────────
|
─────────────────────── ────────────────────────
|
||||||
Tot / % measured: 108s / 89.1% 11.3GiB / 90.7%
|
Tot / % measured: 86.1s / 76.0% 7.23GiB / 76.4%
|
||||||
|
|
||||||
Section ncalls time %tot avg alloc %tot avg
|
Section ncalls time %tot avg alloc %tot avg
|
||||||
────────────────────────────────────────────────────────────────────────────────
|
────────────────────────────────────────────────────────────────────────────────
|
||||||
Propagator computa... 1 52.8s 54.7% 52.8s 5.99GiB 58.3% 5.99GiB
|
Propagator computa... 1 27.1s 41.4% 27.1s 2.09GiB 37.9% 2.09GiB
|
||||||
g5Dw 12 2.80s 2.9% 233ms 220MiB 2.1% 18.3MiB
|
g5Dw 12 2.27s 3.5% 189ms 138MiB 2.4% 11.5MiB
|
||||||
DwdagDw 600 732ms 0.8% 1.22ms 14.7MiB 0.1% 25.1KiB
|
DwdagDw 600 743ms 1.1% 1.24ms 11.0MiB 0.2% 18.7KiB
|
||||||
g5Dw 1.20k 694ms 0.7% 578μs 11.1MiB 0.1% 9.46KiB
|
g5Dw 1.20k 705ms 1.1% 587μs 8.58MiB 0.2% 7.32KiB
|
||||||
SF boundary fix 1.20k 37.2ms 0.0% 31.0μs 3.41MiB 0.0% 2.91KiB
|
SF boundary fix 1.20k 37.7ms 0.1% 31.4μs 2.16MiB 0.0% 1.84KiB
|
||||||
SF boundary fix 12 702ms 0.7% 58.5ms 29.9MiB 0.3% 2.49MiB
|
SF boundary fix 12 391ms 0.6% 32.6ms 61.6MiB 1.1% 5.13MiB
|
||||||
Gauge field reading 1 19.7s 20.4% 19.7s 1.52GiB 14.8% 1.52GiB
|
Gauge field reading 1 19.9s 30.5% 19.9s 1.60GiB 29.0% 1.60GiB
|
||||||
Csw computation 1 16.5s 17.1% 16.5s 1.14GiB 11.1% 1.14GiB
|
Csw computation 1 16.4s 25.0% 16.4s 1.08GiB 19.5% 1.08GiB
|
||||||
Gauge measurements 1 18.2s 18.8% 18.2s 1.83GiB 17.8% 1.83GiB
|
Gauge measurements 1 12.4s 18.9% 12.4s 1.04GiB 18.8% 1.04GiB
|
||||||
Qtop measurement 1 8.15s 8.5% 8.15s 691MiB 6.6% 691MiB
|
Qtop measurement 1 5.19s 7.9% 5.19s 380MiB 6.7% 380MiB
|
||||||
Plaquette measur... 1 1.96s 2.0% 1.96s 157MiB 1.5% 157MiB
|
Plaquette measur... 1 1.46s 2.2% 1.46s 71.4MiB 1.3% 71.4MiB
|
||||||
SF coupling meas... 1 1.82s 1.9% 1.82s 137MiB 1.3% 137MiB
|
SF coupling meas... 1 1.39s 2.1% 1.39s 44.1MiB 0.8% 44.1MiB
|
||||||
Imput reading and ... 1 4.46s 4.6% 4.46s 529MiB 5.0% 529MiB
|
Imput reading and ... 1 4.26s 6.5% 4.26s 332MiB 5.9% 332MiB
|
||||||
Allocating Dirac... 1 157ms 0.2% 157ms 1.15MiB 0.0% 1.15MiB
|
Allocating Dirac... 1 124ms 0.2% 124ms 154KiB 0.0% 154KiB
|
||||||
Allocating YMWor... 1 46.5μs 0.0% 46.5μs 1.03KiB 0.0% 1.03KiB
|
Allocating YMWor... 1 50.5μs 0.0% 50.5μs 1.50KiB 0.0% 1.50KiB
|
||||||
Correlators contra... 1 1.39s 1.4% 1.39s 428MiB 4.1% 428MiB
|
Correlators contra... 1 1.32s 2.0% 1.32s 429MiB 7.6% 429MiB
|
||||||
|
Correlators saving 1 505ms 0.8% 505ms 46.2MiB 0.8% 46.2MiB
|
||||||
────────────────────────────────────────────────────────────────────────────────
|
────────────────────────────────────────────────────────────────────────────────
|
||||||
|
|
|
||||||
|
|
@ -123,7 +123,7 @@ function lT_fun()
|
||||||
return cf
|
return cf
|
||||||
end
|
end
|
||||||
|
|
||||||
f1_fun() = (norm2(C11) + norm2(C21) + norm2(C31) + norm2(C12) + norm2(C22) + norm2(C32))
|
f1_fun() = complex(norm2(C11) + norm2(C21) + norm2(C31) + norm2(C12) + norm2(C22) + norm2(C32))
|
||||||
|
|
||||||
k1_fun() = (dot(C11,dmul(Gamma{6},imm(C12))) + dot(C11,dmul(Gamma{7},-C12)) + dot(C11,dmul(Gamma{8}, imm(C11)))
|
k1_fun() = (dot(C11,dmul(Gamma{6},imm(C12))) + dot(C11,dmul(Gamma{7},-C12)) + dot(C11,dmul(Gamma{8}, imm(C11)))
|
||||||
+dot(C12,dmul(Gamma{6},imm(C11))) + dot(C12,dmul(Gamma{7}, C11)) + dot(C12,dmul(Gamma{8},mimm(C12)))
|
+dot(C12,dmul(Gamma{6},imm(C11))) + dot(C12,dmul(Gamma{7}, C11)) + dot(C12,dmul(Gamma{8},mimm(C12)))
|
||||||
|
|
|
||||||
|
|
@ -92,7 +92,9 @@ function gauge_meas()
|
||||||
println(log_file,"\nGauge meassurements:\n")
|
println(log_file,"\nGauge meassurements:\n")
|
||||||
println(log_file,"Plaquette: ", plaquette(U, lp, gp, ymws))
|
println(log_file,"Plaquette: ", plaquette(U, lp, gp, ymws))
|
||||||
println(log_file,"Topological charge: ", Qtop(U, gp, lp, ymws))
|
println(log_file,"Topological charge: ", Qtop(U, gp, lp, ymws))
|
||||||
println(log_file,"SF coupling (dS/deta): ",sfcoupling(U, lp, gp, ymws)[1])
|
sfc = sfcoupling(U, lp, gp, ymws)
|
||||||
|
println(log_file,"SF coupling (dS / deta): ",sfc[1])
|
||||||
|
println(log_file,"SF coupling (d^2S / deta dnu): ",sfc[2])
|
||||||
flush(log_file)
|
flush(log_file)
|
||||||
|
|
||||||
return nothing
|
return nothing
|
||||||
|
|
|
||||||
Loading…
Add table
Add a link
Reference in a new issue