For each option, we give a summary line, a description of the optional parameter and details of constraints.
The summary line contains:
- the keywords, where the minimum abbreviation of each keyword is underlined (if no characters of an optional qualifier are underlined, the qualifier may be omitted);
- a parameter value,
where the letters , denote options that take character, integer and real values respectively;
- the default value, where the symbol is a generic notation for machine precision (see x02aj).
Keywords and character values are case and white space insensitive.
This option specifies how the initial working set is chosen. With a
Cold Start,
e04nc chooses the initial working set based on the values of the variables and constraints at the initial point. Broadly speaking, the initial working set will include equality constraints and bounds or inequality constraints that violate or ‘nearly’ satisfy their bounds (to within
Crash Tolerance).
With a
Warm Start, you must provide a valid definition of every element of the array
istate.
e04nc will override your specification of
istate if necessary, so that a poor choice of the working set will not cause a fatal error. For instance, any elements of
istate which are set to
,
will be reset to zero, as will any elements which are set to
when the corresponding elements of
bl and
bu are not equal. A warm start will be advantageous if a good estimate of the initial working set is available – for example, when
e04nc is called repeatedly to solve related problems.
This value is used in conjunction with the optional parameter
Cold Start (the default value) when
e04nc selects an initial working set. If
, the initial working set will include (if possible) bounds or general inequality constraints that lie within
of their bounds. In particular, a constraint of the form
will be included in the initial working set if
. If
or
, the default value is used.
This special keyword may be used to reset all optional parameters to their default values.
Feasibility Phase Iteration Limit | | |
Optimality Phase Iteration Limit | | |
The scalars
and
specify the maximum number of iterations allowed in the feasibility and optimality phases. Optional parameter
Optimality Phase Iteration Limit is equivalent to optional parameter
Iteration Limit. Setting
and
means that the workspace needed will be computed and printed, but no iterations will be performed. If
or
, the default value is used.
If , defines the maximum acceptable absolute violation in each constraint at a ‘feasible’ point. For example, if the variables and the coefficients in the general constaints are of order unity, and the latter are correct to about decimal digits, it would be appropriate to specify as . If , the default value is used.
Note that a ‘feasible solution’ is a solution that satisfies the current constraints to within the tolerance .
This option controls the contents of the upper triangular matrix
(see the description of
a in
[Parameters]).
e04nc works exclusively with the transformed and reordered matrix
(8), and hence extra computation is required to form the Hessian itself. If
,
a contains the Cholesky factor of the matrix
with columns ordered as indicated by
kx (see
[Parameters]). If
,
a contains the Cholesky factor of the matrix
, with columns ordered as indicated by
kx.
If , defines the ‘infinite’ bound in the definition of the problem constraints. Any upper bound greater than or equal to will be regarded as (and similarly any lower bound less than or equal to will be regarded as ). If , the default value is used.
If , specifies the magnitude of the change in variables that will be considered a step to an unbounded solution. (Note that an unbounded solution can occur only when the Hessian is singular and the objective contains an explicit linear term.) If the change in during an iteration would exceed the value of , the objective function is considered to be unbounded below in the feasible region. If , the default value is used.
Normally each optional parameter specification is printed as it is supplied. Optional parameter
Nolist may be used to suppress the printing and optional parameter
List may be used to restore printing.
If
and
, monitoring information produced by
e04nc at every iteration is sent to a file with logical unit number
. If
and/or
, no monitoring information is produced.
The value of
controls the amount of printout produced by
e04nc, as indicated below. A detailed description of the printed output is given in
[Description of the Printed Output] (summary output at each iteration and the final solution) and
[Description of Monitoring Information] (monitoring information at each iteration).
The following printout is sent to the current advisory message unit (as defined by
(X04ABF not in this release)):
| Output |
| No output. |
| The final solution only. |
| One line of summary output ( characters; see [Description of the Printed Output]) for each iteration (no printout of the final solution). |
| The final solution and one line of summary output for each iteration. |
The following printout is sent to the logical unit number defined by the optional parameter
Monitoring File:
| Output |
| No output. |
| One long line of output ( characters; see [Description of Monitoring Information]) for each iteration (no printout of the final solution). |
| At each iteration, the Lagrange multipliers, the variables , the constraint values and the constraint status. |
| At each iteration, the diagonal elements of the matrix associated with the factorization (4) (see [Definition of Search Direction]) of the working set, and the diagonal elements of the upper triangular matrix . |
If
and the unit number defined by the optional parameter
Monitoring File is the same as that defined by
(X04ABF not in this release), then the summary output is suppressed.
This option specifies the type of objective function to be minimized during the optimality phase. The following are the nine optional keywords and the dimensions of the arrays that must be specified in order to define the objective function:
LP | a and b not referenced, ; |
QP1 | symmetric, b and cvec not referenced; |
QP2 | symmetric, b not referenced, ; |
QP3 | upper trapezoidal, , b and cvec not referenced; |
QP4 | upper trapezoidal, , b not referenced, ; |
LS1 | , , cvec not referenced; |
LS2 | , , ; |
LS3 | upper trapezoidal, , , cvec not referenced; |
LS4 | upper trapezoidal, , , . |
For problems of type FP, the objective function is omitted and
a,
b and
cvec are not referenced.
The following keywords are also acceptable. The minimum abbreviation of each keyword is underlined.
| Option |
Least | LS1 |
Quadratic | QP2 |
Linear | LP |
In addition, the keywords LS and LSQ are equivalent to the default option LS1, and the keyword QP is equivalent to the option QP2.
If
, i.e., the objective function is purely linear, the efficiency of
e04nc may be increased by specifying
as LP.
Note that this option does not apply to problems of type FP or LP.
The default value of depends on the problem type. If occurs as a least squares matrix, as it does in problem types QP1, LS1 and LS3, then the default value of is . In all other cases, is treated as the ‘square root’ of the Hessian matrix and has the default value .
This parameter enables you to control the estimate of the triangular factor
(see
[Main Iteration]). If
denotes the function
, the rank of
is defined to be smallest index
i such that
. If
, the default value is used.