diff --git a/docs/make.jl b/docs/make.jl
index c528133c..0304c1a8 100644
--- a/docs/make.jl
+++ b/docs/make.jl
@@ -51,6 +51,9 @@ examples = [
     "Implicit -- Heat 1D DG" => "heat_1D_DG",
     "Implicit -- Heat 2D" => "heat_2D",
     "Trixi" => "trixi",
+    "Trixi Rosenbrock" => "trixi_rosenbrock",
+    "Trixi Rosenbrock Warm Bubble" => "trixi_rosenbrock_warm_bubble",
+    "Trixi Rosenbrock Blast Wave" => "trixi_rosenbrock_blast_wave",
     "Trixi IMEX SSP" => "trixi_imex_ssp",
     "Trixi IMEX ARS" => "trixi_imex_ars",
     "Trixi IMEX Von Karman street" => "trixi_imex_von_karman_street",
diff --git a/docs/src/refs.bib b/docs/src/refs.bib
index 49abf10e..f7a58db3 100644
--- a/docs/src/refs.bib
+++ b/docs/src/refs.bib
@@ -2,6 +2,7 @@ @BOOK{Kelley2022
   title     = {Solving nonlinear equations with iterative methods: Solvers and
                examples in Julia},
   author    = {Kelley, C T},
+  year      = {2022},
   publisher = {Society for Industrial and Applied Mathematics},
   year      = {2022},
   doi       = {10.1137/1.9781611977271},
@@ -11,6 +12,7 @@ @BOOK{Kelley2022
 @BOOK{Kelley2003,
   title     = {Solving nonlinear equations with Newton's method},
   author    = {Kelley, C T},
+  year      = {2003},
   publisher = {Society for Industrial and Applied Mathematics},
   year      = {2003},
   doi       = {10.1137/1.9780898718898},
diff --git a/examples/trixi_rosenbrock_blast_wave.jl b/examples/trixi_rosenbrock_blast_wave.jl
index 5a7a1646..2e50c7dc 100644
--- a/examples/trixi_rosenbrock_blast_wave.jl
+++ b/examples/trixi_rosenbrock_blast_wave.jl
@@ -26,7 +26,6 @@ sol = solve(
     ode, Theseus.ROS2();
     dt = 1.0, # solve needs some value here but it will be overwritten by the stepsize_callback
     ode_default_options()..., callback = callbacks,
-    # verbose=1,
     krylov_algo = :gmres,
     assume_p_const = false,
 );
diff --git a/examples/trixi_rosenbrock_warm_bubble.jl b/examples/trixi_rosenbrock_warm_bubble.jl
index 7490c679..7e26b548 100644
--- a/examples/trixi_rosenbrock_warm_bubble.jl
+++ b/examples/trixi_rosenbrock_warm_bubble.jl
@@ -26,7 +26,6 @@ sol = solve(
     ode, Theseus.ROS2();
     dt = 1.0, # solve needs some value here but it will be overwritten by the stepsize_callback
     ode_default_options()..., callback = callbacks,
-    # verbose=1,
     krylov_algo = :gmres,
     assume_p_const = false,
 );
diff --git a/src/Ariadne.jl b/src/Ariadne.jl
index 14628d22..18abf57a 100644
--- a/src/Ariadne.jl
+++ b/src/Ariadne.jl
@@ -495,6 +495,8 @@ function newton_krylov!(
 
     if forcing !== nothing
         η = initial(forcing)
+    else
+        η = missing
     end
 
     verbose > 0 && @info "Jacobian-Free Newton-Krylov" res₀ = norm_res tol tol_rel tol_abs η