NATURE_MI2020/make_graphs.jl at master · jpp46/NATURE_MI2020 · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
using Random, Statistics, DataFrames, CSV
using BSON: @load

using Plots
pyplot(size=(500,400))

const len = 0.22
const time = 19.796841683855746

function normalize(x, minx, maxx)
    return (x .- minx)./(maxx.-minx)
end

function get_data(folder)
	data = Matrix{Float64}(undef, 60, 201)
	for i in 1:60
		File = "$folder/worker_$(lpad(i, 2, "0"))_results.bson"
		@load File results
		for j in 0:200
			data[i, j+1] = results[j+1][3]
		end
	end
	return data
end

function write_data(path, data)
	df = DataFrame([:g => 0:200]..., [Symbol("f_$i") => data[i, :] for i in 1:60]...)
	CSV.write(path, df)
end

data_0 = get_data("genomes_0")
println("0 Open => Max: $((maximum(data_0)/len)/time), Mean: $((mean(data_0, dims=1)[1, end]/len)/time)")
data_semi_0 = get_data("genomes_semi_0")
println("0 Semi => Max: $((maximum(data_semi_0)/len)/time), Mean: $((mean(data_semi_0, dims=1)[1, end]/len)/time)")
data_closed_0 = get_data("genomes_closed_0")
println("0 Closed => Max: $((maximum(data_closed_0)/len)/time), Mean: $((mean(data_closed_0, dims=1)[1, end]/len)/time)")

data_0 = ((data_0./len)./time)
data_semi_0 = ((data_semi_0./len)./time)
data_closed_0 = ((data_closed_0./len)./time)
write_data("extra/flat_open.csv", data_0)
write_data("extra/flat_semi.csv", data_semi_0)
write_data("extra/flat_closed.csv", data_closed_0)

#maxx = maximum([maximum(data_0), maximum(data_semi_0), maximum(data_closed_0)])
#minx = minimum([minimum(data_0), minimum(data_semi_0), minimum(data_closed_0)])
#data_0 = normalize(data_0, minx, maxx)
#data_semi_0 = normalize(data_semi_0, minx, maxx)
#data_closed_0 = normalize(data_closed_0, minx, maxx)
#write_data("extra/flat_open_n.csv", data_0)
#write_data("extra/flat_semi_n.csv", data_semi_0)
#write_data("extra/flat_closed_n.csv", data_closed_0)

mx = maximum(data_0, dims=1)
m = mean(data_0, dims=1)
v = std(data_0, dims=1)
plt = plot(0:200, m[:], ribbon=v[:],
	title="Flat Surface", xlabel="Generations", ylabel="BL/s",
	label="Orientation,Shape,Control", legend=:none, grid=false,
	linecolor=:blue, fill=:blue, fillalpha=0.15, ylims=(0.0,0.201)
)
plot!(plt, mx[:], linecolor=:blue, linestyle=:dash, label=nothing)

mx = maximum(data_semi_0, dims=1)
m = mean(data_semi_0, dims=1)
v = std(data_semi_0, dims=1)
plot!(plt, 0:200, m[:], ribbon=v[:],
	title="Flat Surface", xlabel="Generations", ylabel="BL/s",
	label="Shape,Control", legend=:none, grid=false,
	linecolor=:green, fill=:green, fillalpha=0.15, ylims=(0.0,0.201)
)
plot!(plt, mx[:], linecolor=:green, linestyle=:dash, label=nothing)

mx = maximum(data_closed_0, dims=1)
m = mean(data_closed_0, dims=1)
v = std(data_closed_0, dims=1)
plot!(plt, 0:200, m[:], ribbon=v[:],
	title="Flat Surface", xlabel="Generations", ylabel="BL/s",
	label="Control", legend=:none, grid=false,
	linecolor=:red, fill=:red, fillalpha=0.15, ylims=(0.0,0.201)
)
plot!(plt, mx[:], linecolor=:red, linestyle=:dash, label=nothing)
savefig(plt, "graphs/flat.png")


data_15 = get_data("genomes_15")
println("15 Open => Max: $((maximum(data_15)/len)/time), Mean: $((mean(data_15, dims=1)[1, end]/len)/time)")
data_semi_15 = get_data("genomes_semi_15")
println("15 Semi => Max: $((maximum(data_semi_15)/len)/time), Mean: $((mean(data_semi_15, dims=1)[1, end]/len)/time)")
data_closed_15 = get_data("genomes_closed_15")
println("15 Closed => Max: $((maximum(data_closed_15)/len)/time), Mean: $((mean(data_closed_15, dims=1)[1, end]/len)/time)")

data_15 = ((data_15./len)./time)
data_semi_15 = ((data_semi_15./len)./time)
data_closed_15 = ((data_closed_15./len)./time)
write_data("extra/hill_open.csv", data_15)
write_data("extra/hill_semi.csv", data_semi_15)
write_data("extra/hill_closed.csv", data_closed_15)

#maxx = maximum([maximum(data_15), maximum(data_semi_15), maximum(data_closed_15)])
#minx = minimum([minimum(data_15), minimum(data_semi_15), minimum(data_closed_15)])
#data_15 = normalize(data_15, minx, maxx)
#data_semi_15 = normalize(data_semi_15, minx, maxx)
#data_closed_15 = normalize(data_closed_15, minx, maxx)
#write_data("extra/hill_open_n.csv", data_15)
#write_data("extra/hill_semi_n.csv", data_semi_15)
#write_data("extra/hill_closed_n.csv", data_closed_15)

mx = maximum(data_15, dims=1)
m = mean(data_15, dims=1)
v = std(data_15, dims=1)
plt = plot(0:200, m[:], ribbon=v[:],
	title="Inclined Surface", xlabel="Generations", ylabel="BL/s",
	label="Orientation,Shape,Control", legend=:bottomright, grid=false,
	linecolor=:blue, fill=:blue, fillalpha=0.15, ylims=(-0.075,0.016)
)
plot!(plt, mx[:], linecolor=:blue, linestyle=:dash, label=nothing)

mx = maximum(data_semi_15, dims=1)
m = mean(data_semi_15, dims=1)
v = std(data_semi_15, dims=1)
plot!(plt, 0:200, m[:], ribbon=v[:],
	title="Inclined Surface", xlabel="Generations", ylabel="BL/s",
	label="Shape,Control", legend=:bottomright, grid=false,
	linecolor=:green, fill=:green, fillalpha=0.15, ylims=(-0.075,0.016)
)
plot!(plt, mx[:], linecolor=:green, linestyle=:dash, label=nothing)

mx = maximum(data_closed_15, dims=1)
m = mean(data_closed_15, dims=1)
v = std(data_closed_15, dims=1)
plot!(plt, 0:200, m[:], ribbon=v[:],
	title="Inclined Surface", xlabel="Generations", ylabel="BL/s",
	label="Control", legend=:bottomright, grid=false,
	linecolor=:red, fill=:red, fillalpha=0.15, ylims=(-0.075,0.016)
)
plot!(plt, mx[:], linecolor=:red, linestyle=:dash, label=nothing)
savefig(plt, "graphs/hill.png")