use Δ to name the laplacian

bifurcationkit · Jun 25, 2023 · 848ac29 · 848ac29
1 parent a31f2dc
commit 848ac29
Show file tree

Hide file tree

Showing 4 changed files with 7 additions and 7 deletions.
diff --git a/docs/src/tutorials/dde/HutchinsonDiff.md b/docs/src/tutorials/dde/HutchinsonDiff.md
@@ -46,15 +46,15 @@ X = -Lx .+ 2Lx/Nx*(0:Nx-1) |> collect
 
 # boundary condition
 Q = Neumann0BC(X[2]-X[1])
-Dxx = sparse(CenteredDifference(2, 2, X[2]-X[1], Nx) * Q)[1]
+Δ = sparse(CenteredDifference(2, 2, X[2]-X[1], Nx) * Q)[1]
 nothing #hide
 ```
 
 We are now ready to compute the bifurcation of the trivial (constant in space) solution:
 
 ```@example TUTHut
 # bifurcation problem
-pars = (a = 0.5, d = 1, τ = 1.0, Δ = Dxx, N = Nx)
+pars = (a = 0.5, d = 1, τ = 1.0, Δ = Δ, N = Nx)
 x0 = zeros(Nx)
 
 prob = ConstantDDEBifProblem(Hutchinson, delaysF, x0, pars, (@lens _.a))

diff --git a/examples/HutchinsonDiffusion.jl b/examples/HutchinsonDiffusion.jl
@@ -26,9 +26,9 @@ X = -Lx .+ 2Lx/Nx*(0:Nx-1) |> collect
 
 # boundary condition
 Q = Neumann0BC(X[2]-X[1])
-Dxx = sparse(CenteredDifference(2, 2, X[2]-X[1], Nx) * Q)[1]
+Δ = sparse(CenteredDifference(2, 2, X[2]-X[1], Nx) * Q)[1]
 
-pars = (a = 0.5, d = 1, τ = 1.0, Δ = Dxx, N = Nx)
+pars = (a = 0.5, d = 1, τ = 1.0, Δ = Δ, N = Nx)
 x0 = zeros(Nx)
 
 prob = ConstantDDEBifProblem(Hutchinson, delaysF, x0, pars, (@lens _.a))

diff --git a/examples/TravisWebb.jl b/examples/TravisWebb.jl
@@ -26,9 +26,9 @@ X = -Lx .+ 2Lx/Nx*(0:Nx-1) |> collect
 
 # boundary condition
 Q = Dirichlet0BC(X[2]-X[1]|>typeof)
-Dxx = sparse(CenteredDifference(2, 2, X[2]-X[1], Nx) * Q)[1]
+Δ = sparse(CenteredDifference(2, 2, X[2]-X[1], Nx) * Q)[1]
 
-pars = (a = 0.5, τ = 1.0, Δ = Dxx, N = Nx)
+pars = (a = 0.5, τ = 1.0, Δ = Δ, N = Nx)
 x0 = zeros(Nx)
 
 prob = ConstantDDEBifProblem(TW, delaysF, x0, pars, (@lens _.a))

diff --git a/examples/humpries.jl b/examples/humpries.jl
@@ -92,7 +92,7 @@ probpo = PeriodicOrbitOCollProblem(200, 2; N = 1)
 	end
 	)
 
-plot(br, br_pocoll)
+plot(br);plot!(br_pocoll, plotfold=false, ylabel = "amplitude")
 ################################################################################
 using  DifferentialEquations