diff --git a/.github/workflows/TagBot.yml b/.github/workflows/TagBot.yml
old mode 100644
new mode 100755
diff --git a/.github/workflows/documentation.yml b/.github/workflows/documentation.yml
old mode 100644
new mode 100755
diff --git a/.github/workflows/draft-pdf.yml b/.github/workflows/draft-pdf.yml
old mode 100644
new mode 100755
diff --git a/.gitignore b/.gitignore
old mode 100644
new mode 100755
diff --git a/.travis.yml b/.travis.yml
old mode 100644
new mode 100755
diff --git a/CITATION.bib b/CITATION.bib
old mode 100644
new mode 100755
diff --git a/LICENSE.md b/LICENSE.md
old mode 100644
new mode 100755
diff --git a/Manifest.toml b/Manifest.toml
old mode 100644
new mode 100755
index 64444aac..5168e07d
--- a/Manifest.toml
+++ b/Manifest.toml
@@ -13,14 +13,20 @@ git-tree-sha1 = "03e0550477d86222521d254b741d470ba17ea0b5"
 uuid = "1520ce14-60c1-5f80-bbc7-55ef81b5835c"
 version = "0.3.4"
 
+[[deps.Adapt]]
+deps = ["LinearAlgebra"]
+git-tree-sha1 = "195c5505521008abea5aee4f96930717958eac6f"
+uuid = "79e6a3ab-5dfb-504d-930d-738a2a938a0e"
+version = "3.4.0"
+
 [[deps.ArgTools]]
 uuid = "0dad84c5-d112-42e6-8d28-ef12dabb789f"
 
 [[deps.ArrayInterfaceCore]]
 deps = ["LinearAlgebra", "SparseArrays", "SuiteSparse"]
-git-tree-sha1 = "7d255eb1d2e409335835dc8624c35d97453011eb"
+git-tree-sha1 = "40debc9f72d0511e12d817c7ca06a721b6423ba3"
 uuid = "30b0a656-2188-435a-8636-2ec0e6a096e2"
-version = "0.1.14"
+version = "0.1.17"
 
 [[deps.Artifacts]]
 uuid = "56f22d72-fd6d-98f1-02f0-08ddc0907c33"
@@ -46,6 +52,18 @@ git-tree-sha1 = "91d6baa911283650df649d0aea7c28639273ae7b"
 uuid = "0b7ba130-8d10-5ba8-a3d6-c5182647fed9"
 version = "1.21.1+0"
 
+[[deps.Bzip2_jll]]
+deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
+git-tree-sha1 = "19a35467a82e236ff51bc17a3a44b69ef35185a2"
+uuid = "6e34b625-4abd-537c-b88f-471c36dfa7a0"
+version = "1.0.8+0"
+
+[[deps.Cairo_jll]]
+deps = ["Artifacts", "Bzip2_jll", "Fontconfig_jll", "FreeType2_jll", "Glib_jll", "JLLWrappers", "LZO_jll", "Libdl", "Pixman_jll", "Pkg", "Xorg_libXext_jll", "Xorg_libXrender_jll", "Zlib_jll", "libpng_jll"]
+git-tree-sha1 = "4b859a208b2397a7a623a03449e4636bdb17bcf2"
+uuid = "83423d85-b0ee-5818-9007-b63ccbeb887a"
+version = "1.16.1+1"
+
 [[deps.ChainRulesCore]]
 deps = ["Compat", "LinearAlgebra", "SparseArrays"]
 git-tree-sha1 = "80ca332f6dcb2508adba68f22f551adb2d00a624"
@@ -58,6 +76,42 @@ git-tree-sha1 = "38f7a08f19d8810338d4f5085211c7dfa5d5bdd8"
 uuid = "9e997f8a-9a97-42d5-a9f1-ce6bfc15e2c0"
 version = "0.1.4"
 
+[[deps.CodeTracking]]
+deps = ["InteractiveUtils", "UUIDs"]
+git-tree-sha1 = "1833bda4a027f4b2a1c984baddcf755d77266818"
+uuid = "da1fd8a2-8d9e-5ec2-8556-3022fb5608a2"
+version = "1.1.0"
+
+[[deps.CodecZlib]]
+deps = ["TranscodingStreams", "Zlib_jll"]
+git-tree-sha1 = "ded953804d019afa9a3f98981d99b33e3db7b6da"
+uuid = "944b1d66-785c-5afd-91f1-9de20f533193"
+version = "0.7.0"
+
+[[deps.ColorSchemes]]
+deps = ["ColorTypes", "ColorVectorSpace", "Colors", "FixedPointNumbers", "Random"]
+git-tree-sha1 = "1fd869cc3875b57347f7027521f561cf46d1fcd8"
+uuid = "35d6a980-a343-548e-a6ea-1d62b119f2f4"
+version = "3.19.0"
+
+[[deps.ColorTypes]]
+deps = ["FixedPointNumbers", "Random"]
+git-tree-sha1 = "eb7f0f8307f71fac7c606984ea5fb2817275d6e4"
+uuid = "3da002f7-5984-5a60-b8a6-cbb66c0b333f"
+version = "0.11.4"
+
+[[deps.ColorVectorSpace]]
+deps = ["ColorTypes", "FixedPointNumbers", "LinearAlgebra", "SpecialFunctions", "Statistics", "TensorCore"]
+git-tree-sha1 = "d08c20eef1f2cbc6e60fd3612ac4340b89fea322"
+uuid = "c3611d14-8923-5661-9e6a-0046d554d3a4"
+version = "0.9.9"
+
+[[deps.Colors]]
+deps = ["ColorTypes", "FixedPointNumbers", "Reexport"]
+git-tree-sha1 = "417b0ed7b8b838aa6ca0a87aadf1bb9eb111ce40"
+uuid = "5ae59095-9a9b-59fe-a467-6f913c188581"
+version = "0.12.8"
+
 [[deps.CommonSubexpressions]]
 deps = ["MacroTools", "Test"]
 git-tree-sha1 = "7b8a93dba8af7e3b42fecabf646260105ac373f7"
@@ -80,6 +134,11 @@ git-tree-sha1 = "59d00b3139a9de4eb961057eabb65ac6522be954"
 uuid = "187b0558-2788-49d3-abe0-74a17ed4e7c9"
 version = "1.4.0"
 
+[[deps.Contour]]
+git-tree-sha1 = "d05d9e7b7aedff4e5b51a029dced05cfb6125781"
+uuid = "d38c429a-6771-53c6-b99e-75d170b6e991"
+version = "0.6.2"
+
 [[deps.DataAPI]]
 git-tree-sha1 = "fb5f5316dd3fd4c5e7c30a24d50643b73e37cd40"
 uuid = "9a962f9c-6df0-11e9-0e5d-c546b8b5ee8a"
@@ -91,6 +150,11 @@ git-tree-sha1 = "d1fff3a548102f48987a52a2e0d114fa97d730f0"
 uuid = "864edb3b-99cc-5e75-8d2d-829cb0a9cfe8"
 version = "0.18.13"
 
+[[deps.DataValueInterfaces]]
+git-tree-sha1 = "bfc1187b79289637fa0ef6d4436ebdfe6905cbd6"
+uuid = "e2d170a0-9d28-54be-80f0-106bbe20a464"
+version = "1.0.0"
+
 [[deps.Dates]]
 deps = ["Printf"]
 uuid = "ade2ca70-3891-5945-98fb-dc099432e06a"
@@ -149,11 +213,35 @@ version = "0.1.14"
 deps = ["ArgTools", "FileWatching", "LibCURL", "NetworkOptions"]
 uuid = "f43a241f-c20a-4ad4-852c-f6b1247861c6"
 
+[[deps.EarCut_jll]]
+deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
+git-tree-sha1 = "3f3a2501fa7236e9b911e0f7a588c657e822bb6d"
+uuid = "5ae413db-bbd1-5e63-b57d-d24a61df00f5"
+version = "2.2.3+0"
+
+[[deps.Expat_jll]]
+deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
+git-tree-sha1 = "bad72f730e9e91c08d9427d5e8db95478a3c323d"
+uuid = "2e619515-83b5-522b-bb60-26c02a35a201"
+version = "2.4.8+0"
+
+[[deps.FFMPEG]]
+deps = ["FFMPEG_jll"]
+git-tree-sha1 = "b57e3acbe22f8484b4b5ff66a7499717fe1a9cc8"
+uuid = "c87230d0-a227-11e9-1b43-d7ebe4e7570a"
+version = "0.4.1"
+
+[[deps.FFMPEG_jll]]
+deps = ["Artifacts", "Bzip2_jll", "FreeType2_jll", "FriBidi_jll", "JLLWrappers", "LAME_jll", "Libdl", "Ogg_jll", "OpenSSL_jll", "Opus_jll", "Pkg", "Zlib_jll", "libaom_jll", "libass_jll", "libfdk_aac_jll", "libvorbis_jll", "x264_jll", "x265_jll"]
+git-tree-sha1 = "ccd479984c7838684b3ac204b716c89955c76623"
+uuid = "b22a6f82-2f65-5046-a5b2-351ab43fb4e5"
+version = "4.4.2+0"
+
 [[deps.FileIO]]
 deps = ["Pkg", "Requires", "UUIDs"]
-git-tree-sha1 = "9267e5f50b0e12fdfd5a2455534345c4cf2c7f7a"
+git-tree-sha1 = "94f5101b96d2d968ace56f7f2db19d0a5f592e28"
 uuid = "5789e2e9-d7fb-5bc7-8068-2c6fae9b9549"
-version = "1.14.0"
+version = "1.15.0"
 
 [[deps.FileWatching]]
 uuid = "7b1f6079-737a-58dc-b8bc-7a2ca5c1b5ee"
@@ -170,16 +258,93 @@ git-tree-sha1 = "cb8c5f0074153ace28ce5100714df4378ad885e0"
 uuid = "6a86dc24-6348-571c-b903-95158fe2bd41"
 version = "2.14.0"
 
+[[deps.FixedPointNumbers]]
+deps = ["Statistics"]
+git-tree-sha1 = "335bfdceacc84c5cdf16aadc768aa5ddfc5383cc"
+uuid = "53c48c17-4a7d-5ca2-90c5-79b7896eea93"
+version = "0.8.4"
+
+[[deps.Fontconfig_jll]]
+deps = ["Artifacts", "Bzip2_jll", "Expat_jll", "FreeType2_jll", "JLLWrappers", "Libdl", "Libuuid_jll", "Pkg", "Zlib_jll"]
+git-tree-sha1 = "21efd19106a55620a188615da6d3d06cd7f6ee03"
+uuid = "a3f928ae-7b40-5064-980b-68af3947d34b"
+version = "2.13.93+0"
+
+[[deps.Formatting]]
+deps = ["Printf"]
+git-tree-sha1 = "8339d61043228fdd3eb658d86c926cb282ae72a8"
+uuid = "59287772-0a20-5a39-b81b-1366585eb4c0"
+version = "0.4.2"
+
 [[deps.ForwardDiff]]
 deps = ["CommonSubexpressions", "DiffResults", "DiffRules", "LinearAlgebra", "LogExpFunctions", "NaNMath", "Preferences", "Printf", "Random", "SpecialFunctions", "StaticArrays"]
 git-tree-sha1 = "2f18915445b248731ec5db4e4a17e451020bf21e"
 uuid = "f6369f11-7733-5829-9624-2563aa707210"
 version = "0.10.30"
 
+[[deps.FreeType2_jll]]
+deps = ["Artifacts", "Bzip2_jll", "JLLWrappers", "Libdl", "Pkg", "Zlib_jll"]
+git-tree-sha1 = "87eb71354d8ec1a96d4a7636bd57a7347dde3ef9"
+uuid = "d7e528f0-a631-5988-bf34-fe36492bcfd7"
+version = "2.10.4+0"
+
+[[deps.FriBidi_jll]]
+deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
+git-tree-sha1 = "aa31987c2ba8704e23c6c8ba8a4f769d5d7e4f91"
+uuid = "559328eb-81f9-559d-9380-de523a88c83c"
+version = "1.0.10+0"
+
 [[deps.Future]]
 deps = ["Random"]
 uuid = "9fa8497b-333b-5362-9e8d-4d0656e87820"
 
+[[deps.GLFW_jll]]
+deps = ["Artifacts", "JLLWrappers", "Libdl", "Libglvnd_jll", "Pkg", "Xorg_libXcursor_jll", "Xorg_libXi_jll", "Xorg_libXinerama_jll", "Xorg_libXrandr_jll"]
+git-tree-sha1 = "51d2dfe8e590fbd74e7a842cf6d13d8a2f45dc01"
+uuid = "0656b61e-2033-5cc2-a64a-77c0f6c09b89"
+version = "3.3.6+0"
+
+[[deps.GR]]
+deps = ["Base64", "DelimitedFiles", "GR_jll", "HTTP", "JSON", "Libdl", "LinearAlgebra", "Pkg", "Printf", "Random", "RelocatableFolders", "Serialization", "Sockets", "Test", "UUIDs"]
+git-tree-sha1 = "037a1ca47e8a5989cc07d19729567bb71bfabd0c"
+uuid = "28b8d3ca-fb5f-59d9-8090-bfdbd6d07a71"
+version = "0.66.0"
+
+[[deps.GR_jll]]
+deps = ["Artifacts", "Bzip2_jll", "Cairo_jll", "FFMPEG_jll", "Fontconfig_jll", "GLFW_jll", "JLLWrappers", "JpegTurbo_jll", "Libdl", "Libtiff_jll", "Pixman_jll", "Pkg", "Qt5Base_jll", "Zlib_jll", "libpng_jll"]
+git-tree-sha1 = "c8ab731c9127cd931c93221f65d6a1008dad7256"
+uuid = "d2c73de3-f751-5644-a686-071e5b155ba9"
+version = "0.66.0+0"
+
+[[deps.GeometryBasics]]
+deps = ["EarCut_jll", "IterTools", "LinearAlgebra", "StaticArrays", "StructArrays", "Tables"]
+git-tree-sha1 = "83ea630384a13fc4f002b77690bc0afeb4255ac9"
+uuid = "5c1252a2-5f33-56bf-86c9-59e7332b4326"
+version = "0.4.2"
+
+[[deps.Gettext_jll]]
+deps = ["Artifacts", "CompilerSupportLibraries_jll", "JLLWrappers", "Libdl", "Libiconv_jll", "Pkg", "XML2_jll"]
+git-tree-sha1 = "9b02998aba7bf074d14de89f9d37ca24a1a0b046"
+uuid = "78b55507-aeef-58d4-861c-77aaff3498b1"
+version = "0.21.0+0"
+
+[[deps.Glib_jll]]
+deps = ["Artifacts", "Gettext_jll", "JLLWrappers", "Libdl", "Libffi_jll", "Libiconv_jll", "Libmount_jll", "PCRE_jll", "Pkg", "Zlib_jll"]
+git-tree-sha1 = "a32d672ac2c967f3deb8a81d828afc739c838a06"
+uuid = "7746bdde-850d-59dc-9ae8-88ece973131d"
+version = "2.68.3+2"
+
+[[deps.Graphite2_jll]]
+deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
+git-tree-sha1 = "344bf40dcab1073aca04aa0df4fb092f920e4011"
+uuid = "3b182d85-2403-5c21-9c21-1e1f0cc25472"
+version = "1.3.14+0"
+
+[[deps.Grisu]]
+git-tree-sha1 = "53bb909d1151e57e2484c3d1b53e19552b887fb2"
+uuid = "42e2da0e-8278-4e71-bc24-59509adca0fe"
+version = "1.0.2"
+
 [[deps.Gumbo]]
 deps = ["AbstractTrees", "Gumbo_jll", "Libdl"]
 git-tree-sha1 = "a4b03190a8d79e95b2f418df44f74cbc34289920"
@@ -200,9 +365,21 @@ version = "0.15.7"
 
 [[deps.HDF5_jll]]
 deps = ["Artifacts", "JLLWrappers", "LibCURL_jll", "Libdl", "OpenSSL_jll", "Pkg", "Zlib_jll"]
-git-tree-sha1 = "bab67c0d1c4662d2c4be8c6007751b0b6111de5c"
+git-tree-sha1 = "c003b31e2e818bc512b0ff99d7dce03b0c1359f5"
 uuid = "0234f1f7-429e-5d53-9886-15a909be8d59"
-version = "1.12.1+0"
+version = "1.12.2+1"
+
+[[deps.HTTP]]
+deps = ["Base64", "CodecZlib", "Dates", "IniFile", "Logging", "LoggingExtras", "MbedTLS", "NetworkOptions", "Random", "SimpleBufferStream", "Sockets", "URIs", "UUIDs"]
+git-tree-sha1 = "ed47af35905b7cc8f1a522ca684b35a212269bd8"
+uuid = "cd3eb016-35fb-5094-929b-558a96fad6f3"
+version = "1.2.0"
+
+[[deps.HarfBuzz_jll]]
+deps = ["Artifacts", "Cairo_jll", "Fontconfig_jll", "FreeType2_jll", "Glib_jll", "Graphite2_jll", "JLLWrappers", "Libdl", "Libffi_jll", "Pkg"]
+git-tree-sha1 = "129acf094d168394e80ee1dc4bc06ec835e510a3"
+uuid = "2e76f6c2-a576-52d4-95c1-20adfe4de566"
+version = "2.8.1+1"
 
 [[deps.IOCapture]]
 deps = ["Logging", "Random"]
@@ -210,6 +387,11 @@ git-tree-sha1 = "f7be53659ab06ddc986428d3a9dcc95f6fa6705a"
 uuid = "b5f81e59-6552-4d32-b1f0-c071b021bf89"
 version = "0.2.2"
 
+[[deps.IniFile]]
+git-tree-sha1 = "f550e6e32074c939295eb5ea6de31849ac2c9625"
+uuid = "83e8ac13-25f8-5344-8a64-a9f2b223428f"
+version = "0.5.1"
+
 [[deps.InteractiveUtils]]
 deps = ["Markdown"]
 uuid = "b77e0a4c-d291-57a0-90e8-8db25a27a240"
@@ -225,6 +407,16 @@ git-tree-sha1 = "7fd44fd4ff43fc60815f8e764c0f352b83c49151"
 uuid = "92d709cd-6900-40b7-9082-c6be49f344b6"
 version = "0.1.1"
 
+[[deps.IterTools]]
+git-tree-sha1 = "fa6287a4469f5e048d763df38279ee729fbd44e5"
+uuid = "c8e1da08-722c-5040-9ed9-7db0dc04731e"
+version = "1.4.0"
+
+[[deps.IteratorInterfaceExtensions]]
+git-tree-sha1 = "a3f24677c21f5bbe9d2a714f95dcd58337fb2856"
+uuid = "82899510-4779-5014-852e-03e436cf321d"
+version = "1.0.0"
+
 [[deps.JLD]]
 deps = ["Compat", "FileIO", "HDF5", "Printf"]
 git-tree-sha1 = "390ed210fcbdaffde2efe7890f39bd1e681a62d7"
@@ -249,6 +441,47 @@ git-tree-sha1 = "3c837543ddb02250ef42f4738347454f95079d4e"
 uuid = "682c06a0-de6a-54ab-a142-c8b1cf79cde6"
 version = "0.21.3"
 
+[[deps.JpegTurbo_jll]]
+deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
+git-tree-sha1 = "b53380851c6e6664204efb2e62cd24fa5c47e4ba"
+uuid = "aacddb02-875f-59d6-b918-886e6ef4fbf8"
+version = "2.1.2+0"
+
+[[deps.JuliaInterpreter]]
+deps = ["CodeTracking", "InteractiveUtils", "Random", "UUIDs"]
+git-tree-sha1 = "0f960b1404abb0b244c1ece579a0ec78d056a5d1"
+uuid = "aa1ae85d-cabe-5617-a682-6adf51b2e16a"
+version = "0.9.15"
+
+[[deps.LAME_jll]]
+deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
+git-tree-sha1 = "f6250b16881adf048549549fba48b1161acdac8c"
+uuid = "c1c5ebd0-6772-5130-a774-d5fcae4a789d"
+version = "3.100.1+0"
+
+[[deps.LERC_jll]]
+deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
+git-tree-sha1 = "bf36f528eec6634efc60d7ec062008f171071434"
+uuid = "88015f11-f218-50d7-93a8-a6af411a945d"
+version = "3.0.0+1"
+
+[[deps.LZO_jll]]
+deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
+git-tree-sha1 = "e5b909bcf985c5e2605737d2ce278ed791b89be6"
+uuid = "dd4b983a-f0e5-5f8d-a1b7-129d4a5fb1ac"
+version = "2.10.1+0"
+
+[[deps.LaTeXStrings]]
+git-tree-sha1 = "f2355693d6778a178ade15952b7ac47a4ff97996"
+uuid = "b964fa9f-0449-5b57-a5c2-d3ea65f4040f"
+version = "1.3.0"
+
+[[deps.Latexify]]
+deps = ["Formatting", "InteractiveUtils", "LaTeXStrings", "MacroTools", "Markdown", "Printf", "Requires"]
+git-tree-sha1 = "1a43be956d433b5d0321197150c2f94e16c0aaa0"
+uuid = "23fbe1c1-3f47-55db-b15f-69d7ec21a316"
+version = "0.15.16"
+
 [[deps.LibCURL]]
 deps = ["LibCURL_jll", "MozillaCACerts_jll"]
 uuid = "b27032c2-a3e7-50c8-80cd-2d36dbcbfd21"
@@ -268,6 +501,54 @@ uuid = "29816b5a-b9ab-546f-933c-edad1886dfa8"
 [[deps.Libdl]]
 uuid = "8f399da3-3557-5675-b5ff-fb832c97cbdb"
 
+[[deps.Libffi_jll]]
+deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
+git-tree-sha1 = "0b4a5d71f3e5200a7dff793393e09dfc2d874290"
+uuid = "e9f186c6-92d2-5b65-8a66-fee21dc1b490"
+version = "3.2.2+1"
+
+[[deps.Libgcrypt_jll]]
+deps = ["Artifacts", "JLLWrappers", "Libdl", "Libgpg_error_jll", "Pkg"]
+git-tree-sha1 = "64613c82a59c120435c067c2b809fc61cf5166ae"
+uuid = "d4300ac3-e22c-5743-9152-c294e39db1e4"
+version = "1.8.7+0"
+
+[[deps.Libglvnd_jll]]
+deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg", "Xorg_libX11_jll", "Xorg_libXext_jll"]
+git-tree-sha1 = "7739f837d6447403596a75d19ed01fd08d6f56bf"
+uuid = "7e76a0d4-f3c7-5321-8279-8d96eeed0f29"
+version = "1.3.0+3"
+
+[[deps.Libgpg_error_jll]]
+deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
+git-tree-sha1 = "c333716e46366857753e273ce6a69ee0945a6db9"
+uuid = "7add5ba3-2f88-524e-9cd5-f83b8a55f7b8"
+version = "1.42.0+0"
+
+[[deps.Libiconv_jll]]
+deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
+git-tree-sha1 = "42b62845d70a619f063a7da093d995ec8e15e778"
+uuid = "94ce4f54-9a6c-5748-9c1c-f9c7231a4531"
+version = "1.16.1+1"
+
+[[deps.Libmount_jll]]
+deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
+git-tree-sha1 = "9c30530bf0effd46e15e0fdcf2b8636e78cbbd73"
+uuid = "4b2f31a3-9ecc-558c-b454-b3730dcb73e9"
+version = "2.35.0+0"
+
+[[deps.Libtiff_jll]]
+deps = ["Artifacts", "JLLWrappers", "JpegTurbo_jll", "LERC_jll", "Libdl", "Pkg", "Zlib_jll", "Zstd_jll"]
+git-tree-sha1 = "3eb79b0ca5764d4799c06699573fd8f533259713"
+uuid = "89763e89-9b03-5906-acba-b20f662cd828"
+version = "4.4.0+0"
+
+[[deps.Libuuid_jll]]
+deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
+git-tree-sha1 = "7f3efec06033682db852f8b3bc3c1d2b0a0ab066"
+uuid = "38a345b3-de98-5d2b-a5d3-14cd9215e700"
+version = "2.36.0+0"
+
 [[deps.LineSearches]]
 deps = ["LinearAlgebra", "NLSolversBase", "NaNMath", "Parameters", "Printf"]
 git-tree-sha1 = "f27132e551e959b3667d8c93eae90973225032dd"
@@ -280,13 +561,25 @@ uuid = "37e2e46d-f89d-539d-b4ee-838fcccc9c8e"
 
 [[deps.LogExpFunctions]]
 deps = ["ChainRulesCore", "ChangesOfVariables", "DocStringExtensions", "InverseFunctions", "IrrationalConstants", "LinearAlgebra"]
-git-tree-sha1 = "09e4b894ce6a976c354a69041a04748180d43637"
+git-tree-sha1 = "361c2b088575b07946508f135ac556751240091c"
 uuid = "2ab3a3ac-af41-5b50-aa03-7779005ae688"
-version = "0.3.15"
+version = "0.3.17"
 
 [[deps.Logging]]
 uuid = "56ddb016-857b-54e1-b83d-db4d58db5568"
 
+[[deps.LoggingExtras]]
+deps = ["Dates", "Logging"]
+git-tree-sha1 = "5d4d2d9904227b8bd66386c1138cf4d5ffa826bf"
+uuid = "e6f89c97-d47a-5376-807f-9c37f3926c36"
+version = "0.4.9"
+
+[[deps.LoweredCodeUtils]]
+deps = ["JuliaInterpreter"]
+git-tree-sha1 = "dedbebe234e06e1ddad435f5c6f4b85cd8ce55f7"
+uuid = "6f1432cf-f94c-5a45-995e-cdbf5db27b0b"
+version = "2.2.2"
+
 [[deps.LsqFit]]
 deps = ["Distributions", "ForwardDiff", "LinearAlgebra", "NLSolversBase", "OptimBase", "Random", "StatsBase"]
 git-tree-sha1 = "91aa1442e63a77f101aff01dec5a821a17f43922"
@@ -309,10 +602,21 @@ version = "0.5.9"
 deps = ["Base64"]
 uuid = "d6f4376e-aef5-505a-96c1-9c027394607a"
 
+[[deps.MbedTLS]]
+deps = ["Dates", "MbedTLS_jll", "MozillaCACerts_jll", "Random", "Sockets"]
+git-tree-sha1 = "14cb991ee7ccc6dabda93d310400575c3cae435b"
+uuid = "739be429-bea8-5141-9913-cc70e7f3736d"
+version = "1.1.2"
+
 [[deps.MbedTLS_jll]]
 deps = ["Artifacts", "Libdl"]
 uuid = "c8ffd9c3-330d-5841-b78e-0817d7145fa1"
 
+[[deps.Measures]]
+git-tree-sha1 = "e498ddeee6f9fdb4551ce855a46f54dbd900245f"
+uuid = "442fdcdd-2543-5da2-b0f3-8c86c306513e"
+version = "0.3.1"
+
 [[deps.Missings]]
 deps = ["DataAPI"]
 git-tree-sha1 = "bf210ce90b6c9eed32d25dbcae1ebc565df2687f"
@@ -339,6 +643,12 @@ version = "0.3.7"
 [[deps.NetworkOptions]]
 uuid = "ca575930-c2e3-43a9-ace4-1e988b2c1908"
 
+[[deps.Ogg_jll]]
+deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
+git-tree-sha1 = "887579a3eb005446d514ab7aeac5d1d027658b8f"
+uuid = "e7412a2a-1a6e-54c0-be00-318e2571c051"
+version = "1.3.5+1"
+
 [[deps.OpenBLAS_jll]]
 deps = ["Artifacts", "CompilerSupportLibraries_jll", "Libdl"]
 uuid = "4536629a-c528-5b80-bd46-f80d51c5b363"
@@ -377,11 +687,23 @@ git-tree-sha1 = "9cb1fee807b599b5f803809e85c81b582d2009d6"
 uuid = "87e2bd06-a317-5318-96d9-3ecbac512eee"
 version = "2.0.2"
 
+[[deps.Opus_jll]]
+deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
+git-tree-sha1 = "51a08fb14ec28da2ec7a927c4337e4332c2a4720"
+uuid = "91d4177d-7536-5919-b921-800302f37372"
+version = "1.3.2+0"
+
 [[deps.OrderedCollections]]
 git-tree-sha1 = "85f8e6578bf1f9ee0d11e7bb1b1456435479d47c"
 uuid = "bac558e1-5e72-5ebc-8fee-abe8a469f55d"
 version = "1.4.1"
 
+[[deps.PCRE_jll]]
+deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
+git-tree-sha1 = "b2a7af664e098055a7529ad1a900ded962bca488"
+uuid = "2f80f16e-611a-54ab-bc61-aa92de5b98fc"
+version = "8.44.0+0"
+
 [[deps.PDMats]]
 deps = ["LinearAlgebra", "SparseArrays", "SuiteSparse"]
 git-tree-sha1 = "cf494dca75a69712a72b80bc48f59dcf3dea63ec"
@@ -400,10 +722,34 @@ git-tree-sha1 = "0044b23da09b5608b4ecacb4e5e6c6332f833a7e"
 uuid = "69de0a69-1ddd-5017-9359-2bf0b02dc9f0"
 version = "2.3.2"
 
+[[deps.Pixman_jll]]
+deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
+git-tree-sha1 = "b4f5d02549a10e20780a24fce72bea96b6329e29"
+uuid = "30392449-352a-5448-841d-b1acce4e97dc"
+version = "0.40.1+0"
+
 [[deps.Pkg]]
 deps = ["Artifacts", "Dates", "Downloads", "LibGit2", "Libdl", "Logging", "Markdown", "Printf", "REPL", "Random", "SHA", "Serialization", "TOML", "Tar", "UUIDs", "p7zip_jll"]
 uuid = "44cfe95a-1eb2-52ea-b672-e2afdf69b78f"
 
+[[deps.PlotThemes]]
+deps = ["PlotUtils", "Statistics"]
+git-tree-sha1 = "8162b2f8547bc23876edd0c5181b27702ae58dce"
+uuid = "ccf2f8ad-2431-5c83-bf29-c5338b663b6a"
+version = "3.0.0"
+
+[[deps.PlotUtils]]
+deps = ["ColorSchemes", "Colors", "Dates", "Printf", "Random", "Reexport", "Statistics"]
+git-tree-sha1 = "9888e59493658e476d3073f1ce24348bdc086660"
+uuid = "995b91a9-d308-5afd-9ec6-746e21dbc043"
+version = "1.3.0"
+
+[[deps.Plots]]
+deps = ["Base64", "Contour", "Dates", "Downloads", "FFMPEG", "FixedPointNumbers", "GR", "GeometryBasics", "JSON", "Latexify", "LinearAlgebra", "Measures", "NaNMath", "Pkg", "PlotThemes", "PlotUtils", "Printf", "REPL", "Random", "RecipesBase", "RecipesPipeline", "Reexport", "Requires", "Scratch", "Showoff", "SparseArrays", "Statistics", "StatsBase", "UUIDs", "UnicodeFun", "Unzip"]
+git-tree-sha1 = "0a0da27969e8b6b2ee67c112dcf7001a659049a0"
+uuid = "91a5bcdd-55d7-5caf-9e0b-520d859cae80"
+version = "1.31.4"
+
 [[deps.PositiveFactorizations]]
 deps = ["LinearAlgebra"]
 git-tree-sha1 = "17275485f373e6673f7e7f97051f703ed5b15b20"
@@ -424,6 +770,12 @@ uuid = "de0858da-6303-5e67-8744-51eddeeeb8d7"
 deps = ["Printf"]
 uuid = "9abbd945-dff8-562f-b5e8-e1ebf5ef1b79"
 
+[[deps.Qt5Base_jll]]
+deps = ["Artifacts", "CompilerSupportLibraries_jll", "Fontconfig_jll", "Glib_jll", "JLLWrappers", "Libdl", "Libglvnd_jll", "OpenSSL_jll", "Pkg", "Xorg_libXext_jll", "Xorg_libxcb_jll", "Xorg_xcb_util_image_jll", "Xorg_xcb_util_keysyms_jll", "Xorg_xcb_util_renderutil_jll", "Xorg_xcb_util_wm_jll", "Zlib_jll", "xkbcommon_jll"]
+git-tree-sha1 = "c6c0f690d0cc7caddb74cef7aa847b824a16b256"
+uuid = "ea2cea3b-5b76-57ae-a6ef-0a8af62496e1"
+version = "5.15.3+1"
+
 [[deps.QuadGK]]
 deps = ["DataStructures", "LinearAlgebra"]
 git-tree-sha1 = "78aadffb3efd2155af139781b8a8df1ef279ea39"
@@ -438,17 +790,40 @@ uuid = "3fa0cd96-eef1-5676-8a61-b3b8758bbffb"
 deps = ["SHA", "Serialization"]
 uuid = "9a3f8284-a2c9-5f02-9a11-845980a1fd5c"
 
+[[deps.RecipesBase]]
+git-tree-sha1 = "6bf3f380ff52ce0832ddd3a2a7b9538ed1bcca7d"
+uuid = "3cdcf5f2-1ef4-517c-9805-6587b60abb01"
+version = "1.2.1"
+
+[[deps.RecipesPipeline]]
+deps = ["Dates", "NaNMath", "PlotUtils", "RecipesBase"]
+git-tree-sha1 = "e7eac76a958f8664f2718508435d058168c7953d"
+uuid = "01d81517-befc-4cb6-b9ec-a95719d0359c"
+version = "0.6.3"
+
 [[deps.Reexport]]
 git-tree-sha1 = "45e428421666073eab6f2da5c9d310d99bb12f9b"
 uuid = "189a3867-3050-52da-a836-e630ba90ab69"
 version = "1.2.2"
 
+[[deps.RelocatableFolders]]
+deps = ["SHA", "Scratch"]
+git-tree-sha1 = "22c5201127d7b243b9ee1de3b43c408879dff60f"
+uuid = "05181044-ff0b-4ac5-8273-598c1e38db00"
+version = "0.3.0"
+
 [[deps.Requires]]
 deps = ["UUIDs"]
 git-tree-sha1 = "838a3a4188e2ded87a4f9f184b4b0d78a1e91cb7"
 uuid = "ae029012-a4dd-5104-9daa-d747884805df"
 version = "1.3.0"
 
+[[deps.Revise]]
+deps = ["CodeTracking", "Distributed", "FileWatching", "JuliaInterpreter", "LibGit2", "LoweredCodeUtils", "OrderedCollections", "Pkg", "REPL", "Requires", "UUIDs", "Unicode"]
+git-tree-sha1 = "c73149ff75d4efb19b6d77411d293ae8fb55c58e"
+uuid = "295af30f-e4ad-537b-8983-00126c2a3abe"
+version = "3.3.4"
+
 [[deps.Rmath]]
 deps = ["Random", "Rmath_jll"]
 git-tree-sha1 = "bf3188feca147ce108c76ad82c2792c57abe7b1f"
@@ -470,6 +845,12 @@ git-tree-sha1 = "aa841c3738cec78b5dbccd56dda332710f35f6a5"
 uuid = "322a6be2-4ae8-5d68-aaf1-3e960788d1d9"
 version = "0.2.0"
 
+[[deps.Scratch]]
+deps = ["Dates"]
+git-tree-sha1 = "f94f779c94e58bf9ea243e77a37e16d9de9126bd"
+uuid = "6c6a2e73-6563-6170-7368-637461726353"
+version = "1.1.1"
+
 [[deps.Serialization]]
 uuid = "9e88b42a-f829-5b0c-bbe9-9e923198166b"
 
@@ -483,6 +864,17 @@ version = "1.0.0"
 deps = ["Distributed", "Mmap", "Random", "Serialization"]
 uuid = "1a1011a3-84de-559e-8e89-a11a2f7dc383"
 
+[[deps.Showoff]]
+deps = ["Dates", "Grisu"]
+git-tree-sha1 = "91eddf657aca81df9ae6ceb20b959ae5653ad1de"
+uuid = "992d4aef-0814-514b-bc4d-f2e9a6c4116f"
+version = "1.0.3"
+
+[[deps.SimpleBufferStream]]
+git-tree-sha1 = "874e8867b33a00e784c8a7e4b60afe9e037b74e1"
+uuid = "777ac1f9-54b0-4bf8-805c-2214025038e7"
+version = "1.1.0"
+
 [[deps.Sockets]]
 uuid = "6462fe0b-24de-5631-8697-dd941f90decc"
 
@@ -504,9 +896,9 @@ version = "2.1.7"
 
 [[deps.StaticArrays]]
 deps = ["LinearAlgebra", "Random", "StaticArraysCore", "Statistics"]
-git-tree-sha1 = "e972716025466461a3dc1588d9168334b71aafff"
+git-tree-sha1 = "23368a3313d12a2326ad0035f0db0c0966f438ef"
 uuid = "90137ffa-7385-5640-81b9-e52037218182"
-version = "1.5.1"
+version = "1.5.2"
 
 [[deps.StaticArraysCore]]
 git-tree-sha1 = "66fe9eb253f910fe8cf161953880cfdaef01cdf0"
@@ -525,9 +917,9 @@ version = "1.4.0"
 
 [[deps.StatsBase]]
 deps = ["DataAPI", "DataStructures", "LinearAlgebra", "LogExpFunctions", "Missings", "Printf", "Random", "SortingAlgorithms", "SparseArrays", "Statistics", "StatsAPI"]
-git-tree-sha1 = "472d044a1c8df2b062b23f222573ad6837a615ba"
+git-tree-sha1 = "0005d75f43ff23688914536c5e9d5ac94f8077f7"
 uuid = "2913bbd2-ae8a-5f71-8c99-4fb6c76f3a91"
-version = "0.33.19"
+version = "0.33.20"
 
 [[deps.StatsFuns]]
 deps = ["ChainRulesCore", "InverseFunctions", "IrrationalConstants", "LogExpFunctions", "Reexport", "Rmath", "SpecialFunctions"]
@@ -535,6 +927,12 @@ git-tree-sha1 = "5950925ff997ed6fb3e985dcce8eb1ba42a0bbe7"
 uuid = "4c63d2b9-4356-54db-8cca-17b64c39e42c"
 version = "0.9.18"
 
+[[deps.StructArrays]]
+deps = ["Adapt", "DataAPI", "StaticArrays", "Tables"]
+git-tree-sha1 = "ec47fb6069c57f1cee2f67541bf8f23415146de7"
+uuid = "09ab397b-f2b6-538f-b94a-2f83cf4a842a"
+version = "0.6.11"
+
 [[deps.SuiteSparse]]
 deps = ["Libdl", "LinearAlgebra", "Serialization", "SparseArrays"]
 uuid = "4607b0f0-06f3-5cda-b6b1-a6196a1729e9"
@@ -543,10 +941,28 @@ uuid = "4607b0f0-06f3-5cda-b6b1-a6196a1729e9"
 deps = ["Dates"]
 uuid = "fa267f1f-6049-4f14-aa54-33bafae1ed76"
 
+[[deps.TableTraits]]
+deps = ["IteratorInterfaceExtensions"]
+git-tree-sha1 = "c06b2f539df1c6efa794486abfb6ed2022561a39"
+uuid = "3783bdb8-4a98-5b6b-af9a-565f29a5fe9c"
+version = "1.0.1"
+
+[[deps.Tables]]
+deps = ["DataAPI", "DataValueInterfaces", "IteratorInterfaceExtensions", "LinearAlgebra", "OrderedCollections", "TableTraits", "Test"]
+git-tree-sha1 = "5ce79ce186cc678bbb5c5681ca3379d1ddae11a1"
+uuid = "bd369af6-aec1-5ad0-b16a-f7cc5008161c"
+version = "1.7.0"
+
 [[deps.Tar]]
 deps = ["ArgTools", "SHA"]
 uuid = "a4e569a6-e804-4fa4-b0f3-eef7a1d5b13e"
 
+[[deps.TensorCore]]
+deps = ["LinearAlgebra"]
+git-tree-sha1 = "1feb45f88d133a655e001435632f019a9a1bcdb6"
+uuid = "62fd8b95-f654-4bbd-a8a5-9c27f68ccd50"
+version = "0.1.1"
+
 [[deps.Test]]
 deps = ["InteractiveUtils", "Logging", "Random", "Serialization"]
 uuid = "8dfed614-e22c-5e08-85e1-65c5234f0b40"
@@ -557,6 +973,11 @@ git-tree-sha1 = "216b95ea110b5972db65aa90f88d8d89dcb8851c"
 uuid = "3bb67fe8-82b1-5028-8e26-92a6c54297fa"
 version = "0.9.6"
 
+[[deps.URIs]]
+git-tree-sha1 = "e59ecc5a41b000fa94423a578d29290c7266fc10"
+uuid = "5c2747f8-b7ea-4ff2-ba2e-563bfd36b1d4"
+version = "1.4.0"
+
 [[deps.UUIDs]]
 deps = ["Random", "SHA"]
 uuid = "cf7118a7-6976-5b1a-9a39-7adc72f591a4"
@@ -569,6 +990,167 @@ version = "1.0.2"
 [[deps.Unicode]]
 uuid = "4ec0a83e-493e-50e2-b9ac-8f72acf5a8f5"
 
+[[deps.UnicodeFun]]
+deps = ["REPL"]
+git-tree-sha1 = "53915e50200959667e78a92a418594b428dffddf"
+uuid = "1cfade01-22cf-5700-b092-accc4b62d6e1"
+version = "0.4.1"
+
+[[deps.Unzip]]
+git-tree-sha1 = "34db80951901073501137bdbc3d5a8e7bbd06670"
+uuid = "41fe7b60-77ed-43a1-b4f0-825fd5a5650d"
+version = "0.1.2"
+
+[[deps.Wayland_jll]]
+deps = ["Artifacts", "Expat_jll", "JLLWrappers", "Libdl", "Libffi_jll", "Pkg", "XML2_jll"]
+git-tree-sha1 = "3e61f0b86f90dacb0bc0e73a0c5a83f6a8636e23"
+uuid = "a2964d1f-97da-50d4-b82a-358c7fce9d89"
+version = "1.19.0+0"
+
+[[deps.Wayland_protocols_jll]]
+deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
+git-tree-sha1 = "4528479aa01ee1b3b4cd0e6faef0e04cf16466da"
+uuid = "2381bf8a-dfd0-557d-9999-79630e7b1b91"
+version = "1.25.0+0"
+
+[[deps.XML2_jll]]
+deps = ["Artifacts", "JLLWrappers", "Libdl", "Libiconv_jll", "Pkg", "Zlib_jll"]
+git-tree-sha1 = "58443b63fb7e465a8a7210828c91c08b92132dff"
+uuid = "02c8fc9c-b97f-50b9-bbe4-9be30ff0a78a"
+version = "2.9.14+0"
+
+[[deps.XSLT_jll]]
+deps = ["Artifacts", "JLLWrappers", "Libdl", "Libgcrypt_jll", "Libgpg_error_jll", "Libiconv_jll", "Pkg", "XML2_jll", "Zlib_jll"]
+git-tree-sha1 = "91844873c4085240b95e795f692c4cec4d805f8a"
+uuid = "aed1982a-8fda-507f-9586-7b0439959a61"
+version = "1.1.34+0"
+
+[[deps.Xorg_libX11_jll]]
+deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg", "Xorg_libxcb_jll", "Xorg_xtrans_jll"]
+git-tree-sha1 = "5be649d550f3f4b95308bf0183b82e2582876527"
+uuid = "4f6342f7-b3d2-589e-9d20-edeb45f2b2bc"
+version = "1.6.9+4"
+
+[[deps.Xorg_libXau_jll]]
+deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
+git-tree-sha1 = "4e490d5c960c314f33885790ed410ff3a94ce67e"
+uuid = "0c0b7dd1-d40b-584c-a123-a41640f87eec"
+version = "1.0.9+4"
+
+[[deps.Xorg_libXcursor_jll]]
+deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg", "Xorg_libXfixes_jll", "Xorg_libXrender_jll"]
+git-tree-sha1 = "12e0eb3bc634fa2080c1c37fccf56f7c22989afd"
+uuid = "935fb764-8cf2-53bf-bb30-45bb1f8bf724"
+version = "1.2.0+4"
+
+[[deps.Xorg_libXdmcp_jll]]
+deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
+git-tree-sha1 = "4fe47bd2247248125c428978740e18a681372dd4"
+uuid = "a3789734-cfe1-5b06-b2d0-1dd0d9d62d05"
+version = "1.1.3+4"
+
+[[deps.Xorg_libXext_jll]]
+deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg", "Xorg_libX11_jll"]
+git-tree-sha1 = "b7c0aa8c376b31e4852b360222848637f481f8c3"
+uuid = "1082639a-0dae-5f34-9b06-72781eeb8cb3"
+version = "1.3.4+4"
+
+[[deps.Xorg_libXfixes_jll]]
+deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg", "Xorg_libX11_jll"]
+git-tree-sha1 = "0e0dc7431e7a0587559f9294aeec269471c991a4"
+uuid = "d091e8ba-531a-589c-9de9-94069b037ed8"
+version = "5.0.3+4"
+
+[[deps.Xorg_libXi_jll]]
+deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg", "Xorg_libXext_jll", "Xorg_libXfixes_jll"]
+git-tree-sha1 = "89b52bc2160aadc84d707093930ef0bffa641246"
+uuid = "a51aa0fd-4e3c-5386-b890-e753decda492"
+version = "1.7.10+4"
+
+[[deps.Xorg_libXinerama_jll]]
+deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg", "Xorg_libXext_jll"]
+git-tree-sha1 = "26be8b1c342929259317d8b9f7b53bf2bb73b123"
+uuid = "d1454406-59df-5ea1-beac-c340f2130bc3"
+version = "1.1.4+4"
+
+[[deps.Xorg_libXrandr_jll]]
+deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg", "Xorg_libXext_jll", "Xorg_libXrender_jll"]
+git-tree-sha1 = "34cea83cb726fb58f325887bf0612c6b3fb17631"
+uuid = "ec84b674-ba8e-5d96-8ba1-2a689ba10484"
+version = "1.5.2+4"
+
+[[deps.Xorg_libXrender_jll]]
+deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg", "Xorg_libX11_jll"]
+git-tree-sha1 = "19560f30fd49f4d4efbe7002a1037f8c43d43b96"
+uuid = "ea2f1a96-1ddc-540d-b46f-429655e07cfa"
+version = "0.9.10+4"
+
+[[deps.Xorg_libpthread_stubs_jll]]
+deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
+git-tree-sha1 = "6783737e45d3c59a4a4c4091f5f88cdcf0908cbb"
+uuid = "14d82f49-176c-5ed1-bb49-ad3f5cbd8c74"
+version = "0.1.0+3"
+
+[[deps.Xorg_libxcb_jll]]
+deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg", "XSLT_jll", "Xorg_libXau_jll", "Xorg_libXdmcp_jll", "Xorg_libpthread_stubs_jll"]
+git-tree-sha1 = "daf17f441228e7a3833846cd048892861cff16d6"
+uuid = "c7cfdc94-dc32-55de-ac96-5a1b8d977c5b"
+version = "1.13.0+3"
+
+[[deps.Xorg_libxkbfile_jll]]
+deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg", "Xorg_libX11_jll"]
+git-tree-sha1 = "926af861744212db0eb001d9e40b5d16292080b2"
+uuid = "cc61e674-0454-545c-8b26-ed2c68acab7a"
+version = "1.1.0+4"
+
+[[deps.Xorg_xcb_util_image_jll]]
+deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg", "Xorg_xcb_util_jll"]
+git-tree-sha1 = "0fab0a40349ba1cba2c1da699243396ff8e94b97"
+uuid = "12413925-8142-5f55-bb0e-6d7ca50bb09b"
+version = "0.4.0+1"
+
+[[deps.Xorg_xcb_util_jll]]
+deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg", "Xorg_libxcb_jll"]
+git-tree-sha1 = "e7fd7b2881fa2eaa72717420894d3938177862d1"
+uuid = "2def613f-5ad1-5310-b15b-b15d46f528f5"
+version = "0.4.0+1"
+
+[[deps.Xorg_xcb_util_keysyms_jll]]
+deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg", "Xorg_xcb_util_jll"]
+git-tree-sha1 = "d1151e2c45a544f32441a567d1690e701ec89b00"
+uuid = "975044d2-76e6-5fbe-bf08-97ce7c6574c7"
+version = "0.4.0+1"
+
+[[deps.Xorg_xcb_util_renderutil_jll]]
+deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg", "Xorg_xcb_util_jll"]
+git-tree-sha1 = "dfd7a8f38d4613b6a575253b3174dd991ca6183e"
+uuid = "0d47668e-0667-5a69-a72c-f761630bfb7e"
+version = "0.3.9+1"
+
+[[deps.Xorg_xcb_util_wm_jll]]
+deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg", "Xorg_xcb_util_jll"]
+git-tree-sha1 = "e78d10aab01a4a154142c5006ed44fd9e8e31b67"
+uuid = "c22f9ab0-d5fe-5066-847c-f4bb1cd4e361"
+version = "0.4.1+1"
+
+[[deps.Xorg_xkbcomp_jll]]
+deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg", "Xorg_libxkbfile_jll"]
+git-tree-sha1 = "4bcbf660f6c2e714f87e960a171b119d06ee163b"
+uuid = "35661453-b289-5fab-8a00-3d9160c6a3a4"
+version = "1.4.2+4"
+
+[[deps.Xorg_xkeyboard_config_jll]]
+deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg", "Xorg_xkbcomp_jll"]
+git-tree-sha1 = "5c8424f8a67c3f2209646d4425f3d415fee5931d"
+uuid = "33bec58e-1273-512f-9401-5d533626f822"
+version = "2.27.0+4"
+
+[[deps.Xorg_xtrans_jll]]
+deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
+git-tree-sha1 = "79c31e7844f6ecf779705fbc12146eb190b7d845"
+uuid = "c5fb5394-a638-5e4d-96e5-b29de1b5cf10"
+version = "1.4.0+3"
+
 [[deps.Zlib_jll]]
 deps = ["Libdl"]
 uuid = "83775a58-1f1d-513f-b197-d71354ab007a"
@@ -579,16 +1161,46 @@ git-tree-sha1 = "e45044cd873ded54b6a5bac0eb5c971392cf1927"
 uuid = "3161d3a3-bdf6-5164-811a-617609db77b4"
 version = "1.5.2+0"
 
+[[deps.libaom_jll]]
+deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
+git-tree-sha1 = "3a2ea60308f0996d26f1e5354e10c24e9ef905d4"
+uuid = "a4ae2306-e953-59d6-aa16-d00cac43593b"
+version = "3.4.0+0"
+
+[[deps.libass_jll]]
+deps = ["Artifacts", "Bzip2_jll", "FreeType2_jll", "FriBidi_jll", "HarfBuzz_jll", "JLLWrappers", "Libdl", "Pkg", "Zlib_jll"]
+git-tree-sha1 = "5982a94fcba20f02f42ace44b9894ee2b140fe47"
+uuid = "0ac62f75-1d6f-5e53-bd7c-93b484bb37c0"
+version = "0.15.1+0"
+
 [[deps.libblastrampoline_jll]]
 deps = ["Artifacts", "Libdl", "OpenBLAS_jll"]
 uuid = "8e850b90-86db-534c-a0d3-1478176c7d93"
 
+[[deps.libfdk_aac_jll]]
+deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
+git-tree-sha1 = "daacc84a041563f965be61859a36e17c4e4fcd55"
+uuid = "f638f0a6-7fb0-5443-88ba-1cc74229b280"
+version = "2.0.2+0"
+
+[[deps.libpng_jll]]
+deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg", "Zlib_jll"]
+git-tree-sha1 = "94d180a6d2b5e55e447e2d27a29ed04fe79eb30c"
+uuid = "b53b4c65-9356-5827-b1ea-8c7a1a84506f"
+version = "1.6.38+0"
+
 [[deps.libsass_jll]]
 deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
 git-tree-sha1 = "941afb93587dcec07f89e511057f5efc0bec6f0d"
 uuid = "47bcb7c8-5119-555a-9eeb-0afcc36cd728"
 version = "3.6.4+0"
 
+[[deps.libvorbis_jll]]
+deps = ["Artifacts", "JLLWrappers", "Libdl", "Ogg_jll", "Pkg"]
+git-tree-sha1 = "b910cb81ef3fe6e78bf6acee440bda86fd6ae00c"
+uuid = "f27f6e37-5d2b-51aa-960f-b287f2bc3b7a"
+version = "1.3.7+1"
+
 [[deps.nghttp2_jll]]
 deps = ["Artifacts", "Libdl"]
 uuid = "8e850ede-7688-5339-a07c-302acd2aaf8d"
@@ -596,3 +1208,21 @@ uuid = "8e850ede-7688-5339-a07c-302acd2aaf8d"
 [[deps.p7zip_jll]]
 deps = ["Artifacts", "Libdl"]
 uuid = "3f19e933-33d8-53b3-aaab-bd5110c3b7a0"
+
+[[deps.x264_jll]]
+deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
+git-tree-sha1 = "4fea590b89e6ec504593146bf8b988b2c00922b2"
+uuid = "1270edf5-f2f9-52d2-97e9-ab00b5d0237a"
+version = "2021.5.5+0"
+
+[[deps.x265_jll]]
+deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
+git-tree-sha1 = "ee567a171cce03570d77ad3a43e90218e38937a9"
+uuid = "dfaa095f-4041-5dcd-9319-2fabd8486b76"
+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 = "ece2350174195bb31de1a63bea3a41ae1aa593b6"
+uuid = "d8fb68d0-12a3-5cfd-a85a-d49703b185fd"
+version = "0.9.1+5"
diff --git a/Project.toml b/Project.toml
old mode 100644
new mode 100755
index 30c628df..9957804e
--- a/Project.toml
+++ b/Project.toml
@@ -9,6 +9,7 @@ Distributed = "8ba89e20-285c-5b6f-9357-94700520ee1b"
 Distributions = "31c24e10-a181-5473-b8eb-7969acd0382f"
 Documenter = "e30172f5-a6a5-5a46-863b-614d45cd2de4"
 DocumenterTools = "35a29f4d-8980-5a13-9543-d66fff28ecb8"
+FileIO = "5789e2e9-d7fb-5bc7-8068-2c6fae9b9549"
 ForwardDiff = "f6369f11-7733-5829-9624-2563aa707210"
 HDF5 = "f67ccb44-e63f-5c2f-98bd-6dc0ccc4ba2f"
 JLD = "4138dd39-2aa7-5051-a626-17a0bb65d9c8"
@@ -17,9 +18,12 @@ LineSearches = "d3d80556-e9d4-5f37-9878-2ab0fcc64255"
 LinearAlgebra = "37e2e46d-f89d-539d-b4ee-838fcccc9c8e"
 LsqFit = "2fda8390-95c7-5789-9bda-21331edee243"
 Optim = "429524aa-4258-5aef-a3af-852621145aeb"
+Plots = "91a5bcdd-55d7-5caf-9e0b-520d859cae80"
 Random = "9a3f8284-a2c9-5f02-9a11-845980a1fd5c"
+Revise = "295af30f-e4ad-537b-8983-00126c2a3abe"
 SparseArrays = "2f01184e-e22b-5df5-ae63-d93ebab69eaf"
 Statistics = "10745b16-79ce-11e8-11f9-7d13ad32a3b2"
+StatsBase = "2913bbd2-ae8a-5f71-8c99-4fb6c76f3a91"
 StatsFuns = "4c63d2b9-4356-54db-8cca-17b64c39e42c"
 
 [compat]
diff --git a/README.md b/README.md
old mode 100644
new mode 100755
diff --git a/assets/dabenchmarks.png b/assets/dabenchmarks.png
old mode 100644
new mode 100755
diff --git a/docs/.gitignore b/docs/.gitignore
old mode 100644
new mode 100755
diff --git a/docs/Project.toml b/docs/Project.toml
old mode 100644
new mode 100755
diff --git a/docs/make.jl b/docs/make.jl
old mode 100644
new mode 100755
diff --git a/docs/src/home/DataAssimilationBenchmarks.md b/docs/src/home/DataAssimilationBenchmarks.md
old mode 100644
new mode 100755
diff --git a/docs/src/home/Getting Started.md b/docs/src/home/Getting Started.md
old mode 100644
new mode 100755
diff --git a/docs/src/home/Introduction.md b/docs/src/home/Introduction.md
old mode 100644
new mode 100755
diff --git a/docs/src/index.md b/docs/src/index.md
old mode 100644
new mode 100755
diff --git a/docs/src/submodules/analysis/PlotExperimentData.md b/docs/src/submodules/analysis/PlotExperimentData.md
old mode 100644
new mode 100755
diff --git a/docs/src/submodules/analysis/ProcessExperimentData.md b/docs/src/submodules/analysis/ProcessExperimentData.md
old mode 100644
new mode 100755
diff --git a/docs/src/submodules/experiments/D3VARExps.md b/docs/src/submodules/experiments/D3VARExps.md
new file mode 100755
index 00000000..b6465a2d
--- /dev/null
+++ b/docs/src/submodules/experiments/D3VARExps.md
@@ -0,0 +1,22 @@
+# D3VARExps
+
+This module defines methods for experiments with classical variational data assimilation with
+3D-VAR.  Primal cost functions are defined, with their implicit differentiation
+performed with automatic differentiation with [JuliaDiff](https://github.com/JuliaDiff)
+methods. Development of gradient-based optimization schemes using automatic
+differentiation is ongoing, with future development planned to integrate variational
+benchmark experiments.
+
+The basic 3D-VAR cost function API is defined as follows
+```{julia}
+    D3_var_cost(x::VecA(T), obs::VecA(T), x_background::VecA(T), state_cov::CovM(T),
+    obs_cov::CovM(T), kwargs::StepKwargs) where T <: Real
+```
+where the control variable `x` is optimized, with fixed hyper-parameters defined in a
+wrapping function passed to auto-differentiation.
+
+## Methods
+
+```@autodocs
+Modules = [DataAssimilationBenchmarks.D3VARExps]
+```
\ No newline at end of file
diff --git a/docs/src/submodules/experiments/FilterExps.md b/docs/src/submodules/experiments/FilterExps.md
old mode 100644
new mode 100755
diff --git a/docs/src/submodules/experiments/GenerateTimeSeries.md b/docs/src/submodules/experiments/GenerateTimeSeries.md
old mode 100644
new mode 100755
diff --git a/docs/src/submodules/experiments/ParallelExperimentDriver.md b/docs/src/submodules/experiments/ParallelExperimentDriver.md
old mode 100644
new mode 100755
diff --git a/docs/src/submodules/experiments/SingleExperimentDriver.md b/docs/src/submodules/experiments/SingleExperimentDriver.md
old mode 100644
new mode 100755
diff --git a/docs/src/submodules/experiments/Slurm.md b/docs/src/submodules/experiments/Slurm.md
old mode 100644
new mode 100755
diff --git a/docs/src/submodules/experiments/SmootherExps.md b/docs/src/submodules/experiments/SmootherExps.md
old mode 100644
new mode 100755
diff --git a/docs/src/submodules/experiments/VarAnalysisExperimentDriver.md b/docs/src/submodules/experiments/VarAnalysisExperimentDriver.md
new file mode 100755
index 00000000..5b721e80
--- /dev/null
+++ b/docs/src/submodules/experiments/VarAnalysisExperimentDriver.md
@@ -0,0 +1,8 @@
+# VarAnalysisExperimentDriver
+
+
+## Methods
+
+```@autodocs
+Modules = [DataAssimilationBenchmarks.VarAnalysisExperimentDriver]
+```
diff --git a/docs/src/submodules/methods/DeSolvers.md b/docs/src/submodules/methods/DeSolvers.md
old mode 100644
new mode 100755
diff --git a/docs/src/submodules/methods/EnsembleKalmanSchemes.md b/docs/src/submodules/methods/EnsembleKalmanSchemes.md
old mode 100644
new mode 100755
diff --git a/docs/src/submodules/methods/XdVAR.md b/docs/src/submodules/methods/XdVAR.md
old mode 100644
new mode 100755
diff --git a/docs/src/submodules/models/IEEE39bus.md b/docs/src/submodules/models/IEEE39bus.md
old mode 100644
new mode 100755
diff --git a/docs/src/submodules/models/L96.md b/docs/src/submodules/models/L96.md
old mode 100644
new mode 100755
diff --git a/docs/src/submodules/models/ObsOperators.md b/docs/src/submodules/models/ObsOperators.md
old mode 100644
new mode 100755
diff --git a/paper.bib b/paper.bib
old mode 100644
new mode 100755
diff --git a/paper.md b/paper.md
old mode 100644
new mode 100755
diff --git a/src/DataAssimilationBenchmarks.jl b/src/DataAssimilationBenchmarks.jl
old mode 100644
new mode 100755
diff --git a/src/analysis/ProcessExperimentData.jl b/src/analysis/ProcessExperimentData.jl
old mode 100644
new mode 100755
diff --git a/src/analysis/param_filter/plt_param_rmse_spread_4_pane.py b/src/analysis/param_filter/plt_param_rmse_spread_4_pane.py
old mode 100644
new mode 100755
diff --git a/src/analysis/param_smoother/lorenz-96/plt_param_rmse_spread_4_pane.py b/src/analysis/param_smoother/lorenz-96/plt_param_rmse_spread_4_pane.py
old mode 100644
new mode 100755
diff --git a/src/analysis/param_smoother/plt_smoother_rmse_spread_all_stats_4_pane.py b/src/analysis/param_smoother/plt_smoother_rmse_spread_all_stats_4_pane.py
old mode 100644
new mode 100755
diff --git a/src/analysis/state_filter/plt_filter_rmse_spread_all_stats.py b/src/analysis/state_filter/plt_filter_rmse_spread_all_stats.py
old mode 100644
new mode 100755
diff --git a/src/analysis/state_filter/plt_filter_rmse_spread_v_ens.py b/src/analysis/state_filter/plt_filter_rmse_spread_v_ens.py
old mode 100644
new mode 100755
diff --git a/src/analysis/state_filter/plt_state_rmse_spread_2_pane.py b/src/analysis/state_filter/plt_state_rmse_spread_2_pane.py
old mode 100644
new mode 100755
diff --git a/src/analysis/state_smoother/heat_plots/all_ens_mda_tuned.pdf b/src/analysis/state_smoother/heat_plots/all_ens_mda_tuned.pdf
old mode 100644
new mode 100755
diff --git a/src/analysis/state_smoother/heat_plots/all_ens_sda_adaptive.pdf b/src/analysis/state_smoother/heat_plots/all_ens_sda_adaptive.pdf
old mode 100644
new mode 100755
diff --git a/src/analysis/state_smoother/heat_plots/all_ens_sda_tuned.pdf b/src/analysis/state_smoother/heat_plots/all_ens_sda_tuned.pdf
old mode 100644
new mode 100755
diff --git a/src/analysis/state_smoother/heat_plots/iteration_statistics_tanl_005_shift_1.pdf b/src/analysis/state_smoother/heat_plots/iteration_statistics_tanl_005_shift_1.pdf
old mode 100644
new mode 100755
diff --git a/src/analysis/state_smoother/heat_plots/versus_obs_mda_tuned.pdf b/src/analysis/state_smoother/heat_plots/versus_obs_mda_tuned.pdf
old mode 100644
new mode 100755
diff --git a/src/analysis/state_smoother/heat_plots/versus_obs_sda_tuned.pdf b/src/analysis/state_smoother/heat_plots/versus_obs_sda_tuned.pdf
old mode 100644
new mode 100755
diff --git a/src/analysis/state_smoother/heat_plots/versus_shift_mda_tuned.pdf b/src/analysis/state_smoother/heat_plots/versus_shift_mda_tuned.pdf
old mode 100644
new mode 100755
diff --git a/src/analysis/state_smoother/heat_plots/versus_shift_sda_adaptive.pdf b/src/analysis/state_smoother/heat_plots/versus_shift_sda_adaptive.pdf
old mode 100644
new mode 100755
diff --git a/src/analysis/state_smoother/heat_plots/versus_shift_sda_tuned.pdf b/src/analysis/state_smoother/heat_plots/versus_shift_sda_tuned.pdf
old mode 100644
new mode 100755
diff --git a/src/analysis/state_smoother/heat_plots/versus_tanl_mda_tuned.pdf b/src/analysis/state_smoother/heat_plots/versus_tanl_mda_tuned.pdf
old mode 100644
new mode 100755
diff --git a/src/analysis/state_smoother/heat_plots/versus_tanl_sda_adaptive.pdf b/src/analysis/state_smoother/heat_plots/versus_tanl_sda_adaptive.pdf
old mode 100644
new mode 100755
diff --git a/src/analysis/state_smoother/heat_plots/versus_tanl_sda_tuned.pdf b/src/analysis/state_smoother/heat_plots/versus_tanl_sda_tuned.pdf
old mode 100644
new mode 100755
diff --git a/src/analysis/state_smoother/line_plots/rmse_spread_optimium_forecast_lag_v_nens_tanl_005_mda_tuned.pdf b/src/analysis/state_smoother/line_plots/rmse_spread_optimium_forecast_lag_v_nens_tanl_005_mda_tuned.pdf
old mode 100644
new mode 100755
diff --git a/src/analysis/state_smoother/line_plots/rmse_spread_optimium_smoother_lag_v_nens_tanl_005_mda_tuned.pdf b/src/analysis/state_smoother/line_plots/rmse_spread_optimium_smoother_lag_v_nens_tanl_005_mda_tuned.pdf
old mode 100644
new mode 100755
diff --git a/src/analysis/state_smoother/line_plots/rmse_spread_v_gamma_mda_N_21.pdf b/src/analysis/state_smoother/line_plots/rmse_spread_v_gamma_mda_N_21.pdf
old mode 100644
new mode 100755
diff --git a/src/analysis/state_smoother/line_plots/rmse_spread_v_gamma_sda_N_21.pdf b/src/analysis/state_smoother/line_plots/rmse_spread_v_gamma_sda_N_21.pdf
old mode 100644
new mode 100755
diff --git a/src/analysis/state_smoother/line_plots/rmse_spread_v_lag_nens_21_tanl_005_mda_tuned.pdf b/src/analysis/state_smoother/line_plots/rmse_spread_v_lag_nens_21_tanl_005_mda_tuned.pdf
old mode 100644
new mode 100755
diff --git a/src/analysis/state_smoother/line_plots/rmse_spread_v_lag_nens_21_tanl_005_sda_adaptive.pdf b/src/analysis/state_smoother/line_plots/rmse_spread_v_lag_nens_21_tanl_005_sda_adaptive.pdf
old mode 100644
new mode 100755
diff --git a/src/analysis/state_smoother/line_plots/rmse_spread_v_lag_nens_21_tanl_005_sda_tuned.pdf b/src/analysis/state_smoother/line_plots/rmse_spread_v_lag_nens_21_tanl_005_sda_tuned.pdf
old mode 100644
new mode 100755
diff --git a/src/analysis/state_smoother/line_plots/rmse_spread_v_shift_nens_21_tanl_005_mda_tuned.pdf b/src/analysis/state_smoother/line_plots/rmse_spread_v_shift_nens_21_tanl_005_mda_tuned.pdf
old mode 100644
new mode 100755
diff --git a/src/analysis/state_smoother/line_plots/rmse_spread_v_shift_nens_21_tanl_005_sda_tuned.pdf b/src/analysis/state_smoother/line_plots/rmse_spread_v_shift_nens_21_tanl_005_sda_tuned.pdf
old mode 100644
new mode 100755
diff --git a/src/analysis/state_smoother/plt_number_of_iterations_6_pane.py b/src/analysis/state_smoother/plt_number_of_iterations_6_pane.py
old mode 100644
new mode 100755
diff --git a/src/analysis/state_smoother/plt_optimal_inflation_heat_plot.py b/src/analysis/state_smoother/plt_optimal_inflation_heat_plot.py
old mode 100644
new mode 100755
diff --git a/src/analysis/state_smoother/plt_smoother_rmse_spread_all_stats.py b/src/analysis/state_smoother/plt_smoother_rmse_spread_all_stats.py
old mode 100644
new mode 100755
diff --git a/src/analysis/state_smoother/plt_smoother_rmse_spread_all_stats_4_pane.py b/src/analysis/state_smoother/plt_smoother_rmse_spread_all_stats_4_pane.py
old mode 100644
new mode 100755
diff --git a/src/analysis/state_smoother/plt_smoother_rmse_spread_all_stats_4_pane_verus_obs.py b/src/analysis/state_smoother/plt_smoother_rmse_spread_all_stats_4_pane_verus_obs.py
old mode 100644
new mode 100755
diff --git a/src/analysis/state_smoother/plt_smoother_rmse_spread_all_stats_4_pane_verus_shift.py b/src/analysis/state_smoother/plt_smoother_rmse_spread_all_stats_4_pane_verus_shift.py
old mode 100644
new mode 100755
diff --git a/src/analysis/state_smoother/plt_smoother_rmse_spread_all_stats_4_pane_verus_tanl.py b/src/analysis/state_smoother/plt_smoother_rmse_spread_all_stats_4_pane_verus_tanl.py
old mode 100644
new mode 100755
diff --git a/src/analysis/state_smoother/plt_smoother_rmse_spread_fixed_lag_v_shift.py b/src/analysis/state_smoother/plt_smoother_rmse_spread_fixed_lag_v_shift.py
old mode 100644
new mode 100755
diff --git a/src/analysis/state_smoother/plt_smoother_rmse_spread_obs_v_gamma.py b/src/analysis/state_smoother/plt_smoother_rmse_spread_obs_v_gamma.py
old mode 100644
new mode 100755
diff --git a/src/analysis/state_smoother/plt_smoother_rmse_spread_obs_v_lag.py b/src/analysis/state_smoother/plt_smoother_rmse_spread_obs_v_lag.py
old mode 100644
new mode 100755
diff --git a/src/analysis/state_smoother/plt_smoother_rmse_spread_optimum_forecast_lag_versus_nens.py b/src/analysis/state_smoother/plt_smoother_rmse_spread_optimum_forecast_lag_versus_nens.py
old mode 100644
new mode 100755
diff --git a/src/analysis/state_smoother/plt_smoother_rmse_spread_optimum_lag_versus_nens.py b/src/analysis/state_smoother/plt_smoother_rmse_spread_optimum_lag_versus_nens.py
old mode 100644
new mode 100755
diff --git a/src/analysis/state_smoother/plt_smoother_rmse_spread_optimum_lag_versus_tanl.py b/src/analysis/state_smoother/plt_smoother_rmse_spread_optimum_lag_versus_tanl.py
old mode 100644
new mode 100755
diff --git a/src/analysis/state_smoother/plt_smoother_rmse_spread_optimum_smoother_lag_versus_nens.py b/src/analysis/state_smoother/plt_smoother_rmse_spread_optimum_smoother_lag_versus_nens.py
old mode 100644
new mode 100755
diff --git a/src/analysis/state_smoother/plt_smoother_rmse_spread_v_lag.py b/src/analysis/state_smoother/plt_smoother_rmse_spread_v_lag.py
old mode 100644
new mode 100755
diff --git a/src/analysis/state_smoother/plt_smoother_rmse_spread_v_lag_mda_v_sda.py b/src/analysis/state_smoother/plt_smoother_rmse_spread_v_lag_mda_v_sda.py
old mode 100644
new mode 100755
diff --git a/src/analysis/state_smoother/plt_smoother_rmse_spread_v_shift.py b/src/analysis/state_smoother/plt_smoother_rmse_spread_v_shift.py
old mode 100644
new mode 100755
diff --git a/src/analysis/state_smoother/plt_smoother_state_rmse_spread_2_pane.py b/src/analysis/state_smoother/plt_smoother_state_rmse_spread_2_pane.py
old mode 100644
new mode 100755
diff --git a/src/analysis/state_smoother/plt_smoother_state_rmse_spread_6_pane.py b/src/analysis/state_smoother/plt_smoother_state_rmse_spread_6_pane.py
old mode 100644
new mode 100755
diff --git a/src/analysis/state_smoother/plt_tuned_inflation_3_pane.py b/src/analysis/state_smoother/plt_tuned_inflation_3_pane.py
old mode 100644
new mode 100755
diff --git a/src/data/Rename.jl b/src/data/Rename.jl
old mode 100644
new mode 100755
diff --git a/src/data/rm_data.py b/src/data/rm_data.py
old mode 100644
new mode 100755
diff --git a/src/experiments/D3VARExps.jl b/src/experiments/D3VARExps.jl
new file mode 100755
index 00000000..d2bc5a95
--- /dev/null
+++ b/src/experiments/D3VARExps.jl
@@ -0,0 +1,315 @@
+##############################################################################################
+module D3VARExps
+##############################################################################################
+# imports and exports
+using Random, Distributions, LinearAlgebra, StatsBase, Statistics, Measures
+using JLD2, HDF5, Plots
+using ..DataAssimilationBenchmarks, ..ObsOperators, ..DeSolvers, ..XdVAR
+##############################################################################################
+# Main 3DVAR experiments
+##############################################################################################
+"""
+    D3_var_filter_analysis_simple()
+"""
+function D3_var_filter_analysis_simple()
+    # time the experiment
+    t1 = time()
+
+    # Define experiment parameters
+    # number of cycles in experiment
+    nanl = 40
+    diffusion = 0.0
+    tanl = 0.05
+    γ = [8.0]
+
+    # define the observation operator HARD-CODED in this line
+    H_obs = alternating_obs_operator
+
+    # set the integration step size for the ensemble at 0.01 - we are assuming SDE
+    h = 0.01
+
+    # define derivative parameter
+    dx_params = Dict{String, Vector{Float64}}("F" => [8.0])
+
+    # define the dynamical model derivative for this experiment - we are assuming 
+    # Lorenz-96 model
+    dx_dt = L96.dx_dt
+
+    # define integration method
+    step_model! = rk4_step!
+
+    # number of discrete forecast steps
+    f_steps = convert(Int64, tanl / h)
+
+    # set seed
+    seed = 234
+    Random.seed!(seed)
+
+    # define the initialization
+    # observation noise
+    v = rand(Normal(0, 1), 40)
+    # define the initial observation range and truth reference solution
+    x_b = zeros(40)
+    x_t = x_b + v
+    
+    # define kwargs for the analysis method
+    # and the underlying dynamical model
+    kwargs = Dict{String,Any}(
+                              "dx_dt" => dx_dt,
+                              "f_steps" => f_steps,
+                              "step_model" => step_model!,
+                              "dx_params" => dx_params,
+                              "h" => h,
+                              "diffusion" => diffusion,
+                              "gamma" => γ,
+                             )
+
+    # create storage for the forecast and analysis statistics
+    fore_rmse = Vector{Float64}(undef, nanl)
+    filt_rmse = Vector{Float64}(undef, nanl)
+    
+    for i in 1:nanl
+        #print("Iteration: ")
+        #display(i)
+        for j in 1:f_steps
+            # M(x^b)
+            step_model!(x_b, 0.0, kwargs)
+            # M(x^t)
+            step_model!(x_t, 0.0, kwargs)
+        end
+
+        # multivariate - rand(MvNormal(zeros(40), I))
+        w = rand(Normal(0, 1), 40)
+        obs = x_t + w
+
+        state_cov = I
+        obs_cov = I
+
+        # optimized cost function input and value
+        x_opt = XdVAR.D3_var_NewtonOp(x_b, obs, x_b, state_cov, H_obs, obs_cov, kwargs)
+
+        # compare model forecast and truth twin via RMSE
+        rmse_forecast = sqrt(msd(x_b, x_t))
+        fore_rmse[i] = rmse_forecast
+
+        # compare optimal forecast and truth twin via RMSE
+        rmse_filter = sqrt(msd(x_opt, x_t))
+        filt_rmse[i] = rmse_filter
+
+        # reinitializing x_b and x_t for next cycle
+        x_b = x_opt
+    end
+
+    data = Dict{String,Any}(
+                            "seed" => seed,
+                            "diffusion" => diffusion,
+                            "dx_params" => dx_params,
+                            "gamma" => γ,
+                            "tanl" => tanl,
+                            "nanl" => nanl,
+                            "h" =>  h,
+                            "fore_rmse" => fore_rmse,
+                            "filt_rmse" => filt_rmse
+                           )
+
+    path = pkgdir(DataAssimilationBenchmarks) * "/src/data/time_series/"
+    name = "L96_3DVAR_time_series_seed_" * lpad(seed, 4, "0") *
+           "_diff_" * rpad(diffusion, 5, "0") *
+           #"_F_" * lpad(, 4, "0") *
+           "_tanl_" * rpad(tanl, 4, "0") *
+           "_nanl_" * lpad(nanl, 5, "0") *
+           "_h_" * rpad(h, 5, "0") *
+           ".jld2"
+
+    save(path * name, data)
+
+    # output time
+    print("Runtime " * string(round((time() - t1)  / 60.0, digits=4))  * " minutes\n")
+
+    # make plot
+    path = pkgdir(DataAssimilationBenchmarks) * "/src/analysis/var_exp/"
+    t = 1:nanl
+    plot(t, fore_rmse, marker=(:circle,5), label = "Forecast", 
+        title="Update: Root-Mean-Square Error vs. Time", 
+        legend_position = :outertopright, 
+        margin=15mm, size=(800,500), dpi = 600)
+    plot!(t, filt_rmse, marker=(:circle,5), label = "Filter")
+    xlabel!("Time [Cycles]")
+    ylabel!("Root-Mean-Square Error [RMSE]")
+    savefig(path * "I_Update_SIMPLE")
+end
+
+
+##############################################################################################
+"""
+function D3_var_filter_analysis((time_series, γ, is_informed, tuning_factor, is_updated)::NamedTuple{
+    (:time_series,:γ,:is_informed,:tuning_factor,:is_updated),<:Tuple{String,
+        Float64,Bool,Float64,Bool}})
+    Plotting capabilities are commented out for parallel experiment.
+"""
+function D3_var_filter_analysis((time_series, γ, is_informed, tuning_factor, is_updated)::NamedTuple{
+                        (:time_series,:γ,:is_informed,:tuning_factor,:is_updated),<:Tuple{String,
+                            Float64,Bool,Float64,Bool}})
+    
+    # time the experiment
+    t1 = time()
+
+    # load the path, timeseries, and associated parameters
+    path = pkgdir(DataAssimilationBenchmarks) * "/src/data/"
+    ts = load(path * time_series)::Dict{String,Any}
+    diffusion = ts["diffusion"]::Float64
+    dx_params = ts["dx_params"]::ParamDict(Float64)
+    tanl = ts["tanl"]::Float64
+    nanl = ts["nanl"]::Int64
+    # set the integration step size for the ensemble 
+    h = ts["h"]::Float64
+
+    # define the observation operator HARD-CODED in this line
+    H_obs = alternating_obs_operator
+
+    # define the dynamical model derivative for this experiment - we are assuming 
+    # Lorenz-96 model
+    dx_dt = L96.dx_dt
+
+    # define integration method
+    step_model! = rk4_step!
+
+    # number of discrete forecast steps
+    f_steps = convert(Int64, tanl / h)
+
+    # set seed
+    seed = ts["seed"]::Int64
+    Random.seed!(seed)
+
+    # define the initialization
+    o = ts["obs"]::Array{Float64, 2}
+    obs_un = 1
+    obs_cov = obs_un^2.0 * I
+
+    # define state covaraince based on input
+    if is_informed == true
+        c = cov(o, dims = 2)
+        state_cov = tuning_factor*Symmetric(c)
+    else
+        state_cov = tuning_factor*I
+    end
+
+    x_t = o[:,1]
+    # observation noise
+    v = rand(MvNormal(zeros(40), I))
+    # define the initial background state
+    x_b = x_t + v;
+    
+    # define kwargs for the analysis method
+    # and the underlying dynamical model
+    kwargs = Dict{String,Any}(
+                              "dx_dt" => dx_dt,
+                              "f_steps" => f_steps,
+                              "step_model" => step_model!,
+                              "dx_params" => dx_params,
+                              "h" => h,
+                              "diffusion" => diffusion,
+                              "γ" => γ,
+                              "gamma" => γ,
+                              "obs_un" => obs_un,
+                              "obs_cov" => obs_cov,
+                              "state_cov" => state_cov
+                             )
+    
+    # create storage for the forecast and analysis statistics
+    fore_rmse = Vector{Float64}(undef, nanl)
+    filt_rmse = Vector{Float64}(undef, nanl)
+    
+    for i in 1:(nanl-1)
+        for j in 1:f_steps
+            # M(x^b)
+            step_model!(x_b, 0.0, kwargs)
+        end
+
+        w = rand(MvNormal(zeros(40), I))
+        obs = o[:, i+1] + w
+    
+        # optimized cost function input and value 
+        x_opt = XdVAR.D3_var_NewtonOp(x_b, obs, x_b, state_cov, H_obs, obs_cov, kwargs)
+        
+        # generate actual observation value
+        x_t = o[:, i+1]
+
+        # compare model forecast and filter via RMSE
+        rmse_forecast = sqrt(msd(x_b, x_t))
+        fore_rmse[i] = rmse_forecast
+        rmse_filter = sqrt(msd(x_opt, x_t))
+        filt_rmse[i] = rmse_filter
+
+        # reinitializing x_b for next cycle if updated
+        if is_updated == true
+            x_b = x_opt
+        end
+    end
+
+    data = Dict{String,Any}(
+                            "seed" => seed,
+                            "diffusion" => diffusion,
+                            "dx_params" => dx_params,
+                            "gamma" => γ,
+                            "γ" => γ,
+                            "tanl" => tanl,
+                            "nanl" => nanl,
+                            "h" =>  h,
+                            "fore_rmse" => fore_rmse,
+                            "filt_rmse" => filt_rmse
+                           )
+
+    if is_informed == true
+        inf = "true"
+    else
+        inf = "false"
+    end
+
+    if is_updated == true
+        upd = "true"
+    else
+        upd = "false"
+    end
+
+    path = pkgdir(DataAssimilationBenchmarks) * "/src/data/d3_var_exp/"
+    name = "D3_var_filter_analysis_" * "L96_time_series_seed_" * lpad(seed, 4, "0") *
+           "_gam_" * rpad(γ, 5, "0") *
+           "_Informed_" * lpad(inf, 4, "0") *
+           "_Updated_" * lpad(upd, 4, "0") *
+           "_Tuned_" * rpad(tuning_factor, 5, "0") *
+           ".jld2"
+
+    save(path * name, data)
+    # output time
+    print("Runtime " * string(round((time() - t1)  / 60.0, digits=4))  * " minutes\n")
+
+    #= path = pkgdir(DataAssimilationBenchmarks) * "/src/analysis/var_exp/"
+    name = "D3_var_filter_analysis_" * "L96_time_series_seed_" * lpad(seed, 4, "0") *
+           "_gam_" * rpad(γ, 5, "0") *
+           "_Informed_" * lpad(is_informed, 4, "0") *
+           "_Updated_" * rpad(is_informed, 4, "0") *
+           "_Tuned_" * lpad(tuning_factor, 5, "0")
+
+    # make plot
+    t = 1:nanl
+    fore_rmse_ra = Vector{Float64}(undef, nanl)
+    filt_rmse_ra = Vector{Float64}(undef, nanl)
+
+    for i in 1:nanl
+        fore_rmse_ra[i] = sum(fore_rmse[1:i])/i
+        filt_rmse_ra[i] = sum(filt_rmse[1:i])/i
+    end
+
+    plot(t, fore_rmse_ra, label = "Forecast", title="Average Analysis RMSE vs. Time")
+    plot!(t, filt_rmse_ra, label = "Filter")
+    plot!([0, 5000], [1, 1])
+    xlabel!("Time [Cycles]")
+    ylabel!("Average Analysis RMSE")
+    savefig(path * name)=#
+end
+
+##############################################################################################
+# end module
+
+end
\ No newline at end of file
diff --git a/src/experiments/FilterExps.jl b/src/experiments/FilterExps.jl
old mode 100644
new mode 100755
diff --git a/src/experiments/GenerateTimeSeries.jl b/src/experiments/GenerateTimeSeries.jl
old mode 100644
new mode 100755
diff --git a/src/experiments/ParallelExperimentDriver.jl b/src/experiments/ParallelExperimentDriver.jl
old mode 100644
new mode 100755
diff --git a/src/experiments/SingleExperimentDriver.jl b/src/experiments/SingleExperimentDriver.jl
old mode 100644
new mode 100755
diff --git a/src/experiments/SmootherExps.jl b/src/experiments/SmootherExps.jl
old mode 100644
new mode 100755
diff --git a/src/experiments/VarAnalysisExperimentDriver.jl b/src/experiments/VarAnalysisExperimentDriver.jl
new file mode 100755
index 00000000..aeae340e
--- /dev/null
+++ b/src/experiments/VarAnalysisExperimentDriver.jl
@@ -0,0 +1,62 @@
+##############################################################################################
+module VarAnalysisExperimentDriver
+##############################################################################################
+# imports and exports
+using ..DataAssimilationBenchmarks
+export D3_var_filter_tuning
+
+##############################################################################################
+# Utility methods and definitions
+##############################################################################################
+
+path = pkgdir(DataAssimilationBenchmarks) * "/src/analysis/var_exp/"
+
+##############################################################################################
+# Filters
+##############################################################################################
+"""
+    args, exp = D3_var_filter_tuning()
+Constucts a parameter map and experiment wrapper for sensitivity test of covariance tuning.
+"""
+function D3_var_filter_tuning()
+
+    exp = DataAssimilationBenchmarks.D3VARExps.D3_var_filter_analysis
+    function wrap_exp(arguments)
+        try
+            exp(arguments)
+        catch
+            print("Error on " * string(arguments) * "\n")
+        end
+    end
+
+    # set filter parameters
+    time_series = "L96_time_series_seed_0000_dim_40_diff_0.000_F_08.0_tanl_0.05_nanl_05000_spin_1500_h_0.050.jld2"
+    γ = 1.0
+    is_informed = true
+    is_updated = true
+
+    # define tuning range
+    tuning_min = 0.01
+    tuning_step = 0.001
+    tuning_max = 0.03
+    
+    # load the experiments
+    args = Vector{Any}()
+    for t in tuning_min:tuning_step:tuning_max
+        tmp = (
+                time_series = time_series,
+                γ = γ,
+                is_informed = is_informed,
+                tuning_factor = t,
+                is_updated = is_updated
+            )
+        push!(args, tmp)
+    end
+    return args, wrap_exp
+end
+
+
+##############################################################################################
+# end module
+
+end
\ No newline at end of file
diff --git a/src/experiments/plot_var_analysis.jl b/src/experiments/plot_var_analysis.jl
new file mode 100755
index 00000000..4368a358
--- /dev/null
+++ b/src/experiments/plot_var_analysis.jl
@@ -0,0 +1,48 @@
+##############################################################################################
+module plot_var_analysis
+##############################################################################################
+# imports and exports
+using Random, Distributions, LinearAlgebra, StatsBase, Statistics
+using JLD2, HDF5, Plots
+using Measures
+using DataAssimilationBenchmarks
+
+##############################################################################################
+# define tuning range
+tuning_min = 0.01
+tuning_step = 0.01
+tuning_max = 1.0
+
+
+t = collect(tuning_min:tuning_step:tuning_max)
+stab_fore_rmse = Vector{Float64}(undef, length(t))
+stab_filt_rmse = Vector{Float64}(undef, length(t))
+
+loadpath = pkgdir(DataAssimilationBenchmarks) * "/src/data/d3_var_exp/"
+
+for i in 1:length(t)
+      time_series = "D3_var_filter_analysis_L96_time_series_seed_0000_gam_1.000_Informed_false" * 
+      "_Updated_true_Tuned_" * rpad(t[i], 5, "0") * ".jld2"
+
+      ts = load(loadpath * time_series)::Dict{String,Any}
+      nanl = ts["nanl"]::Int64
+      fore_rmse = ts["fore_rmse"]::Array{Float64, 1}
+      filt_rmse = ts["filt_rmse"]::Array{Float64, 1}
+
+      stab_fore_rmse[i] = sum(fore_rmse[1:(nanl-1)])/(nanl-1)
+      stab_filt_rmse[i] = sum(filt_rmse[1:(nanl-1)])/(nanl-1)
+end
+
+print("Minimum: ")
+print(findmin(stab_filt_rmse))
+plot(t, stab_fore_rmse, label = "Forecast", title="Uninformative Prior: Stabilized Analysis RMSE vs. Tuning Parameter", legend_position = :topright, margin=15mm, size=(800,500), dpi = 600)
+plot!(t, stab_filt_rmse, label = "Filter", )
+plot!([tuning_min, tuning_max], [1, 1], label = "Observation Error Std")
+xlabel!("Tuning Parameter")
+ylabel!("Stabilized Analysis RMSE")
+savefig("stab_I")
+
+##############################################################################################
+# end module
+
+end
diff --git a/src/experiments/run_var_analysis_test.jl b/src/experiments/run_var_analysis_test.jl
new file mode 100755
index 00000000..1d97536b
--- /dev/null
+++ b/src/experiments/run_var_analysis_test.jl
@@ -0,0 +1,29 @@
+##############################################################################################
+module run_var_analysis_test
+##############################################################################################
+# imports and exports
+using Distributed
+@everywhere using DataAssimilationBenchmarks
+##############################################################################################
+
+config = VarAnalysisExperimentDriver.D3_var_filter_tuning
+
+print("Generating experiment configurations from " * string(config) * "\n")
+print("Generate truth twin\n")
+
+args, exp = config()
+num_exps = length(args)
+
+print("Configuration ready\n")
+print("\n")
+print("Running " * string(num_exps) * " configurations on " * string(nworkers()) *
+      " total workers\n")
+print("Begin pmap\n")
+pmap(exp, args)
+print("Experiments completed, verify outputs in the appropriate directory under:\n")
+print(pkgdir(DataAssimilationBenchmarks) * "/src/data/d3_var_exp/\n")
+
+##############################################################################################
+# end module
+
+end
diff --git a/src/methods/DeSolvers.jl b/src/methods/DeSolvers.jl
old mode 100644
new mode 100755
diff --git a/src/methods/EnsembleKalmanSchemes.jl b/src/methods/EnsembleKalmanSchemes.jl
old mode 100644
new mode 100755
diff --git a/src/methods/XdVAR.jl b/src/methods/XdVAR.jl
old mode 100644
new mode 100755
index 05996ddf..063d53f5
--- a/src/methods/XdVAR.jl
+++ b/src/methods/XdVAR.jl
@@ -91,7 +91,7 @@ function D3_var_NewtonOp(x::VecA(T), obs::VecA(T), x_background::VecA(T), state_
     hess_x = Array{Float64}(undef, (sys_dim, sys_dim))
 
     while j <= j_max
-        print("Iteration: " * string(j) * "\n")
+        #print("Iteration: " * string(j) * "\n")
         # compute the gradient and hessian
         grad!(grad_x, x)
         hess!(hess_x, x)
diff --git a/src/models/IEEE39bus.jl b/src/models/IEEE39bus.jl
old mode 100644
new mode 100755
diff --git a/src/models/IEEE39bus_inputs/NE_EffectiveNetworkParams.jld2 b/src/models/IEEE39bus_inputs/NE_EffectiveNetworkParams.jld2
old mode 100644
new mode 100755
diff --git a/src/models/L96.jl b/src/models/L96.jl
old mode 100644
new mode 100755
diff --git a/src/models/ObsOperators.jl b/src/models/ObsOperators.jl
old mode 100644
new mode 100755
diff --git a/test/Test3dVAR.jl b/test/Test3dVAR.jl
old mode 100644
new mode 100755
diff --git a/test/TestClassicSmootherExps.jl b/test/TestClassicSmootherExps.jl
old mode 100644
new mode 100755
diff --git a/test/TestDeSolvers.jl b/test/TestDeSolvers.jl
old mode 100644
new mode 100755
diff --git a/test/TestFilterExps.jl b/test/TestFilterExps.jl
old mode 100644
new mode 100755
diff --git a/test/TestGenerateTimeSeries.jl b/test/TestGenerateTimeSeries.jl
old mode 100644
new mode 100755
diff --git a/test/TestIEEE39bus.jl b/test/TestIEEE39bus.jl
old mode 100644
new mode 100755
diff --git a/test/TestIterativeSmootherExps.jl b/test/TestIterativeSmootherExps.jl
old mode 100644
new mode 100755
diff --git a/test/TestL96.jl b/test/TestL96.jl
old mode 100644
new mode 100755
diff --git a/test/TestObsOperators.jl b/test/TestObsOperators.jl
old mode 100644
new mode 100755
diff --git a/test/TestSingleIterationSmootherExps.jl b/test/TestSingleIterationSmootherExps.jl
old mode 100644
new mode 100755
diff --git a/test/runtests.jl b/test/runtests.jl
old mode 100644
new mode 100755
index 3c1d1fef..54795774
--- a/test/runtests.jl
+++ b/test/runtests.jl
@@ -19,14 +19,6 @@ include("TestSingleIterationSmootherExps.jl")
 ##############################################################################################
 # Run tests
 
-# test set 9: test 3D-VAR
-@testset "3DVAR" begin
-    @test Test3dVAR.testCost()
-    @test Test3dVAR.testGrad()
-    @test Test3dVAR.testNewton()
-    @test Test3dVAR.testNewtonNoise()
-end
-
 # test set 0: test Observation Operators jacobian
 @testset "Observation Operators" begin
     @test TestObsOperators.alternating_obs_jacobian_pos()