Make examples runnable

Author: ArnoStrouwen

docs/common_first_steps.jl

@@ -1,6 +1,25 @@
 using Markdown
-function first_steps(name, solver)
-    Markdown.parse("""## Installation
+function first_steps(name, solver; fixed_timesteps = false)
+    sym_name =  Symbol(name)
+    sym_solver = Symbol(solver)
+    @eval begin
+    using $sym_name
+    function lorenz!(du, u, p, t)
+        du[1] = 10.0 * (u[2] - u[1])
+        du[2] = u[1] * (28.0 - u[3]) - u[2]
+        du[3] = u[1] * u[2] - (8 / 3) * u[3]
+    end
+    u0 = [1.0; 0.0; 0.0]
+    tspan = (0.0, 100.0)
+    prob = ODEProblem(lorenz!, u0, tspan)
+    if !($fixed_timesteps)
+        sol = solve(prob, ($sym_solver)())
+    else
+        sol = solve(prob, ($sym_solver)(), dt = 1/100)
+    end
+    sol_end = sol[end]
+    end
+    installation = """## Installation
     To be able to access the solvers in `$name`, you must first install them use the Julia package manager:
 
     ```julia
@@ -10,20 +29,53 @@ function first_steps(name, solver)
     This will only install the solvers listed at the bottom of this page.
     If you want to explore other solvers for your problem,
     you will need to install some of the other libraries listed in the navigation bar on the left.
+    """
+    example = if !fixed_timesteps
+        """
 
-    ## Example usage
+        ## Example usage
 
-    ```julia
-    using $name
+        ```julia
+        using $name
 
-    function lorenz!(du, u, p, t)
-        du[1] = 10.0 * (u[2] - u[1])
-        u[2] = u[1] * (28.0 - u[3]) - u[2]
-        du[3] = u[1] * u[2] - (8 / 3) * u[3]
+        function lorenz!(du, u, p, t)
+            du[1] = 10.0 * (u[2] - u[1])
+            du[2] = u[1] * (28.0 - u[3]) - u[2]
+            du[3] = u[1] * u[2] - (8 / 3) * u[3]
+        end
+        u0 = [1.0; 0.0; 0.0]
+        tspan = (0.0, 100.0)
+        prob = ODEProblem(lorenz!, u0, tspan)
+        sol = solve(prob, $solver())
+        sol[end]
+        ```
+        ```julia
+        $sol_end
+        ```
+        """
+    else
+        """
+
+        ## Example usage
+
+        ```julia
+        using $name
+
+        function lorenz!(du, u, p, t)
+            du[1] = 10.0 * (u[2] - u[1])
+            du[2] = u[1] * (28.0 - u[3]) - u[2]
+            du[3] = u[1] * u[2] - (8 / 3) * u[3]
+        end
+        u0 = [1.0; 0.0; 0.0]
+        tspan = (0.0, 100.0)
+        prob = ODEProblem(lorenz!, u0, tspan)
+        sol = solve(prob, $solver(), dt=1/100)
+        sol[end]
+        ```
+        ```julia
+        $sol_end
+        ```
+        """
     end
-    u0 = [1.0; 0.0; 0.0]
-    tspan = (0.0, 100.0)
-    prob = ODEProblem(lorenz!, u0, tspan)
-    sol = solve(prob, $solver())
-    ```""")
+    Markdown.parse(installation*example)
 end

docs/src/explicit/LowStorageRK.md

@@ -12,7 +12,7 @@ are generally only recommended in situations which are RAM bound, like large-sca
 
 ```@eval
 first_steps = evalfile("./common_first_steps.jl")
-first_steps("OrdinaryDiffEqLowStorageRK", "CarpenterKennedy2N54")
+first_steps("OrdinaryDiffEqLowStorageRK", "CKLLSRK43_2")
 ```
 
 ## Full list of solvers

docs/src/explicit/PRK.md

@@ -8,7 +8,7 @@ Explicit solvers optimized for a certain number of parallel calls of the system
 
 ```@eval
 first_steps = evalfile("./common_first_steps.jl")
-first_steps("OrdinaryDiffEqPRK", "KuttaPRK2p5")
+first_steps("OrdinaryDiffEqPRK", "KuttaPRK2p5", fixed_timesteps = true)
 ```
 
 ## Full list of solvers

docs/src/explicit/SSPRK.md

@@ -14,7 +14,7 @@ Note that for SSPRK methods, a algorithm utility `OrdinaryDiffEqCore.ssp_coeffic
 
 ```@eval
 first_steps = evalfile("./common_first_steps.jl")
-first_steps("OrdinaryDiffEqSSPRK", "SSPRK22")
+first_steps("OrdinaryDiffEqSSPRK", "SSPRK432")
 ```
 
 ## Full list of solvers

docs/src/implicit/PDIRK.md

@@ -24,7 +24,7 @@ which is why these are the parallel DIRK methods.
 
 ```@eval
 first_steps = evalfile("./common_first_steps.jl")
-first_steps("OrdinaryDiffEqPDIRK", "PDIRK44")
+first_steps("OrdinaryDiffEqPDIRK", "PDIRK44", fixed_timesteps = true)
 ```
 
 ## Full list of solvers

docs/src/implicit/SDIRK.md

@@ -8,7 +8,7 @@ This article is a stub.
 
 ```@eval
 first_steps = evalfile("./common_first_steps.jl")
-first_steps("OrdinaryDiffEqSDIRK", "PDIRK44")
+first_steps("OrdinaryDiffEqSDIRK", "KenCarp3")
 ```
 
 ## Full list of solvers

docs/src/semiimplicit/ExponentialRK.md

@@ -23,7 +23,7 @@ you will need to install some of the other libraries listed in the navigation ba
 
 ```@eval
 first_steps = evalfile("./common_first_steps.jl")
-first_steps("OrdinaryDiffEqExponentialRK", "EPIRK5s3")
+first_steps("OrdinaryDiffEqExponentialRK", "Exprb43", fixed_timesteps = true)
 ```
 
 ## Full list of solvers