text:
### A Pluto.jl notebook ###
# v0.18.0

using Markdown
using InteractiveUtils

# This Pluto notebook uses @bind for interactivity. When running this notebook outside of Pluto, the following 'mock version' of @bind gives bound variables a default value (instead of an error).
macro bind(def, element)
    quote
        local iv = try Base.loaded_modules[Base.PkgId(Base.UUID("6e696c72-6542-2067-7265-42206c756150"), "AbstractPlutoDingetjes")].Bonds.initial_value catch; b -> missing; end
        local el = $(esc(element))
        global $(esc(def)) = Core.applicable(Base.get, el) ? Base.get(el) : iv(el)
        el
    end
end

# ╔═╡ 4478d0a8-360a-4792-8859-244d80d2902b
using Colors, FileIO, ImageIO, ImageShow, QuartzImageIO

# ╔═╡ 0e4ef764-37f3-4db8-b7d5-2d40d81e2797
using Statistics

# ╔═╡ 1750a808-0c77-466c-a7e8-e569082c60af
using PlutoUI

# ╔═╡ a3966b72-930d-4485-94f7-177c3c5a7391
using Plots

# ╔═╡ c4c0dd00-655c-419a-9b8d-3afb78f139b1
using LsqFit

# ╔═╡ 670a7c9a-7dd9-11ec-2f9e-5daa27e49f3c
md"""
# Dveco, non-jpeg image compression
For Pablo
"""

# ╔═╡ 8e642f27-ea51-4cf1-b798-d8a586770b57
md"""
## Libraries and stuff to get up and running
"""

# ╔═╡ c6e8d6fb-7993-4d9e-bf0a-c4dc1d1e4b4a
PlutoUI.TableOfContents()

# ╔═╡ 419dfb6d-f7ca-46eb-984a-99d10b804e89
Image = Matrix{<:Color}

# ╔═╡ b389b734-e61a-4348-9720-3475b654ad8d
Channel = Matrix{<:Number}

# ╔═╡ 386615cc-6bd9-4a49-ba70-39f277c0e1db
struct FancyPixel
	value::Number
	size::Tuple{Int64, Int64}
end

# ╔═╡ 6b6bb9b9-b2c1-409e-9fc1-60285c1d2177
function luminance(pixel::Color)
	0.2126*pixel.r + 0.7152*pixel.g + 0.0722*pixel.b
end

# ╔═╡ 85f6c36d-a395-442f-8ba7-f2a8506d1544
luminance(image::Image) = luminance.(image)

# ╔═╡ 394fb978-f350-4ddc-9de5-d5114e9a87ed
md"""
## Main implementation
"""

# ╔═╡ 4377d345-71b9-4aff-b8d9-a45c4e9eb05e
md"""
### Compression
"""

# ╔═╡ 9bbc975a-d0d8-4f8e-b6b6-85fc2a9823ca
md"""
First, we make a function to mortar a matrix (divide in 4)
"""

# ╔═╡ b68f03ea-e884-4614-b1ff-c1d92118c30b
function mortar(matrix::Matrix{T})::Matrix{Matrix{T}} where {T}
	(fx, fy) = size(matrix)

	x = fx÷2
	y = fy÷2

	brick = [
		[matrix[1:x, 1:y]] [matrix[1:x, y+1:end]]
		[matrix[x+1:end, 1:y]] [matrix[x+1:end, y+1:end]]
	]

	brick
end

# ╔═╡ b4ad67b5-5b71-444b-ad68-5064587905d4
md"""
Then, we make a recursive compression algorithm that takes the average if the standard deviation of a region is above a certain threshold
"""

# ╔═╡ d5623e00-571f-407f-a88c-f1a020102e1d
function compress(channel::Channel, threashold::Number)
	dev = length(channel) == 1 ? 0 : std(channel)

	if dev > threashold
		compress.(mortar(channel), threashold)
	else
		s = size(channel)
		if s == (1, 1)
			channel[1, 1]
		else
			FancyPixel(mean(channel), s)
		end
	end
end

# ╔═╡ dccb5b97-4865-436a-b2cb-e774b2b1d79b
md"""
### Decompression
"""

# ╔═╡ daf3699e-267c-4b65-a39a-af0c1978430c
md"""
To decompress the image we need to repeat the value stored
"""

# ╔═╡ 280de3c0-6329-4a74-b684-3adebefe033e
function reconstruct(pixel::Real)
	pixel
end

# ╔═╡ 7bcb0d03-a859-496f-b5f4-4d6ba17ab07a
function reconstruct(pixel::FancyPixel)
	fill(pixel.value, pixel.size)
end

# ╔═╡ 5e91243f-8f61-4076-a1bc-b60b5ce8fb4c
md"""
And then we need to flatten the mortified (xd) matrix
"""

# ╔═╡ 2fe9e1dc-81ab-450c-9544-6429c2a80d6b
function flatten(matrix::Matrix{<:Any})
	caca = nothing
	for i in 1:size(matrix)[1]
		caca = matrix[i, :]
	end
	caca

	vcat([
		hcat(matrix[i, :]...) for i in 1:size(matrix)[1]
	]...)
end

# ╔═╡ 83769dd6-6b73-4821-b87d-f9b3ae12efe6
function reconstruct(compressed::Array{<:Any})
	compressed = map(reconstruct, compressed)
	flatten(compressed)
end

# ╔═╡ c7cbae3a-248e-4f9a-ba42-b2c94a570665
md"""
### Demo
"""

# ╔═╡ 98eb608b-c7f1-45a0-911c-4d5ef1a1fdc5
function show_compressed(channel::Channel, threshold::Number)
	RGB.(reconstruct(compress(channel, threshold)))
end

# ╔═╡ 347f8c01-96d6-41d0-a785-15e51886f89d
begin
	url = "https://upload.wikimedia.org/wikipedia/commons/3/39/GodfreyKneller-IsaacNewton-1689.jpg"
	philip_filename = download(url)
	test = load(philip_filename)
end

# ╔═╡ c32aeff2-72d4-4869-aaa7-e9bb26d59ab8
@bind threshold_bnw Slider(0.3:-0.001:0)

# ╔═╡ 751e0e94-7899-4c25-8519-6ce98b7b5f24
md"""
## Adding color

Kinda rudimentary
"""

# ╔═╡ cf2d80f0-9b2a-4162-a1b5-acc41a3dedcc
function channels(image::Image)
	[red.(image), green.(image), blue.(image)]
end

# ╔═╡ 0320566c-c8ca-4696-9aff-4d737be2ea97
function compress_by_channel(image::Image, threshold::Number)
	compress.(channels(image), threshold)
end

# ╔═╡ 8802eb09-e619-4d04-ae11-37fb19649fa1
function show_compressed(image::Image, threshold::Number)
	comp = compress_by_channel(image, threshold)
	comp = reconstruct.(comp)
	[RGB(r, g, b) for (r, g, b) in zip(comp...)]
end

# ╔═╡ 7cfd7bc3-fcd6-4034-9d8d-36a87eac1627
show_compressed(luminance(test), threshold_bnw)

# ╔═╡ c4d5c2f5-6603-4309-8fbf-ed41a8e90675
md"""
## File size

Rudimentary.
"""

# ╔═╡ e0452dc4-8bd1-470a-9cfd-782dfc4e301e
function filesize(pixel::Real, size_per_pixel = 3)
	return 1
end

# ╔═╡ 4d7c4da8-58aa-42af-a8a1-89160e73e08f
function filesize(compressed::Matrix, size_per_pixel = 3)
	tally = 0
	for i in compressed
		tally += filesize(i, size_per_pixel)
	end
	tally
end

# ╔═╡ 414a7f26-6c02-4deb-888f-62acbf52238c
function filesize(pixel::FancyPixel, size_per_pixel = 3)
	size_per_pixel
end

# ╔═╡ bf58ce82-c5ec-416b-8d42-d5673c2d156d
filesize(compress(luminance(test), 0.01))

# ╔═╡ 3ce732b7-d7fd-434b-95f9-f37c33192934
function compression_quality(original::Channel, compressed::Channel)
	sum([abs(o - c) for (o, c) in zip(original, compressed)])
end

# ╔═╡ 2db8f5bc-0b2a-4a1e-9cdb-125af1b001c4
function compression_stats(image::Image, threshold::Number)
	s = size(image)
	original = s[1] * s[2] * 3
	compressed = compress_by_channel(image, threshold)
	reconstructed = [RGB(r,g,b) for (r,g,b) in zip(reconstruct.(compressed)...)]
	csize = sum(filesize.(compressed))

	diff = 0
	for (o, p) in zip(image, reconstructed)
		diff += colordiff(o, p)
	end

	(original_in_kb = original/1000,
	comp_size_in_kb = csize/1000,
	ratio = csize/original,
	diff = diff)
end

# ╔═╡ 449a2e13-fe1f-49f8-a16d-7d160d396aa2
md"""
### Compression given a threshold
"""

# ╔═╡ 92c2398d-e393-48b3-a795-2a4afc799613
@bind threshold_color Slider(0.3:-0.001:0)

# ╔═╡ cb3e9ead-2e22-4b94-9ac2-8b63808df445
show_compressed(test, threshold_color)

# ╔═╡ 02839b7d-9944-4430-909a-cdba85cc64c2
compression_stats(test, threshold_color)

# ╔═╡ 703cf884-aa9a-4387-9a1b-afebdae5530e
md"Threshold is: $threshold_color"

# ╔═╡ 8ab423a8-4b33-477c-b445-1a10f86c260e
show_compressed(test, threshold_color)

# ╔═╡ 4b0750aa-70df-4bbe-8d90-fafea639f9af
md"""
### Average RGB
"""

# ╔═╡ c7404c7e-f847-4c27-be83-d8798a3c5dde
let
	c = mean([p for p in test])
	(r=c.r, g=c.g, b=c.b)
end

# ╔═╡ c3221ca6-7023-4bcc-865a-8c754ff45182
md"""
# Plotting!
"""

# ╔═╡ dab7a6d7-a679-41f6-ab05-cf0a6cc5418c
plotly()

# ╔═╡ 220b0f16-4569-4187-a4e9-badf4dca946e
md"""
Collecting the data takes a good while (~20s)
"""

# ╔═╡ 8db87c95-64d6-4176-830a-c906c21125da
compression_stats(test, 0.3)

# ╔═╡ cc3a57d7-ee3b-4975-9973-bad0bb9f2f53
md"""
# Regression
"""

# ╔═╡ d3d3b840-70ef-4c4d-932f-0c259682d141
model(t, p) = p[1] * exp.(-p[2] * t)

# ╔═╡ c0e07072-6c6e-4ec5-9795-bc5d6bb74fe7
tdata = 0.3:-0.01:0

# ╔═╡ 01b82335-b1bd-457a-8bf8-9ec3de0bd375
test_data = [compression_stats(test, x) for x in tdata]

# ╔═╡ 0fc5c49f-2f2d-419e-a0cb-76a53b641fbe
let
	range = tdata
	canvas = plot(; xaxis="std threshold")
	plot!(range, [d.ratio for d in test_data], label="compression ratio")
	plot!(range, [1 / (log(100, d.diff + 1) + 1) for d in test_data], label="compression quality")
end

# ╔═╡ 4973ccbb-ea9b-49c3-b549-7b53c14dd685
begin
	fit1 = curve_fit(model, tdata, [d.ratio for d in test_data], [0.5, 0.5])
	p1 = fit1.param

	fit2 = curve_fit(model, tdata, [d.diff for d in test_data], [0.5, 0.5])
	p2 = fit2.param
	(ratio = p1, quality = p2)
end

# ╔═╡ 2b53da0c-c5dc-4788-9f28-36f5bb822559
plot(tdata, [model(t, p1) for t in tdata], label="ratio")

# ╔═╡ a772b2ce-5ef8-47a4-952b-683796fe3959
plot(tdata, [model(t, p2) for t in tdata], label="difference")

# ╔═╡ 133d23fa-d7d3-4dc8-b1f4-c8f068a8c07f
function camera_input(;max_size=500, default_url="https://i.imgur.com/SUmi94P.png")
"""

			<button id="stop" title="Stop video">✖</button>
			<button id="pop-out" title="Pop out/pop in">⏏</button>

		<video playsinline autoplay></video>

		<button id="shutter" title="Click to take a picture">📷</button>

		Enable webcam

""" |> HTML
end

# ╔═╡ 5adf6861-9a4e-4347-8770-47df895038a6
function process_raw_camera_data(raw_camera_data)
	# the raw image data is a long byte array, we need to transform it into something
	# more "Julian" - something with more _structure_.

	# The encoding of the raw byte stream is:
	# every 4 bytes is a single pixel
	# every pixel has 4 values: Red, Green, Blue, Alpha
	# (we ignore alpha for this notebook)

	# So to get the red values for each pixel, we take every 4th value, starting at
	# the 1st:
	reds_flat = UInt8.(raw_camera_data["data"][1:4:end])
	greens_flat = UInt8.(raw_camera_data["data"][2:4:end])
	blues_flat = UInt8.(raw_camera_data["data"][3:4:end])

	# but these are still 1-dimensional arrays, nicknamed 'flat' arrays
	# We will 'reshape' this into 2D arrays:

	width = raw_camera_data["width"]
	height = raw_camera_data["height"]

	# shuffle and flip to get it in the right shape
	reds = reshape(reds_flat, (width, height))' / 255.0
	greens = reshape(greens_flat, (width, height))' / 255.0
	blues = reshape(blues_flat, (width, height))' / 255.0

	# we have our 2D array for each color
	# Let's create a single 2D array, where each value contains the R, G and B value of
	# that pixel

	return RGB.(reds, greens, blues)
end

# ╔═╡ 00000000-0000-0000-0000-000000000001
PLUTO_PROJECT_TOML_CONTENTS = """
[deps]
Colors = "5ae59095-9a9b-59fe-a467-6f913c188581"
FileIO = "5789e2e9-d7fb-5bc7-8068-2c6fae9b9549"
ImageIO = "82e4d734-157c-48bb-816b-45c225c6df19"
ImageShow = "4e3cecfd-b093-5904-9786-8bbb286a6a31"
LsqFit = "2fda8390-95c7-5789-9bda-21331edee243"
Plots = "91a5bcdd-55d7-5caf-9e0b-520d859cae80"
PlutoUI = "7f904dfe-b85e-4ff6-b463-dae2292396a8"
QuartzImageIO = "dca85d43-d64c-5e67-8c65-017450d5d020"
Statistics = "10745b16-79ce-11e8-11f9-7d13ad32a3b2"

[compat]
Colors = "~0.12.8"
FileIO = "~1.12.0"
ImageIO = "~0.6.0"
ImageShow = "~0.3.3"
LsqFit = "~0.12.1"
Plots = "~1.25.8"
PlutoUI = "~0.7.34"
QuartzImageIO = "~0.7.4"
"""

# ╔═╡ 00000000-0000-0000-0000-000000000002
PLUTO_MANIFEST_TOML_CONTENTS = """
# This file is machine-generated - editing it directly is not advised

[[AbstractFFTs]]
deps = ["ChainRulesCore", "LinearAlgebra"]
git-tree-sha1 = "6f1d9bc1c08f9f4a8fa92e3ea3cb50153a1b40d4"
uuid = "621f4979-c628-5d54-868e-fcf4e3e8185c"
version = "1.1.0"

[[AbstractPlutoDingetjes]]
deps = ["Pkg"]
git-tree-sha1 = "8eaf9f1b4921132a4cff3f36a1d9ba923b14a481"
uuid = "6e696c72-6542-2067-7265-42206c756150"
version = "1.1.4"

[[Adapt]]
deps = ["LinearAlgebra"]
git-tree-sha1 = "af92965fb30777147966f58acb05da51c5616b5f"
uuid = "79e6a3ab-5dfb-504d-930d-738a2a938a0e"
version = "3.3.3"

[[ArgTools]]
uuid = "0dad84c5-d112-42e6-8d28-ef12dabb789f"
version = "1.1.1"

[[ArrayInterface]]
deps = ["Compat", "IfElse", "LinearAlgebra", "Requires", "SparseArrays", "Static"]
git-tree-sha1 = "745233d77146ad221629590b6d82fe7f1ddb478f"
uuid = "4fba245c-0d91-5ea0-9b3e-6abc04ee57a9"
version = "4.0.3"

[[Artifacts]]
uuid = "56f22d72-fd6d-98f1-02f0-08ddc0907c33"

[[Base64]]
uuid = "2a0f44e3-6c83-55bd-87e4-b1978d98bd5f"

[[Bzip2_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
git-tree-sha1 = "19a35467a82e236ff51bc17a3a44b69ef35185a2"
uuid = "6e34b625-4abd-537c-b88f-471c36dfa7a0"
version = "1.0.8+0"

[[CEnum]]
git-tree-sha1 = "215a9aa4a1f23fbd05b92769fdd62559488d70e9"
uuid = "fa961155-64e5-5f13-b03f-caf6b980ea82"
version = "0.4.1"

[[Cairo]]
deps = ["Cairo_jll", "Colors", "Glib_jll", "Graphics", "Libdl", "Pango_jll"]
git-tree-sha1 = "d0b3f8b4ad16cb0a2988c6788646a5e6a17b6b1b"
uuid = "159f3aea-2a34-519c-b102-8c37f9878175"
version = "1.0.5"

[[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"

[[ChainRulesCore]]
deps = ["Compat", "LinearAlgebra", "SparseArrays"]
git-tree-sha1 = "54fc4400de6e5c3e27be6047da2ef6ba355511f8"
uuid = "d360d2e6-b24c-11e9-a2a3-2a2ae2dbcce4"
version = "1.11.6"

[[ChangesOfVariables]]
deps = ["ChainRulesCore", "LinearAlgebra", "Test"]
git-tree-sha1 = "bf98fa45a0a4cee295de98d4c1462be26345b9a1"
uuid = "9e997f8a-9a97-42d5-a9f1-ce6bfc15e2c0"
version = "0.1.2"

[[ColorSchemes]]
deps = ["ColorTypes", "Colors", "FixedPointNumbers", "Luxor", "Random"]
git-tree-sha1 = "5b7d2a8b53c68dfdbce545e957a3b5cc4da80b01"
uuid = "35d6a980-a343-548e-a6ea-1d62b119f2f4"
version = "3.17.0"

[[ColorTypes]]
deps = ["FixedPointNumbers", "Random"]
git-tree-sha1 = "024fe24d83e4a5bf5fc80501a314ce0d1aa35597"
uuid = "3da002f7-5984-5a60-b8a6-cbb66c0b333f"
version = "0.11.0"

[[ColorVectorSpace]]
deps = ["ColorTypes", "FixedPointNumbers", "LinearAlgebra", "SpecialFunctions", "Statistics", "TensorCore"]
git-tree-sha1 = "3f1f500312161f1ae067abe07d13b40f78f32e07"
uuid = "c3611d14-8923-5661-9e6a-0046d554d3a4"
version = "0.9.8"

[[Colors]]
deps = ["ColorTypes", "FixedPointNumbers", "Reexport"]
git-tree-sha1 = "417b0ed7b8b838aa6ca0a87aadf1bb9eb111ce40"
uuid = "5ae59095-9a9b-59fe-a467-6f913c188581"
version = "0.12.8"

[[CommonSubexpressions]]
deps = ["MacroTools", "Test"]
git-tree-sha1 = "7b8a93dba8af7e3b42fecabf646260105ac373f7"
uuid = "bbf7d656-a473-5ed7-a52c-81e309532950"
version = "0.3.0"

[[Compat]]
deps = ["Base64", "Dates", "DelimitedFiles", "Distributed", "InteractiveUtils", "LibGit2", "Libdl", "LinearAlgebra", "Markdown", "Mmap", "Pkg", "Printf", "REPL", "Random", "SHA", "Serialization", "SharedArrays", "Sockets", "SparseArrays", "Statistics", "Test", "UUIDs", "Unicode"]
git-tree-sha1 = "44c37b4636bc54afac5c574d2d02b625349d6582"
uuid = "34da2185-b29b-5c13-b0c7-acf172513d20"
version = "3.41.0"

[[CompilerSupportLibraries_jll]]
deps = ["Artifacts", "Libdl"]
uuid = "e66e0078-7015-5450-92f7-15fbd957f2ae"
version = "0.5.0+0"

[[Contour]]
deps = ["StaticArrays"]
git-tree-sha1 = "9f02045d934dc030edad45944ea80dbd1f0ebea7"
uuid = "d38c429a-6771-53c6-b99e-75d170b6e991"
version = "0.5.7"

[[DataAPI]]
git-tree-sha1 = "cc70b17275652eb47bc9e5f81635981f13cea5c8"
uuid = "9a962f9c-6df0-11e9-0e5d-c546b8b5ee8a"
version = "1.9.0"

[[DataStructures]]
deps = ["Compat", "InteractiveUtils", "OrderedCollections"]
git-tree-sha1 = "3daef5523dd2e769dad2365274f760ff5f282c7d"
uuid = "864edb3b-99cc-5e75-8d2d-829cb0a9cfe8"
version = "0.18.11"

[[DataValueInterfaces]]
git-tree-sha1 = "bfc1187b79289637fa0ef6d4436ebdfe6905cbd6"
uuid = "e2d170a0-9d28-54be-80f0-106bbe20a464"
version = "1.0.0"

[[Dates]]
deps = ["Printf"]
uuid = "ade2ca70-3891-5945-98fb-dc099432e06a"

[[DelimitedFiles]]
deps = ["Mmap"]
uuid = "8bb1440f-4735-579b-a4ab-409b98df4dab"

[[DensityInterface]]
deps = ["InverseFunctions", "Test"]
git-tree-sha1 = "80c3e8639e3353e5d2912fb3a1916b8455e2494b"
uuid = "b429d917-457f-4dbc-8f4c-0cc954292b1d"
version = "0.4.0"

[[DiffResults]]
deps = ["StaticArrays"]
git-tree-sha1 = "c18e98cba888c6c25d1c3b048e4b3380ca956805"
uuid = "163ba53b-c6d8-5494-b064-1a9d43ac40c5"
version = "1.0.3"

[[DiffRules]]
deps = ["IrrationalConstants", "LogExpFunctions", "NaNMath", "Random", "SpecialFunctions"]
git-tree-sha1 = "84083a5136b6abf426174a58325ffd159dd6d94f"
uuid = "b552c78f-8df3-52c6-915a-8e097449b14b"
version = "1.9.1"

[[Distributed]]
deps = ["Random", "Serialization", "Sockets"]
uuid = "8ba89e20-285c-5b6f-9357-94700520ee1b"

[[Distributions]]
deps = ["ChainRulesCore", "DensityInterface", "FillArrays", "LinearAlgebra", "PDMats", "Printf", "QuadGK", "Random", "SparseArrays", "SpecialFunctions", "Statistics", "StatsBase", "StatsFuns", "Test"]
git-tree-sha1 = "2e97190dfd4382499a4ac349e8d316491c9db341"
uuid = "31c24e10-a181-5473-b8eb-7969acd0382f"
version = "0.25.46"

[[DocStringExtensions]]
deps = ["LibGit2"]
git-tree-sha1 = "b19534d1895d702889b219c382a6e18010797f0b"
uuid = "ffbed154-4ef7-542d-bbb7-c09d3a79fcae"
version = "0.8.6"

[[Downloads]]
deps = ["ArgTools", "FileWatching", "LibCURL", "NetworkOptions"]
uuid = "f43a241f-c20a-4ad4-852c-f6b1247861c6"
version = "1.6.0"

[[EarCut_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
git-tree-sha1 = "3f3a2501fa7236e9b911e0f7a588c657e822bb6d"
uuid = "5ae413db-bbd1-5e63-b57d-d24a61df00f5"
version = "2.2.3+0"

[[Expat_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
git-tree-sha1 = "ae13fcbc7ab8f16b0856729b050ef0c446aa3492"
uuid = "2e619515-83b5-522b-bb60-26c02a35a201"
version = "2.4.4+0"

[[FFMPEG]]
deps = ["FFMPEG_jll"]
git-tree-sha1 = "b57e3acbe22f8484b4b5ff66a7499717fe1a9cc8"
uuid = "c87230d0-a227-11e9-1b43-d7ebe4e7570a"
version = "0.4.1"

[[FFMPEG_jll]]
deps = ["Artifacts", "Bzip2_jll", "FreeType2_jll", "FriBidi_jll", "JLLWrappers", "LAME_jll", "Libdl", "Ogg_jll", "OpenSSL_jll", "Opus_jll", "Pkg", "Zlib_jll", "libass_jll", "libfdk_aac_jll", "libvorbis_jll", "x264_jll", "x265_jll"]
git-tree-sha1 = "d8a578692e3077ac998b50c0217dfd67f21d1e5f"
uuid = "b22a6f82-2f65-5046-a5b2-351ab43fb4e5"
version = "4.4.0+0"

[[FileIO]]
deps = ["Pkg", "Requires", "UUIDs"]
git-tree-sha1 = "67551df041955cc6ee2ed098718c8fcd7fc7aebe"
uuid = "5789e2e9-d7fb-5bc7-8068-2c6fae9b9549"
version = "1.12.0"

[[FileWatching]]
uuid = "7b1f6079-737a-58dc-b8bc-7a2ca5c1b5ee"

[[FillArrays]]
deps = ["LinearAlgebra", "Random", "SparseArrays", "Statistics"]
git-tree-sha1 = "deed294cde3de20ae0b2e0355a6c4e1c6a5ceffc"
uuid = "1a297f60-69ca-5386-bcde-b61e274b549b"
version = "0.12.8"

[[FiniteDiff]]
deps = ["ArrayInterface", "LinearAlgebra", "Requires", "SparseArrays", "StaticArrays"]
git-tree-sha1 = "6eae72e9943d8992d14359c32aed5f892bda1569"
uuid = "6a86dc24-6348-571c-b903-95158fe2bd41"
version = "2.10.0"

[[FixedPointNumbers]]
deps = ["Statistics"]
git-tree-sha1 = "335bfdceacc84c5cdf16aadc768aa5ddfc5383cc"
uuid = "53c48c17-4a7d-5ca2-90c5-79b7896eea93"
version = "0.8.4"

[[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"

[[Formatting]]
deps = ["Printf"]
git-tree-sha1 = "8339d61043228fdd3eb658d86c926cb282ae72a8"
uuid = "59287772-0a20-5a39-b81b-1366585eb4c0"
version = "0.4.2"

[[ForwardDiff]]
deps = ["CommonSubexpressions", "DiffResults", "DiffRules", "LinearAlgebra", "LogExpFunctions", "NaNMath", "Preferences", "Printf", "Random", "SpecialFunctions", "StaticArrays"]
git-tree-sha1 = "1bd6fc0c344fc0cbee1f42f8d2e7ec8253dda2d2"
uuid = "f6369f11-7733-5829-9624-2563aa707210"
version = "0.10.25"

[[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"

[[FriBidi_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
git-tree-sha1 = "aa31987c2ba8704e23c6c8ba8a4f769d5d7e4f91"
uuid = "559328eb-81f9-559d-9380-de523a88c83c"
version = "1.0.10+0"

[[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"

[[GR]]
deps = ["Base64", "DelimitedFiles", "GR_jll", "HTTP", "JSON", "Libdl", "LinearAlgebra", "Pkg", "Printf", "Random", "RelocatableFolders", "Serialization", "Sockets", "Test", "UUIDs"]
git-tree-sha1 = "4a740db447aae0fbeb3ee730de1afbb14ac798a1"
uuid = "28b8d3ca-fb5f-59d9-8090-bfdbd6d07a71"
version = "0.63.1"

[[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 = "aa22e1ee9e722f1da183eb33370df4c1aeb6c2cd"
uuid = "d2c73de3-f751-5644-a686-071e5b155ba9"
version = "0.63.1+0"

[[GeometryBasics]]
deps = ["EarCut_jll", "IterTools", "LinearAlgebra", "StaticArrays", "StructArrays", "Tables"]
git-tree-sha1 = "58bcdf5ebc057b085e58d95c138725628dd7453c"
uuid = "5c1252a2-5f33-56bf-86c9-59e7332b4326"
version = "0.4.1"

[[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"

[[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"

[[Graphics]]
deps = ["Colors", "LinearAlgebra", "NaNMath"]
git-tree-sha1 = "1c5a84319923bea76fa145d49e93aa4394c73fc2"
uuid = "a2bd30eb-e257-5431-a919-1863eab51364"
version = "1.1.1"

[[Graphite2_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
git-tree-sha1 = "344bf40dcab1073aca04aa0df4fb092f920e4011"
uuid = "3b182d85-2403-5c21-9c21-1e1f0cc25472"
version = "1.3.14+0"

[[Grisu]]
git-tree-sha1 = "53bb909d1151e57e2484c3d1b53e19552b887fb2"
uuid = "42e2da0e-8278-4e71-bc24-59509adca0fe"
version = "1.0.2"

[[HTTP]]
deps = ["Base64", "Dates", "IniFile", "Logging", "MbedTLS", "NetworkOptions", "Sockets", "URIs"]
git-tree-sha1 = "0fa77022fe4b511826b39c894c90daf5fce3334a"
uuid = "cd3eb016-35fb-5094-929b-558a96fad6f3"
version = "0.9.17"

[[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"

[[Hyperscript]]
deps = ["Test"]
git-tree-sha1 = "8d511d5b81240fc8e6802386302675bdf47737b9"
uuid = "47d2ed2b-36de-50cf-bf87-49c2cf4b8b91"
version = "0.0.4"

[[HypertextLiteral]]
git-tree-sha1 = "2b078b5a615c6c0396c77810d92ee8c6f470d238"
uuid = "ac1192a8-f4b3-4bfe-ba22-af5b92cd3ab2"
version = "0.9.3"

[[IOCapture]]
deps = ["Logging", "Random"]
git-tree-sha1 = "f7be53659ab06ddc986428d3a9dcc95f6fa6705a"
uuid = "b5f81e59-6552-4d32-b1f0-c071b021bf89"
version = "0.2.2"

[[IfElse]]
git-tree-sha1 = "debdd00ffef04665ccbb3e150747a77560e8fad1"
uuid = "615f187c-cbe4-4ef1-ba3b-2fcf58d6d173"
version = "0.1.1"

[[ImageBase]]
deps = ["ImageCore", "Reexport"]
git-tree-sha1 = "b51bb8cae22c66d0f6357e3bcb6363145ef20835"
uuid = "c817782e-172a-44cc-b673-b171935fbb9e"
version = "0.1.5"

[[ImageCore]]
deps = ["AbstractFFTs", "ColorVectorSpace", "Colors", "FixedPointNumbers", "Graphics", "MappedArrays", "MosaicViews", "OffsetArrays", "PaddedViews", "Reexport"]
git-tree-sha1 = "9a5c62f231e5bba35695a20988fc7cd6de7eeb5a"
uuid = "a09fc81d-aa75-5fe9-8630-4744c3626534"
version = "0.9.3"

[[ImageIO]]
deps = ["FileIO", "Netpbm", "OpenEXR", "PNGFiles", "QOI", "Sixel", "TiffImages", "UUIDs"]
git-tree-sha1 = "816fc866edd8307a6e79a575e6585bfab8cef27f"
uuid = "82e4d734-157c-48bb-816b-45c225c6df19"
version = "0.6.0"

[[ImageShow]]
deps = ["Base64", "FileIO", "ImageBase", "ImageCore", "OffsetArrays", "StackViews"]
git-tree-sha1 = "d0ac64c9bee0aed6fdbb2bc0e5dfa9a3a78e3acc"
uuid = "4e3cecfd-b093-5904-9786-8bbb286a6a31"
version = "0.3.3"

[[Imath_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
git-tree-sha1 = "87f7662e03a649cffa2e05bf19c303e168732d3e"
uuid = "905a6f67-0a94-5f89-b386-d35d92009cd1"
version = "3.1.2+0"

[[IndirectArrays]]
git-tree-sha1 = "012e604e1c7458645cb8b436f8fba789a51b257f"
uuid = "9b13fd28-a010-5f03-acff-a1bbcff69959"
version = "1.0.0"

[[Inflate]]
git-tree-sha1 = "f5fc07d4e706b84f72d54eedcc1c13d92fb0871c"
uuid = "d25df0c9-e2be-5dd7-82c8-3ad0b3e990b9"
version = "0.1.2"

[[IniFile]]
deps = ["Test"]
git-tree-sha1 = "098e4d2c533924c921f9f9847274f2ad89e018b8"
uuid = "83e8ac13-25f8-5344-8a64-a9f2b223428f"
version = "0.5.0"

[[InteractiveUtils]]
deps = ["Markdown"]
uuid = "b77e0a4c-d291-57a0-90e8-8db25a27a240"

[[InverseFunctions]]
deps = ["Test"]
git-tree-sha1 = "a7254c0acd8e62f1ac75ad24d5db43f5f19f3c65"
uuid = "3587e190-3f89-42d0-90ee-14403ec27112"
version = "0.1.2"

[[IrrationalConstants]]
git-tree-sha1 = "7fd44fd4ff43fc60815f8e764c0f352b83c49151"
uuid = "92d709cd-6900-40b7-9082-c6be49f344b6"
version = "0.1.1"

[[IterTools]]
git-tree-sha1 = "fa6287a4469f5e048d763df38279ee729fbd44e5"
uuid = "c8e1da08-722c-5040-9ed9-7db0dc04731e"
version = "1.4.0"

[[IteratorInterfaceExtensions]]
git-tree-sha1 = "a3f24677c21f5bbe9d2a714f95dcd58337fb2856"
uuid = "82899510-4779-5014-852e-03e436cf321d"
version = "1.0.0"

[[JLLWrappers]]
deps = ["Preferences"]
git-tree-sha1 = "22df5b96feef82434b07327e2d3c770a9b21e023"
uuid = "692b3bcd-3c85-4b1f-b108-f13ce0eb3210"
version = "1.4.0"

[[JSON]]
deps = ["Dates", "Mmap", "Parsers", "Unicode"]
git-tree-sha1 = "8076680b162ada2a031f707ac7b4953e30667a37"
uuid = "682c06a0-de6a-54ab-a142-c8b1cf79cde6"
version = "0.21.2"

[[JpegTurbo_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
git-tree-sha1 = "b53380851c6e6664204efb2e62cd24fa5c47e4ba"
uuid = "aacddb02-875f-59d6-b918-886e6ef4fbf8"
version = "2.1.2+0"

[[Juno]]
deps = ["Base64", "Logging", "Media", "Profile"]
git-tree-sha1 = "07cb43290a840908a771552911a6274bc6c072c7"
uuid = "e5e0dc1b-0480-54bc-9374-aad01c23163d"
version = "0.8.4"

[[LAME_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
git-tree-sha1 = "f6250b16881adf048549549fba48b1161acdac8c"
uuid = "c1c5ebd0-6772-5130-a774-d5fcae4a789d"
version = "3.100.1+0"

[[LZO_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
git-tree-sha1 = "e5b909bcf985c5e2605737d2ce278ed791b89be6"
uuid = "dd4b983a-f0e5-5f8d-a1b7-129d4a5fb1ac"
version = "2.10.1+0"

[[LaTeXStrings]]
git-tree-sha1 = "f2355693d6778a178ade15952b7ac47a4ff97996"
uuid = "b964fa9f-0449-5b57-a5c2-d3ea65f4040f"
version = "1.3.0"

[[Latexify]]
deps = ["Formatting", "InteractiveUtils", "LaTeXStrings", "MacroTools", "Markdown", "Printf", "Requires"]
git-tree-sha1 = "a8f4f279b6fa3c3c4f1adadd78a621b13a506bce"
uuid = "23fbe1c1-3f47-55db-b15f-69d7ec21a316"
version = "0.15.9"

[[LibCURL]]
deps = ["LibCURL_jll", "MozillaCACerts_jll"]
uuid = "b27032c2-a3e7-50c8-80cd-2d36dbcbfd21"
version = "0.6.3"

[[LibCURL_jll]]
deps = ["Artifacts", "LibSSH2_jll", "Libdl", "MbedTLS_jll", "Zlib_jll", "nghttp2_jll"]
uuid = "deac9b47-8bc7-5906-a0fe-35ac56dc84c0"
version = "7.73.0+4"

[[LibGit2]]
deps = ["Base64", "NetworkOptions", "Printf", "SHA"]
uuid = "76f85450-5226-5b5a-8eaa-529ad045b433"

[[LibSSH2_jll]]
deps = ["Artifacts", "Libdl", "MbedTLS_jll"]
uuid = "29816b5a-b9ab-546f-933c-edad1886dfa8"
version = "1.9.1+2"

[[Libdl]]
uuid = "8f399da3-3557-5675-b5ff-fb832c97cbdb"

[[Libffi_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
git-tree-sha1 = "0b4a5d71f3e5200a7dff793393e09dfc2d874290"
uuid = "e9f186c6-92d2-5b65-8a66-fee21dc1b490"
version = "3.2.2+1"

[[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"

[[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"

[[Libgpg_error_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
git-tree-sha1 = "c333716e46366857753e273ce6a69ee0945a6db9"
uuid = "7add5ba3-2f88-524e-9cd5-f83b8a55f7b8"
version = "1.42.0+0"

[[Libiconv_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
git-tree-sha1 = "42b62845d70a619f063a7da093d995ec8e15e778"
uuid = "94ce4f54-9a6c-5748-9c1c-f9c7231a4531"
version = "1.16.1+1"

[[Libmount_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
git-tree-sha1 = "9c30530bf0effd46e15e0fdcf2b8636e78cbbd73"
uuid = "4b2f31a3-9ecc-558c-b454-b3730dcb73e9"
version = "2.35.0+0"

[[Librsvg_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl", "Pango_jll", "Pkg", "gdk_pixbuf_jll"]
git-tree-sha1 = "25d5e6b4eb3558613ace1c67d6a871420bfca527"
uuid = "925c91fb-5dd6-59dd-8e8c-345e74382d89"
version = "2.52.4+0"

[[Libtiff_jll]]
deps = ["Artifacts", "JLLWrappers", "JpegTurbo_jll", "Libdl", "Pkg", "Zlib_jll", "Zstd_jll"]
git-tree-sha1 = "340e257aada13f95f98ee352d316c3bed37c8ab9"
uuid = "89763e89-9b03-5906-acba-b20f662cd828"
version = "4.3.0+0"

[[Libuuid_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
git-tree-sha1 = "7f3efec06033682db852f8b3bc3c1d2b0a0ab066"
uuid = "38a345b3-de98-5d2b-a5d3-14cd9215e700"
version = "2.36.0+0"

[[LinearAlgebra]]
deps = ["Libdl", "libblastrampoline_jll"]
uuid = "37e2e46d-f89d-539d-b4ee-838fcccc9c8e"

[[LogExpFunctions]]
deps = ["ChainRulesCore", "ChangesOfVariables", "DocStringExtensions", "InverseFunctions", "IrrationalConstants", "LinearAlgebra"]
git-tree-sha1 = "e5718a00af0ab9756305a0392832c8952c7426c1"
uuid = "2ab3a3ac-af41-5b50-aa03-7779005ae688"
version = "0.3.6"

[[Logging]]
uuid = "56ddb016-857b-54e1-b83d-db4d58db5568"

[[LsqFit]]
deps = ["Distributions", "ForwardDiff", "LinearAlgebra", "NLSolversBase", "OptimBase", "Random", "StatsBase"]
git-tree-sha1 = "91aa1442e63a77f101aff01dec5a821a17f43922"
uuid = "2fda8390-95c7-5789-9bda-21331edee243"
version = "0.12.1"

[[Luxor]]
deps = ["Base64", "Cairo", "Colors", "Dates", "FFMPEG", "FileIO", "Juno", "LaTeXStrings", "Random", "Requires", "Rsvg"]
git-tree-sha1 = "81a4fd2c618ba952feec85e4236f36c7a5660393"
uuid = "ae8d54c2-7ccd-5906-9d76-62fc9837b5bc"
version = "3.0.0"

[[MacroTools]]
deps = ["Markdown", "Random"]
git-tree-sha1 = "3d3e902b31198a27340d0bf00d6ac452866021cf"
uuid = "1914dd2f-81c6-5fcd-8719-6d5c9610ff09"
version = "0.5.9"

[[MappedArrays]]
git-tree-sha1 = "e8b359ef06ec72e8c030463fe02efe5527ee5142"
uuid = "dbb5928d-eab1-5f90-85c2-b9b0edb7c900"
version = "0.4.1"

[[Markdown]]
deps = ["Base64"]
uuid = "d6f4376e-aef5-505a-96c1-9c027394607a"

[[MbedTLS]]
deps = ["Dates", "MbedTLS_jll", "Random", "Sockets"]
git-tree-sha1 = "1c38e51c3d08ef2278062ebceade0e46cefc96fe"
uuid = "739be429-bea8-5141-9913-cc70e7f3736d"
version = "1.0.3"

[[MbedTLS_jll]]
deps = ["Artifacts", "Libdl"]
uuid = "c8ffd9c3-330d-5841-b78e-0817d7145fa1"
version = "2.24.0+2"

[[Measures]]
git-tree-sha1 = "e498ddeee6f9fdb4551ce855a46f54dbd900245f"
uuid = "442fdcdd-2543-5da2-b0f3-8c86c306513e"
version = "0.3.1"

[[Media]]
deps = ["MacroTools", "Test"]
git-tree-sha1 = "75a54abd10709c01f1b86b84ec225d26e840ed58"
uuid = "e89f7d12-3494-54d1-8411-f7d8b9ae1f27"
version = "0.5.0"

[[Missings]]
deps = ["DataAPI"]
git-tree-sha1 = "bf210ce90b6c9eed32d25dbcae1ebc565df2687f"
uuid = "e1d29d7a-bbdc-5cf2-9ac0-f12de2c33e28"
version = "1.0.2"

[[Mmap]]
uuid = "a63ad114-7e13-5084-954f-fe012c677804"

[[MosaicViews]]
deps = ["MappedArrays", "OffsetArrays", "PaddedViews", "StackViews"]
git-tree-sha1 = "b34e3bc3ca7c94914418637cb10cc4d1d80d877d"
uuid = "e94cdb99-869f-56ef-bcf0-1ae2bcbe0389"
version = "0.3.3"

[[MozillaCACerts_jll]]
uuid = "14a3606d-f60d-562e-9121-12d972cd8159"
version = "2020.7.22"

[[NLSolversBase]]
deps = ["DiffResults", "Distributed", "FiniteDiff", "ForwardDiff"]
git-tree-sha1 = "50310f934e55e5ca3912fb941dec199b49ca9b68"
uuid = "d41bc354-129a-5804-8e4c-c37616107c6c"
version = "7.8.2"

[[NaNMath]]
git-tree-sha1 = "f755f36b19a5116bb580de457cda0c140153f283"
uuid = "77ba4419-2d1f-58cd-9bb1-8ffee604a2e3"
version = "0.3.6"

[[Netpbm]]
deps = ["FileIO", "ImageCore"]
git-tree-sha1 = "18efc06f6ec36a8b801b23f076e3c6ac7c3bf153"
uuid = "f09324ee-3d7c-5217-9330-fc30815ba969"
version = "1.0.2"

[[NetworkOptions]]
uuid = "ca575930-c2e3-43a9-ace4-1e988b2c1908"
version = "1.2.0"

[[OffsetArrays]]
deps = ["Adapt"]
git-tree-sha1 = "043017e0bdeff61cfbb7afeb558ab29536bbb5ed"
uuid = "6fe1bfb0-de20-5000-8ca7-80f57d26f881"
version = "1.10.8"

[[Ogg_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
git-tree-sha1 = "887579a3eb005446d514ab7aeac5d1d027658b8f"
uuid = "e7412a2a-1a6e-54c0-be00-318e2571c051"
version = "1.3.5+1"

[[OpenBLAS_jll]]
deps = ["Artifacts", "CompilerSupportLibraries_jll", "Libdl"]
uuid = "4536629a-c528-5b80-bd46-f80d51c5b363"
version = "0.3.17+2"

[[OpenEXR]]
deps = ["Colors", "FileIO", "OpenEXR_jll"]
git-tree-sha1 = "327f53360fdb54df7ecd01e96ef1983536d1e633"
uuid = "52e1d378-f018-4a11-a4be-720524705ac7"
version = "0.3.2"

[[OpenEXR_jll]]
deps = ["Artifacts", "Imath_jll", "JLLWrappers", "Libdl", "Pkg", "Zlib_jll"]
git-tree-sha1 = "923319661e9a22712f24596ce81c54fc0366f304"
uuid = "18a262bb-aa17-5467-a713-aee519bc75cb"
version = "3.1.1+0"

[[OpenLibm_jll]]
deps = ["Artifacts", "Libdl"]
uuid = "05823500-19ac-5b8b-9628-191a04bc5112"
version = "0.8.1+0"

[[OpenSSL_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
git-tree-sha1 = "648107615c15d4e09f7eca16307bc821c1f718d8"
uuid = "458c3c95-2e84-50aa-8efc-19380b2a3a95"
version = "1.1.13+0"

[[OpenSpecFun_jll]]
deps = ["Artifacts", "CompilerSupportLibraries_jll", "JLLWrappers", "Libdl", "Pkg"]
git-tree-sha1 = "13652491f6856acfd2db29360e1bbcd4565d04f1"
uuid = "efe28fd5-8261-553b-a9e1-b2916fc3738e"
version = "0.5.5+0"

[[OptimBase]]
deps = ["NLSolversBase", "Printf", "Reexport"]
git-tree-sha1 = "9cb1fee807b599b5f803809e85c81b582d2009d6"
uuid = "87e2bd06-a317-5318-96d9-3ecbac512eee"
version = "2.0.2"

[[Opus_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
git-tree-sha1 = "51a08fb14ec28da2ec7a927c4337e4332c2a4720"
uuid = "91d4177d-7536-5919-b921-800302f37372"
version = "1.3.2+0"

[[OrderedCollections]]
git-tree-sha1 = "85f8e6578bf1f9ee0d11e7bb1b1456435479d47c"
uuid = "bac558e1-5e72-5ebc-8fee-abe8a469f55d"
version = "1.4.1"

[[PCRE_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
git-tree-sha1 = "b2a7af664e098055a7529ad1a900ded962bca488"
uuid = "2f80f16e-611a-54ab-bc61-aa92de5b98fc"
version = "8.44.0+0"

[[PDMats]]
deps = ["LinearAlgebra", "SparseArrays", "SuiteSparse"]
git-tree-sha1 = "ee26b350276c51697c9c2d88a072b339f9f03d73"
uuid = "90014a1f-27ba-587c-ab20-58faa44d9150"
version = "0.11.5"

[[PNGFiles]]
deps = ["Base64", "CEnum", "ImageCore", "IndirectArrays", "OffsetArrays", "libpng_jll"]
git-tree-sha1 = "6d105d40e30b635cfed9d52ec29cf456e27d38f8"
uuid = "f57f5aa1-a3ce-4bc8-8ab9-96f992907883"
version = "0.3.12"

[[PaddedViews]]
deps = ["OffsetArrays"]
git-tree-sha1 = "03a7a85b76381a3d04c7a1656039197e70eda03d"
uuid = "5432bcbf-9aad-5242-b902-cca2824c8663"
version = "0.5.11"

[[Pango_jll]]
deps = ["Artifacts", "Cairo_jll", "Fontconfig_jll", "FreeType2_jll", "FriBidi_jll", "Glib_jll", "HarfBuzz_jll", "JLLWrappers", "Libdl", "Pkg"]
git-tree-sha1 = "3a121dfbba67c94a5bec9dde613c3d0cbcf3a12b"
uuid = "36c8627f-9965-5494-a995-c6b170f724f3"
version = "1.50.3+0"

[[Parsers]]
deps = ["Dates"]
git-tree-sha1 = "0b5cfbb704034b5b4c1869e36634438a047df065"
uuid = "69de0a69-1ddd-5017-9359-2bf0b02dc9f0"
version = "2.2.1"

[[Pixman_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
git-tree-sha1 = "b4f5d02549a10e20780a24fce72bea96b6329e29"
uuid = "30392449-352a-5448-841d-b1acce4e97dc"
version = "0.40.1+0"

[[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"
version = "1.8.0"

[[PkgVersion]]
deps = ["Pkg"]
git-tree-sha1 = "a7a7e1a88853564e551e4eba8650f8c38df79b37"
uuid = "eebad327-c553-4316-9ea0-9fa01ccd7688"
version = "0.1.1"

[[PlotThemes]]
deps = ["PlotUtils", "Requires", "Statistics"]
git-tree-sha1 = "a3a964ce9dc7898193536002a6dd892b1b5a6f1d"
uuid = "ccf2f8ad-2431-5c83-bf29-c5338b663b6a"
version = "2.0.1"

[[PlotUtils]]
deps = ["ColorSchemes", "Colors", "Dates", "Printf", "Random", "Reexport", "Statistics"]
git-tree-sha1 = "6f1b25e8ea06279b5689263cc538f51331d7ca17"
uuid = "995b91a9-d308-5afd-9ec6-746e21dbc043"
version = "1.1.3"

[[Plots]]
deps = ["Base64", "Contour", "Dates", "Downloads", "FFMPEG", "FixedPointNumbers", "GR", "GeometryBasics", "JSON", "Latexify", "LinearAlgebra", "Measures", "NaNMath", "PlotThemes", "PlotUtils", "Printf", "REPL", "Random", "RecipesBase", "RecipesPipeline", "Reexport", "Requires", "Scratch", "Showoff", "SparseArrays", "Statistics", "StatsBase", "UUIDs", "UnicodeFun", "Unzip"]
git-tree-sha1 = "eb1432ec2b781f70ce2126c277d120554605669a"
uuid = "91a5bcdd-55d7-5caf-9e0b-520d859cae80"
version = "1.25.8"

[[PlutoUI]]
deps = ["AbstractPlutoDingetjes", "Base64", "ColorTypes", "Dates", "Hyperscript", "HypertextLiteral", "IOCapture", "InteractiveUtils", "JSON", "Logging", "Markdown", "Random", "Reexport", "UUIDs"]
git-tree-sha1 = "8979e9802b4ac3d58c503a20f2824ad67f9074dd"
uuid = "7f904dfe-b85e-4ff6-b463-dae2292396a8"
version = "0.7.34"

[[Preferences]]
deps = ["TOML"]
git-tree-sha1 = "2cf929d64681236a2e074ffafb8d568733d2e6af"
uuid = "21216c6a-2e73-6563-6e65-726566657250"
version = "1.2.3"

[[Printf]]
deps = ["Unicode"]
uuid = "de0858da-6303-5e67-8744-51eddeeeb8d7"

[[Profile]]
deps = ["Printf"]
uuid = "9abbd945-dff8-562f-b5e8-e1ebf5ef1b79"

[[ProgressMeter]]
deps = ["Distributed", "Printf"]
git-tree-sha1 = "afadeba63d90ff223a6a48d2009434ecee2ec9e8"
uuid = "92933f4c-e287-5a05-a399-4b506db050ca"
version = "1.7.1"

[[QOI]]
deps = ["ColorTypes", "FileIO", "FixedPointNumbers"]
git-tree-sha1 = "18e8f4d1426e965c7b532ddd260599e1510d26ce"
uuid = "4b34888f-f399-49d4-9bb3-47ed5cae4e65"
version = "1.0.0"

[[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 = "ad368663a5e20dbb8d6dc2fddeefe4dae0781ae8"
uuid = "ea2cea3b-5b76-57ae-a6ef-0a8af62496e1"
version = "5.15.3+0"

[[QuadGK]]
deps = ["DataStructures", "LinearAlgebra"]
git-tree-sha1 = "78aadffb3efd2155af139781b8a8df1ef279ea39"
uuid = "1fd47b50-473d-5c70-9696-f719f8f3bcdc"
version = "2.4.2"

[[QuartzImageIO]]
deps = ["FileIO", "ImageCore", "Libdl"]
git-tree-sha1 = "16de3b880ffdfbc8fc6707383c00a2e076bb0221"
uuid = "dca85d43-d64c-5e67-8c65-017450d5d020"
version = "0.7.4"

[[REPL]]
deps = ["InteractiveUtils", "Markdown", "Sockets", "Unicode"]
uuid = "3fa0cd96-eef1-5676-8a61-b3b8758bbffb"

[[Random]]
deps = ["SHA", "Serialization"]
uuid = "9a3f8284-a2c9-5f02-9a11-845980a1fd5c"

[[RecipesBase]]
git-tree-sha1 = "6bf3f380ff52ce0832ddd3a2a7b9538ed1bcca7d"
uuid = "3cdcf5f2-1ef4-517c-9805-6587b60abb01"
version = "1.2.1"

[[RecipesPipeline]]
deps = ["Dates", "NaNMath", "PlotUtils", "RecipesBase"]
git-tree-sha1 = "37c1631cb3cc36a535105e6d5557864c82cd8c2b"
uuid = "01d81517-befc-4cb6-b9ec-a95719d0359c"
version = "0.5.0"

[[Reexport]]
git-tree-sha1 = "45e428421666073eab6f2da5c9d310d99bb12f9b"
uuid = "189a3867-3050-52da-a836-e630ba90ab69"
version = "1.2.2"

[[RelocatableFolders]]
deps = ["SHA", "Scratch"]
git-tree-sha1 = "cdbd3b1338c72ce29d9584fdbe9e9b70eeb5adca"
uuid = "05181044-ff0b-4ac5-8273-598c1e38db00"
version = "0.1.3"

[[Requires]]
deps = ["UUIDs"]
git-tree-sha1 = "838a3a4188e2ded87a4f9f184b4b0d78a1e91cb7"
uuid = "ae029012-a4dd-5104-9daa-d747884805df"
version = "1.3.0"

[[Rmath]]
deps = ["Random", "Rmath_jll"]
git-tree-sha1 = "bf3188feca147ce108c76ad82c2792c57abe7b1f"
uuid = "79098fc4-a85e-5d69-aa6a-4863f24498fa"
version = "0.7.0"

[[Rmath_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
git-tree-sha1 = "68db32dff12bb6127bac73c209881191bf0efbb7"
uuid = "f50d1b31-88e8-58de-be2c-1cc44531875f"
version = "0.3.0+0"

[[Rsvg]]
deps = ["Cairo", "Glib_jll", "Librsvg_jll"]
git-tree-sha1 = "3d3dc66eb46568fb3a5259034bfc752a0eb0c686"
uuid = "c4c386cf-5103-5370-be45-f3a111cca3b8"
version = "1.0.0"

[[SHA]]
uuid = "ea8e919c-243c-51af-8825-aaa63cd721ce"
version = "0.7.0"

[[Scratch]]
deps = ["Dates"]
git-tree-sha1 = "0b4b7f1393cff97c33891da2a0bf69c6ed241fda"
uuid = "6c6a2e73-6563-6170-7368-637461726353"
version = "1.1.0"

[[Serialization]]
uuid = "9e88b42a-f829-5b0c-bbe9-9e923198166b"

[[SharedArrays]]
deps = ["Distributed", "Mmap", "Random", "Serialization"]
uuid = "1a1011a3-84de-559e-8e89-a11a2f7dc383"

[[Showoff]]
deps = ["Dates", "Grisu"]
git-tree-sha1 = "91eddf657aca81df9ae6ceb20b959ae5653ad1de"
uuid = "992d4aef-0814-514b-bc4d-f2e9a6c4116f"
version = "1.0.3"

[[Sixel]]
deps = ["Dates", "FileIO", "ImageCore", "IndirectArrays", "OffsetArrays", "REPL", "libsixel_jll"]
git-tree-sha1 = "8fb59825be681d451c246a795117f317ecbcaa28"
uuid = "45858cf5-a6b0-47a3-bbea-62219f50df47"
version = "0.1.2"

[[Sockets]]
uuid = "6462fe0b-24de-5631-8697-dd941f90decc"

[[SortingAlgorithms]]
deps = ["DataStructures"]
git-tree-sha1 = "b3363d7460f7d098ca0912c69b082f75625d7508"
uuid = "a2af1166-a08f-5f64-846c-94a0d3cef48c"
version = "1.0.1"

[[SparseArrays]]
deps = ["LinearAlgebra", "Random"]
uuid = "2f01184e-e22b-5df5-ae63-d93ebab69eaf"

[[SpecialFunctions]]
deps = ["ChainRulesCore", "IrrationalConstants", "LogExpFunctions", "OpenLibm_jll", "OpenSpecFun_jll"]
git-tree-sha1 = "e08890d19787ec25029113e88c34ec20cac1c91e"
uuid = "276daf66-3868-5448-9aa4-cd146d93841b"
version = "2.0.0"

[[StackViews]]
deps = ["OffsetArrays"]
git-tree-sha1 = "46e589465204cd0c08b4bd97385e4fa79a0c770c"
uuid = "cae243ae-269e-4f55-b966-ac2d0dc13c15"
version = "0.1.1"

[[Static]]
deps = ["IfElse"]
git-tree-sha1 = "d4da8b728580709d736704764e55d6ef38cb7c87"
uuid = "aedffcd0-7271-4cad-89d0-dc628f76c6d3"
version = "0.5.3"

[[StaticArrays]]
deps = ["LinearAlgebra", "Random", "Statistics"]
git-tree-sha1 = "a635a9333989a094bddc9f940c04c549cd66afcf"
uuid = "90137ffa-7385-5640-81b9-e52037218182"
version = "1.3.4"

[[Statistics]]
deps = ["LinearAlgebra", "SparseArrays"]
uuid = "10745b16-79ce-11e8-11f9-7d13ad32a3b2"

[[StatsAPI]]
git-tree-sha1 = "d88665adc9bcf45903013af0982e2fd05ae3d0a6"
uuid = "82ae8749-77ed-4fe6-ae5f-f523153014b0"
version = "1.2.0"

[[StatsBase]]
deps = ["DataAPI", "DataStructures", "LinearAlgebra", "LogExpFunctions", "Missings", "Printf", "Random", "SortingAlgorithms", "SparseArrays", "Statistics", "StatsAPI"]
git-tree-sha1 = "51383f2d367eb3b444c961d485c565e4c0cf4ba0"
uuid = "2913bbd2-ae8a-5f71-8c99-4fb6c76f3a91"
version = "0.33.14"

[[StatsFuns]]
deps = ["ChainRulesCore", "InverseFunctions", "IrrationalConstants", "LogExpFunctions", "Reexport", "Rmath", "SpecialFunctions"]
git-tree-sha1 = "f35e1879a71cca95f4826a14cdbf0b9e253ed918"
uuid = "4c63d2b9-4356-54db-8cca-17b64c39e42c"
version = "0.9.15"

[[StructArrays]]
deps = ["Adapt", "DataAPI", "StaticArrays", "Tables"]
git-tree-sha1 = "d21f2c564b21a202f4677c0fba5b5ee431058544"
uuid = "09ab397b-f2b6-538f-b94a-2f83cf4a842a"
version = "0.6.4"

[[SuiteSparse]]
deps = ["Libdl", "LinearAlgebra", "Serialization", "SparseArrays"]
uuid = "4607b0f0-06f3-5cda-b6b1-a6196a1729e9"

[[TOML]]
deps = ["Dates"]
uuid = "fa267f1f-6049-4f14-aa54-33bafae1ed76"
version = "1.0.0"

[[TableTraits]]
deps = ["IteratorInterfaceExtensions"]
git-tree-sha1 = "c06b2f539df1c6efa794486abfb6ed2022561a39"
uuid = "3783bdb8-4a98-5b6b-af9a-565f29a5fe9c"
version = "1.0.1"

[[Tables]]
deps = ["DataAPI", "DataValueInterfaces", "IteratorInterfaceExtensions", "LinearAlgebra", "TableTraits", "Test"]
git-tree-sha1 = "bb1064c9a84c52e277f1096cf41434b675cd368b"
uuid = "bd369af6-aec1-5ad0-b16a-f7cc5008161c"
version = "1.6.1"

[[Tar]]
deps = ["ArgTools", "SHA"]
uuid = "a4e569a6-e804-4fa4-b0f3-eef7a1d5b13e"
version = "1.10.0"

[[TensorCore]]
deps = ["LinearAlgebra"]
git-tree-sha1 = "1feb45f88d133a655e001435632f019a9a1bcdb6"
uuid = "62fd8b95-f654-4bbd-a8a5-9c27f68ccd50"
version = "0.1.1"

[[Test]]
deps = ["InteractiveUtils", "Logging", "Random", "Serialization"]
uuid = "8dfed614-e22c-5e08-85e1-65c5234f0b40"

[[TiffImages]]
deps = ["ColorTypes", "DataStructures", "DocStringExtensions", "FileIO", "FixedPointNumbers", "IndirectArrays", "Inflate", "OffsetArrays", "PkgVersion", "ProgressMeter", "UUIDs"]
git-tree-sha1 = "991d34bbff0d9125d93ba15887d6594e8e84b305"
uuid = "731e570b-9d59-4bfa-96dc-6df516fadf69"
version = "0.5.3"

[[URIs]]
git-tree-sha1 = "97bbe755a53fe859669cd907f2d96aee8d2c1355"
uuid = "5c2747f8-b7ea-4ff2-ba2e-563bfd36b1d4"
version = "1.3.0"

[[UUIDs]]
deps = ["Random", "SHA"]
uuid = "cf7118a7-6976-5b1a-9a39-7adc72f591a4"

[[Unicode]]
uuid = "4ec0a83e-493e-50e2-b9ac-8f72acf5a8f5"

[[UnicodeFun]]
deps = ["REPL"]
git-tree-sha1 = "53915e50200959667e78a92a418594b428dffddf"
uuid = "1cfade01-22cf-5700-b092-accc4b62d6e1"
version = "0.4.1"

[[Unzip]]
git-tree-sha1 = "34db80951901073501137bdbc3d5a8e7bbd06670"
uuid = "41fe7b60-77ed-43a1-b4f0-825fd5a5650d"
version = "0.1.2"

[[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"

[[Wayland_protocols_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
git-tree-sha1 = "66d72dc6fcc86352f01676e8f0f698562e60510f"
uuid = "2381bf8a-dfd0-557d-9999-79630e7b1b91"
version = "1.23.0+0"

[[XML2_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl", "Libiconv_jll", "Pkg", "Zlib_jll"]
git-tree-sha1 = "1acf5bdf07aa0907e0a37d3718bb88d4b687b74a"
uuid = "02c8fc9c-b97f-50b9-bbe4-9be30ff0a78a"
version = "2.9.12+0"

[[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"

[[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"

[[Xorg_libXau_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
git-tree-sha1 = "4e490d5c960c314f33885790ed410ff3a94ce67e"
uuid = "0c0b7dd1-d40b-584c-a123-a41640f87eec"
version = "1.0.9+4"

[[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"

[[Xorg_libXdmcp_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
git-tree-sha1 = "4fe47bd2247248125c428978740e18a681372dd4"
uuid = "a3789734-cfe1-5b06-b2d0-1dd0d9d62d05"
version = "1.1.3+4"

[[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"

[[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"

[[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"

[[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"

[[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"

[[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"

[[Xorg_libpthread_stubs_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
git-tree-sha1 = "6783737e45d3c59a4a4c4091f5f88cdcf0908cbb"
uuid = "14d82f49-176c-5ed1-bb49-ad3f5cbd8c74"
version = "0.1.0+3"

[[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"

[[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"

[[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"

[[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"

[[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"

[[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"

[[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"

[[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"

[[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"

[[Xorg_xtrans_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
git-tree-sha1 = "79c31e7844f6ecf779705fbc12146eb190b7d845"
uuid = "c5fb5394-a638-5e4d-96e5-b29de1b5cf10"
version = "1.4.0+3"

[[Zlib_jll]]
deps = ["Libdl"]
uuid = "83775a58-1f1d-513f-b197-d71354ab007a"
version = "1.2.12+1"

[[Zstd_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
git-tree-sha1 = "e45044cd873ded54b6a5bac0eb5c971392cf1927"
uuid = "3161d3a3-bdf6-5164-811a-617609db77b4"
version = "1.5.2+0"

[[gdk_pixbuf_jll]]
deps = ["Artifacts", "Glib_jll", "JLLWrappers", "JpegTurbo_jll", "Libdl", "Libtiff_jll", "Pkg", "Xorg_libX11_jll", "libpng_jll"]
git-tree-sha1 = "c23323cd30d60941f8c68419a70905d9bdd92808"
uuid = "da03df04-f53b-5353-a52f-6a8b0620ced0"
version = "2.42.6+1"

[[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"

[[libblastrampoline_jll]]
deps = ["Artifacts", "Libdl", "OpenBLAS_jll"]
uuid = "8e850b90-86db-534c-a0d3-1478176c7d93"
version = "5.0.1+0"

[[libfdk_aac_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
git-tree-sha1 = "daacc84a041563f965be61859a36e17c4e4fcd55"
uuid = "f638f0a6-7fb0-5443-88ba-1cc74229b280"
version = "2.0.2+0"

[[libpng_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg", "Zlib_jll"]
git-tree-sha1 = "94d180a6d2b5e55e447e2d27a29ed04fe79eb30c"
uuid = "b53b4c65-9356-5827-b1ea-8c7a1a84506f"
version = "1.6.38+0"

[[libsixel_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
git-tree-sha1 = "78736dab31ae7a53540a6b752efc61f77b304c5b"
uuid = "075b6546-f08a-558a-be8f-8157d0f608a5"
version = "1.8.6+1"

[[libvorbis_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl", "Ogg_jll", "Pkg"]
git-tree-sha1 = "b910cb81ef3fe6e78bf6acee440bda86fd6ae00c"
uuid = "f27f6e37-5d2b-51aa-960f-b287f2bc3b7a"
version = "1.3.7+1"

[[nghttp2_jll]]
deps = ["Artifacts", "Libdl"]
uuid = "8e850ede-7688-5339-a07c-302acd2aaf8d"
version = "1.41.0+1"

[[p7zip_jll]]
deps = ["Artifacts", "Libdl"]
uuid = "3f19e933-33d8-53b3-aaab-bd5110c3b7a0"
version = "16.2.1+1"

[[x264_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
git-tree-sha1 = "4fea590b89e6ec504593146bf8b988b2c00922b2"
uuid = "1270edf5-f2f9-52d2-97e9-ab00b5d0237a"
version = "2021.5.5+0"

[[x265_jll]]
deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
git-tree-sha1 = "ee567a171cce03570d77ad3a43e90218e38937a9"
uuid = "dfaa095f-4041-5dcd-9319-2fabd8486b76"
version = "3.5.0+0"

[[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"
"""

# ╔═╡ Cell order:
# ╠═670a7c9a-7dd9-11ec-2f9e-5daa27e49f3c
# ╟─8e642f27-ea51-4cf1-b798-d8a586770b57
# ╠═4478d0a8-360a-4792-8859-244d80d2902b
# ╠═0e4ef764-37f3-4db8-b7d5-2d40d81e2797
# ╠═1750a808-0c77-466c-a7e8-e569082c60af
# ╟─c6e8d6fb-7993-4d9e-bf0a-c4dc1d1e4b4a
# ╠═419dfb6d-f7ca-46eb-984a-99d10b804e89
# ╠═b389b734-e61a-4348-9720-3475b654ad8d
# ╠═386615cc-6bd9-4a49-ba70-39f277c0e1db
# ╠═6b6bb9b9-b2c1-409e-9fc1-60285c1d2177
# ╠═85f6c36d-a395-442f-8ba7-f2a8506d1544
# ╟─394fb978-f350-4ddc-9de5-d5114e9a87ed
# ╟─4377d345-71b9-4aff-b8d9-a45c4e9eb05e
# ╟─9bbc975a-d0d8-4f8e-b6b6-85fc2a9823ca
# ╠═b68f03ea-e884-4614-b1ff-c1d92118c30b
# ╟─b4ad67b5-5b71-444b-ad68-5064587905d4
# ╠═d5623e00-571f-407f-a88c-f1a020102e1d
# ╟─dccb5b97-4865-436a-b2cb-e774b2b1d79b
# ╟─daf3699e-267c-4b65-a39a-af0c1978430c
# ╠═280de3c0-6329-4a74-b684-3adebefe033e
# ╠═7bcb0d03-a859-496f-b5f4-4d6ba17ab07a
# ╠═83769dd6-6b73-4821-b87d-f9b3ae12efe6
# ╟─5e91243f-8f61-4076-a1bc-b60b5ce8fb4c
# ╠═2fe9e1dc-81ab-450c-9544-6429c2a80d6b
# ╟─c7cbae3a-248e-4f9a-ba42-b2c94a570665
# ╠═98eb608b-c7f1-45a0-911c-4d5ef1a1fdc5
# ╠═347f8c01-96d6-41d0-a785-15e51886f89d
# ╠═c32aeff2-72d4-4869-aaa7-e9bb26d59ab8
# ╠═7cfd7bc3-fcd6-4034-9d8d-36a87eac1627
# ╟─751e0e94-7899-4c25-8519-6ce98b7b5f24
# ╠═cf2d80f0-9b2a-4162-a1b5-acc41a3dedcc
# ╠═0320566c-c8ca-4696-9aff-4d737be2ea97
# ╠═8802eb09-e619-4d04-ae11-37fb19649fa1
# ╠═cb3e9ead-2e22-4b94-9ac2-8b63808df445
# ╟─c4d5c2f5-6603-4309-8fbf-ed41a8e90675
# ╠═e0452dc4-8bd1-470a-9cfd-782dfc4e301e
# ╠═4d7c4da8-58aa-42af-a8a1-89160e73e08f
# ╠═414a7f26-6c02-4deb-888f-62acbf52238c
# ╠═bf58ce82-c5ec-416b-8d42-d5673c2d156d
# ╠═3ce732b7-d7fd-434b-95f9-f37c33192934
# ╠═2db8f5bc-0b2a-4a1e-9cdb-125af1b001c4
# ╟─449a2e13-fe1f-49f8-a16d-7d160d396aa2
# ╠═02839b7d-9944-4430-909a-cdba85cc64c2
# ╟─703cf884-aa9a-4387-9a1b-afebdae5530e
# ╟─92c2398d-e393-48b3-a795-2a4afc799613
# ╠═8ab423a8-4b33-477c-b445-1a10f86c260e
# ╟─4b0750aa-70df-4bbe-8d90-fafea639f9af
# ╠═c7404c7e-f847-4c27-be83-d8798a3c5dde
# ╟─c3221ca6-7023-4bcc-865a-8c754ff45182
# ╠═a3966b72-930d-4485-94f7-177c3c5a7391
# ╠═dab7a6d7-a679-41f6-ab05-cf0a6cc5418c
# ╟─220b0f16-4569-4187-a4e9-badf4dca946e
# ╠═01b82335-b1bd-457a-8bf8-9ec3de0bd375
# ╠═0fc5c49f-2f2d-419e-a0cb-76a53b641fbe
# ╠═8db87c95-64d6-4176-830a-c906c21125da
# ╟─cc3a57d7-ee3b-4975-9973-bad0bb9f2f53
# ╠═c4c0dd00-655c-419a-9b8d-3afb78f139b1
# ╠═d3d3b840-70ef-4c4d-932f-0c259682d141
# ╠═c0e07072-6c6e-4ec5-9795-bc5d6bb74fe7
# ╠═4973ccbb-ea9b-49c3-b549-7b53c14dd685
# ╠═2b53da0c-c5dc-4788-9f28-36f5bb822559
# ╠═a772b2ce-5ef8-47a4-952b-683796fe3959
# ╟─133d23fa-d7d3-4dc8-b1f4-c8f068a8c07f
# ╟─5adf6861-9a4e-4347-8770-47df895038a6
# ╟─00000000-0000-0000-0000-000000000001
# ╟─00000000-0000-0000-0000-000000000002

max_stars_repo_name: Odilf/Dveco

max_issues_repo_name: Odilf/Dveco

max_forks_repo_name: Odilf/Dveco


text:
using Flux

# throttle, steering, bVx, bVy
Chain(
  Dense(25,8,relu),
  MaxPool(8,4),
  # flatten
  Dense(4,8,relu),
  Dense(8,2)
)

ext: jl


text:
abstract type AbstractRectangle <: AbstractShape end

struct Rectangle{I <: Integer} <: AbstractRectangle
    position::Point{I}
    height::I
    width::I
end

struct FilledRectangle{I <: Integer} <: AbstractRectangle
    position::Point{I}
    height::I
    width::I
end

struct ThickRectangle{I <: Integer} <: AbstractRectangle
    position::Point{I}
    height::I
    width::I
    thickness::I
end

#####
##### AbstractRectangle
#####

is_valid(shape::AbstractRectangle) = shape.height > zero(shape.height) && shape.width > zero(shape.width)

get_i_min(shape::AbstractRectangle) = shape.position.i
get_i_max(shape::AbstractRectangle) = shape.position.i + shape.height - one(shape.height)

get_j_min(shape::AbstractRectangle) = shape.position.j
get_j_max(shape::AbstractRectangle) = shape.position.j + shape.width - one(shape.width)

get_drawing_optimization_style(::AbstractRectangle) = CHECK_BOUNDS

#####
##### Rectangle
#####

function draw!(f::F, image::AbstractMatrix, shape::Rectangle, color) where {F <: Function}
    @assert is_valid(shape) "Cannot draw invalid shape $(shape)"

    i_min, i_max = get_i_extrema(shape)
    j_min, j_max = get_j_extrema(shape)

    draw!(f, image, VerticalLine(i_min, i_max, j_min), color)
    draw!(f, image, HorizontalLine(i_min, j_min, j_max), color)
    draw!(f, image, HorizontalLine(i_max, j_min, j_max), color)
    draw!(f, image, VerticalLine(i_min, i_max, j_max), color)

    return nothing
end

#####
##### FilledRectangle
#####

function clip(image::AbstractMatrix, shape::FilledRectangle)
    i_min_shape, i_max_shape = get_i_extrema(shape)
    i_min_image, i_max_image = get_i_extrema(image)

    I = typeof(i_min_shape)

    j_min_shape, j_max_shape = get_j_extrema(shape)
    j_min_image, j_max_image = get_j_extrema(image)

    i_min = max(i_min_shape, i_min_image)
    i_max = min(i_max_shape, i_max_image)

    j_min = max(j_min_shape, j_min_image)
    j_max = min(j_max_shape, j_max_image)

    return FilledRectangle(Point(i_min, j_min), i_max - i_min + one(I), j_max - j_min + one(I))
end

get_drawing_optimization_style(::FilledRectangle) = CLIP

function draw!(f::F, image::AbstractMatrix, shape::FilledRectangle, color) where {F <: Function}
    @assert is_valid(shape) "Cannot draw invalid shape $(shape)"

    i_min, i_max = get_i_extrema(shape)
    j_min, j_max = get_j_extrema(shape)

    for j in j_min:j_max
        for i in i_min:i_max
            f(image, i, j, color)
        end
    end

    return nothing
end

#####
##### ThickRectangle
#####

function is_valid(shape::ThickRectangle)
    height = shape.height
    width = shape.width
    thickness = shape.thickness

    I = typeof(height)

    return height > zero(I) && width > zero(I) && thickness > zero(I) && thickness <= min(height, width)
end

function draw!(f::F, image::AbstractMatrix, shape::ThickRectangle, color) where {F <: Function}
    @assert is_valid(shape) "Cannot draw invalid shape $(shape)"

    position = shape.position
    height = shape.height
    width = shape.width
    thickness = shape.thickness

    i_min = position.i
    j_min = position.j

    draw!(f, image, FilledRectangle(position, height, thickness), color)
    draw!(f, image, FilledRectangle(position, thickness, width), color)
    draw!(f, image, FilledRectangle(Point(i_min + height - thickness, j_min), thickness, width), color)
    draw!(f, image, FilledRectangle(Point(i_min, j_min + width - thickness), height, thickness), color)

    return nothing
end


text:
# Section 4.4 of
# A. A. Ahmadi, and A. Majumdar
# DSOS and SDSOS Optimization: More Tractable Alternatives to Sum of Squares and Semidefinite Optimization
# 2017

using MultivariateMoments

function options_pricing_test(optimizer, config::MOIT.Config,
                              cone::SumOfSquares.PolyJuMP.PolynomialSet,
                              K::Int, expected::Float64)
    VERSION < v"1.0.3" && return # see https://github.com/jump-dev/SumOfSquares.jl/issues/48
    atol = 100config.atol
    rtol = 100config.rtol

    @polyvar x y z
    σ = [184.04, 164.88, 164.88, 184.04, 164.88, 184.04]
    X = [x^2, x*y, x*z, y^2, y*z, z^2, x, y, z, 1]
    μ = measure([σ .+ 44.21^2; 44.21 * ones(3); 1], X)

    cocone = CopositiveInner(cone)

    model = _model(optimizer)
    @variable(model, p, Poly(X))
    @constraint(model, p in cocone)
    @constraint(model, p - (x - K) in cocone)
    @constraint(model, p - (y - K) in cocone)
    @constraint(model, p - (z - K) in cocone)
    @objective(model, Min, dot(μ, p))
    if MOI.Test._supports(config, MOI.optimize!)
        JuMP.optimize!(model)
        @test JuMP.primal_status(model) == MOI.FEASIBLE_POINT
        @test JuMP.objective_value(model) ≈ expected atol=atol rtol=rtol
    end
end

function options_pricing_test(optimizer, config::MOIT.Config,
                              cone::SumOfSquares.PolyJuMP.PolynomialSet,
                              Ks::Vector{Int},
                              expected::Vector{Float64})
    @testset "K = $K" for (K, exp) in zip(Ks, expected)
        options_pricing_test(optimizer, config, cone, K, exp)
    end
end

const K = [30, 35, 40, 45, 50]
const dsos_codsos_exp   = [132.63, 132.63, 132.63, 132.63, 132.63]
const sdsos_cosdsos_exp = [ 21.51,  17.17,  13.20,   9.85,   7.30]

sos_options_pricing_test(optimizer, config)   = options_pricing_test(optimizer, config, SOSCone(), K, sdsos_cosdsos_exp)
sd_tests["sos_options_pricing"] = sos_options_pricing_test
sdsos_options_pricing_test(optimizer, config) = options_pricing_test(optimizer, config, SDSOSCone(), K, sdsos_cosdsos_exp)
soc_tests["sdsos_options_pricing"] = sdsos_options_pricing_test
dsos_options_pricing_test(optimizer, config)  = options_pricing_test(optimizer, config, DSOSCone(), K, dsos_codsos_exp)
linear_tests["dsos_options_pricing"] = dsos_options_pricing_test

max_stars_repo_name: blegat/SumOfSquares

max_issues_repo_name: blegat/SumOfSquares

max_forks_repo_name: blegat/SumOfSquares


text:
###################################################
## TimingsInference.jl
## Package usings, includes and exports
###################################################

module TimingsInference

using RoutingNetworks, LightGraphs, NearestNeighbors
using Distributions, DataStructures
using JuMP, Gurobi, Mosek, MathProgBase
using Dates, Printf, SparseArrays, Distributed, Random, LinearAlgebra, Statistics
using PyPlot, Colors

import RoutingNetworks: inPolygon
import Base.<

#network
export NetworkTimings, NetworkTrip, NetworkData
export noisyVirtualData, perfectVirtualData, simpleVirtualData
export findNetworkDependence, simplifyPath, evaluateTimes, getFullPathEdges
export findNearEdges, findNearEdgesSameType, clusterEdges
export centerCongestion, twoCongestions, squareCongestion, gradientCongestion
export typicalSpeeds
export yenKSP

#iterative
export IterativeState, doIteration!, LimitedPaths, heuristicPaths, GreedyEdges, updateState!
export lpTimes, fraclpTimes, mipTimes, heuristicTimes, socpTimes, constantSpeedTimes
export socpTimesCont, socpTimesContNbhd

#geo data
export GeoTrip, GeoData, fromNYCTaxiCSV, removeOutliers, isRegular,  tripDistance, stats
export inTimeWindow, onlyWeekdays, inPolygon, getPolygon, DataSplit, RandomSplit, LocSplit, trainSet
export testSet, GeoStats, RealGeoStats, NetworkStats, RealNetworkStats, VirtNetworkStats
export printStats, plotStats
#projectors
export NetworkProjector, NearestNode, AvgRadius, NearestRoad
#geo-timings
export GeoTimings, TimingsFromNetwork, estimateTime, KnnTimings
#visualization
export CompareTimes

include("network/networktypes.jl")

include("network/tools/edges.jl")
include("network/tools/graphdep.jl")

include("network/synthetic/virtualdata.jl")
include("network/synthetic/congestion.jl")

include("network/iterative/state.jl")
include("network/iterative/limitedpaths.jl")
include("network/iterative/greedy.jl")

include("network/stats/statistics.jl")
include("network/stats/stats.jl")

include("network/solvers/lp.jl")
include("network/solvers/socp.jl")
include("network/solvers/fraclp.jl")
include("network/solvers/mip.jl")
include("network/solvers/heuristic.jl")
include("network/solvers/roadtypespeed.jl")
include("network/solvers/pathchoice.jl")

include("network/tools/speeds.jl")
include("network/tools/yenksp.jl")

include("geo/trips.jl")
include("geo/data/nyctaxi.jl")
include("geo/data/filters.jl")

include("geo/projectors/networkprojector.jl")
include("geo/projectors/nearestnode.jl")
include("geo/projectors/avgradius.jl")
include("geo/projectors/nearestroad.jl")

include("geo/timings/geotimings.jl")

include("geo/stats/datasplit.jl")
include("geo/stats/statistics.jl")
include("geo/stats/stats.jl")

include("geo/timings/fromnetwork.jl")
include("geo/timings/knn.jl")

end # module


text:
using CardinalDicts
using Base.Test

#=
  create an CardinalDict with indices 1:20 that holds Int64 values
  check length
  confirm keymax
  populate it
  use it
  unset an index
  check it
  reassign an indexable value
=#

factorials = CardinalDict{Int64}(20);

@test length(factorials) == 0
@test keymax(factorials) == 20

for i in 1:20
    setindex!(factorials, factorial(i), i)
end

@test length(factorials) == 20
@test keymax(factorials) == 20

@test haskey(factorials, 17) == true
@test factorials[17] == factorial(17)

delete!(factorials, 17)
@test haskey(factorials, 17) == false
@test get(factorials, 17, 0) == 0

factorials[17] = factorial(17)
@test factorials[17] == factorial(17)

vec = [1.0, 3.0, 2.0]
fromvec = CardinalDict(vec)
@test fromvec[2] == 3.0
@test keymax(fromvec) == 3
@test isfull(fromvec) == true

# Container, Associative interface

tenfold = CardinalDict{String}(40)
@test length(tenfold) == 0
@test endof(tenfold) == 0
@test keymax(tenfold) == 40
@test keys(tenfold) == []
@test values(tenfold) == []

tenfold[20] = "200"
tenfold[25] = "250"
tenfold[26] = "260"

@test length(tenfold) == 3
@test endof(tenfold) == 26
@test keymax(tenfold) == 40
@test keys(tenfold) == Int8[20, 25, 26]
@test values(tenfold) == String["200", "250", "260"]
@test eltype(tenfold) == Pair{Int8, String}

delete!(tenfold, 20)
@test haskey(tenfold, 20) == false
@test get(tenfold, 20, "0") == "0"

tenfold[20] = "200"
@test haskey(tenfold, 20) == true
@test get(tenfold, 20, "0") == "200"

# etc

tenfold2 = eval(parse(string(tenfold)))
@test tenfold == tenfold2


text:
function compute()

	TOTAL=200

	ways=fill(0,TOTAL+1)
	ways[1]=1

	for coin in [1,2,5,10,20,50,100,200]
		for i =1: length(ways)-coin
			ways[i+coin]+=ways[i]
		end
	end

	return ways[length(ways)]

end

println(compute())


text:
using ContinuousTimeMarkov, Test, LinearAlgebra, SparseArrays, Distributions, StatsBase

import Random

Random.seed!(1)

include("test_transition_matrices.jl")
include("test_competing_poisson.jl")


text:
module TuringGLM

import Base: ==, length
using LinearAlgebra: I
using Statistics: mean, median, std
using StatsBase: mad, coefnames
using StatsModels: ConstantTerm, FormulaTerm, FunctionTerm, Term
using StatsModels: hasintercept, response, term
using MixedModels: MixedModel, RandomEffectsTerm

# Different modelmatrix
using MixedModels: MixedModels
using StatsModels: StatsModels

# Tables API stuff
using Tables: Tables
using TableOperations: TableOperations

# LazyArrays
using LazyArrays: @~, LazyArray

# Distributions
using Distributions: UnivariateDistribution

using Reexport: @reexport

@reexport begin
    using Turing
    using StatsModels: @formula
    using Distributions:
        Bernoulli,
        Binomial,
        Gamma,
        InverseGaussian,
        NegativeBinomial,
        Normal,
        Poisson,
        TDist
end

include("utils.jl")
include("data_constructors.jl")
include("priors.jl")
include("turing_model.jl")

export turing_model
export CustomPrior, DefaultPrior

end # module


text:
using DashBootstrapComponents, DashHtmlComponents

toast = dbc_toast(
    [html_p("This is the content of the toast", className = "mb-0")],
    header = "This is the header",
);

max_stars_repo_name: glsdown/dash-bootstrap-components

max_issues_repo_name: glsdown/dash-bootstrap-components

max_forks_repo_name: glsdown/dash-bootstrap-components


text:
module CoinGeckoAPI

# dependacies

using HTTP
using JSON3

# exports
export ping
export get_price
export get_token_price
export get_supported_vs_currencies
export get_coins_list
export get_coins_markets
export get_coin_by_id
export get_coin_ticker_by_id
export get_coin_history_by_id
export get_coin_market_chart_by_id
export get_coin_market_chart_range_by_id
export get_coin_status_updates_by_id
export get_coin_ohlc_by_id
export get_coin_info_from_contract_address_by_id
export get_coin_market_chart_from_contract_address_by_id
export get_coin_market_chart_range_from_contract_address_by_id
export get_asset_platforms
export get_coins_categories_list
export get_coins_categories
export get_exchanges_list
export get_exchanges_id_name_list
export get_exchanges_by_id
export get_exchanges_tickers_by_id
export get_exchanges_status_updates_by_id
export get_exchanges_volume_chart_by_id
export get_finance_platforms
export get_finance_products

# globals constants
const global url_base = "https://api.coingecko.com/api/v3/"

"""
    ping()
Function to check if the API is up and running

# Arguments
    None

# Returns
"""
function ping()
    apiurl ="ping";
    r =  HTTP.request("GET", url_base * apiurl);
    r = String(r.body);
    return JSON3.read(r);

end

"""
    _api_url_params(api_url, params, api_url_has_params = false)

Function to build the API url with the parameters

# Arguments
    `api_url: string` : the API url
    `params: dict` : the parameters to be added to the API url
    `api_url_has_params: bool` : if the API url has parameters

# Returns
    `api_url: String` : the API url with the parameters

"""
function _api_url_params(api_url, params, api_url_has_params = false)

    if api_url_has_params == true
        api_url = api_url * "&";
    else
        api_url = api_url * "?";
    end

    for a_pair in params
        api_url *=  a_pair.first * "=" * a_pair.second * "&";
    end

    return api_url[1:end-1];
end

"""
    get_price(ids, vs_currencies, kargs...)

Get the price of a cryptocurrency in any other supported currencies needed

# Arguments
    `ids: list` : the cryptocurrency ids
    `vs_currencies: list` : the currencies to get the price in
    `kargs: dict` : the parameters to be added to the API url

# Returns
    `prices: dict` : the prices of the cryptocurrencies in the currencies

# Example:

```julia
r = get_price("bitcoin", "usd", "include_market_cap" => "true", "include_24hr_vol" => "true", "include_24hr_change" => "true")
```

"""
function get_price(ids, vs_currencies, kargs...)
    apiurl = "simple/price";
    kwards= Dict(kargs);
    kwards["ids"] = ids;
    kwards["vs_currencies"] = vs_currencies;

    api_url = _api_url_params(apiurl, kwards);
    r =  HTTP.request("GET", url_base * api_url);
    r = String(r.body);
    return JSON3.read(r);

end

function get_token_price(id, contract_addresses, vs_currencies, kargs...)
    apiurl = "simple/token_price";
    kwards= Dict(kargs);
    kwards["id"] = id;
    kwards["contract_addresses"] = contract_addresses;
    kwards["vs_currencies"] = vs_currencies;

    api_url = _api_url_params(apiurl, kwards);
    r =  HTTP.request("GET", url_base * api_url);
    r = String(r.body);
    return JSON3.read(r);

end

"""
    get_supported_vs_currencies()

Get list of supported_vs_currencies

# Arguments
    None

# Returns
    `supported_vs_currencies: Dict` : the list of supported_vs_currencies
"""
function get_supported_vs_currencies()
    apiurl = "simple/supported_vs_currencies";
    r =  HTTP.request("GET", url_base * apiurl);
    r = String(r.body);
    return JSON3.read(r);

end

"""
    get_coins(kargs...)

List all coins with data (name, price, market, developer, community, etc)

# Arguments
    `kargs: dict` : the parameters to be added to the API url

# Returns
    `coins: Dict` : the list of coins with id, name and symbol

"""
function get_coins_list(kargs...)
    apiurl = "coins/list";
    kwards= Dict(kargs);

    api_url = _api_url_params(apiurl, kwards);
    r =  HTTP.request("GET", url_base * api_url);
    r = String(r.body);
    return JSON3.read(r);

end

"""
    get_coins_markets(vs_currencies, kargs...)

Use this to obtain all the coins market data (price, market cap, volume)

# Arguments
    `vs_currencies: list` : the currencies to get the price in
    `kargs: dict` : the parameters to be added to the API url

# Returns
    `coins_markets: Dict` : the list of coins with market data

# Example:

```julia
r = get_coins_markets("usd")

````
"""
function get_coins_markets(vs_currency, kargs...)
    apiurl = "coins/markets";
    kwards= Dict(kargs);
    kwards["vs_currency"] = vs_currency;

    api_url = _api_url_params(apiurl, kwards);
    r =  HTTP.request("GET", url_base * api_url);
    r = String(r.body);
    return JSON3.read(r);

end

"""
    get_coin_by_id(id, kargs...)

Get current data (name, price, market, including exchange tickers for a coin)

# Arguments
    `id: string` : the coin id
    `kargs: dict` : the parameters to be added to the API url

# Returns
    `coin: Dict` : the current data for a coin

"""
function get_coin_by_id(id, kargs...)
    apiurl = "coins";
    kwards= Dict(kargs);
    kwards["id"] = id;

    api_url = _api_url_params(apiurl, kwards);
    r =  HTTP.request("GET", url_base * api_url);
    r = String(r.body);
    return JSON3.read(r);

end

"""
    get_coin_ticker_by_id(id, kargs...)

Get coin tickers (paginated to 100 items)

# Arguments
    `id: string` : the coin id
    `kargs: dict` : the parameters to be added to the API url

# Returns
    `coin_ticker: Dict` : the coin tickers

"""
function get_coin_ticker_by_id(id, kargs...)
    apiurl = "coins/$id/tickers";
    kwards= Dict(kargs);
    kwards["id"] = id;

    api_url = _api_url_params(apiurl, kwards);
    r =  HTTP.request("GET", url_base * api_url);
    r = String(r.body);
    return JSON3.read(r);

end

"""
    get_coin_history_by_id(id, date, kargs...)

Get Historical data (name, price, market, stats ) at a given date for a coin.

# Arguments
    `id: string` : the coin id
    `date: string` : the date to get the data for
    `kargs: dict` : the parameters to be added to the API url

# Returns
    `coin_history: Dict` : the historical data for a coin

# Example:
```julia
r = get_coin_history_by_id("bitcoin", "30-12-2017")

````
"""
function get_coin_history_by_id(id, date, kargs...)
    apiurl = "coins/$id/history";
    kwards= Dict(kargs);
    kwards["id"] = id;
    kwards["date"] = date;

    api_url = _api_url_params(apiurl, kwards);
    r =  HTTP.request("GET", url_base * api_url);
    r = String(r.body);
    return JSON3.read(r);

end

"""
    get_coin_market_chart_by_id(id, vs_currency, days, kargs...)

Get historical market data include price, market cap, and 24h volume (granularity auto)
Minutely data will be used for duration within 1 day, Hourly data will be used for duration between 1 day and 90 days,
    Daily data will be used for duration above 90 days.

# Arguments
    `id: string` : the coin id
    `vs_currency: string` : the currency to get the data in
    `days: String` : the number of days to get the data for
    `kargs: dict` : the parameters to be added to the API url

# Returns
    `coin_market_chart: Dict` : the historical market data for a coin

# Example:
```julia
r = get_coin_market_chart_by_id("bitcoin", "usd", "30")
```

"""
function get_coin_market_chart_by_id(id, vs_currency, days, kargs...)
    apiurl = "coins/$id/market_chart";
    kwards= Dict(kargs);
    kwards["id"] = id;
    kwards["vs_currency"] = vs_currency;
    kwards["days"] = days;

    api_url = _api_url_params(apiurl, kwards);
    r =  HTTP.request("GET", url_base * api_url);
    r = String(r.body);
    return JSON3.read(r);

end

"""
    get_coin_market_chart_range_by_id(id, vs_currency, from_timestamp, to_timestamp, kargs...)

Get historical market data include price, market cap, and 24h volume within a range of timestamp (granularity auto)
Data granularity is automatic (cannot be adjusted)
1 day from query time = 5 minute interval data
1 - 90 days from query time = hourly data
above 90 days from query time = daily data (00:00 UTC)

# Arguments
    `id: string` : the coin id
    `vs_currency: string` : the currency to get the data in
    `from_timestamp: String` : the timestamp to start the data from (UNIX Timestamp (eg 1392577232))
    `to_timestamp: String` : the timestamp to end the data from (UNIX Timestamp (eg 1422577232))
    `kargs: dict` : the parameters to be added to the API url

# Returns
    `coin_market_chart: Dict` : the historical market data for a coin

# Example:
```julia
r = get_coin_market_chart_range_by_id("bitcoin", "usd", "1392577232", "1422577232")
```
"""
function get_coin_market_chart_range_by_id(id, vs_currency, from_timestamp, to_timestamp, kargs...)
    apiurl = "coins/$id/market_chart/range";
    kwards= Dict(kargs);

    kwards["vs_currency"] = vs_currency;
    kwards["from"] = from_timestamp;
    kwards["to"] = to_timestamp;

    api_url = _api_url_params(apiurl, kwards);

    r =  HTTP.request("GET", url_base * api_url);
    r = String(r.body);
    return JSON3.read(r);

end

"""
    get_coin_status_updates_by_id(id, kargs...)

Get status updates for a given coin.

# Arguments
    `id: string` : the coin id  (eg bitcoin)
    `kargs: dict` : the parameters to be added to the API url

# Returns
    `coin_status_updates: Dict` : the status updates for a coin

"""
function get_coin_status_updates_by_id(id, kargs...)
    apiurl = "coins/$id/status_updates";
    kwards= Dict(kargs);
    kwards["id"] = id;

    api_url = _api_url_params(apiurl, kwards);
    r =  HTTP.request("GET", url_base * api_url);
    r = String(r.body);
    return JSON3.read(r);

end

"""
    get_coin_ohlc_by_id(id, vs_currency, days)

Get coin's ohlc
Candles's body:
    1-2 days: 30 minutes
    3-30 days: 4 hours
    31 and before: 4 days

# Arguments
    `id: string` : the coin id  (eg bitcoin)
    `vs_currency: string` : the currency to get the data in
    `days: String` : the number of days to get the data for
    `kargs: dict` : the parameters to be added to the API url

# Returns
    `coin_ohlc: Dict` : the ohlc data for a coin

# Example:
```julia
r = get_coin_ohlc_by_id("bitcoin", "usd", "30")
```
"""
function get_coin_ohlc_by_id(id, vs_currency, days)
    apiurl = "coins/$id/ohlc";
    kwards= Dict();
    kwards["vs_currency"] = vs_currency;
    kwards["days"] = days;

    api_url = _api_url_params(apiurl, kwards);
    r =  HTTP.request("GET", url_base * api_url);
    r = String(r.body);
    return JSON3.read(r);

end

"""
    get_coin_info_from_contract_address_by_id(id, contract_address)

Get coin info from contract address

# Arguments
    `id: string` : Asset platform (See asset_platforms endpoint for list of options)
    `contract_address: string` : the token's contract address

# Returns
    `coin_info: Dict` : the coin info from contract address

# Example:
```julia
get_coin_info_from_contract_address_by_id("ethereum", "0xc00e94cb662c3520282e6f5717214004a7f26888")
```

"""
function get_coin_info_from_contract_address_by_id(id, contract_address)
    apiurl = "coins/$id/contract/$contract_address";
    kwards= Dict();

    api_url = _api_url_params(apiurl, kwards);
    r =  HTTP.request("GET", url_base * api_url);
    r = String(r.body);
    return JSON3.read(r);

end

"""
    get_coin_market_chart_from_contract_address_by_id(id, contract_address, vs_currency, days)

Get historical market data include price, market cap and 24h volume (granularity auto)

# Arguments
    `id: string` : The id of the platform issuing tokens (See asset_platforms endpoint for list of options)
    `contract_address: string` : Token's contract address
    `vs_currency: string` : the currency to get the data in
    `days: String` : Data up to number of days ago (eg. 1,14,30,max)

# Returns
    `coin_market_chart: Dict` : the historical market data for a coin

# Example:
```julia
r = get_coin_market_chart_from_contract_address_by_id("ethereum", "0xc00e94cb662c3520282e6f5717214004a7f26888", "usd", "1")
```
"""
function get_coin_market_chart_from_contract_address_by_id(id, contract_address, vs_currency, days)
    apiurl = "coins/$id/contract/$contract_address/market_chart";
    kwards= Dict();
    kwards["vs_currency"] = vs_currency;
    kwards["days"] = days;

    api_url = _api_url_params(apiurl, kwards);
    r =  HTTP.request("GET", url_base * api_url);
    r = String(r.body);
    return JSON3.read(r);

end

"""
    get_coin_market_chart_range_from_contract_address_by_id(id, contract_address, vs_currency, from_timestamp, to_timestamp)

Get historical market data include price, market cap, and 24h volume within a range of timestamp (granularity auto)

# Arguments
    `id: string` : The id of the platform issuing tokens (See asset_platforms endpoint for list of options)
    `contract_address: string` : Token's contract address
    `vs_currency: string` : the currency to get the data in
    `from_timestamp: String` : From date in UNIX Timestamp (eg. 1392577232)
    `to_timestamp: String` : To date in UNIX Timestamp (eg. 1422577232)

# Returns
    `coin_market_chart: Dict` : the historical market data for a coin

"""
function get_coin_market_chart_range_from_contract_address_by_id(id, contract_address, vs_currency, from_timestamp, to_timestamp)
    apiurl = "coins/$id/contract/$contract_address/market_chart/range";
    kwards= Dict();
    kwards["vs_currency"] = vs_currency;
    kwards["from"] = from_timestamp;
    kwards["to"] = to_timestamp;

    api_url = _api_url_params(apiurl, kwards);
    r =  HTTP.request("GET", url_base * api_url);
    r = String(r.body);
    return JSON3.read(r);

end

"""
    get_asset_platforms()

List all the asset platforms (blockchain networks)

# Arguments
    None

# Returns
    `asset_platforms: Dict` : the asset platforms
"""
function get_asset_platforms()
    apiurl = "asset_platforms";

    r =  HTTP.request("GET", url_base * apiurl);
    r = String(r.body);
    return JSON3.read(r);

end

"""
    get_coins_categories_list()

List all the coins categories

# Arguments
    None

# Returns
    `coins_categories: Dict` : the coins categories
"""
function get_coins_categories_list()
    apiurl = "coins/categories/list";

    r =  HTTP.request("GET", url_base * apiurl);
    r = String(r.body);
    return JSON3.read(r);

end

"""
    get_coins_categories(kargs...)

List all categories with market data.

# Arguments
    `kargs: Dict` : the parameters to be added to the API url
    valid values: market_cap_desc (default), market_cap_asc,
    name_desc, name_asc, market_cap_change_24h_desc and market_cap_change_24h_asc

# Returns
    `coins_categories: Dict` : the coins categories

# Example:
```julia
get_coins_categories()
get_coins_categories("market_cap_desc" => "true")

"""
function get_coins_categories(kargs...)
    apiurl = "coins/categories";
    kwards= Dict(kargs);
    api_url = _api_url_params(apiurl, kwards);
    r =  HTTP.request("GET", url_base * api_url);
    r = String(r.body);
    return JSON3.read(r);
end

"""
    function get_exchanges_list(kargs...)

List all exchanges.

# Arguments
    `kargs: Dict` : the parameters to be added to the API url
    valid values:
    per_page : number of results per page 1..250. Total results per page. Default value is 100
    page : Page through results.

# Returns
    `exchanges_list: Dict` : the exchanges list

# Example:

```julia
get_exchanges_list("per_page" => "10")
```
"""
function get_exchanges_list(kargs...)
    apiurl = "exchanges";
    kwards= Dict(kargs);
    api_url = _api_url_params(apiurl, kwards);
    r =  HTTP.request("GET", url_base * api_url);
    r = String(r.body);
    return JSON3.read(r);
end

"""
    get_exchanges_id_name_list()

Use this to obtain all the markets' id in order to make API calls

# Arguments
    None

# Returns
    `exchanges_id_name_list: Dict` : the exchanges id name list

"""
function get_exchanges_id_name_list()
    apiurl = "exchanges/list";

    r =  HTTP.request("GET", url_base * apiurl);
    r = String(r.body);
    return JSON3.read(r);

end

"""
    get_exchanges_by_id(id)

Get exchange volume in BTC and top 100 tickers

# Arguments
    `id: string` : the id of the exchange

# Returns
    `exchanges_by_id: Dict` : the exchange data
"""
function get_exchanges_by_id(id)
    apiurl = "exchanges/$id";

    r =  HTTP.request("GET", url_base * apiurl);
    r = String(r.body);
    return JSON3.read(r);

end

"""
    get_exchanges_tickers_by_id(id, kargs... )

Get exchange tickers (paginated)

# Arguments
    `id: string` : the id of the exchange (can be obtained from exchanges/list eg.binance)
    `kargs: Dict` : the parameters to be added to the API url

# Returns
    `exchanges_tickers_by_id: Dict` : the exchange tickers

"""
function get_exchanges_tickers_by_id(id, kargs... )
    apiurl = "exchanges/$id/tickers";
    kwards= Dict(kargs);
    api_url = _api_url_params(apiurl, kwards);
    r =  HTTP.request("GET", url_base * api_url);
    r = String(r.body);
    return JSON3.read(r);

end

"""
    get_exchanges_status_updates_by_id(id, kargs... )

Get status updates for a given exchange

# Arguments
    `id: string` : the id of the exchange (can be obtained from exchanges/list eg.binance)
    `kargs: Dict` : the parameters to be added to the API url

# Returns
    `exchanges_status_updates_by_id: Dict` : the exchange status updates

# Example:
```julia
get_exchanges_status_updates_by_id("binance")
```
"""
function get_exchanges_status_updates_by_id(id, kargs... )
    apiurl = "exchanges/$id/status_updates";
    kwards= Dict(kargs);
    api_url = _api_url_params(apiurl, kwards);
    r =  HTTP.request("GET", url_base * api_url);
    r = String(r.body);
    return JSON3.read(r);

end

"""
    get_exchanges_volume_chart_by_id(id, days)

Get volume_chart data for a given exchange

# Arguments
    `id: string` : the id of the exchange (can be obtained from exchanges/list eg.binance)
    `days: String` : data up to number of days ago

# Returns
    `exchanges_volume_chart_by_id: Dict` : the exchange volume chart
"""
function get_exchanges_volume_chart_by_id(id, days)
    apiurl = "exchanges/$id/volume_chart";
    kwards= Dict();
    kwards["days"] = days;
    api_url = _api_url_params(apiurl, kwards);
    r =  HTTP.request("GET", url_base * api_url);
    r = String(r.body);
    return JSON3.read(r);

end

"""
    function get_finance_platforms(kargs...)

List all finance platforms

# Arguments
    `kargs: Dict` : the parameters to be added to the API url
    valid values:
    per_page : number of results per page 1..250. Total results per page. Default value is 100
    page : Page of results (paginated to 100 by default)

# Returns
    `finance_platforms: Dict` : the finance platforms

"""
function get_finance_platforms(kargs...)
    apiurl = "finance_platforms";
    kwards= Dict(kargs);
    api_url = _api_url_params(apiurl, kwards);
    r =  HTTP.request("GET", url_base * api_url);
    r = String(r.body);
    return JSON3.read(r);

end

"""
    get_finance_products(kargs...)

List all finance finance_products

# Arguments
    `kargs: Dict` : the parameters to be added to the API url

# Returns
    `finance_products: Dict` : the finance products
"""
function get_finance_products(kargs...)
    apiurl = "finance_products";
    kwards= Dict(kargs);
    api_url = _api_url_params(apiurl, kwards);
    r =  HTTP.request("GET", url_base * api_url);
    r = String(r.body);
    return JSON3.read(r);

end

end


text:
export sum_autobatch

"""
    sum_autobatch(v)

it's equivalent to `sum(v, dims=1)` except for the case of a vector, where it
will return a scalar instead of a 1-element vector.
"""
sum_autobatch(v) = sum(v, dims=1)
sum_autobatch(v::AbstractVector) = sum(v)

"""
    _batched_outer_prod!(R, vb, wb)

Efficiently performs the outer product R[i] .= vb[i] .* wb[i]'
along the batch dimension i, assuming that the batch dimension
is the last.

R[i] is a matrix, vb[i] and wb[i] are vectors.

Internally uses the fact that R is a StridedView.
"""
@inline function _batched_outer_prod!(R::StridedView, vb, wb)
    @inbounds for i=1:size(R, 3)
        @inbounds for j=1:size(wb, 1)
            @simd for k= @inbounds 1:size(vb, 1)
                @inbounds R[k,j,i] = vb[k,i]*conj(wb[j,i])
            end
        end
    end
    return R
end

@inline function _batched_outer_prod_noconj!(R::StridedView, vb, wb)
    @inbounds for i=1:size(R, 3)
        @inbounds for j=1:size(wb, 1)
            @simd for k= @inbounds 1:size(vb, 1)
                @inbounds R[k,j,i] = vb[k,i]*wb[j,i]
            end
        end
    end
    return R
end

@inline function _batched_outer_prod!(R::StridedView, α, vb, wb)
    @inbounds for i=1:size(R, 3)
        @inbounds for j=1:size(wb, 1)
            @simd for k= @inbounds 1:size(vb, 1)
                @inbounds R[k,j,i] = α * vb[k,i]*conj(wb[j,i])
            end
        end
    end
    return R
end

@inline function _batched_outer_prod_∑!(R::StridedView, α, vb, wb, vb2, wb2)
    @inbounds for i=1:size(R, 3)
        @inbounds for j=1:size(wb, 1)
            @simd for k= @inbounds 1:size(vb, 1)
                @inbounds R[k,j,i] = α * (vb[k,i]*conj(wb[j,i]) + vb2[k,i]*conj(wb2[j,i]))
            end
        end
    end
    return R
end

@inline function _batched_outer_prod_Δ!(R::StridedView, α, vb, wb, vb2, wb2)
    @inbounds for i=1:size(R, 3)
        @inbounds for j=1:size(wb, 1)
            @simd for k= @inbounds 1:size(vb, 1)
                @inbounds R[k,j,i] = α * (vb[k,i]*conj(wb[j,i]) - vb2[k,i]*conj(wb2[j,i]))
            end
        end
    end
    return R
end


text:
module Roaring

import Base.insert!

include("roaring_containers.jl")

end


text:
using CGP
using Logging
using ArgParse

include("../graph_utils.jl")

function read_data(dfile::String="data/lif.data")
    df = open(dfile, "r")
    data = readdlm(dfile, ' ')
    ntrain = Int64(round(size(data, 1) * 3.0 / 4.0))
    ntest = size(data, 1) - ntrain
    training = data[1:ntrain, :]
    test = data[ntrain+(1:ntest), :]
    5, 5, training', test'
end

function regression(c::Chromosome, data::Array{Float64}, nin::Int64, nout::Int64)
    error = 0
    nsamples = size(data, 2)
    for d in 1:nsamples
        outputs = process(c, data[1:nin, d])
        for p in eachindex(outputs)
            error += (outputs[p] - data[nin+p, d])^2
        end
    end
    error /= nsamples
    -error
end

s = ArgParseSettings()
@add_arg_table s begin
    "--seed"
    arg_type = Int
    default = 0
    "--log"
    arg_type = String
    required = true
end

args = parse_args(s)

CGP.Config.init("cfg/base.yaml")
CGP.Config.init("cfg/functions.yaml")

srand(args["seed"])
Logging.configure(filename=args["log"], level=INFO)
nin, nout, train, test = read_data()

fit = x->regression(x, train, nin, nout)
refit = x->regression(x, test, nin, nout)
maxfit, best = oneplus(CGPChromo, nin, nout, fit; seed=args["seed"],
                       record_best=true, record_fitness=refit)


text:
"""
Simply-supported circular plate with uniform distributed load
"""
module ss_circular_plate_udl_examples

using FinEtools
using FinEtoolsDeforLinear
using FinEtoolsFlexStructures.FESetShellT3Module: FESetShellT3
using FinEtoolsFlexStructures.FEMMShellT3FFModule
using FinEtoolsFlexStructures.RotUtilModule: initial_Rfield, update_rotation_field!
using VisualStructures: plot_nodes, plot_midline, render, plot_space_box, plot_midsurface, space_aspectratio, save_to_json

function _execute(mesh_procedure = :q4_t3, n = 2, t_radius_ratio = 0.01, visualize = true)
    E = 200*phun("GPa");
    nu = 0.3;
    a = 1.0*phun("m");
    thickness = a * t_radius_ratio;
    D = E*thickness^3/12/(1-nu^2)
    q = 1.0e10*t_radius_ratio^3
    # analytical solution for the vertical deflection under the load
    # Includes the effect of shear deformation.
    analyt_sol = -q*a^4/64/D*((5+nu)/(1+nu)+4/(1-nu)*thickness^2/a^2);
    formul = FEMMShellT3FFModule

    tolerance = a/n/1000
    if mesh_procedure == :q4_t3
        fens, fes = Q4circlen(a, n);
        fens, fes = Q4toT3(fens, fes)
    elseif mesh_procedure == :t3_nice
        fens, fes = T3circlen(a, n);
    elseif mesh_procedure == :t3_poor
        fens, fes = T3circleseg(pi/2, a, n, n);
    else
        @error "Unknown mesh procedure"
    end
    fens.xyz = xyz3(fens)

    mater = MatDeforElastIso(DeforModelRed3D, E, nu)

    @info "Mesh: $n elements per side"

    sfes = FESetShellT3()
    accepttodelegate(fes, sfes)
    femm = formul.make(IntegDomain(fes, TriRule(1), thickness), mater)
    stiffness = formul.stiffness
    associategeometry! = formul.associategeometry!

    # Construct the requisite fields, geometry and displacement
    # Initialize configuration variables
    geom0 = NodalField(fens.xyz)
    u0 = NodalField(zeros(size(fens.xyz,1), 3))
    Rfield0 = initial_Rfield(fens)
    dchi = NodalField(zeros(size(fens.xyz,1), 6))

    # Apply EBC's
    bfes = meshboundary(fes)
    lx = selectelem(fens, bfes; box = Float64[0 0 -Inf Inf -Inf Inf], inflate = tolerance)
    ly = selectelem(fens, bfes; box = Float64[-Inf Inf 0 0 -Inf Inf], inflate = tolerance)
    lc = setdiff(1:count(bfes), vcat(lx, ly))
    # plane of symmetry perpendicular to X
    l1 = connectednodes(subset(bfes, lx))
    for i in [1,5,6]
        setebc!(dchi, l1, true, i)
    end
    # plane of symmetry perpendicular to Y
    l1 = connectednodes(subset(bfes, ly))
    for i in [2,4,6]
        setebc!(dchi, l1, true, i)
    end
    # ss support
    l1 = connectednodes(subset(bfes, lc))
    for i in [1,2,3,]
        setebc!(dchi, l1, true, i)
    end
    applyebc!(dchi)
    numberdofs!(dchi);

    # Assemble the system matrix
    associategeometry!(femm, geom0)
    K = stiffness(femm, geom0, u0, Rfield0, dchi);

    # Load
    nl = selectnode(fens; box = Float64[0 0 0 0 -Inf Inf], tolerance = tolerance)
    lfemm = FEMMBase(IntegDomain(fes, TriRule(3)))
    fi = ForceIntensity(FFlt[0, 0, -q, 0, 0, 0]);
    F = distribloads(lfemm, geom0, dchi, fi, 2);

    # @infiltrate
    # Solve
    U = K\F
    scattersysvec!(dchi, U[:])
    targetu =  dchi.values[nl, 3][1]
    @info "Target: $(round(targetu, digits=8)),  $(round(targetu/analyt_sol, digits = 4)*100)%"

    # Visualization
    if !visualize
        return true
    end
    scattersysvec!(dchi, (a/4)/maximum(abs.(U)).*U)
    update_rotation_field!(Rfield0, dchi)
    plots = cat(plot_space_box([[0 0 -a/2]; [a a a/2]]),
        #plot_nodes(fens),
        plot_midsurface(fens, fes; x = geom0.values, u = dchi.values[:, 1:3], R = Rfield0.values);
    dims = 1)
    pl = render(plots)
    return true
end

function test_convergence()
    t_radius_ratio = 0.05
    @info "Simply supported square plate with uniform load,"
    @info "thickness/length = $t_radius_ratio "
    for n in [2, 4, 8, 16, 32, 64]
        _execute(:t3_poor, n, t_radius_ratio, false)
    end
    return true
end

function allrun()
    println("#####################################################")
    println("# test_convergence ")
    test_convergence()
    return true
end # function allrun

@info "All examples may be executed with "
println("using .$(@__MODULE__); $(@__MODULE__).allrun()")

end # module
nothing