Skip to content

Commit dc4f0f4

Browse files
Merge pull request #325 from avik-pal/ap/doc
Clarify unstable equilibrium for SteadyStateProblems
2 parents 1d2ef5c + 0e03d5a commit dc4f0f4

File tree

1 file changed

+14
-7
lines changed

1 file changed

+14
-7
lines changed

docs/src/solvers/SteadyStateSolvers.md

Lines changed: 14 additions & 7 deletions
Original file line numberDiff line numberDiff line change
@@ -1,27 +1,34 @@
11
# [Steady State Solvers](@id ss_solvers)
22

3-
`solve(prob::SteadyStateProblem,alg;kwargs)`
3+
`solve(prob::SteadyStateProblem, alg; kwargs)`
44

55
Solves for the steady states in the problem defined by `prob` using the algorithm
66
`alg`. If no algorithm is given, a default algorithm will be chosen.
77

88
## Recommended Methods
99

1010
Conversion to a NonlinearProblem is generally the fastest method. However, this will not
11-
guarantee the preferred root, and thus if the preferred root is required, then it's
12-
recommended that one uses `DynamicSS`. For `DynamicSS`, often an adaptive stiff solver,
13-
like a Rosenbrock or BDF method (`Rodas5` or `QNDF`), is a good way to allow for very large
14-
time steps as the steady state approaches.
11+
guarantee the preferred root (the stable equilibrium), and thus if the preferred root is
12+
required, then it's recommended that one uses `DynamicSS`. For `DynamicSS`, often an
13+
adaptive stiff solver, like a Rosenbrock or BDF method (`Rodas5` or `QNDF`), is a good way
14+
to allow for very large time steps as the steady state approaches.
1515

1616
!!! note
1717

1818
The SteadyStateDiffEq.jl methods on a `SteadyStateProblem` respect the time definition
19-
in the nonlinear definition, i.e., `u' = f(u,t)` uses the correct values for `t` as the
19+
in the nonlinear definition, i.e., `u' = f(u, t)` uses the correct values for `t` as the
2020
solution evolves. A conversion of a `SteadyStateProblem` to a `NonlinearProblem`
21-
replaces this with the nonlinear system `u' = f(u,∞)`, and thus the direct
21+
replaces this with the nonlinear system `u' = f(u, ∞)`, and thus the direct
2222
`SteadyStateProblem` approach can give different answers (i.e., the correct unique
2323
fixed point) on ODEs with non-autonomous dynamics.
2424

25+
!!! note
26+
27+
If you have an unstable equilibrium and you want to solve for the unstable equilibrium,
28+
then `DynamicSS` might converge to the equilibrium based on the initial condition.
29+
However, Nonlinear Solvers don't suffer from this issue, and thus it's recommended to
30+
use a nonlinear solver if you want to solve for the unstable equilibrium.
31+
2532
## Full List of Methods
2633

2734
### Conversion to NonlinearProblem

0 commit comments

Comments
 (0)