Merge pull request #114 from alexfikl/fix-brusselator

Fix Brusselator right-hand side in some examples
SciFracX · Aug 22, 2024 · eccc025 · eccc025 · github-actions · Aug 22, 2024
2 parents 8576eb0 + 491b643
commit eccc025
Show file tree

Hide file tree

Showing 3 changed files with 5 additions and 5 deletions.
diff --git a/bench/runbenchmarks.jl b/bench/runbenchmarks.jl
@@ -10,7 +10,7 @@ u0 = [0.2, 0.03]
 function Brusselator!(du, u, p, t)
     a, μ = 1, 4
     du[1] = a - (μ + 1) * u[1] + (u[1])^2 * u[2]
-    return du[2] = μ * u[1] - (u[1])^2 * u[1]
+    return du[2] = μ * u[1] - (u[1])^2 * u[2]
 end
 
 prob = FODEProblem(Brusselator!, α, u0, (0, 20))

diff --git a/docs/ChaosGallery/Gallery/Brusselator.md b/docs/ChaosGallery/Gallery/Brusselator.md
@@ -3,7 +3,7 @@
 ```math
 \begin{cases}
 {_{t_0}D_t^\alpha}y_1(t)=a-(\mu+1)y_1(t)+(y_1(t))^2y_2(t)\\
-{_{t_0}D_t^\alpha}y_2(t)=\mu y_1(t)-(y_1(t))^2y_1(t)\\
+{_{t_0}D_t^\alpha}y_2(t)=\mu y_1(t)-(y_1(t))^2y_2(t)\\
 y_1(t_0)=y_{1,0}\\
 y_2(t_0)=y_{2,0}
 \end{cases}
@@ -20,7 +20,7 @@ tspan = (0, 100)
 function Brusselator!(du, u, p, t)
     a, μ = 1, 4
     du[1] = a-(μ+1)*u[1]+(u[1])^2*u[2]
-    du[2] = μ*u[1]-(u[1])^2*u[1]
+    du[2] = μ*u[1]-(u[1])^2*u[2]
 end
 
 prob = FODEProblem(Brusselator!, α, u0, tspan)
@@ -35,4 +35,4 @@ plot(sol, vars=(0,1,2))
 
 ![BrusselatorPhase](./assets/Brusselator.png)
 
-![BrusselatorTime](./assets/BrusselatorTime.png)
+![BrusselatorTime](./assets/BrusselatorTime.png)
diff --git a/examples/Brusselator.jl b/examples/Brusselator.jl
@@ -7,7 +7,7 @@ h = 0.001
 function Brusselator!(du, u, p, t)
     a, μ = 1, 4
     du[1] = a - (μ + 1) * u[1] + (u[1])^2 * u[2]
-    du[2] = μ * u[1] - (u[1])^2 * u[1]
+    du[2] = μ * u[1] - (u[1])^2 * u[2]
 end
 
 prob = FODESystem(Brusselator!, α, u0, (0, 100))
Benchmark suite	Current: `eccc025`	Previous: `8576eb0`	Ratio
`FLMM/Trapezoid`	`19841908` ns	`20342051` ns	`0.98`
`FLMM/NewtonGregory`	`18515591.5` ns	`16888141` ns	`1.10`
`FLMM/BDF`	`18477555` ns	`16795517` ns	`1.10`