A better way to chose a random neighbor

[stevengj]

@lsaravia

x=[i,j] .+ rand(neighborhood)

This allocates an array [i,j] on every iteration. Better to use a tuple (i,j) .+ rand(neighborhood).

@lsaravia

# check boundaries
if !(x[1] > fil || x[1] <1 || x[2] > col || x[2]<1)

I would usually recommend employing ghost cells to avoid these kinds of boundary checks in inner loops.