TuringPatterns.jl

using TuringPatterns

sz = (200, 200)

patterns = [
    SimplePattern(Params(2,   4,   0.01), sz),
    SimplePattern(Params(5,   10,  0.02), sz),
    SimplePattern(Params(10,  20,  0.03), sz),
    SimplePattern(Params(20,  40,  0.04), sz),
    SimplePattern(Params(50,  100, 0.05), sz),
]

simulate(
    rand(sz...), # initial conditions
    patterns,    # patterns
    BoxBlur(sz), # blur
    100          # iterations
)

using GeoStats
using TuringPatterns

using Plots, GeoStatsPlots

# define simulation problem for a variable "z"
# request 3 realizations (i.e. images)
problem = SimulationProblem(CartesianGrid(200,200), :z => Float64, 3)

# define Turing patterns solver
# see docstring for options
solver  = TPS()

# solve problem over any domain
# that has grid topology
solution = solve(problem, solver)

plot(solution)