This object serves for storing parameters and data used in numerical continuation (for more details about the continuation see the Getfem++ user documentation).
S = gf_cont_struct(model md, string dataname_parameter[,string dataname_init, string dataname_final, string dataname_current], scalar sc_fac[, ...])
General constructor for cont_struct objects.
This object serves for storing parameters and data used in numerical continuation (for more details about the continuation see the Getfem++ user documentation).
S = gf_cont_struct(model md, string dataname_parameter[,string dataname_init, string dataname_final, string dataname_current], scalar sc_fac[, ...])
The variable dataname_parameter
should parametrise the model given by
md
. If the parametrisation is done via some vector datum,
dataname_init
and dataname_final
should store two given values of
this datum determining the parametrisation, and dataname_current
serves for actual values of this datum. sc_fac
is a scale factor
involved in the norm used in the continuation.
Additional options:
- 'lsolver', string SOLVER_NAME
name of the solver to be used for the incorporated linear systems
(the default value is 'auto', which lets getfem choose itself);
possible values are 'superlu', 'mumps' (if supported), 'cg/ildlt',
'gmres/ilu' and 'gmres/ilut';
- 'max_iter', int NIT
maximum number of iterations allowed in the correction (the default
value is 10);
- 'thr_iter', int TIT
threshold number of iterations of the correction for enlarging the
step size (the default value is 8);
- 'max_res', scalar RES
target residual value of the new point (the default value is 1e-6);
- 'max_diff', scalar DIFF
determines a convergence criterion to the new tangent vector (the
default value is 1e-9);
- 'min_ang', scalar ANG
minimal value of the cosine of the angle between tangents to the
solution curve at the old point and the new one (the default value
is 0.9);
- 'h_init', scalar HIN
initial step size (the default value is 1e-2);
- 'h_max', scalar HMAX
maximal step size (the default value is 1e-1);
- 'h_min', scalar HMIN
minimal step size (the default value is 1e-5);
- 'h_inc', scalar HINC
factor for enlarging the step size (the default value is 1.3);
- 'h_dec', scalar HDEC
factor for diminishing the step size (the default value is 0.5);
- 'epsilon', scalar EPS
increment to be used to compute the incorporated finite
differences (the default value is 1e-8);
- 'max_res_solve', scalar RES_SOLVE
target residual value for the linear systems to be solved (the
default value is 1e-7);
- 'noisy' or 'very_noisy'
determines how detailed information has to be displayed during the
process (residual values etc.).
Y. Collette