NAG Library Routine Document

E04UHF/E04UHA

+− Contents

1 Purpose

+− 2 Specification

2.1 Specification for E04UHF

2.2 Specification for E04UHA

3 Description

4 References

5 Parameters

6 Error Indicators and Warnings

7 Accuracy

8 Further Comments

+− 9 Example

9.1 Program Text

9.2 Program Data

9.3 Program Results

1 Purpose

To supply optional parameters to E04UGF/E04UGA from an external file. More precisely, E04UHF must be used to supply optional parameters to E04UGF and E04UHA must be used to supply optional parameters to E04UGA.

E04UHA is a version of E04UHF that has additional parameters in order to make it safe for use in multithreaded applications (see Section 5). The initialization routine E04WBF must have been called before calling E04UHA.

2 Specification

2.1 Specification for E04UHF

SUBROUTINE E04UHF (

IOPTNS, INFORM)

INTEGER

IOPTNS, INFORM

2.2 Specification for E04UHA

SUBROUTINE E04UHA (

IOPTNS, LWSAV, IWSAV, RWSAV, INFORM)

INTEGER	IOPTNS, IWSAV(550), INFORM
REAL (KIND=nag_wp)	RWSAV(550)
LOGICAL	LWSAV(20)

3 Description

E04UHF/E04UHA may be used to supply values for optional parameters to E04UGF/E04UGA. E04UHF/E04UHA reads an external file and each line of the file defines a single optional parameter. It is only necessary to supply values for those parameters whose values are to be different from their default values.

Each optional parameter is defined by a single character string, of up to

72

characters, consisting of one or more items. The items associated with a given option must be separated by spaces, or equals signs

[=]

. Alphabetic characters may be upper or lower case. The string

Print Level = 1

is an example of a string used to set an optional parameter. For each option the string contains one or more of the following items:

–	a mandatory keyword;
–	a phrase that qualifies the keyword;
–	a number that specifies an integer or real value. Such numbers may be up to $16$ contiguous characters in Fortran's I, F, E or D formats, terminated by a space if this is not the last item on the line.

Blank strings and comments are ignored. A comment begins with an asterisk (*) and all subsequent characters in the string are regarded as part of the comment.

The file containing the options must start with Begin and must finish with End. An example of a valid options file is:

Begin * Example options file 
   Print level = 5 
End

For E04UHF each line of the file is normally printed as it is read, on the current advisory message unit (see X04ABF), but printing may be suppressed using the keyword Nolist. To suppress printing of Begin, Nolist must be the first option supplied as in the file:

Begin 
   Nolist 
   Print level = 5 
End

Printing will automatically be turned on again after a call to E04UGF or E04UHF and may be turned on again at any time using the keyword List.

For E04UHA printing is turned off by default, but may be turned on at any time using the keyword List.

Optional parameter settings are preserved following a call to E04UGF/E04UGA and so the keyword Defaults is provided to allow you to reset all the optional parameters to their default values before a subsequent call to E04UGF/E04UGA.

A complete list of optional parameters, their abbreviations, synonyms and default values is given in Section 11 in E04UGF/E04UGA.

4 References

Hock W and Schittkowski K (1981) Test Examples for Nonlinear Programming Codes. Lecture Notes in Economics and Mathematical Systems 187 Springer–Verlag

5 Parameters

1: IOPTNS – INTEGERInput: On entry: the unit number of the options file to be read.

Constraint: $0 \leq IOPTNS \leq 99$ .
2: INFORM – INTEGEROutput: Note: for E04UHA, INFORM does not occur in this position in the parameter list. See the additional parameters described below.
On exit: contains zero if the options file has been successfully read and a $value > 0$ otherwise (see Section 6).; Note: the following are additional parameters for specific use with E04UHA. Users of E04UHF therefore need not read the remainder of this description.
2: LWSAV( $20$ ) – LOGICAL arrayCommunication Array
3: IWSAV( $550$ ) – INTEGER arrayCommunication Array
4: RWSAV( $550$ ) – REAL (KIND=nag_wp) arrayCommunication Array: The arrays LWSAV, IWSAV and RWSAV must not be altered between calls to any of the routines E04UHA, E04UGA, E04UJA or E04WBF.
5: INFORM – INTEGEROutput: Note: see the parameter description for INFORM above.

6 Error Indicators and Warnings

Errors or warnings detected by the routine:

$INFORM = 1$: IOPTNS is not in the range $[0, 99]$ .

$INFORM = 2$: Begin was found, but end-of-file was found before End was found.

$INFORM = 3$: end-of-file was found before Begin was found.

$INFORM = 4$: Not used.

$INFORM = 5$: One or more lines of the options file is invalid. Check that all keywords are neither ambiguous nor misspelt.

7 Accuracy

Not applicable.

8 Further Comments

E04UJF/E04UJA may also be used to supply optional parameters to E04UGF/E04UGA.

9 Example

This is Problem 45 from Hock and Schittkowski (1981) and involves the minimization of the nonlinear function

f (x) = 2 - \frac{1}{120} \times x_{1} x_{2} x_{3} x_{4} x_{5}

subject to the bounds

\begin{matrix} 0 \leq x_{1} \leq 1, \\ 0 \leq x_{2} \leq 2, \\ 0 \leq x_{3} \leq 3, \\ 0 \leq x_{4} \leq 4, \\ 0 \leq x_{5} \leq 5 . \end{matrix}

The initial point, which is infeasible, is

x_{0} = {(2, 2, 2, 2, 2)}^{T},

and

f (x_{0}) = 1.7333

(to five figures).

The optimal solution is

x^{*} = {(1, 2, 3, 4, 5)}^{T},

and

f (x^{*}) = 1

. All the bounds are active at the solution.

In this example the options file read by E04UHF/E04UHA is appended to the data file for the program (see Section 9.2). It would usually be more convenient in practice to keep the data file and the options file separate.