diff --git a/docs/make.jl b/docs/make.jl
index e1156af..cf8a7ce 100644
--- a/docs/make.jl
+++ b/docs/make.jl
@@ -50,6 +50,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",
 ]
 
 for (_, name) in examples
diff --git a/docs/src/refs.bib b/docs/src/refs.bib
index 5abda41..563900f 100644
--- a/docs/src/refs.bib
+++ b/docs/src/refs.bib
@@ -1,7 +1,8 @@
-@BOOK{Kelley2022
+@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},
   location  = {Philadelphia, PA},
   date      = {2022-01},
@@ -14,6 +15,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},
   location  = {Philadelphia, Pa.},
   date      = {2003},
@@ -40,7 +42,8 @@ @ARTICLE{Kan2022-ko
                   Jacobian-free {Newton–Krylov​} method},
   author       = {Kan, Ziyun and Song, Ningning and Peng, Haijun and Chen,
                   Biaosong},
-  journaltitle = {Journal of computational and applied mathematics},
+  year         = {2022},
+  journal      = {Journal of Computational and Applied Mathematics},
   publisher    = {Elsevier BV},
   volume       = {399},
   issue        = {113732},
@@ -55,7 +58,8 @@ @ARTICLE{Kan2022-ko
 @ARTICLE{Hosea1996-xv,
   title        = {Analysis and implementation of {TR}-{BDF2}},
   author       = {Hosea, M E and Shampine, L F},
-  journaltitle = {Applied numerical mathematics: transactions of IMACS},
+  year         = {1996},
+  journal      = {Applied Numerical Mathematics: Transactions of IMACS},
   publisher    = {Elsevier BV},
   volume       = {20},
   issue        = {1-2},
@@ -71,7 +75,8 @@ @ARTICLE{Bank1985-gh
   title        = {Transient simulation of silicon devices and circuits},
   author       = {Bank, R E and Coughran, W M and Fichtner, W and Grosse, E H
                   and Rose, D J and Smith, R K},
-  journaltitle = {IEEE transactions on electron devices},
+  year         = {1985},
+  journal      = {IEEE Transactions on Electron Devices},
   publisher    = {Institute of Electrical and Electronics Engineers (IEEE)},
   volume       = {32},
   issue        = {10},
@@ -87,7 +92,8 @@ @ARTICLE{Bonaventura2021-za
   title        = {The {TR}-{BDF2} method for second order problems in structural
                   mechanics},
   author       = {Bonaventura, Luca and Mármol, Macarena Gómez},
-  journaltitle = {Computers \& mathematics with applications (Oxford, England:
+  year         = {2021},
+  journal      = {Computers \& mathematics with applications (Oxford, England:
                   1987)},
   publisher    = {Elsevier BV},
   volume       = {92},
diff --git a/examples/trixi_rosenbrock.jl b/examples/trixi_rosenbrock.jl
index 39bd63c..0c27a4c 100644
--- a/examples/trixi_rosenbrock.jl
+++ b/examples/trixi_rosenbrock.jl
@@ -26,6 +26,5 @@ sol = solve(
     ode, Theseus.SSPKnoth();
     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
 );
diff --git a/examples/trixi_rosenbrock_blast_wave.jl b/examples/trixi_rosenbrock_blast_wave.jl
index 5a7a164..2e50c7d 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 7490c67..7e26b54 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 c126847..397f3d3 100644
--- a/src/Ariadne.jl
+++ b/src/Ariadne.jl
@@ -373,6 +373,8 @@ function newton_krylov!(
 
     if forcing !== nothing
         η = initial(forcing)
+    else
+        η = missing
     end
 
     verbose > 0 && @info "Jacobian-Free Newton-Krylov" algo res₀ = n_res tol tol_rel tol_abs η