Fix threshold detection for reaction substepper (#131)

termi-official · Aug 8, 2024 · 56945d6 · 56945d6
1 parent b087c9c
commit 56945d6
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/solver/time/partitioned_solver.jl b/src/solver/time/partitioned_solver.jl
@@ -111,7 +111,7 @@ Adapt.@adapt_structure AdaptiveForwardEulerSubstepperCache
     # TODO get Cₘ
     cell_rhs!(du_local, u_local, x, t, cell_model)
 
-    if du_local[φₘidx] < cache.reaction_threshold
+    if abs(du_local[φₘidx]) < cache.reaction_threshold
         for j in 1:length(u_local)
             u_local[j] += Δt*du_local[j]
         end