From 2e56390c76ec24de94c097986d82d69371d977b1 Mon Sep 17 00:00:00 2001 From: Torkel Loman Date: Thu, 23 May 2024 15:16:17 -0400 Subject: [PATCH] Update src/steady_state_stability.jl Co-authored-by: Sam Isaacson --- src/steady_state_stability.jl | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/steady_state_stability.jl b/src/steady_state_stability.jl index e65e4ad921..8ac4e27a83 100644 --- a/src/steady_state_stability.jl +++ b/src/steady_state_stability.jl @@ -67,7 +67,7 @@ function steady_state_stability(u::Vector, rs::ReactionSystem, ps; tol = 10*sqrt # Computes stability (by checking that the real part of all eigenvalues is negative). max_eig = maximum(real(ev) for ev in eigvals(J)) if abs(max_eig) < tol - error("The system Jacobian's maximum eigenvalue at the steady state is within the tolerance range (abs($max_eig) < $tol). Hence, stability could not be reliably determined. If you still wish to compute the stability, reduce the `tol` argument.") + error("The system Jacobian's maximum eigenvalue at the steady state is within the tolerance range (abs($max_eig) < $tol). Hence, stability could not be reliably determined. If you still wish to compute the stability, reduce the `tol` argument, but note that floating point error in the eigenvalue calculation could lead to incorrect results.") end return max_eig < 0 end