Multi-scale Turing pattern simulation solver for the GeoStats.jl framework.
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
)