Skip to content

F
FerFlow.jl
Commit dc2fcd73 authored by Fernando Pérez Panadero's avatar Fernando Pérez Panadero
Browse files
Options

One derivative with AD.

parent 7b2d14d5
Hide whitespace changes
Inline Side-by-side
Showing with 610 additions and 3244 deletions
Manifest.toml
View file @ dc2fcd73
......@@ -2,7 +2,7 @@
julia_version = "1.8.5"
manifest_format = "2.0"
project_hash = "7a31d20d9e16a1e751c8e1b248aa5f00852fd2de"
project_hash = "d7d01fc8c20a58b96112d1ce64384f703066537c"
[[deps.ADerrors]]
deps = ["BDIO", "FFTW", "FastGaussQuadrature", "ForwardDiff", "LaTeXStrings", "LinearAlgebra", "Nettle", "PGFPlotsX", "Plots", "Printf", "QuadGK", "Roots", "Statistics", "Test", "UnicodePlots"]
......@@ -88,9 +88,9 @@ version = "0.1.9"
[[deps.Bzip2_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
git-tree-sha1 = "9e2a6b69137e6969bab0152632dcb3bc108c8bdd"
git-tree-sha1 = "8873e196c2eb87962a2048b3b8e08946535864a1"
uuid = "6e34b625-4abd-537c-b88f-471c36dfa7a0"
version = "1.0.8+1"
version = "1.0.8+4"
[[deps.CEnum]]
git-tree-sha1 = "389ad5c84de1ae7cf0e28e381131c98ea87d54fc"
......@@ -139,6 +139,12 @@ git-tree-sha1 = "799b25ca3a8a24936ae7b5c52ad194685fc3e6ef"
uuid = "9e997f8a-9a97-42d5-a9f1-ce6bfc15e2c0"
version = "0.1.9"
[[deps.CodeTracking]]
deps = ["InteractiveUtils", "UUIDs"]
git-tree-sha1 = "7eee164f122511d3e4e1ebadb7956939ea7e1c77"
uuid = "da1fd8a2-8d9e-5ec2-8556-3022fb5608a2"
version = "1.3.6"
[[deps.CodecZlib]]
deps = ["TranscodingStreams", "Zlib_jll"]
git-tree-sha1 = "bce6804e5e6044c6daab27bb533d1295e4a2e759"
......@@ -218,6 +224,12 @@ git-tree-sha1 = "249fe38abf76d48563e2f4556bebd215aa317e15"
uuid = "a8cc5b0e-0ffa-5ad4-8c14-923d3ee1735f"
version = "4.1.1"
[[deps.Cthulhu]]
deps = ["CodeTracking", "FoldingTrees", "InteractiveUtils", "JuliaSyntax", "PrecompileTools", "Preferences", "REPL", "TypedSyntax", "UUIDs", "Unicode", "WidthLimitedIO"]
git-tree-sha1 = "42df01693cfacf5599819a1f0bcc767c7759d0ab"
uuid = "f68482b8-f384-11e8-15f7-abe071a5a75f"
version = "2.9.6"
[[deps.DataAPI]]
git-tree-sha1 = "abe83f3a2f1b857aac70ef8b269080af17764bbe"
uuid = "9a962f9c-6df0-11e9-0e5d-c546b8b5ee8a"
......@@ -297,9 +309,9 @@ version = "0.7.8"
[[deps.EpollShim_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl"]
git-tree-sha1 = "8e9441ee83492030ace98f9789a654a6d0b1f643"
git-tree-sha1 = "8a4be429317c42cfae6a7fc03c31bad1970c310d"
uuid = "2702e6a9-849d-5ed8-8c21-79e8b8f9ee43"
version = "0.0.20230411+0"
version = "0.0.20230411+1"
[[deps.ExceptionUnwrapping]]
deps = ["Test"]
......@@ -338,9 +350,9 @@ version = "1.8.0"
[[deps.FFTW_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
git-tree-sha1 = "c6033cc3892d0ef5bb9cd29b7f2f0331ea5184ea"
git-tree-sha1 = "4d81ed14783ec49ce9f2e168208a12ce1815aa25"
uuid = "f5851436-0d7a-5f13-b9de-f02708fd171a"
version = "3.3.10+0"
version = "3.3.10+3"
[[deps.FastGaussQuadrature]]
deps = ["LinearAlgebra", "SpecialFunctions", "StaticArrays"]
......@@ -357,12 +369,26 @@ git-tree-sha1 = "05882d6995ae5c12bb5f36dd2ed3f61c98cbb172"
uuid = "53c48c17-4a7d-5ca2-90c5-79b7896eea93"
version = "0.8.5"
[[deps.FoldingTrees]]
deps = ["AbstractTrees", "REPL"]
git-tree-sha1 = "c1b0164369256b26f71d9830df9000a9c39757fc"
uuid = "1eca21be-9b9b-4ed8-839a-6d8ae26b1781"
version = "1.2.2"
[[deps.Fontconfig_jll]]
deps = ["Artifacts", "Bzip2_jll", "Expat_jll", "FreeType2_jll", "JLLWrappers", "Libdl", "Libuuid_jll", "Zlib_jll"]
git-tree-sha1 = "db16beca600632c95fc8aca29890d83788dd8b23"
uuid = "a3f928ae-7b40-5064-980b-68af3947d34b"
version = "2.13.96+0"
[[deps.FormalSeries]]
deps = ["LinearAlgebra", "PrecompileTools", "StaticArrays"]
git-tree-sha1 = "42e910d9de5af437f282e8561e6a386d951fc4c1"
repo-rev = "master"
repo-url = "https://igit.ific.uv.es/alramos/formalseries.jl.git"
uuid = "a55ae9b5-0907-4af6-81ab-c6ed2f8e0777"
version = "1.5.0"
[[deps.Format]]
git-tree-sha1 = "9c68794ef81b08086aeb32eeaf33531668d5f5fc"
uuid = "1fa38f19-a742-5d3f-a2b9-30dd87b9d5f8"
......@@ -392,9 +418,9 @@ uuid = "9fa8497b-333b-5362-9e8d-4d0656e87820"
[[deps.GLFW_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl", "Libglvnd_jll", "Xorg_libXcursor_jll", "Xorg_libXi_jll", "Xorg_libXinerama_jll", "Xorg_libXrandr_jll", "libdecor_jll", "xkbcommon_jll"]
git-tree-sha1 = "532f9126ad901533af1d4f5c198867227a7bb077"
git-tree-sha1 = "fcb0584ff34e25155876418979d4c8971243bb89"
uuid = "0656b61e-2033-5cc2-a64a-77c0f6c09b89"
version = "3.4.0+1"
version = "3.4.0+2"
[[deps.GMP_jll]]
deps = ["Artifacts", "Libdl"]
......@@ -457,9 +483,9 @@ version = "2.80.2+0"
[[deps.Graphite2_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
git-tree-sha1 = "344bf40dcab1073aca04aa0df4fb092f920e4011"
git-tree-sha1 = "01979f9b37367603e2848ea225918a3b3861b606"
uuid = "3b182d85-2403-5c21-9c21-1e1f0cc25472"
version = "1.3.14+0"
version = "1.3.14+1"
[[deps.Grisu]]
git-tree-sha1 = "53bb909d1151e57e2484c3d1b53e19552b887fb2"
......@@ -551,6 +577,11 @@ git-tree-sha1 = "af433a10f3942e882d3c671aacb203e006a5808f"
uuid = "9c1d0b0a-7046-5b2e-a33f-ea22f176ac7e"
version = "0.2.1+0"
[[deps.JuliaSyntax]]
git-tree-sha1 = "937da4713526b96ac9a178e2035019d3b78ead4a"
uuid = "70703baa-626e-46a2-a12c-08ffd08c73b4"
version = "0.4.10"
[[deps.KernelAbstractions]]
deps = ["Adapt", "Atomix", "EnzymeCore", "InteractiveUtils", "LinearAlgebra", "MacroTools", "PrecompileTools", "Requires", "SparseArrays", "StaticArrays", "UUIDs", "UnsafeAtomics", "UnsafeAtomicsLLVM"]
git-tree-sha1 = "cb1cff88ef2f3a157cbad75bbe6b229e1975e498"
......@@ -594,9 +625,9 @@ version = "17.0.6+0"
[[deps.LZO_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl"]
git-tree-sha1 = "70c5da094887fd2cae843b8db33920bac4b6f07d"
git-tree-sha1 = "854a9c268c43b77b0a27f22d7fab8d33cdb3a731"
uuid = "dd4b983a-f0e5-5f8d-a1b7-129d4a5fb1ac"
version = "2.10.2+0"
version = "2.10.2+3"
[[deps.LaTeXStrings]]
git-tree-sha1 = "50901ebc375ed41dbf8058da26f9de442febbbec"
......@@ -610,10 +641,8 @@ uuid = "23fbe1c1-3f47-55db-b15f-69d7ec21a316"
version = "0.16.5"
[[deps.LatticeGPU]]
deps = ["BDIO", "CUDA", "Random", "TimerOutputs"]
git-tree-sha1 = "4527979a0779235ac1a16cbf6b0592f9ae3352f9"
repo-rev = "master"
repo-url = "https://igit.ific.uv.es/alramos/latticegpu.jl.git"
deps = ["BDIO", "CUDA", "FormalSeries", "Random", "TimerOutputs"]
path = "/home/fperez/.julia/dev/LatticeGPU"
uuid = "958c3683-801a-4582-9cfa-2d6e2ae1763b"
version = "0.2.0"
......@@ -650,9 +679,9 @@ uuid = "8f399da3-3557-5675-b5ff-fb832c97cbdb"
[[deps.Libffi_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
git-tree-sha1 = "0b4a5d71f3e5200a7dff793393e09dfc2d874290"
git-tree-sha1 = "27ecae93dd25ee0909666e6835051dd684cc035e"
uuid = "e9f186c6-92d2-5b65-8a66-fee21dc1b490"
version = "3.2.2+1"
version = "3.2.2+2"
[[deps.Libgcrypt_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl", "Libgpg_error_jll"]
......@@ -674,9 +703,9 @@ version = "1.49.0+0"
[[deps.Libiconv_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl"]
git-tree-sha1 = "f9557a255370125b405568f9767d6d195822a175"
git-tree-sha1 = "61dfdba58e585066d8bce214c5a51eaa0539f269"
uuid = "94ce4f54-9a6c-5748-9c1c-f9c7231a4531"
version = "1.17.0+0"
version = "1.17.0+1"
[[deps.Libmount_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl"]
......@@ -836,15 +865,15 @@ version = "1.4.3"
[[deps.OpenSSL_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl"]
git-tree-sha1 = "1b35263570443fdd9e76c76b7062116e2f374ab8"
git-tree-sha1 = "7493f61f55a6cce7325f197443aa80d32554ba10"
uuid = "458c3c95-2e84-50aa-8efc-19380b2a3a95"
version = "3.0.15+0"
version = "3.0.15+3"
[[deps.OpenSpecFun_jll]]
deps = ["Artifacts", "CompilerSupportLibraries_jll", "JLLWrappers", "Libdl", "Pkg"]
git-tree-sha1 = "13652491f6856acfd2db29360e1bbcd4565d04f1"
uuid = "efe28fd5-8261-553b-a9e1-b2916fc3738e"
version = "0.5.5+0"
version = "0.5.5+2"
[[deps.Opus_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl"]
......@@ -1169,6 +1198,12 @@ git-tree-sha1 = "e84b3a11b9bece70d14cce63406bbc79ed3464d2"
uuid = "3bb67fe8-82b1-5028-8e26-92a6c54297fa"
version = "0.11.2"
[[deps.TypedSyntax]]
deps = ["CodeTracking", "JuliaSyntax"]
git-tree-sha1 = "d07928dd489d9ea2d6527e3acf9cb401f5bf782a"
uuid = "d265eb64-f81a-44ad-a842-4247ee1503de"
version = "1.2.4"
[[deps.URIs]]
git-tree-sha1 = "67db6cc7b3821e19ebe75791a9dd19c9b1188f2b"
uuid = "5c2747f8-b7ea-4ff2-ba2e-563bfd36b1d4"
......@@ -1234,9 +1269,9 @@ version = "1.3.243+0"
[[deps.Wayland_jll]]
deps = ["Artifacts", "EpollShim_jll", "Expat_jll", "JLLWrappers", "Libdl", "Libffi_jll", "Pkg", "XML2_jll"]
git-tree-sha1 = "7558e29847e99bc3f04d6569e82d0f5c54460703"
git-tree-sha1 = "85c7811eddec9e7f22615371c3cc81a504c508ee"
uuid = "a2964d1f-97da-50d4-b82a-358c7fce9d89"
version = "1.21.0+1"
version = "1.21.0+2"
[[deps.Wayland_protocols_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
......@@ -1244,6 +1279,12 @@ git-tree-sha1 = "93f43ab61b16ddfb2fd3bb13b3ce241cafb0e6c9"
uuid = "2381bf8a-dfd0-557d-9999-79630e7b1b91"
version = "1.31.0+0"
[[deps.WidthLimitedIO]]
deps = ["Unicode"]
git-tree-sha1 = "71142739e695823729a335e9bc124ef41ec1433a"
uuid = "b8c1c048-cf81-46c6-9da0-18c1d99e41f2"
version = "1.0.1"
[[deps.XML2_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl", "Libiconv_jll", "Zlib_jll"]
git-tree-sha1 = "1165b0443d0eca63ac1e32b8c0eb69ed2f4f8127"
......@@ -1276,15 +1317,15 @@ version = "1.2.4+0"
[[deps.Xorg_libX11_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl", "Xorg_libxcb_jll", "Xorg_xtrans_jll"]
git-tree-sha1 = "afead5aba5aa507ad5a3bf01f58f82c8d1403495"
git-tree-sha1 = "9dafcee1d24c4f024e7edc92603cedba72118283"
uuid = "4f6342f7-b3d2-589e-9d20-edeb45f2b2bc"
version = "1.8.6+0"
version = "1.8.6+3"
[[deps.Xorg_libXau_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl"]
git-tree-sha1 = "6035850dcc70518ca32f012e46015b9beeda49d8"
git-tree-sha1 = "2b0e27d52ec9d8d483e2ca0b72b3cb1a8df5c27a"
uuid = "0c0b7dd1-d40b-584c-a123-a41640f87eec"
version = "1.0.11+0"
version = "1.0.11+3"
[[deps.Xorg_libXcursor_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg", "Xorg_libXfixes_jll", "Xorg_libXrender_jll"]
......@@ -1294,15 +1335,15 @@ version = "1.2.0+4"
[[deps.Xorg_libXdmcp_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl"]
git-tree-sha1 = "34d526d318358a859d7de23da945578e8e8727b7"
git-tree-sha1 = "02054ee01980c90297412e4c809c8694d7323af3"
uuid = "a3789734-cfe1-5b06-b2d0-1dd0d9d62d05"
version = "1.1.4+0"
version = "1.1.4+3"
[[deps.Xorg_libXext_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl", "Xorg_libX11_jll"]
git-tree-sha1 = "d2d1a5c49fae4ba39983f63de6afcbea47194e85"
git-tree-sha1 = "d7155fea91a4123ef59f42c4afb5ab3b4ca95058"
uuid = "1082639a-0dae-5f34-9b06-72781eeb8cb3"
version = "1.3.6+0"
version = "1.3.6+3"
[[deps.Xorg_libXfixes_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg", "Xorg_libX11_jll"]
......@@ -1330,27 +1371,27 @@ version = "1.5.2+4"
[[deps.Xorg_libXrender_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl", "Xorg_libX11_jll"]
git-tree-sha1 = "47e45cd78224c53109495b3e324df0c37bb61fbe"
git-tree-sha1 = "a490c6212a0e90d2d55111ac956f7c4fa9c277a6"
uuid = "ea2f1a96-1ddc-540d-b46f-429655e07cfa"
version = "0.9.11+0"
version = "0.9.11+1"
[[deps.Xorg_libpthread_stubs_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl"]
git-tree-sha1 = "8fdda4c692503d44d04a0603d9ac0982054635f9"
git-tree-sha1 = "fee57a273563e273f0f53275101cd41a8153517a"
uuid = "14d82f49-176c-5ed1-bb49-ad3f5cbd8c74"
version = "0.1.1+0"
version = "0.1.1+3"
[[deps.Xorg_libxcb_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl", "XSLT_jll", "Xorg_libXau_jll", "Xorg_libXdmcp_jll", "Xorg_libpthread_stubs_jll"]
git-tree-sha1 = "bcd466676fef0878338c61e655629fa7bbc69d8e"
git-tree-sha1 = "1a74296303b6524a0472a8cb12d3d87a78eb3612"
uuid = "c7cfdc94-dc32-55de-ac96-5a1b8d977c5b"
version = "1.17.0+0"
version = "1.17.0+3"
[[deps.Xorg_libxkbfile_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl", "Xorg_libX11_jll"]
git-tree-sha1 = "730eeca102434283c50ccf7d1ecdadf521a765a4"
git-tree-sha1 = "dbc53e4cf7701c6c7047c51e17d6e64df55dca94"
uuid = "cc61e674-0454-545c-8b26-ed2c68acab7a"
version = "1.1.2+0"
version = "1.1.2+1"
[[deps.Xorg_xcb_util_cursor_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl", "Xorg_xcb_util_image_jll", "Xorg_xcb_util_jll", "Xorg_xcb_util_renderutil_jll"]
......@@ -1390,9 +1431,9 @@ version = "0.4.1+1"
[[deps.Xorg_xkbcomp_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl", "Xorg_libxkbfile_jll"]
git-tree-sha1 = "330f955bc41bb8f5270a369c473fc4a5a4e4d3cb"
git-tree-sha1 = "ab2221d309eda71020cdda67a973aa582aa85d69"
uuid = "35661453-b289-5fab-8a00-3d9160c6a3a4"
version = "1.4.6+0"
version = "1.4.6+1"
[[deps.Xorg_xkeyboard_config_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl", "Xorg_xkbcomp_jll"]
......@@ -1402,9 +1443,9 @@ version = "2.39.0+0"
[[deps.Xorg_xtrans_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl"]
git-tree-sha1 = "e92a1a012a10506618f10b7047e478403a046c77"
git-tree-sha1 = "b9ead2d2bdb27330545eb14234a2e300da61232e"
uuid = "c5fb5394-a638-5e4d-96e5-b29de1b5cf10"
version = "1.5.0+0"
version = "1.5.0+3"
[[deps.Zlib_jll]]
deps = ["Libdl"]
......@@ -1413,9 +1454,9 @@ version = "1.2.12+3"
[[deps.Zstd_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl"]
git-tree-sha1 = "e678132f07ddb5bfa46857f0d7620fb9be675d3b"
git-tree-sha1 = "555d1076590a6cc2fdee2ef1469451f872d8b41b"
uuid = "3161d3a3-bdf6-5164-811a-617609db77b4"
version = "1.5.6+0"
version = "1.5.6+3"
[[deps.eudev_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg", "gperf_jll"]
......@@ -1431,9 +1472,9 @@ version = "0.53.0+0"
[[deps.gperf_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
git-tree-sha1 = "3516a5630f741c9eecb3720b1ec9d8edc3ecc033"
git-tree-sha1 = "0ba42241cb6809f1a278d0bcb976e0483c3f1f2d"
uuid = "1a1c6b14-54f6-533d-8383-74cd7377aa70"
version = "3.1.1+0"
version = "3.1.1+1"
[[deps.libaom_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl"]
......@@ -1524,6 +1565,6 @@ version = "3.5.0+0"
[[deps.xkbcommon_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg", "Wayland_jll", "Wayland_protocols_jll", "Xorg_libxcb_jll", "Xorg_xkeyboard_config_jll"]
git-tree-sha1 = "9c304562909ab2bab0262639bd4f444d7bc2be37"
git-tree-sha1 = "63406453ed9b33a0df95d570816d5366c92b7809"
uuid = "d8fb68d0-12a3-5cfd-a85a-d49703b185fd"
version = "1.4.1+1"
version = "1.4.1+2"
Project.toml
View file @ dc2fcd73
......@@ -3,6 +3,8 @@ ADerrors = "5e92007d-7bf1-471c-8ceb-4591b8b567a9"
ArgParse = "c7e460c6-2fb9-53a9-8c5b-16f535851c63"
BDIO = "375f315e-f2c4-11e9-2ef9-134f02f79e27"
CUDA = "052768ef-5323-5732-b1bb-66c8b64840ba"
Cthulhu = "f68482b8-f384-11e8-15f7-abe071a5a75f"
FormalSeries = "a55ae9b5-0907-4af6-81ab-c6ed2f8e0777"
LatticeGPU = "958c3683-801a-4582-9cfa-2d6e2ae1763b"
MD5 = "6ac74813-4b46-53a4-afec-0b5dc9d7885c"
TimerOutputs = "a759f4b9-e2f1-59dc-863e-4aeb61b1ea8f"
main.jl
View file @ dc2fcd73
......@@ -9,6 +9,7 @@ using CUDA
using InteractiveUtils
using BDIO
using MD5
using FormalSeries
reset_timer!()
......
nohup.out
View file @ dc2fcd73
WARNING: ignoring conflicting import of FormalSeries.norm into Spinors
Activating project at `~/Git/ferflow_quenched_AD_test`
ERROR: LoadError: MethodError: no method matching +(::Series{ComplexF64, 2}, ::Series{Float64, 2})
Closest candidates are:
+(::Any, ::Any, !Matched::Any, !Matched::Any...) at operators.jl:591
+(::Series{T, N}) where {T, N} at ~/.julia/packages/FormalSeries/VrbHe/src/FormalSeriesTypes.jl:106
+(::Series{T, N}, !Matched::Number) where {T, N} at ~/.julia/packages/FormalSeries/VrbHe/src/FormalSeriesTypes.jl:108
...
Stacktrace:
[1] macro expansion
@ ~/Git/ferflow_quenched_AD_test/src/meas.jl:196 [inlined]
[2] macro expansion
@ ~/.julia/packages/TimerOutputs/Lw5SP/src/TimerOutput.jl:237 [inlined]
[3] Backflow_pt()
@ Main ~/Git/ferflow_quenched_AD_test/src/meas.jl:193
[4] macro expansion
@ ~/.julia/packages/TimerOutputs/Lw5SP/src/TimerOutput.jl:237 [inlined]
[5] top-level scope
@ ~/Git/ferflow_quenched_AD_test/main.jl:38
in expression starting at /home/fperez/Git/ferflow_quenched_AD_test/main.jl:29
output/Template_run.log
View file @ dc2fcd73
This source diff could not be displayed because it is too large. You can view the blob instead.
src/io.jl
View file @ dc2fcd73
......@@ -66,13 +66,18 @@ function load_structs()
CUDA.device!(parsed_args["G"])
global lp = SpaceParm{4}(tuple(params["Space"]["size"]...), tuple(params["Space"]["blocks"]...),params["Space"]["bc"], (0,0,0,0,0,0))
global gp = GaugeParm{Float64}(SU3{Float64},params["HMC"]["beta"],params["HMC"]["c0"],(params["HMC"]["cG"],0.0),(0.0,0.0),lp.iL);
global dws = DiracWorkspace(SU3fund,Float64,lp);
global dws = DiracWorkspace(SU3fund,2,Float64,lp);
global ymws = YMworkspace(SU3,Float64,lp);
global int = wfl_rk3(Float64, params["Frontflow"]["epsilon"], params["Backflow"]["tol"])
global intsch = omf4(Float64,params["HMC"]["eps"], params["HMC"]["ns"]);
global fernames = [replace(k, "Fermion" => "") for k in keys(params) if startswith(k, "Fermion")]
global dpar = [DiracParam{Float64}(SU3fund,(1/(2*params["Fermion"*f]["kappa"])) - 4,params["Fermion"*f]["csw"],ntuple(i -> exp(((i!=4)*im*params["Fermion"*f]["theta"]/lp.iL[i]) + ((i==4)*im*params["Fermion"*f]["theta_t"]/lp.iL[i])),4),0.0,params["Fermion"*f]["ct"]) for f in fernames];
# SU3 assumed in the constructor
global dpar = [DiracParam{Float64}(
Series{Float64,2}(( (1/(2*params["Fermion"*f]["kappa"])) - 4, 1.0 )),
Series{Float64,2}(( params["Fermion"*f]["csw"], 0.0 )),
ntuple(i -> exp(((i!=4)*im*params["Fermion"*f]["theta"]/lp.iL[i]) + ((i==4)*im*params["Fermion"*f]["theta_t"]/lp.iL[i])), 4),
Series{Float64,2}(( 0.0, 0.0 )),
Series{Float64,2}(( params["Fermion"*f]["ct"], 0.0 )) ) for f in fernames];
global flow_times = params["Backflow"]["Flow_times"]
global MCid = 0
......@@ -89,7 +94,7 @@ end
function write_log()
println(log_file,"Running Ferflow.jl by ", params["Run"]["user"],". Name of the run: ",params["Run"]["name"])
println(log_file,"Running Ferflow.jl adfer v1.0 by ", params["Run"]["user"],". Name of the run: ",params["Run"]["name"])
println(log_file,"")
print(log_file,"Calling: ")
......@@ -123,6 +128,8 @@ function write_log()
println(log_file,"t_source backflow = ", params["Backflow"]["tsource"])
println(log_file,"Backflow tolerance = ", params["Backflow"]["tol"])
println(log_file,"Dirac Param struct")
println(log_file,dpar)
flush(log_file)
return nothing
end
......@@ -130,14 +137,14 @@ end
function save_data()
ihdr = [convert(Int32,1730280201)]
ihdr = [convert(Int32,1740655571)]
fname = "./output/"*params["Run"]["name"]*".bdio"
if isfile(fname)
fb = BDIO_open(fname, "a")
println(log_file,"\nAppending output to "*fname*"\n")
else
fb = BDIO_open(fname, "w", "BDIO output from ferflow.jl")
fb = BDIO_open(fname, "w", "BDIO output from ferflow.jl adfer v1.0")
println(log_file,"Creating new BDIO output file "*fname)
BDIO_start_record!(fb, BDIO_BIN_GENERIC, 14)
......@@ -148,9 +155,9 @@ function save_data()
BDIO_write!(fb, [convert(Int32, 4)])
BDIO_write!(fb, [convert(Int32, lp.iL[i]) for i in 1:4])
BDIO_write!(fb, [convert(Int32, length(fernames))])
for k in 1:length(fernames)
BDIO_write!(fb, [dpar[k].m0, dpar[k].csw])
BDIO_write!(fb, [dpar[k].th[i] for i in 1:4])
for f in fernames
BDIO_write!(fb, [params["Fermion"*f]["kappa"], params["Fermion"*f]["csw"]])
BDIO_write!(fb, [params["Fermion"*f]["theta"],])
end
BDIO_write!(fb, [convert(Int32,params["Frontflow"]["N_noise"])])
BDIO_write!(fb, [convert(Int32,params["Frontflow"]["nsteps"])])
......@@ -170,42 +177,65 @@ function save_data()
for f in 1:length(dpar)
for noi in 1:params["Frontflow"]["N_noise"]
BDIO_write!(fb,pp_corr_t0[:,noi,f])
BDIO_write!(fb,ap_corr_t0[:,noi,f])
BDIO_write!(fb,pphat_t0[:,noi,f])
BDIO_write!(fb,pptilde_t0[:,noi,f])
BDIO_write!(fb, getindex.(getfield.(pp_corr_t0[:,noi,f], :c), 1))
BDIO_write!(fb, getindex.(getfield.(ap_corr_t0[:,noi,f], :c), 1))
BDIO_write!(fb, getindex.(getfield.(pphat_t0[:,noi,f], :c), 1))
BDIO_write!(fb, getindex.(getfield.(pptilde_t0[:,noi,f], :c), 1))
for fl in 1:params["Frontflow"]["nsteps"]+1
BDIO_write!(fb,pp_corr_t[:,noi,fl,f])
BDIO_write!(fb,ap_corr_t[:,noi,fl,f])
BDIO_write!(fb,pphat_t[:,noi,fl,f])
BDIO_write!(fb,pptilde_t[:,noi,fl,f])
BDIO_write!(fb, getindex.(getfield.(pp_corr_t[:,noi,fl,f], :c), 1))
BDIO_write!(fb, getindex.(getfield.(ap_corr_t[:,noi,fl,f], :c), 1))
BDIO_write!(fb, getindex.(getfield.(pphat_t[:,noi,fl,f], :c), 1))
BDIO_write!(fb, getindex.(getfield.(pptilde_t[:,noi,fl,f], :c), 1))
end
end
end
for f in 1:length(dpar)
for noi in 1:params["Frontflow"]["N_noise"]
BDIO_write!(fb, getindex.(getfield.(pp_corr_t0[:,noi,f], :c), 2))
BDIO_write!(fb, getindex.(getfield.(ap_corr_t0[:,noi,f], :c), 2))
BDIO_write!(fb, getindex.(getfield.(pphat_t0[:,noi,f], :c), 2))
BDIO_write!(fb, getindex.(getfield.(pptilde_t0[:,noi,f], :c), 2))
for fl in 1:params["Frontflow"]["nsteps"]+1
BDIO_write!(fb, getindex.(getfield.(pp_corr_t[:,noi,fl,f], :c), 2))
BDIO_write!(fb, getindex.(getfield.(ap_corr_t[:,noi,fl,f], :c), 2))
BDIO_write!(fb, getindex.(getfield.(pphat_t[:,noi,fl,f], :c), 2))
BDIO_write!(fb, getindex.(getfield.(pptilde_t[:,noi,fl,f], :c), 2))
end
end
end
BDIO_write!(fb,Eoft0)
BDIO_write!(fb,Qt0)
BDIO_write!(fb, Eoft0)
BDIO_write!(fb, Qt0)
for fl in 1:params["Frontflow"]["nsteps"]+1
BDIO_write!(fb,Eoft[fl,:])
BDIO_write!(fb,Qt[fl,:])
BDIO_write!(fb, Eoft[fl,:])
BDIO_write!(fb, Qt[fl,:])
end
for f in 1:length(dpar)
for noi in 1:params["Backflow"]["N_noise"]
for fl in 1:length(flow_times)
BDIO_write!(fb,pp_corr_tfl[:,noi,fl,f])
BDIO_write!(fb,ap_corr_tfl[:,noi,fl,f])
BDIO_write!(fb,Quark_cond[:,noi,fl,f])
BDIO_write!(fb,Quark_cond_cfl[:,noi,fl,f])
BDIO_write!(fb,Quark_cond2[:,noi,fl,f])
BDIO_write!(fb,Quark_cond2_cfl[:,noi,fl,f])
BDIO_write!(fb,ChiDchi[:,noi,fl,f])
BDIO_write!(fb,ChiDchi_cfl[:,noi,fl,f])
BDIO_write!(fb, getindex.(getfield.(pp_corr_tfl[:,noi,fl,f], :c), 1))
BDIO_write!(fb, getindex.(getfield.(ap_corr_tfl[:,noi,fl,f], :c), 1))
BDIO_write!(fb, getindex.(getfield.(Quark_cond[:,noi,fl,f], :c), 1))
BDIO_write!(fb, getindex.(getfield.(Quark_cond_cfl[:,noi,fl,f], :c), 1))
BDIO_write!(fb, getindex.(getfield.(Quark_cond2[:,noi,fl,f], :c), 1))
BDIO_write!(fb, getindex.(getfield.(Quark_cond2_cfl[:,noi,fl,f], :c), 1))
BDIO_write!(fb, getindex.(getfield.(ChiDchi[:,noi,fl,f], :c), 1))
BDIO_write!(fb, getindex.(getfield.(ChiDchi_cfl[:,noi,fl,f], :c), 1))
end
end
end
for f in 1:length(dpar)
for noi in 1:params["Backflow"]["N_noise"]
for fl in 1:length(flow_times)
BDIO_write!(fb, getindex.(getfield.(pp_corr_tfl[:,noi,fl,f], :c), 2))
BDIO_write!(fb, getindex.(getfield.(ap_corr_tfl[:,noi,fl,f], :c), 2))
BDIO_write!(fb, getindex.(getfield.(Quark_cond[:,noi,fl,f], :c), 2))
BDIO_write!(fb, getindex.(getfield.(Quark_cond_cfl[:,noi,fl,f], :c), 2))
BDIO_write!(fb, getindex.(getfield.(Quark_cond2[:,noi,fl,f], :c), 2))
BDIO_write!(fb, getindex.(getfield.(Quark_cond2_cfl[:,noi,fl,f], :c), 2))
BDIO_write!(fb, getindex.(getfield.(ChiDchi[:,noi,fl,f], :c), 2))
BDIO_write!(fb, getindex.(getfield.(ChiDchi_cfl[:,noi,fl,f], :c), 2))
end
end
end
......
src/meas.jl
View file @ dc2fcd73
......@@ -9,34 +9,38 @@ function load_fields()
fill!(U,one(SU3{Float64}));
global U_CPU = Array(U)
global psi = scalar_field(Spinor{4,SU3fund{Float64}},lp)
global psi = scalar_field(Spinor{4,SU3fund{Float64,2}},lp)
fill!(psi,zero(Spinor{4, SU3fund{Float64, 2}}))
global psi_CPU = Array(psi)
global psit_CPU = Array(psi)
global pp_density = Array{Float64}(undef,lp.bsz,lp.rsz);
global ap_density = Array{Complex{Float64}}(undef,lp.bsz,lp.rsz);
global pp_density = Array{Series{Float64,2}}(undef,lp.bsz,lp.rsz);
global ap_density = Array{Series{Complex{Float64},2}}(undef,lp.bsz,lp.rsz);
global complex_scalar = scalar_field_point(Series{ComplexF64, 2}, lp)
global real_scalar = scalar_field_point(Series{Float64, 2}, lp)
# Obs ending in 0 are only relevant if tzero > 0
global pp_corr_t0 = fill(zero(Float64),(lp.iL[4],params["Frontflow"]["N_noise"],length(dpar)));
global ap_corr_t0 = fill(zero(ComplexF64),(lp.iL[4],params["Frontflow"]["N_noise"],length(dpar)));
global pp_corr_t0 = fill(zero(Series{Float64,2}),(lp.iL[4],params["Frontflow"]["N_noise"],length(dpar)));
global ap_corr_t0 = fill(zero(Series{Complex{Float64},2}),(lp.iL[4],params["Frontflow"]["N_noise"],length(dpar)));
global pphat_t0 = Array{Float64}(undef,lp.iL[4],params["Frontflow"]["N_noise"],length(dpar));
global pptilde_t0 = Array{Complex{Float64}}(undef,lp.iL[4],params["Frontflow"]["N_noise"],length(dpar));
global pphat_t0 = Array{Series{Float64,2}}(undef,lp.iL[4],params["Frontflow"]["N_noise"],length(dpar));
global pptilde_t0 = Array{Series{Complex{Float64},2}}(undef,lp.iL[4],params["Frontflow"]["N_noise"],length(dpar));
global Qt0 = Array{Complex{Float64}}(undef,lp.iL[4]);
global Eoft0 = Array{Complex{Float64}}(undef,lp.iL[4]);
# Frontflow
global pp_corr_t = fill(zero(Float64),(lp.iL[4],params["Frontflow"]["N_noise"],params["Frontflow"]["nsteps"] + 1,length(dpar)));
global ap_corr_t = fill(zero(ComplexF64),(lp.iL[4],params["Frontflow"]["N_noise"],params["Frontflow"]["nsteps"] + 1,length(dpar)));
global pp_corr_t = fill(zero(Series{Float64,2}),(lp.iL[4],params["Frontflow"]["N_noise"],params["Frontflow"]["nsteps"] + 1,length(dpar)));
global ap_corr_t = fill(zero(Series{Complex{Float64},2}),(lp.iL[4],params["Frontflow"]["N_noise"],params["Frontflow"]["nsteps"] + 1,length(dpar)));
global pp_corr_tfl = fill(zero(Float64),(lp.iL[4],params["Backflow"]["N_noise"],length(flow_times),length(dpar)));
global ap_corr_tfl = fill(zero(ComplexF64),(lp.iL[4],params["Backflow"]["N_noise"],length(flow_times),length(dpar)));
global pp_corr_tfl = fill(zero(Series{Float64,2}),(lp.iL[4],params["Backflow"]["N_noise"],length(flow_times),length(dpar)));
global ap_corr_tfl = fill(zero(Series{Complex{Float64},2}),(lp.iL[4],params["Backflow"]["N_noise"],length(flow_times),length(dpar)));
global pphat_t = Array{Float64}(undef,lp.iL[4],params["Frontflow"]["N_noise"],params["Frontflow"]["nsteps"] + 1,length(dpar));
global pptilde_t = Array{Complex{Float64}}(undef,lp.iL[4],params["Frontflow"]["N_noise"],params["Frontflow"]["nsteps"] + 1,length(dpar));
global pphat_t = Array{Series{Float64,2}}(undef,lp.iL[4],params["Frontflow"]["N_noise"],params["Frontflow"]["nsteps"] + 1,length(dpar));
global pptilde_t = Array{Series{Complex{Float64},2}}(undef,lp.iL[4],params["Frontflow"]["N_noise"],params["Frontflow"]["nsteps"] + 1,length(dpar));
global Qteuc = Array{Float64}(undef,lp.iL[4]);
global Qt = Array{Complex{Float64}}(undef,params["Frontflow"]["nsteps"] + 1,lp.iL[4]);
......@@ -45,13 +49,13 @@ function load_fields()
# Backflow
global Quark_cond = Array{Complex{Float64}}(undef,lp.iL[4],params["Backflow"]["N_noise"],length(flow_times),length(dpar));
global Quark_cond_cfl = Array{Complex{Float64}}(undef,lp.iL[4],params["Backflow"]["N_noise"],length(flow_times),length(dpar));
global Quark_cond = Array{Series{Complex{Float64},2}}(undef,lp.iL[4],params["Backflow"]["N_noise"],length(flow_times),length(dpar));
global Quark_cond_cfl = Array{Series{Complex{Float64},2}}(undef,lp.iL[4],params["Backflow"]["N_noise"],length(flow_times),length(dpar));
global Quark_cond2 = Array{Complex{Float64}}(undef,lp.iL[4],params["Backflow"]["N_noise"],length(flow_times),length(dpar));
global Quark_cond2_cfl = Array{Complex{Float64}}(undef,lp.iL[4],params["Backflow"]["N_noise"],length(flow_times),length(dpar));
global ChiDchi = Array{Complex{Float64}}(undef,lp.iL[4],params["Backflow"]["N_noise"],length(flow_times),length(dpar));
global ChiDchi_cfl = Array{Complex{Float64}}(undef,lp.iL[4],params["Backflow"]["N_noise"],length(flow_times),length(dpar));
global Quark_cond2 = Array{Series{Complex{Float64},2}}(undef,lp.iL[4],params["Backflow"]["N_noise"],length(flow_times),length(dpar));
global Quark_cond2_cfl = Array{Series{Complex{Float64},2}}(undef,lp.iL[4],params["Backflow"]["N_noise"],length(flow_times),length(dpar));
global ChiDchi = Array{Series{Complex{Float64},2}}(undef,lp.iL[4],params["Backflow"]["N_noise"],length(flow_times),length(dpar));
global ChiDchi_cfl = Array{Series{Complex{Float64},2}}(undef,lp.iL[4],params["Backflow"]["N_noise"],length(flow_times),length(dpar));
return nothing
end
......@@ -81,8 +85,8 @@ function Frontflow_pt() # Will be Frontflow
println(log_file,"CG converged after ",niter," iterations with residue ",CUDA.mapreduce(x -> norm2(x), +, dws.sr))
flush(log_file)
pp_corr_t0[:,noi,f] .= zero(Float64)
ap_corr_t0[:,noi,f] .= zero(ComplexF64)
pp_corr_t0[:,noi,f] .= zero(Series{Float64,2})
ap_corr_t0[:,noi,f] .= zero(Series{ComplexF64,2})
@timeit "Volume sum" begin
pp_corr_t0[:,noi,f] .= scalar_contraction(psi)
ap_corr_t0[:,noi,f] .= gammazero_contraction(psi)
......@@ -103,8 +107,8 @@ function Frontflow_pt() # Will be Frontflow
Qt[tstep,:] .= Qteuc
end
pp_corr_t[:,noi,tstep,f] .= zero(Float64)
ap_corr_t[:,noi,tstep,f] .= zero(ComplexF64)
pp_corr_t[:,noi,tstep,f] .= zero(Series{Float64,2})
ap_corr_t[:,noi,tstep,f] .= zero(Series{ComplexF64,2})
@timeit "Volume sum" begin
pp_corr_t[:,noi,tstep,f] .= scalar_contraction(psi)
ap_corr_t[:,noi,tstep,f] .= gammazero_contraction(psi)
......@@ -120,8 +124,8 @@ function Frontflow_pt() # Will be Frontflow
Qt[end,:] .= Qteuc
end
pp_corr_t[:,noi,end,f] .= zero(Float64)
ap_corr_t[:,noi,end,f] .= zero(ComplexF64)
pp_corr_t[:,noi,end,f] .= zero(Series{Float64,2})
ap_corr_t[:,noi,end,f] .= zero(Series{ComplexF64,2})
@timeit "Volume sum" begin
pp_corr_t[:,noi,end,f] .= scalar_contraction(psi)
ap_corr_t[:,noi,end,f] .= gammazero_contraction(psi)
......@@ -182,25 +186,23 @@ function Backflow_pt() # Will be backflow
@timeit "CPU to GPU" copyto!(dws.st,psi_CPU)
Quark_cond[:,noi,fl,f] .= zero(ComplexF64)
ap_density .= Array(dot.(dws.st,psi))
Quark_cond_cfl[:,noi,fl,f] .= zero(ComplexF64)
pp_density .= Array(norm2.(dws.st))
Quark_cond_cfl[:,noi,fl,f] .= zero(Series{ComplexF64,2})
ap_density .= Array(norm2.(dws.st))
@timeit "Volume sum" for b in 1:lp.bsz for r in 1:lp.rsz
t = point_time((b,r),lp)
Quark_cond[t,noi,fl,f] += ap_density[b,r]
Quark_cond_cfl[t,noi,fl,f] += -pp_density[b,r]
Quark_cond_cfl[t,noi,fl,f] += -ap_density[b,r]
end end
Quark_cond[:,noi,fl,f] .= zero(Series{ComplexF64,2})
@timeit "Volume sum" begin
Quark_cond[:,noi,fl,f] .= dot_contraction(dws.st,psi)
end
pp_density .= Array(norm2.(psi))
ap_density .= Array(dot.(psi,dmul.(Gamma{4},psi)))
pp_corr_tfl[:,noi,fl,f] .= zero(Float64)
ap_corr_tfl[:,noi,fl,f] .= zero(ComplexF64)
@timeit "Volume sum" for b in 1:lp.bsz for r in 1:lp.rsz
t = point_time((b,r),lp)
pp_corr_tfl[t,noi,fl,f] += pp_density[b,r]
ap_corr_tfl[t,noi,fl,f] += ap_density[b,r]
end end
pp_corr_tfl[:,noi,fl,f] .= zero(Series{Float64,2})
ap_corr_tfl[:,noi,fl,f] .= zero(Series{ComplexF64,2})
@timeit "Volume sum" begin
pp_corr_tfl[:,noi,fl,f] .= scalar_contraction(psi)
ap_corr_tfl[:,noi,fl,f] .= gammazero_contraction(psi)
end
@timeit "CPU to GPU" copyto!(U,U_CPU)
......@@ -211,7 +213,7 @@ function Backflow_pt() # Will be backflow
@timeit "CPU to GPU" copyto!(dws.st,psit_CPU)
ChiDchi[:,noi,fl,f] .= zero(Float64)
ChiDchi[:,noi,fl,f] .= zero(Series{Complex{Float64},2})
ap_density .= Array(dot.(dws.st,dws.sp))
@timeit "Volume sum" for b in 1:lp.bsz for r in 1:lp.rsz
t = point_time((b,r),lp)
......@@ -224,7 +226,7 @@ function Backflow_pt() # Will be backflow
ChiDchi[t,noi,fl,f] += -ap_density[b,r]
end end
Quark_cond2[:,noi,fl,f] .= zero(ComplexF64)
Quark_cond2[:,noi,fl,f] .= zero(Series{ComplexF64,2})
ap_density .= Array(dot.(dws.st,psi))
@timeit "Volume sum" for b in 1:lp.bsz for r in 1:lp.rsz
t = point_time((b,r),lp)
......@@ -237,7 +239,7 @@ function Backflow_pt() # Will be backflow
flw_adapt(U, psi, int, flow_times[fl], gp, dpar[f], lp, ymws, dws)
Dw!(dws.sp,U,psi,dpar[f],dws,lp)
ChiDchi_cfl[:,noi,fl,f] .= zero(ComplexF64)
ChiDchi_cfl[:,noi,fl,f] .= zero(Series{ComplexF64,2})
ap_density .= Array(dot.(dws.st,dws.sp))
@timeit "Volume sum" for b in 1:lp.bsz for r in 1:lp.rsz
t = point_time((b,r),lp)
......@@ -250,7 +252,7 @@ function Backflow_pt() # Will be backflow
ChiDchi_cfl[t,noi,fl,f] += -ap_density[b,r]
end end
Quark_cond2_cfl[:,noi,fl,f] .= zero(ComplexF64)
Quark_cond2_cfl[:,noi,fl,f] .= zero(Series{ComplexF64,2})
ap_density .= Array(dot.(dws.st,psi))
@timeit "Volume sum" for b in 1:lp.bsz for r in 1:lp.rsz
t = point_time((b,r),lp)
......@@ -312,8 +314,8 @@ function Two_pt_lagrange() # Will be 2pt lagrange mult
println(log_file,"CG converged after ",niter," iterations with residue ",CUDA.mapreduce(x -> norm2(x), +, dws.sr))
flush(log_file)
pphat_t0[:,noi,f] .= zero(Float64)
pptilde_t0[:,noi,f] .= zero(ComplexF64)
pphat_t0[:,noi,f] .= zero(Series{Float64,2})
pptilde_t0[:,noi,f] .= zero(Series{ComplexF64,2})
@timeit "Volume sum" begin
pphat_t0[:,noi,f] .= -scalar_contraction(psi)
pptilde_t0[:,noi,f] .= dot_contraction(dws.st,psi)
......@@ -333,8 +335,8 @@ function Two_pt_lagrange() # Will be 2pt lagrange mult
for fl in 1:params["Frontflow"]["nsteps"]
pphat_t[:,noi,fl,f] .= zero(Float64)
pptilde_t[:,noi,fl,f] .= zero(ComplexF64)
pphat_t[:,noi,fl,f] .= zero(Series{Float64,2})
pptilde_t[:,noi,fl,f] .= zero(Series{ComplexF64,2})
@timeit "Volume sum" begin
pphat_t[:,noi,fl,f] .= -scalar_contraction(psi)
pptilde_t[:,noi,fl,f] .= dot_contraction(dws.st,psi)
......@@ -345,8 +347,8 @@ function Two_pt_lagrange() # Will be 2pt lagrange mult
flw(U, dws.st, int, 1, params["Frontflow"]["epsilon"], gp, dpar[f], lp, ymws, dws)
end
pphat_t[:,noi,end,f] .= zero(Float64)
pptilde_t[:,noi,end,f] .= zero(ComplexF64)
pphat_t[:,noi,end,f] .= zero(Series{Float64,2})
pptilde_t[:,noi,end,f] .= zero(Series{ComplexF64,2})
@timeit "Volume sum" begin
pphat_t[:,noi,end,f] .= -scalar_contraction(psi)
pptilde_t[:,noi,end,f] .= dot_contraction(dws.st,psi)
......@@ -376,15 +378,15 @@ function krnl_scalar_contraction!(rm, psi, lp)
b = Int64(CUDA.threadIdx().x)
r = Int64(CUDA.blockIdx().x)
I = point_coord((b,r), lp)
rm[I] = norm2(psi[b,r])
rm[I] = ctoreal(norm2(psi[b,r]))
end
return nothing
end
function scalar_contraction(psi)
CUDA.@sync begin
CUDA.@cuda threads=lp.bsz blocks=lp.rsz krnl_scalar_contraction!(ymws.rm, psi, lp)
CUDA.@cuda threads=lp.bsz blocks=lp.rsz krnl_scalar_contraction!(real_scalar, psi, lp)
end
return reshape(Array(CUDA.reduce(+, ymws.rm; dims=(1,2,3))),lp.iL[end])
return reshape(Array(CUDA.reduce(+, real_scalar; dims=(1,2,3))),lp.iL[end])
end
function krnl_dot_contraction!(cm,st, psi, lp)
......@@ -398,9 +400,9 @@ function krnl_dot_contraction!(cm,st, psi, lp)
end
function dot_contraction(st,psi)
CUDA.@sync begin
CUDA.@cuda threads=lp.bsz blocks=lp.rsz krnl_dot_contraction!(ymws.cm,st, psi, lp)
CUDA.@cuda threads=lp.bsz blocks=lp.rsz krnl_dot_contraction!(complex_scalar,st, psi, lp)
end
return reshape(Array(CUDA.reduce(+, ymws.cm; dims=(1,2,3))),lp.iL[end])
return reshape(Array(CUDA.reduce(+, complex_scalar; dims=(1,2,3))),lp.iL[end])
end
function krnl_gammazero_contraction!(cm, psi, lp)
......@@ -414,7 +416,7 @@ function krnl_gammazero_contraction!(cm, psi, lp)
end
function gammazero_contraction(psi)
CUDA.@sync begin
CUDA.@cuda threads=lp.bsz blocks=lp.rsz krnl_gammazero_contraction!(ymws.cm, psi, lp)
CUDA.@cuda threads=lp.bsz blocks=lp.rsz krnl_gammazero_contraction!(complex_scalar, psi, lp)
end
return reshape(Array(CUDA.reduce(+, ymws.cm; dims=(1,2,3))),lp.iL[end])
return reshape(Array(CUDA.reduce(+, complex_scalar; dims=(1,2,3))),lp.iL[end])
end
src/utils.jl
View file @ dc2fcd73
......@@ -12,7 +12,7 @@ function read_ff(name::String)
end
ihdr = Vector{Int32}(undef,1)
BDIO_read(file, ihdr)
if ihdr[1] == 1730280201
if ihdr[1] == 1740655571
println("Reading file ",name)
else
error("Wrong IHDR in file ", name)
......@@ -40,7 +40,7 @@ function read_ff(name::String)
BDIO_read(file, Nf)
dpars=[Vector{Float64}(undef,2) for _ in 1:Nf[1]] # m0,csw
th=[Vector{ComplexF64}(undef,4) for _ in 1:Nf[1]]
th=[Vector{Float64}(undef,1) for _ in 1:Nf[1]]
for i in 1:Nf[1]
BDIO_read(file, dpars[i])
......@@ -95,6 +95,55 @@ function read_ff(name::String)
ChiDchi = Array{ComplexF64}(undef,iL[4],N_noisebf,nflow,Nf);
ChiDchi_cfl = Array{ComplexF64}(undef,iL[4],N_noisebf,nflow,Nf);
#
## d/dm
#
pp_corr_dm = Array{Float64}(undef,iL[4],N_noiseff,Nf);
ap_corr_dm = Array{ComplexF64}(undef,iL[4],N_noiseff,Nf);
pphat0_dm = Array{Float64}(undef,iL[4],N_noiseff,Nf);
pptilde0_dm = Array{ComplexF64}(undef,iL[4],N_noiseff,Nf);
pp_corr_t_dm = Array{Float64}(undef,iL[4],N_noiseff,nsteps+1,Nf);
ap_corr_t_dm = Array{ComplexF64}(undef,iL[4],N_noiseff,nsteps+1,Nf);
pphat_t_dm = Array{Float64}(undef,iL[4],N_noiseff,nsteps+1,Nf);
pptilde_t_dm = Array{ComplexF64}(undef,iL[4],N_noiseff,nsteps+1,Nf);
Eoft_dm = Array{Complex{Float64}}(undef,1+nsteps,iL[4]);
Eoft0_dm = Array{Complex{Float64}}(undef,iL[4]);
Qt_dm = Array{Complex{Float64}}(undef,1+nsteps,iL[4]);
Qt0_dm = Array{Complex{Float64}}(undef,iL[4]);
pp_corr_tfl_dm = Array{Float64}(undef,iL[4],N_noisebf,nflow,Nf);
ap_corr_tfl_dm = Array{ComplexF64}(undef,iL[4],N_noisebf,nflow,Nf);
Sigma_dm = Array{ComplexF64}(undef,iL[4],N_noisebf,nflow,Nf);
Sigma_cfl_dm = Array{ComplexF64}(undef,iL[4],N_noisebf,nflow,Nf);
Sigma2_dm = Array{ComplexF64}(undef,iL[4],N_noisebf,nflow,Nf);
Sigma2_cfl_dm = Array{ComplexF64}(undef,iL[4],N_noisebf,nflow,Nf);
ChiDchi_dm = Array{ComplexF64}(undef,iL[4],N_noisebf,nflow,Nf);
ChiDchi_cfl_dm = Array{ComplexF64}(undef,iL[4],N_noisebf,nflow,Nf);
pp_corr_mc = [Vector{Float64}() for _ in 1:iL[4], _ in 1:N_noiseff, _ in 1:Nf]
ap_corr_mc = [Vector{ComplexF64}() for _ in 1:iL[4], _ in 1:N_noiseff, _ in 1:Nf]
pphat0_mc = [Vector{Float64}() for _ in 1:iL[4], _ in 1:N_noiseff, _ in 1:Nf]
pptilde0_mc = [Vector{ComplexF64}() for _ in 1:iL[4], _ in 1:N_noiseff, _ in 1:Nf]
pp_corr_t_mc = [Vector{Float64}() for _ in 1:iL[4], _ in 1:N_noiseff, _ in 1:nsteps+1, _ in 1:Nf]
ap_corr_t_mc = [Vector{ComplexF64}() for _ in 1:iL[4], _ in 1:N_noiseff, _ in 1:nsteps+1, _ in 1:Nf]
pphat_t_mc = [Vector{Float64}() for _ in 1:iL[4], _ in 1:N_noiseff, _ in 1:nsteps+1, _ in 1:Nf]
pptilde_t_mc = [Vector{ComplexF64}() for _ in 1:iL[4], _ in 1:N_noiseff, _ in 1:nsteps+1, _ in 1:Nf]
Eoft_mc = [Vector{Complex{Float64}}() for _ in 1:1+nsteps, _ in 1:iL[4]]
Eoft0_mc = [Vector{Complex{Float64}}() for _ in 1:iL[4]]
Qt_mc = [Vector{Complex{Float64}}() for _ in 1:1+nsteps, _ in 1:iL[4]]
Qt0_mc = [Vector{Complex{Float64}}() for _ in 1:iL[4]]
pp_corr_tfl_mc = [Vector{Float64}() for _ in 1:iL[4], _ in 1:N_noisebf, _ in 1:nflow, _ in 1:Nf]
ap_corr_tfl_mc = [Vector{ComplexF64}() for _ in 1:iL[4], _ in 1:N_noisebf, _ in 1:nflow, _ in 1:Nf]
Sigma_mc = [Vector{ComplexF64}() for _ in 1:iL[4], _ in 1:N_noisebf, _ in 1:nflow, _ in 1:Nf]
Sigma_cfl_mc = [Vector{ComplexF64}() for _ in 1:iL[4], _ in 1:N_noisebf, _ in 1:nflow, _ in 1:Nf]
Sigma2_mc = [Vector{ComplexF64}() for _ in 1:iL[4], _ in 1:N_noisebf, _ in 1:nflow, _ in 1:Nf]
Sigma2_cfl_mc = [Vector{ComplexF64}() for _ in 1:iL[4], _ in 1:N_noisebf, _ in 1:nflow, _ in 1:Nf]
ChiDchi_mc = [Vector{ComplexF64}() for _ in 1:iL[4], _ in 1:N_noisebf, _ in 1:nflow, _ in 1:Nf]
ChiDchi_cfl_mc = [Vector{ComplexF64}() for _ in 1:iL[4], _ in 1:N_noisebf, _ in 1:nflow, _ in 1:Nf]
pp_corr_mc = [Vector{Float64}() for _ in 1:iL[4], _ in 1:N_noiseff, _ in 1:Nf]
ap_corr_mc = [Vector{ComplexF64}() for _ in 1:iL[4], _ in 1:N_noiseff, _ in 1:Nf]
pphat0_mc = [Vector{Float64}() for _ in 1:iL[4], _ in 1:N_noiseff, _ in 1:Nf]
......@@ -117,6 +166,45 @@ function read_ff(name::String)
ChiDchi_mc = [Vector{ComplexF64}() for _ in 1:iL[4], _ in 1:N_noisebf, _ in 1:nflow, _ in 1:Nf]
ChiDchi_cfl_mc = [Vector{ComplexF64}() for _ in 1:iL[4], _ in 1:N_noisebf, _ in 1:nflow, _ in 1:Nf]
pp_corr_mc_dm = [Vector{Float64}() for _ in 1:iL[4], _ in 1:N_noiseff, _ in 1:Nf]
ap_corr_mc_dm = [Vector{ComplexF64}() for _ in 1:iL[4], _ in 1:N_noiseff, _ in 1:Nf]
pphat0_mc_dm = [Vector{Float64}() for _ in 1:iL[4], _ in 1:N_noiseff, _ in 1:Nf]
pptilde0_mc_dm = [Vector{ComplexF64}() for _ in 1:iL[4], _ in 1:N_noiseff, _ in 1:Nf]
pp_corr_t_mc_dm = [Vector{Float64}() for _ in 1:iL[4], _ in 1:N_noiseff, _ in 1:nsteps+1, _ in 1:Nf]
ap_corr_t_mc_dm = [Vector{ComplexF64}() for _ in 1:iL[4], _ in 1:N_noiseff, _ in 1:nsteps+1, _ in 1:Nf]
pphat_t_mc_dm = [Vector{Float64}() for _ in 1:iL[4], _ in 1:N_noiseff, _ in 1:nsteps+1, _ in 1:Nf]
pptilde_t_mc_dm = [Vector{ComplexF64}() for _ in 1:iL[4], _ in 1:N_noiseff, _ in 1:nsteps+1, _ in 1:Nf]
pp_corr_tfl_mc_dm = [Vector{Float64}() for _ in 1:iL[4], _ in 1:N_noisebf, _ in 1:nflow, _ in 1:Nf]
ap_corr_tfl_mc_dm = [Vector{ComplexF64}() for _ in 1:iL[4], _ in 1:N_noisebf, _ in 1:nflow, _ in 1:Nf]
Sigma_mc_dm = [Vector{ComplexF64}() for _ in 1:iL[4], _ in 1:N_noisebf, _ in 1:nflow, _ in 1:Nf]
Sigma_cfl_mc_dm = [Vector{ComplexF64}() for _ in 1:iL[4], _ in 1:N_noisebf, _ in 1:nflow, _ in 1:Nf]
Sigma2_mc_dm = [Vector{ComplexF64}() for _ in 1:iL[4], _ in 1:N_noisebf, _ in 1:nflow, _ in 1:Nf]
Sigma2_cfl_mc_dm = [Vector{ComplexF64}() for _ in 1:iL[4], _ in 1:N_noisebf, _ in 1:nflow, _ in 1:Nf]
ChiDchi_mc_dm = [Vector{ComplexF64}() for _ in 1:iL[4], _ in 1:N_noisebf, _ in 1:nflow, _ in 1:Nf]
ChiDchi_cfl_mc_dm = [Vector{ComplexF64}() for _ in 1:iL[4], _ in 1:N_noisebf, _ in 1:nflow, _ in 1:Nf]
pp_corr_mc_dm = [Vector{Float64}() for _ in 1:iL[4], _ in 1:N_noiseff, _ in 1:Nf]
ap_corr_mc_dm = [Vector{ComplexF64}() for _ in 1:iL[4], _ in 1:N_noiseff, _ in 1:Nf]
pphat0_mc_dm = [Vector{Float64}() for _ in 1:iL[4], _ in 1:N_noiseff, _ in 1:Nf]
pptilde0_mc_dm = [Vector{ComplexF64}() for _ in 1:iL[4], _ in 1:N_noiseff, _ in 1:Nf]
pp_corr_t_mc_dm = [Vector{Float64}() for _ in 1:iL[4], _ in 1:N_noiseff, _ in 1:nsteps+1, _ in 1:Nf]
ap_corr_t_mc_dm = [Vector{ComplexF64}() for _ in 1:iL[4], _ in 1:N_noiseff, _ in 1:nsteps+1, _ in 1:Nf]
pphat_t_mc_dm = [Vector{Float64}() for _ in 1:iL[4], _ in 1:N_noiseff, _ in 1:nsteps+1, _ in 1:Nf]
pptilde_t_mc_dm = [Vector{ComplexF64}() for _ in 1:iL[4], _ in 1:N_noiseff, _ in 1:nsteps+1, _ in 1:Nf]
pp_corr_tfl_mc_dm = [Vector{Float64}() for _ in 1:iL[4], _ in 1:N_noisebf, _ in 1:nflow, _ in 1:Nf]
ap_corr_tfl_mc_dm = [Vector{ComplexF64}() for _ in 1:iL[4], _ in 1:N_noisebf, _ in 1:nflow, _ in 1:Nf]
Sigma_mc_dm = [Vector{ComplexF64}() for _ in 1:iL[4], _ in 1:N_noisebf, _ in 1:nflow, _ in 1:Nf]
Sigma_cfl_mc_dm = [Vector{ComplexF64}() for _ in 1:iL[4], _ in 1:N_noisebf, _ in 1:nflow, _ in 1:Nf]
Sigma2_mc_dm = [Vector{ComplexF64}() for _ in 1:iL[4], _ in 1:N_noisebf, _ in 1:nflow, _ in 1:Nf]
Sigma2_cfl_mc_dm = [Vector{ComplexF64}() for _ in 1:iL[4], _ in 1:N_noisebf, _ in 1:nflow, _ in 1:Nf]
ChiDchi_mc_dm = [Vector{ComplexF64}() for _ in 1:iL[4], _ in 1:N_noisebf, _ in 1:nflow, _ in 1:Nf]
ChiDchi_cfl_mc_dm = [Vector{ComplexF64}() for _ in 1:iL[4], _ in 1:N_noisebf, _ in 1:nflow, _ in 1:Nf]
while BDIO_get_uinfo(file) != 8
BDIO_seek!(file)
end
......@@ -152,6 +240,35 @@ function read_ff(name::String)
end
end
end
for f in 1:Nf
for noi in 1:N_noiseff
BDIO_read(file,TvecR)
pp_corr_dm[:,noi,f] .= TvecR
BDIO_read(file,TvecC)
ap_corr_dm[:,noi,f] .= TvecC
BDIO_read(file,TvecR)
pphat0_dm[:,noi,f] .= TvecR
BDIO_read(file,TvecC)
pptilde0_dm[:,noi,f] .= TvecC
for fl in 1:nsteps+1
BDIO_read(file,TvecR)
pp_corr_t_dm[:,noi,fl,f].= TvecR
BDIO_read(file,TvecC)
ap_corr_t_dm[:,noi,fl,f].= TvecC
BDIO_read(file,TvecR)
pphat_t_dm[:,noi,fl,f].= TvecR
BDIO_read(file,TvecC)
pptilde_t_dm[:,noi,fl,f].= TvecC
end
end
end
BDIO_read(file,TvecC)
Eoft0[:] .= TvecC
......@@ -194,6 +311,35 @@ function read_ff(name::String)
end
end
end
for f in 1:Nf
for noi in 1:N_noisebf
for fl in 1:nflow
BDIO_read(file,TvecR)
pp_corr_tfl_dm[:,noi,fl,f].= TvecR
BDIO_read(file,TvecC)
ap_corr_tfl_dm[:,noi,fl,f].= TvecC
BDIO_read(file,TvecC)
Sigma_dm[:,noi,fl,f] .= TvecC
BDIO_read(file,TvecC)
Sigma_cfl_dm[:,noi,fl,f] .= TvecC
BDIO_read(file,TvecC)
Sigma2_dm[:,noi,fl,f] .= TvecC
BDIO_read(file,TvecC)
Sigma2_cfl_dm[:,noi,fl,f] .= TvecC
BDIO_read(file,TvecC)
ChiDchi_dm[:,noi,fl,f] .= TvecC
BDIO_read(file,TvecC)
ChiDchi_cfl_dm[:,noi,fl,f] .= TvecC
end
end
end
push!.(Eoft0_mc, Eoft0)
push!.(Eoft_mc, Eoft)
......@@ -218,13 +364,36 @@ function read_ff(name::String)
push!.(ChiDchi_mc, ChiDchi)
push!.(ChiDchi_cfl_mc, ChiDchi_cfl)
push!.(pp_corr_mc_dm, pp_corr_dm)
push!.(ap_corr_mc_dm, ap_corr_dm)
push!.(pphat0_mc_dm, pphat0_dm)
push!.(pptilde0_mc_dm, pptilde0_dm)
push!.(pp_corr_t_mc_dm, pp_corr_t_dm)
push!.(ap_corr_t_mc_dm, ap_corr_t_dm)
push!.(pphat_t_mc_dm, pphat_t_dm)
push!.(pptilde_t_mc_dm, pptilde_t_dm)
push!.(pp_corr_tfl_mc_dm, pp_corr_tfl_dm)
push!.(ap_corr_tfl_mc_dm, ap_corr_tfl_dm)
push!.(Sigma_mc_dm, Sigma_dm)
push!.(Sigma_cfl_mc_dm, Sigma_cfl_dm)
push!.(Sigma2_mc_dm, Sigma2_dm)
push!.(Sigma2_cfl_mc_dm, Sigma2_cfl_dm)
push!.(ChiDchi_mc_dm, ChiDchi_dm)
push!.(ChiDchi_cfl_mc_dm, ChiDchi_cfl_dm)
BDIO_seek!(file)
BDIO_seek!(file)
end
BDIO_close!(file)
return Eoft0_mc,Eoft_mc,Qt0_mc,Qt_mc,pp_corr_mc,ap_corr_mc,pp_corr_t_mc,ap_corr_t_mc,pp_corr_tfl_mc,ap_corr_tfl_mc,Sigma_mc,Sigma_cfl_mc,Sigma2_mc,Sigma2_cfl_mc,ChiDchi_mc,ChiDchi_cfl_mc,pphat_t_mc,pptilde_t_mc,pphat0_mc,pptilde0_mc
return Eoft0_mc,Eoft_mc,Qt0_mc,Qt_mc,pp_corr_mc,ap_corr_mc,pp_corr_t_mc,ap_corr_t_mc,
pp_corr_tfl_mc,ap_corr_tfl_mc,Sigma_mc,Sigma_cfl_mc,Sigma2_mc,Sigma2_cfl_mc,ChiDchi_mc,
ChiDchi_cfl_mc,pphat_t_mc,pptilde_t_mc,pphat0_mc,pptilde0_mc,
pp_corr_mc_dm,ap_corr_mc_dm,pp_corr_t_mc_dm,ap_corr_t_mc_dm,
pp_corr_tfl_mc_dm,ap_corr_tfl_mc_dm,Sigma_mc_dm,Sigma_cfl_mc_dm,Sigma2_mc_dm,Sigma2_cfl_mc_dm,ChiDchi_mc_dm,
ChiDchi_cfl_mc_dm,pphat_t_mc_dm,pptilde_t_mc_dm,pphat0_mc_dm,pptilde0_mc_dm
end
......@@ -234,13 +403,16 @@ function uwff(file::String)
Loads in the variables 'Eoft','Eoft0','Qt','Qt0', 'pp_corr', 'ap_corr', 'pp_corr_t','ap_corr_t', 'pphat_t', 'pptilde_t','pphat0', 'pptilde0','pp_corr_tfl','ap_corr_tfl' 'Sigma'(v1 and v2), 'Sigma_cfl'(v1 and v2), 'ChiDchi' and 'ChiDchi_cfl' the corresponding
'uwreal' quantities. Noise average has been performed. The indices are, when corresponding, eucliden time and flow times
"""
function uwff(file::String)
Eoft0_mc,Eoft_mc,Qt0_mc,Qt_mc,pp_corr_mc,ap_corr_mc,pp_corr_t_mc,ap_corr_t_mc,pp_corr_tfl_mc,ap_corr_tfl_mc,Sigma_mc,Sigma_cfl_mc,Sigma2_mc,Sigma2_cfl_mc,ChiDchi_mc,ChiDchi_cfl_mc,pphat_t_mc,pptilde_t_mc,pphat0_mc,pptilde0_mc = read_ff(file)
function uwff(file::String; tfl = 0)
Eoft0_mc,Eoft_mc,Qt0_mc,Qt_mc,pp_corr_mc,ap_corr_mc,pp_corr_t_mc,ap_corr_t_mc,pp_corr_tfl_mc,ap_corr_tfl_mc,Sigma_mc,Sigma_cfl_mc,Sigma2_mc,Sigma2_cfl_mc,ChiDchi_mc,ChiDchi_cfl_mc,
pphat_t_mc,pptilde_t_mc,pphat0_mc,pptilde0_mc,pp_corr_mc_dm,ap_corr_mc_dm,pp_corr_t_mc_dm,ap_corr_t_mc_dm,pp_corr_tfl_mc_dm,ap_corr_tfl_mc_dm,Sigma_mc_dm,Sigma_cfl_mc_dm,Sigma2_mc_dm,
Sigma2_cfl_mc_dm,ChiDchi_mc_dm,ChiDchi_cfl_mc_dm,pphat_t_mc_dm,pptilde_t_mc_dm,pphat0_mc_dm,pptilde0_mc_dm = read_ff(file)
runame = String(split(split(file,"/")[end],".")[1])
T = size(pp_corr_mc)[1]
Nsff = size(pp_corr_mc)[2]
nsteps = size(pp_corr_t_mc)[3]-1
tfl == 0 ? nsteps = size(pp_corr_t_mc)[3]-1 : nsteps = tfl
Nsbf = size(pp_corr_tfl_mc)[2]
Nfl = size(pp_corr_tfl_mc)[3]
Nf = size(pp_corr_mc)[3]
......@@ -260,15 +432,56 @@ function uwff(file::String)
global pphat0 = [uwreal(sum(pphat0_mc,dims = 2)[i,1,f]./Nsff ,runame) for i in 1:T, f in 1:Nf]
global pptilde0 = [uwreal(real.(sum(pptilde0_mc,dims = 2)[i,1,f]./Nsff) ,runame) for i in 1:T, f in 1:Nf]
temp_matrix = Array{uwreal}(undef,T,nsteps+1,Nf)
pp_corr_t_mc = sum(pp_corr_t_mc,dims = 2)[:,1,:,:]./Nsff
global pp_corr_t = [uwreal(pp_corr_t_mc[i,k,f] ,runame) for i in 1:T, k in 1:nsteps+1, f in 1:Nf]
ap_corr_t_mc = real.(sum(ap_corr_t_mc,dims = 2)[:,1,:,:]./Nsff)
global ap_corr_t = [uwreal(ap_corr_t_mc[i,k,f] ,runame) for i in 1:T, k in 1:nsteps+1, f in 1:Nf]
for f in 1:Nf
for k in 1:(nsteps+1)
for i in 1:T
temp_matrix[i,k,f] = uwreal(pp_corr_t_mc[i,k,f] ,runame)
end
end
end
global pp_corr_t = copy(temp_matrix)
pp_corr_t_mc = nothing
# global pp_corr_t = [uwreal(pp_corr_t_mc[i,k,f] ,runame) for i in 1:T, k in 1:nsteps+1, f in 1:Nf]
# ap_corr_t_mc = real.(sum(ap_corr_t_mc, dims = 2)[:,1,:,:]./Nsff)
# # ap_corr_t_mc = real.(sum(ap_corr_t_mc,dims = 2)[:,1,:,:]./Nsff)
# for f in 1:Nf
# for k in 1:(nsteps+1)
# for i in 1:T
# temp_matrix[i,k,f] = uwreal(ap_corr_t_mc[i,k,f], runame)
# end
# end
# end
# global ap_corr_t = copy(temp_matrix)
ap_corr_t_mc = nothing
# # global ap_corr_t = [uwreal(ap_corr_t_mc[i,k,f] ,runame) for i in 1:T, k in 1:nsteps+1, f in 1:Nf]
pphat_t_mc = sum(pphat_t_mc,dims = 2)[:,1,:,:]./Nsff
global pphat_t = [uwreal(pphat_t_mc[i,k,f] ,runame) for i in 1:T, k in 1:nsteps+1, f in 1:Nf]
for f in 1:Nf
for k in 1:(nsteps+1)
for i in 1:T
temp_matrix[i,k,f] = uwreal(pphat_t_mc[i,k,f], runame)
end
end
end
global pphat_t = copy(temp_matrix)
pphat_t_mc = nothing
# global pphat_t = [uwreal(pphat_t_mc[i,k,f] ,runame) for i in 1:T, k in 1:nsteps+1, f in 1:Nf]
pptilde_t_mc = real.(sum(pptilde_t_mc,dims = 2)[:,1,:,:]./Nsff)
global pptilde_t = [uwreal(pptilde_t_mc[i,k,f] ,runame) for i in 1:T, k in 1:nsteps+1, f in 1:Nf]
for f in 1:Nf
for k in 1:(nsteps+1)
for i in 1:T
temp_matrix[i,k,f] = uwreal(pptilde_t_mc[i,k,f], runame)
end
end
end
global pptilde_t = copy(temp_matrix)
pptilde_t_mc = nothing
# global pptilde_t = [uwreal(pptilde_t_mc[i,k,f] ,runame) for i in 1:T, k in 1:nsteps+1, f in 1:Nf]
pp_corr_tfl_mc = sum(pp_corr_tfl_mc,dims = 2)[:,1,:,:]./Nsff
global pp_corr_tfl = [uwreal(pp_corr_tfl_mc[i,k,f] ,runame) for i in 1:T, k in 1:Nfl, f in 1:Nf]
......@@ -290,5 +503,68 @@ function uwff(file::String)
ChiDchi_cfl_mc = real.(sum(ChiDchi_cfl_mc,dims = 2)[:,1,:,:]./Nsbf)
global ChiDchi_cfl = [uwreal(ChiDchi_cfl_mc[i,k,f] ,runame) for i in 1:T, k in 1:Nfl, f in 1:Nf]
# d/dm
pp_corr_mc_dm = sum(pp_corr_mc_dm,dims = 2)[:,1,:]./Nsff
global pp_corr_dm = [uwreal(pp_corr_mc_dm[i,f] ,runame) for i in 1:T, f in 1:Nf]
ap_corr_mc_dm = real.(sum(ap_corr_mc_dm,dims = 2)[:,1,:]./Nsff)
global ap_corr_dm = [uwreal(ap_corr_mc_dm[i,f] ,runame) for i in 1:T, f in 1:Nf]
global pphat0_dm = [uwreal(sum(pphat0_mc_dm,dims = 2)[i,1,f]./Nsff ,runame) for i in 1:T, f in 1:Nf]
global pptilde0_dm = [uwreal(real.(sum(pptilde0_mc_dm,dims = 2)[i,1,f]./Nsff) ,runame) for i in 1:T, f in 1:Nf]
pp_corr_t_mc_dm = sum(pp_corr_t_mc_dm,dims = 2)[:,1,:,:]./Nsff
for k in 1:nsteps+1, f in 1:Nf, i in 1:T
temp_matrix[i,k,f] = uwreal(pp_corr_t_mc_dm[i,k,f], runame)
end
global pp_corr_t_dm = copy(temp_matrix)
pp_corr_t_mc_dm = nothing
# global pp_corr_t_dm = [uwreal(pp_corr_t_mc_dm[i,k,f] ,runame) for i in 1:T, k in 1:nsteps+1, f in 1:Nf]
# ap_corr_t_mc_dm = real.(sum(ap_corr_t_mc_dm,dims = 2)[:,1,:,:]./Nsff)
# for k in 1:nsteps+1, f in 1:Nf, i in 1:T
# temp_matrix[i,k,f] = uwreal(ap_corr_t_mc_dm[i,k,f], runame)
# end
# global ap_corr_t_dm = copy(temp_matrix)
ap_corr_t_mc_dm = nothing
# # global ap_corr_t_dm = [uwreal(ap_corr_t_mc_dm[i,k,f] ,runame) for i in 1:T, k in 1:nsteps+1, f in 1:Nf]
pphat_t_mc_dm = sum(pphat_t_mc_dm,dims = 2)[:,1,:,:]./Nsff
for k in 1:nsteps+1, f in 1:Nf, i in 1:T
temp_matrix[i,k,f] = uwreal(pphat_t_mc_dm[i,k,f], runame)
end
global pphat_t_dm = copy(temp_matrix)
pphat_t_mc_dm = nothing
# global pphat_t_dm = [uwreal(pphat_t_mc_dm[i,k,f], runame) for i in 1:T, k in 1:nsteps+1, f in 1:Nf]
pptilde_t_mc_dm = real.(sum(pptilde_t_mc_dm,dims = 2)[:,1,:,:]./Nsff)
for k in 1:nsteps+1, f in 1:Nf, i in 1:T
temp_matrix[i,k,f] = uwreal(pptilde_t_mc_dm[i,k,f], runame)
end
global pptilde_t_dm = copy(temp_matrix)
pptilde_t_mc_dm = nothing
# global pptilde_t_dm = [uwreal(pptilde_t_mc_dm[i,k,f], runame) for i in 1:T, k in 1:nsteps+1, f in 1:Nf]
pp_corr_tfl_mc_dm = sum(pp_corr_tfl_mc_dm,dims = 2)[:,1,:,:]./Nsff
global pp_corr_tfl_dm = [uwreal(pp_corr_tfl_mc_dm[i,k,f] ,runame) for i in 1:T, k in 1:Nfl, f in 1:Nf]
ap_corr_tfl_mc_dm = real.(sum(ap_corr_tfl_mc_dm,dims = 2)[:,1,:,:]./Nsff)
global ap_corr_tfl_dm = [uwreal(ap_corr_tfl_mc_dm[i,k,f] ,runame) for i in 1:T, k in 1:Nfl, f in 1:Nf]
Sigma_mc_dm = real.(sum(Sigma_mc_dm,dims = 2)[:,1,:,:]./Nsbf)
global Sigma_dm = [uwreal(Sigma_mc_dm[i,k,f] ,runame) for i in 1:T, k in 1:Nfl, f in 1:Nf]
Sigma_cfl_mc_dm = real.(sum(Sigma_cfl_mc_dm,dims = 2)[:,1,:,:]./Nsbf)
global Sigma_cfl_dm = [uwreal(Sigma_cfl_mc_dm[i,k,f] ,runame) for i in 1:T, k in 1:Nfl, f in 1:Nf]
Sigma2_mc_dm = real.(sum(Sigma2_mc_dm,dims = 2)[:,1,:,:]./Nsbf)
global Sigma2_dm = [uwreal(Sigma2_mc_dm[i,k,f] ,runame) for i in 1:T, k in 1:Nfl, f in 1:Nf]
Sigma2_cfl_mc_dm = real.(sum(Sigma2_cfl_mc_dm,dims = 2)[:,1,:,:]./Nsbf)
global Sigma2_cfl_dm = [uwreal(Sigma2_cfl_mc_dm[i,k,f] ,runame) for i in 1:T, k in 1:Nfl, f in 1:Nf]
ChiDchi_mc_dm = real.(sum(ChiDchi_mc_dm,dims = 2)[:,1,:,:]./Nsbf)
global ChiDchi_dm = [uwreal(ChiDchi_mc_dm[i,k,f] ,runame) for i in 1:T, k in 1:Nfl, f in 1:Nf]
ChiDchi_cfl_mc_dm = real.(sum(ChiDchi_cfl_mc_dm,dims = 2)[:,1,:,:]./Nsbf)
global ChiDchi_cfl_dm = [uwreal(ChiDchi_cfl_mc_dm[i,k,f] ,runame) for i in 1:T, k in 1:Nfl, f in 1:Nf]
return nothing
end
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment