59 lines
1.2 KiB
Julia
59 lines
1.2 KiB
Julia
include("../src/utils/constants.jl")
|
|
include("../src/fhn_solver.jl")
|
|
#include("../src/gray_scott_solver.jl")
|
|
include("../src/visualization.jl")
|
|
|
|
using Observables
|
|
using GLMakie
|
|
|
|
using .Constants
|
|
using .Visualization
|
|
|
|
"""
|
|
# GSParams
|
|
N = 128
|
|
dx = 1.0
|
|
Du, Dv = Observable(0.16), Observable(0.08)
|
|
F, k = Observable(0.060), Observable(0.062)
|
|
|
|
param_observables = (
|
|
Du=Du,
|
|
Dv=Dv,
|
|
F=F,
|
|
k=k,
|
|
)
|
|
"""
|
|
|
|
# FHNParams
|
|
N = 128
|
|
dx = 1.0
|
|
Du, Dv = Observable(0.016), Observable(0.1)
|
|
ϵ, a, b = Observable(0.1), Observable(0.5), Observable(0.9)
|
|
param_observables = (
|
|
Du=Du,
|
|
Dv=Dv,
|
|
ϵ=ϵ,
|
|
a=a,
|
|
b=b
|
|
)
|
|
|
|
#params_obs = Observable(Constants.GSParams(N, dx, Du[], Dv[], F[], k[]))
|
|
params_obs = Observable(Constants.FHNParams(N=N, dx=dx, Du=Du[], Dv=Dv[], ϵ=ϵ[], a=a[], b=b[]))
|
|
|
|
"""
|
|
lift(Du, Dv, F, k) do u, v, f, ki
|
|
params_obs[] = Constants.GSParams(N, dx, u, v, f, ki)
|
|
end
|
|
"""
|
|
|
|
lift(Du, Dv, ϵ, a, b) do d_u, d_v, eps, aa, bb
|
|
params_obs[] = Constants.FHNParams(N=N, dx=dx, Du=d_u, Dv=d_v, ϵ=eps, a=aa, b=bb)
|
|
end
|
|
|
|
U = ones(N, N)
|
|
V = zeros(N, N)
|
|
heat_obs = Observable(U)
|
|
|
|
fig = Visualization.build_ui(U, V, param_observables, params_obs, heat_obs)
|
|
display(fig)
|