diff --git a/src/multiTracer.jl b/src/multiTracer.jl
index 8466d3a3..44840a5b 100644
--- a/src/multiTracer.jl
+++ b/src/multiTracer.jl
@@ -47,7 +47,7 @@ function split_state_function_and_Jacobian(Ts::Tuple, Ls::Tuple, NLs::Tuple, nb)
     return F, L, NL, ∇ₓF, ∇ₓL, ∇ₓNL, T
 end
 function split_state_function_and_Jacobian(T, L, NL, nb)
-    F(x,p) = NL(x,p) + L(p) * x - T(p) * x                     # full 𝐹(𝑥)
+    F(x,p) = NL(x,p) + L(x,p) - T(p) * x                     # full 𝐹(𝑥)
     ∇ₓNL(x,p) = sparse(Diagonal(𝔇(NL(x .+ ε,p))))     # Jacobian of nonlinear part
     ∇ₓL(p) = sparse(Diagonal(𝔇(L(zeros(nb) .+ ε,p))))     # Jacobian of nonlinear part
     ∇ₓF(x,p) = ∇ₓNL(x,p) + ∇ₓL(p) - T(p)       # full Jacobian ∇ₓ𝐹(𝑥) = -T + ∇ₓ𝐺(𝑥)
diff --git a/src/time_steppers.jl b/src/time_steppers.jl
index c208deb3..f69fe7cd 100644
--- a/src/time_steppers.jl
+++ b/src/time_steppers.jl
@@ -100,8 +100,8 @@ end
 Returns the Crank-Nicolson-leapfrog-step increment (out of place) using stored factors.
 See [`crank_nicolson_leapfrog_step_A⁺_and_A⁻`](@ref) for more information.
 """
-function crank_nicolson_leapfrog_step(xᵢ, xᵢ₋₁, p, A⁺::Factorization, A⁻, G)
-    return A⁺ \ (G(xᵢ,p) + A⁻ * xᵢ₋₁)
+function crank_nicolson_leapfrog_step(xᵢ, xᵢ₋₁, p, A⁺::Factorization, A⁻, NL)
+    return A⁺ \ (NL(xᵢ,p) + A⁻ * xᵢ₋₁)
 end
 """
     crank_nicolson_leapfrog_step_A⁺_and_A⁻(p, δt, T, ∇ₓL)