About BF options
How to create the window for displaying the change of eigenvalues every time BF runs
When calculating a bifurcation set on an arbitrary parameter plane using BF, a figure that shows the change of the positions of eigenvalues in the Gauss plane, i.e., root loci, is created in addition to the figure of a 2parameter plane to show calculation results in real time. The figure is to monitor whether bifurcation points with codimension 2 occur during calculation of the bifurcation set. In the default settings, the previous results are cleared every time BF runs, and the root loci are shown in the same window. You may want to see differences in the root loci under different parameter settings. In this case, change the configuration such that the display window starts every time BF runs.
How to change the convergent precision of the Newton's method
When bifurcation sets are calcurated, information regarding precise coordinates of equilibrium points or periodic points of the Poincare map is necessary. Newton¡Çs method is used to obtain these coordinates. The convergence precision is directly correlated to the accuracy of the calculations to obtain the bifurcation set. During calculations, if the convergence precision of Newton¡Çs method is inaccurate, it may give rise to problems such as the following:
In the default settings, when BF is initially started, the convergence precision is set to 1E6. If a discontinuous point is included the set when a bifurcation set is tracking, users must be .... The behavior of BF may improve by increasing the convergence precision. The following steps are used to change the convergence precision of Newton¡Çs method.
How to change the divergent condition of the Newton's method
Currently, the following are used to determine convergence or divergence as the divergence criterion of Newton¡Çs method*1:
In particular, if the Newton¡Çs method does not converge in the number of iterations iter, Newton¡Çs method considers the system divergent or no solution. Therefore, iter must be appropriately changed depending on the problem. gmax and emax are used to terminate the program when abnormal values are detected internally. In the default settings, when BF is initially started, values that are considered appropriate in general are configured:
The following steps are used to change the divergence criteria for Newton¡Çs method.
How to control the automatic switching of parameters
In calculating a bifurcation set on an arbitrary twoparameter plane, one unknown parameter and one control parameter are used. The bifurcation parameter value is calculated when the control parameter is fixed. Next, similar calculations are conducted when the control parameter is slightly changed. The procedure is repeated to obtain a set of parameters that result in bifurcation in the designated parameter plane. Here, if the curvature of the parameter set is very large, the bifurcation set may be efficiently obtained by switching the configured unknown and control parameters. In the actual calculations, the curvature is calculated every time a bifurcation point is obtained. auto sw gives the maximum curvature. If the curvature of the bifurcation set exceeds this value, the unknown and control parameters are swiched and the calculation of the bifurcation set is continued. The following steps are used to change the criteria for parameter switching.
How to stop the calculation of the inverse set
The bifurcation set is tracked from an initial bifurcation point. However, unless the set is ringlike, the calculation of the bifurcation set fails at some parameter. There are many possible causes, including:
Here, after the calculations of a bifurcation set stopped, the next step is to calculate the bifurcation set from the initial bifurcation point with the parameters changing in the other direction. In the default settings of BF, after the calculation of the bifurcation set stopped for some reason, the bifurcation set is automatically tracked from the initial bifurcation point with the parameters changing in the other direction. In many cases, the user only needs to obtain the bifurcation set in a given range, and therefore in some cases, calculation of the bifurcation set in the other direction is not necessary. Therefore, to stop the automatic tracking of the bifurcation set in the other direction, only the following switch needs to be turned off.
How to change the numerical integration method(Caution !)
It is necessary to numerically integrate the solution when searching for bifurcation points of limit cycles and periodic solutions. Therefore, selecting what solver to use for numerical integration is an important issue. There may be some situations such that a stiff solver must be selected. Basically, to find stable limit cycles or periodic solutions, simulation using PP is conducted to obtain attractor information. Information on the solver used at this stage is passed directly to FIX, therefore the numerical calculation is safer when the solver is not changed at this stage as much as possible. However, if the convergence of the Newton¡Çs method is bad, or does not converge, the match between the solver used and the problem may not be good. Then, it may need to try searching bifurcation points after changing the solver. The following steps are used to change the solver.
How to improve the precision of numerical integration(Caution !)
As seen in How to change the numerical integration method, it is necessary to numerically integrate the solution when searching for bifurcation points of limit cycles and periodic solutions. The precision of the numerical integration is an important issue. Basically, to find for stable limit cycles or periodic solutions, simulation using PP is conducted to obtain attractor information. Information on the solver used at this stage (including precision) is passed directly to FIX, and therefore the numerical calculation is safer when the precision is not changed at this stage as much as possible. However, if the convergence of Newton¡Çs method is bad, or does not converge, the configured precision may be too low. Then it may need to try searching bifurcation points after changing the precision. The following steps are used to change the precision.
