From 5790c8f8639714cfd2dafc0c5a7d19017ebdea96 Mon Sep 17 00:00:00 2001 From: "Fernando P. Panadero" Date: Tue, 27 Feb 2024 09:20:05 +0100 Subject: [PATCH] Minor changes. --- output/Template_run.log | 161 +++++++++++++++++++++------------------- src/corr.jl | 2 +- src/meas.jl | 4 +- 3 files changed, 87 insertions(+), 80 deletions(-) diff --git a/output/Template_run.log b/output/Template_run.log index e3688fe..aa94ff3 100644 --- a/output/Template_run.log +++ b/output/Template_run.log @@ -28,35 +28,38 @@ tolerance = 1.0e-13 maxiter = 1000 Reading gauge field with LGPU native convention from: ./cnfg/randfield... +MD5 checksum of gauge config: 89c5bd25d1cd2c74e59f0ecfcf95ebef + Gauge meassurements: Plaquette: 0.3306405864818687 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 -CG converged for c=1 and s=1 after 50 iterations with absolute residue 0.0009654241012493543 -CG converged for c=2 and s=1 after 50 iterations with absolute residue 0.0009792505060407602 -CG converged for c=3 and s=1 after 50 iterations with absolute residue 0.0009410537206344275 -CG converged for c=1 and s=2 after 50 iterations with absolute residue 0.0009706380660725444 -CG converged for c=2 and s=2 after 50 iterations with absolute residue 0.0009348723257035598 -CG converged for c=3 and s=2 after 50 iterations with absolute residue 0.0009646157878217441 +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.0009792505060407621 +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.000970638066072543 +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.000964615787821745 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=2 and s=1 after 50 iterations with absolute residue 0.0009591325218333536 -CG converged for c=3 and s=1 after 50 iterations with absolute residue 0.000993180324656404 -CG converged for c=1 and s=2 after 50 iterations with absolute residue 0.0009696740933237723 -CG converged for c=2 and s=2 after 50 iterations with absolute residue 0.0009753894671580818 -CG converged for c=3 and s=2 after 50 iterations with absolute residue 0.0009677391305951339 +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.0009591325218333531 +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.0009696740933237713 +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.0009677391305951349 Computing correlators Computing f1... f1 = -7.968153189691683e-10 +7.968153189691686e-10 + 0.0im Computing k1... @@ -70,115 +73,117 @@ fP = 0.39164297606334497 + 0.0im 0.03838253084495811 + 0.0im 0.0036829603591832434 + 0.0im -0.00034990928932360324 + 0.0im +0.0003499092893236032 + 0.0im 3.1848709670960246e-5 + 0.0im -2.83087866678673e-6 + 0.0im +2.8308786667867295e-6 + 0.0im 2.552038846719931e-7 + 0.0im Computing fA... fA = 0.0 + 0.0im --0.28022427481948287 + 3.9466441386130056e-20im --0.028432536079057023 - 9.446402595355576e-22im --0.002717462099554097 + 4.301339185275744e-23im --0.0002553997594898297 - 2.8822074468644558e-24im --2.3190475107607993e-5 - 1.34416849539867e-24im --2.078551613006522e-6 - 3.0342745858421e-25im --1.8801834874904572e-7 - 1.7443292336247013e-26im +-0.28022427481948287 + 2.968585782330306e-20im +-0.028432536079057023 + 1.435985543392056e-21im +-0.002717462099554098 + 1.8611563782443123e-23im +-0.0002553997594898297 + 2.1067256225959924e-23im +-2.3190475107607993e-5 + 3.723928343622517e-24im +-2.0785516130065226e-6 - 1.0905213153775267e-25im +-1.880183487490457e-7 + 3.6085653248429964e-26im Computing gP... gP = 0.0 + 0.0im -2.538597063821986e-7 + 0.0im +2.5385970638219865e-7 + 0.0im 2.8229409631152233e-6 + 0.0im -3.096526353494133e-5 + 0.0im +3.0965263534941325e-5 + 0.0im 0.0003331843094846162 + 0.0im 0.0036254847521787863 + 0.0im 0.03837127818992878 + 0.0im -0.39410791485294716 + 0.0im +0.39410791485294727 + 0.0im Computing gA... gA = 0.0 + 0.0im --1.8696507552685477e-7 - 6.096317075547979e-27im --2.0658936442177055e-6 - 3.842068152788416e-26im --2.2674085372631826e-5 - 8.56261180963095e-25im --0.00024433351353064797 + 9.060212646869881e-24im --0.0026482638745495414 + 2.161009350294785e-22im --0.028194731180271545 + 5.375019620369574e-21im --0.28340818434504744 + 2.511320339710992e-20im +-1.8696507552685472e-7 - 6.203602158658037e-27im +-2.0658936442177055e-6 - 1.7405427012796404e-25im +-2.2674085372631822e-5 + 2.1422685395416303e-24im +-0.000244333513530648 + 5.8458404852771e-23im +-0.0026482638745495414 + 3.7698756417104236e-22im +-0.028194731180271545 + 6.309733712554495e-21im +-0.28340818434504744 + 1.9051623868321333e-20im Computing kV... kV = 0.0 + 0.0im -0.33632172225753315 - 1.0340860564396276e-20im -0.03186802210515694 - 9.926167350636303e-23im -0.0029308316253718663 - 1.1084220208210557e-22im -0.0002663612299706971 + 3.171720661258017e-23im -2.3221979105382962e-5 - 1.3495537858449775e-24im -1.9759435694009697e-6 + 3.989827059409007e-25im -1.7038427603748213e-7 + 2.5655860266867883e-26im +0.33632172225753315 + 1.0940842235368113e-21im +0.03186802210515694 + 1.0742318888355314e-21im +0.0029308316253718663 - 9.760731228125728e-23im +0.00026636122997069706 + 1.672886624241357e-23im +2.3221979105382955e-5 + 1.2580038482577253e-24im +1.9759435694009697e-6 - 2.610519543848177e-25im +1.7038427603748213e-7 + 1.5226067285305813e-26im Computing lV... lV = 0.0 + 0.0im -1.692538762882869e-7 - 1.2470313189733742e-26im -1.972297868357158e-6 + 1.9703431420432253e-25im -2.2655589268307515e-5 + 4.094974855373194e-24im -0.0002539424459369091 - 2.540133797714921e-23im -0.0028788471359388916 - 4.1193594505140788e-22im -0.031803280548268274 + 2.5322755818956662e-21im -0.3390357554359771 - 7.90564084104014e-21im +1.692538762882869e-7 + 3.2328571710497395e-26im +1.9722978683571574e-6 + 2.1049754032009423e-25im +2.2655589268307515e-5 - 3.08684848382419e-24im +0.0002539424459369091 + 3.9704669402545317e-23im +0.002878847135938892 + 5.79026428787119e-23im +0.03180328054826828 + 2.4749243927586577e-21im +0.339035755435977 + 2.883441324611509e-20im Computing kT... kT = 0.0 + 0.0im --0.30173882074287517 + 5.293955920339405e-22im --0.028918463611070797 - 1.548482106699267e-21im --0.002655774263798736 + 6.81045371001993e-23im --0.00024049035597286907 + 3.794691060087014e-23im --2.0961481048454104e-5 - 1.4734154661100807e-24im --1.7864824922599887e-6 + 1.9037001727701538e-25im --1.543158528170733e-7 + 2.862618452866008e-26im +-0.30173882074287517 + 9.03501810404587e-21im +-0.0289184636110708 + 1.6234798155707435e-21im +-0.002655774263798736 + 2.451211881865472e-22im +-0.00024049035597286904 + 5.2388105461691765e-24im +-2.0961481048454104e-5 - 1.861156378244313e-24im +-1.786482492259989e-6 + 1.9925574651342477e-25im +-1.543158528170733e-7 - 4.2577452591128694e-26im Computing lT... lT = 0.0 + 0.0im --1.5314862576847136e-7 + 3.022494262990794e-26im --1.781979951825098e-6 + 3.204247815553729e-26im --2.0486213262577787e-5 - 4.26515003347655e-25im --0.00022976131918157333 - 4.718376077437893e-23im --0.002600005011614236 + 1.9521462456251453e-22im --0.028781733302634345 - 3.211666591672554e-21im --0.3047801284340767 + 6.988021814847987e-21im +-1.531486257684714e-7 + 5.2170001198616245e-27im +-1.781979951825098e-6 - 8.427979548473216e-26im +-2.0486213262577787e-5 + 1.921471631242971e-24im +-0.00022976131918157338 - 2.1437764208665974e-23im +-0.0026000050116142363 - 5.018229049488372e-23im +-0.02878173330263435 - 1.5793635162345812e-21im +-0.3047801284340767 + 1.4470146182260945e-20im +Creating new BDIO output file ./output/Template_run.bdio ──────────────────────────────────────────────────────────────────────────────── 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 ──────────────────────────────────────────────────────────────────────────────── - Propagator computa... 1 52.8s 54.7% 52.8s 5.99GiB 58.3% 5.99GiB - g5Dw 12 2.80s 2.9% 233ms 220MiB 2.1% 18.3MiB - DwdagDw 600 732ms 0.8% 1.22ms 14.7MiB 0.1% 25.1KiB - g5Dw 1.20k 694ms 0.7% 578μs 11.1MiB 0.1% 9.46KiB - SF boundary fix 1.20k 37.2ms 0.0% 31.0μs 3.41MiB 0.0% 2.91KiB - SF boundary fix 12 702ms 0.7% 58.5ms 29.9MiB 0.3% 2.49MiB - Gauge field reading 1 19.7s 20.4% 19.7s 1.52GiB 14.8% 1.52GiB - Csw computation 1 16.5s 17.1% 16.5s 1.14GiB 11.1% 1.14GiB - Gauge measurements 1 18.2s 18.8% 18.2s 1.83GiB 17.8% 1.83GiB - Qtop measurement 1 8.15s 8.5% 8.15s 691MiB 6.6% 691MiB - Plaquette measur... 1 1.96s 2.0% 1.96s 157MiB 1.5% 157MiB - SF coupling meas... 1 1.82s 1.9% 1.82s 137MiB 1.3% 137MiB - Imput reading and ... 1 4.46s 4.6% 4.46s 529MiB 5.0% 529MiB - Allocating Dirac... 1 157ms 0.2% 157ms 1.15MiB 0.0% 1.15MiB - Allocating YMWor... 1 46.5μs 0.0% 46.5μs 1.03KiB 0.0% 1.03KiB - Correlators contra... 1 1.39s 1.4% 1.39s 428MiB 4.1% 428MiB + Propagator computa... 1 27.1s 41.4% 27.1s 2.09GiB 37.9% 2.09GiB + g5Dw 12 2.27s 3.5% 189ms 138MiB 2.4% 11.5MiB + DwdagDw 600 743ms 1.1% 1.24ms 11.0MiB 0.2% 18.7KiB + g5Dw 1.20k 705ms 1.1% 587μs 8.58MiB 0.2% 7.32KiB + SF boundary fix 1.20k 37.7ms 0.1% 31.4μs 2.16MiB 0.0% 1.84KiB + SF boundary fix 12 391ms 0.6% 32.6ms 61.6MiB 1.1% 5.13MiB + Gauge field reading 1 19.9s 30.5% 19.9s 1.60GiB 29.0% 1.60GiB + Csw computation 1 16.4s 25.0% 16.4s 1.08GiB 19.5% 1.08GiB + Gauge measurements 1 12.4s 18.9% 12.4s 1.04GiB 18.8% 1.04GiB + Qtop measurement 1 5.19s 7.9% 5.19s 380MiB 6.7% 380MiB + Plaquette measur... 1 1.46s 2.2% 1.46s 71.4MiB 1.3% 71.4MiB + SF coupling meas... 1 1.39s 2.1% 1.39s 44.1MiB 0.8% 44.1MiB + Imput reading and ... 1 4.26s 6.5% 4.26s 332MiB 5.9% 332MiB + Allocating Dirac... 1 124ms 0.2% 124ms 154KiB 0.0% 154KiB + Allocating YMWor... 1 50.5μs 0.0% 50.5μs 1.50KiB 0.0% 1.50KiB + 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 ──────────────────────────────────────────────────────────────────────────────── diff --git a/src/corr.jl b/src/corr.jl index 462c7f4..4fbd6bf 100644 --- a/src/corr.jl +++ b/src/corr.jl @@ -123,7 +123,7 @@ function lT_fun() return cf 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))) +dot(C12,dmul(Gamma{6},imm(C11))) + dot(C12,dmul(Gamma{7}, C11)) + dot(C12,dmul(Gamma{8},mimm(C12))) diff --git a/src/meas.jl b/src/meas.jl index 84ef6ad..3d8b3c6 100644 --- a/src/meas.jl +++ b/src/meas.jl @@ -92,7 +92,9 @@ function gauge_meas() println(log_file,"\nGauge meassurements:\n") println(log_file,"Plaquette: ", plaquette(U, lp, gp, 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) return nothing