Next: The Continuation State Window
Up: Attributes of Interface Windows
Previous: The Periodic Orbits Window
Contents
The Equilibrium Continuation Window
Figure 2.26:
The Continuation window.
 |
- Window title:
- DsTool: Equilibrium Continuation
- Function:
- The Equilibrium Continuation window allows the user to compute curves of bifurcations.
- Description:
- The Equilibrium Continuation window is opened by selecting
the Equilibrium Continuation option from the Panels menu button
located in the Command window. The window allows the
user to compute curves of equilibrium points with one
varying parameter and curves of saddle-node and Hopf bifurcations
with two varying parameters. A brief tutorial on
continuation calculations is at the end of this section.
The different Hopf bifurcation algorithms are described in
[2] and [3].
Color coding for hyperbolic equlibrium points
is performed according to the
value of the Monitor switch setting. For the
bifurcation types, color coding is performed according to the choice of
colormap, as shown in the Continuation Colors window.
In the default colormap, Hopf bifurcation
points are displayed in magenta, saddle-node bifurcation points are
displayed in green, degenerate Hopf points are displayed in
orange, and resonant saddle-node points are displayed in sea
green.
- Panel items:
-
- Iterations numeric field:
Displays the number of points along
a curve of equilibrium points that will be computed. The
type of curve that will be computed is determined by the
Mode stack setting. The default value is 10.
- Monitor switch exclusive setting:
Allows the user to choose whether or not to color equilibrium points with the
DsTool convention: blue
for sink, green for saddle and red for source.
The default value is on.
- Mode stack setting:
Allows the user to choose the system of equations to be solved with
continuation. The choices are:
- Static Bifurcation:
-
for continuing curves of
equilibrium points with a single active parameter.
- Saddle Node (
) :
-
for computing curves
of saddle-node bifurcations with two active parameters.
The determinant of the Jacobian is computed and used
as an augmenting equation to the equilibrium equations.
- Saddle Node Bifurcation:
- for computing curves
of saddle-node bifurcations with two active parameters.
- Hopf Bif (
) :
-
for computing curves of
Hopf bifurcations with two active parameters, using
the bialternate product.
- Hopf Bif (
) :
- for computing curves of
Hopf bifurcations with two active parameters, using
the Bezoutian.
- Hopf Bif (JGR):
- for computing curves of
Hopf bifurcations with two active parameters, using
the algorithm developed by Jepson, Griewank and Reddien.
- Hopf Bif (kubicek 1):
- for computing curves of
Hopf bifurcations with two active parameters, using
the algorithm developed by Kubicek.
- Aug params listbox:
Allows the user to choose active parameters. Static Bifurcation
calculations require that the user select one active
parameter. Saddle-node and Hopf bifurcation calculations
require the user to select two active parameters. An error message will
be printed if the user tries to select too many active parameters, or
tries to initiate a calculation with too few active parameters.
- Cont. param stack setting:
Allows the user to choose one of the phase space
variables or one of the active parameters to parameterize
the curve of equilibrium points.
- Parameter fix switch exclusive setting:
Allows the user to choose whether the setting in Cont. param stack may be changed
during the calculation. The choice Vary allows the continuation algorithm to change
the variable used as the continuation parameter adaptively.
The choice Fix uses requires the continuation algorithm to
use only the variable selected in the Cont. param
stack setting. The default is the ``Vary'' setting.
- Jacobian update switch exclusive setting:
Allows the user to choose to update the Jacobian of the augmented system
at each step or only at the initial continuation step. The default is to update the
Jacobian at each step.
- Absolute error read-write text field:
Displays the absolute error tolerance used in solving the augmented equations.
The default value is
.
- Relative error read-write text field:
Displays the relative error tolerance used in
solving the augmented equations. The default value is
.
- Minimum step read-write text field:
Displays the minimum step size used by the
continuation algorithm for continuation steps of the
continuation parameter. The default value is
.
- Maximum step read-write text field:
Displays the maximum step size used by the
continuation algorithm for continuation steps of the
continuation parameter. The default value is
.
- Sugg. step read-write text field:
Displays the step size used by the
continuation algorithm for the initial continuation step of the
continuation parameter. The default value is
.
- Debug level numeric field:
Displays an integer for determining the
amount of diagnostic information printed to the terminal
window from which DsTool is run. The default value
is
, which prints no diagnostic information to the terminal.
- Forwards command button:
Initiates a continuation calculation in the direction of
increasing values of the continuation parameter for the number
of steps in the Iterations field.
- Backwards command button:
Initiates a continuation calculation in the direction of
decreasing values of the continuation parameter for the number
of steps in the Iterations field.
- Continue command button:
Continues the last continuation calculation for the number
of steps in the Iterations field.
- Search command button:
Uses Monte Carlo seeds to locate solutions of the augmented
equations. The number of seeds used is determiend by the
value of the Iterations text field.
- State window button:
Opens and brings to the foreground the Continuation State window.
- Dismiss command button:
Closes the Equilibrium Continuation window.
Appendix: Continuation Calculations
One parameter continuation is a systematic strategy for computing
curves of solutions to
equations
in
variables. The mathematical foundation for continuation
algorithms is the implicit function theorem. If
is differentiable
and the Jacobian of
is surjective at a point
where
vanishes, then the solutions form a curve in a neighborhood of
. Moreover, the implicit function theorem gives a formula for the
tangent vector to the solution curve. Continuation algorithms exploit
this information. They use an initial value solver for ordinary
differential equations (often just the Euler method) to step along an
approximation to the solution curve. They then appply a root finding
algorithm to significantly improve the approximation. The alternation
of root finding with numerical integration steps distinguish
continuation methods. Choices of how to parametrize the solution
curve, choose time steps and restrict the equations to a hyperplane to
obtain a square system of equations with a unique solution (locally)
need to be made. Various continuation packages take different
approaches to these matters: the continuation ``engine'' used by DsTool
is PITCON (Rheinboldt ...).
If
is a
parameter family of vector fields,
then we assemble systems of equations
for varied calculations.
To do so we restrict ourselves to a submanifold of dimension
.
Most frequently this submanifold is obtained by fixing
of the parameters of
(called inactive parameters) and varying
active parameters. In the simplest case,
, the number of equations is
,
and the continuation calculation computes curves of equilibria
with a single active parameter. To compute curves of bifurcations, one adds
defining equations to the system of equations
to produce an
augmented system. The number of independent defining equations added
to
is the codimension of the bifurcation.
(In some circumstances, the augmented system uses
additional auxiliary variables such as eigenvalues of the Jacobian
and has a corresponding number of additional equations.)
The continuation window of DsTool includes calculations for saddle-node and
Hopf bifurcations. These are the codimension one bifurcations of
equilibrium points.
For saddle-node bifurcations, the defining equation is
or an algorithm that computes another scalar quantity, such as the
smallest singular value, that vanishes precisely when
is singular.
For Hopf bifurcations, the defining equation(s) compute where the
matrix
has a pair of pure imaginary eigenvalues. Explicit
expressions
in the coefficients of
that compute whether
has a pair of
complex eigenvalues are very complicated. Alternate ways of performing
this computation are discussed in [3].
There are also other algorithms that perform the calculation
by introducing auxiliary variables for the pure imaginary eigenvalues,
and in some cases the eigenvectors of
.
Next: The Continuation State Window
Up: Attributes of Interface Windows
Previous: The Periodic Orbits Window
Contents
2010-01-04