set options for dae solver
Sets a variable
%DAEOPTIONS that is used in function
to set solver options.
This daeoptions function interactively displays a window which sets various options of the dae solver.
%DAEOPTIONS is a list with the following elements:
The default value is:
The meaning of the elements is described below.
if it is 0,
daewill use the dassl/dasrt solver, and if it is 1,
daewill use daskr.
a real scalar or an empty matrix, gives the maximum time for which
gis allowed to be evaluated. An empty matrix means "no limits" imposed for time.
if it is 0,
daereturns only the user specified time point values, and if it is 1,
daereturns its intermediate computed values.
a two components vector which give the definition
[ml,mu]of band matrix computed by
r(i - j + ml + mu + 1,j)=
jacreturns a full matrix set
a scalar or an empty matrix, the maximum step size, empty matrix means "no limitation".
a scalar or an empty matrix, the minimum step size, empty matrix means "not specified".
a scalar, must be set to 0 if the solution is known to be non negative . In the other case it must be set to 1.
a scalar, must be set to 0 if the given initial condition is compatible:
xdot0is just an estimation, then set to 1 if you are using dassl, or set to
[+-1, ..., +-1], with:
1if y(i) is a differential variable and
-1if y(i) is an algebraic variable (if its derivatives do not appear explicitly in the function g(t, y, ydot)).
a scalar, direct / Krylov method. Set to
1and provide a routine
dae()if you want the solver to use Krylov iterations, else (daskr's direct method) set to
a vector, treat as dummy argument if you have set method=0. Otherwise, set to
[maxl kmp nrmax epli], where: -
maxl= maximum number of iterations in the GMRes algorithm (default min(5, neq)), -
kmp= number of vectors on which orthogonalization is done in the GMRes algorithm (default maxl), -
nrmax= maximum number of restarts of the GMRes algorithm per nonlinear iteration (default 5), -
epli= convergence test constant in GMRes algorithm (default 0.05).
a scalar, treat as dummy argument if you have set isest=0. Set to 1 if the solver should stop right after computation of the initial conditions, else set to 0.
a scalar, preconditioner computation and LU-factorization routine for
psol. Treat as dummy argument if method=0. Set to 1 and provide a
daeif the external
psolshould use a specific routine, else set to 0.
a vector, treat as dummy argument if you have set isest=0if you wish to control errors locally on all the variables then set to . Otherwise, set to [+-1, ..., +-1], 1 if
y(i)is a differential variable and -1 if
y(i)is an algebraic variable (if its derivatives do not appear explicitly in the function
g(t, y, ydot)).
a vector, heuristic parameters. Treat as dummy argument if isest=0. Otherwise, set to  or to
[mxnit mxnj mxnh lsoff stptol epinit], where: -
mxnit= maximum number of Newton iterations per Jacobian or preconditioner evaluation (default 5 if info(8)=0, 15 otherwise), -
mxnj= maximum number of Jacobian or preconditioner evaluations (default 6 if info(8)=0, 2 otherwise), -
mxnh= maximum number of values of the artificial stepsize parameter h to be tried if info(7) ≠  (default 5), -
lsoff= flag to turn off the linesearch algorithm (lsoff = 0 means linesearch is on, lsoff = 1 means it is turned off) (default 0), -
stptol= minimum scaled step in linesearch algorithm (default (unit roundoff)^(2/3)), -
epinit= swing factor in the Newton iteration convergence test (default 0.01).
a scalar, verbosity. Set to 1 for minimal extra printing, 2 for full printing, else set to 0.
- dae — Differential algebraic equations solver
|Report an issue|
|<< dae||Differential calculus, Integration||daskr >>|