src/fhn_solver.jl

@@ -16,7 +16,7 @@ using .Laplacian
     # Arguments:
     `U`: activator matrix
     `V`: inhibitor matrix
-    `param_obs`: used parameters from CombinedPDEParams
+    `param_obs`: used parameters from FHNParams
     `dx`: dx
     `dt`: dt
 """

src/gray_scott_solver.jl

@@ -14,7 +14,7 @@ using .Laplacian
     # Arguments:
     `U`: activator matrix
     `V`: inhibitor matrix
-    `param_obs`: used parameters from CombinedPDEParams
+    `param_obs`: used parameters from FHNParams
     `dx`: dx
 """
 function step_gray_scott!(U, V, params_obs::Observable; dx=1)

src/utils/constants.jl

@@ -2,6 +2,25 @@
 
 abstract type PDEParams end
 
+struct FHNParams <: PDEParams
+    N::Int
+    dx::Float64 # grid spacing
+    Du::Float64
+    Dv::Float64
+    ϵ::Float64
+    a::Float64
+    b::Float64
+end
+
+struct GSParams <: PDEParams
+    N::Int       # grid size
+    dx::Float64  # grid spacing
+    Du::Float64  # diffusion rate U
+    Dv::Float64  # diffusion rate V
+    F::Float64   # feed rate
+    k::Float64   # kill rate
+
+end
 
 struct CombinedPDEParams <: PDEParams
     N::Int