From e1f9912b95b5b4d539e27434f55d11c827516d0e Mon Sep 17 00:00:00 2001 From: Gerd Steinebach <64948537+gstein3m@users.noreply.github.com> Date: Sun, 4 Feb 2024 20:09:18 +0100 Subject: [PATCH] Update rosenbrock_perform_step.jl Bug in out-of-place version --- src/perform_step/rosenbrock_perform_step.jl | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/perform_step/rosenbrock_perform_step.jl b/src/perform_step/rosenbrock_perform_step.jl index b932c9c1fb..6114431e03 100644 --- a/src/perform_step/rosenbrock_perform_step.jl +++ b/src/perform_step/rosenbrock_perform_step.jl @@ -1034,7 +1034,7 @@ end integrator.k[3] = h2_21 * k1 + h2_22 * k2 + h2_23 * k3 + h2_24 * k4 + h2_25 * k5 if integrator.opts.adaptive calculate_interpoldiff!(k1, k2, uprev, du, u, integrator.k[1], integrator.k[2], integrator.k[3]) - atm = calculate_residuals!(k2, uprev, k1, integrator.opts.abstol, + atmp = calculate_residuals(k2, uprev, k1, integrator.opts.abstol, integrator.opts.reltol, integrator.opts.internalnorm, t) EEst = max(EEst,integrator.opts.internalnorm(atmp, t)) #-- role of t unclear end