#Pkg.clone("https://github.com/goropikari/Percolation.jl")
using Percolation

N = 50; p = 0.6
sq = Square(N, p, "nn")
sq.lattice_config.lattice_sites

label_components!(sq)
sq.lattice_properties.labeled_lattice_sites

heatmap!(sq)

N = 100; p = 0.45
sq = Square(N, p, "nnn")
label_components!(sq)

heatmap!(sq)

N = 100; p = 0.45
tri = Triangular(N, p, "nnn")
label_components!(tri)
heatmap!(tri)

N = 100; p = 0.45
h = Honeycomb(N, p, "nnn")
label_components!(h)
heatmap!(h)

linsize = 50
ps = 0.1      
pinc = 0.025
pf = 1
nsample = 50
plot_percolation_prob("square", linsize, ps, pinc, pf, nsample)

linsize = 50
ps = 0
pinc = 0.025
pf = 1
nsample = 50
plot_percolation_prob("triangular", linsize, ps, pinc, pf, nsample)

linsize = 50
ps = 0
pinc = 0.025
pf = 1
nsample = 50
plot_percolation_prob("honeycomb", linsize, ps, pinc, pf, nsample)

N = 50; p = 0.6
sq = Square(N, p, "nn")
println(ispercolation!(sq))
heatmap!(sq)

clustersize!(sq)

clustersizefreq!(sq)

?clusternumber!

clusternumber!(sq)

average_clustersize!(sq)

strength!(sq)

linsize = 50
ps = 0
pinc = 0.025
pf = 1
nsample = 50
plot_lifetime(linsize, ps, pinc, pf, nsample)