From 142cd258622ad87318a6a5e743b342e756c7739f Mon Sep 17 00:00:00 2001
From: Gerd Steinebach <64948537+gstein3m@users.noreply.github.com>
Date: Thu, 8 Aug 2024 13:11:39 +0200
Subject: [PATCH] Update rosenbrock_perform_step.jl

Changes according to issue https://github.com/SciML/OrdinaryDiffEq.jl/issues/2315
---
 lib/OrdinaryDiffEqRosenbrock/src/rosenbrock_perform_step.jl | 6 +++---
 1 file changed, 3 insertions(+), 3 deletions(-)

diff --git a/lib/OrdinaryDiffEqRosenbrock/src/rosenbrock_perform_step.jl b/lib/OrdinaryDiffEqRosenbrock/src/rosenbrock_perform_step.jl
index 9393c2c573..c476601573 100644
--- a/lib/OrdinaryDiffEqRosenbrock/src/rosenbrock_perform_step.jl
+++ b/lib/OrdinaryDiffEqRosenbrock/src/rosenbrock_perform_step.jl
@@ -1716,7 +1716,7 @@ end
             else
                 du2 = mass_matrix * du1 - du
             end
-            EEst = norm(du2) / (integrator.opts.abstol + integrator.opts.reltol * norm(k2))
+            EEst = norm(du2) / norm(integrator.opts.abstol .+ integrator.opts.reltol .* k2)
             integrator.EEst = max(EEst, integrator.EEst)
         end
     end
@@ -1978,7 +1978,7 @@ end
                 mul!(_vec(du2), mass_matrix, _vec(du1))
                 du2 = du2 - du
             end
-            EEst = norm(du2) / (integrator.opts.abstol + integrator.opts.reltol * norm(k2))
+            EEst = norm(du2) / norm(integrator.opts.abstol .+ integrator.opts.reltol .* k2)
             integrator.EEst = max(EEst, integrator.EEst)
         end
     end
@@ -1986,4 +1986,4 @@ end
 end
 
 @RosenbrockW6S4OS(:init)
-@RosenbrockW6S4OS(:performstep)
\ No newline at end of file
+@RosenbrockW6S4OS(:performstep)