NAG Library Routine Document

E05JCF

▸▿ Contents

1 Purpose

2 Specification

3 Description

4 References

5 Parameters

6 Error Indicators and Warnings

7 Accuracy

8 Parallelism and Performance

9 Further Comments

▸▿ 10 Example

10.1 Program Text

10.2 Program Data

10.3 Program Results

1 Purpose

E05JCF may be used to supply optional parameters to E05JBF from an external file. The initialization routine E05JAF must have been called before calling E05JCF.

2 Specification

SUBROUTINE E05JCF (

IOPTS, COMM, LCOMM, IFAIL)

INTEGER	IOPTS, LCOMM, IFAIL
REAL (KIND=nag_wp)	COMM(LCOMM)

3 Description

E05JCF may be used to supply values for optional parameters to E05JBF. E05JCF 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 optional parameter must be separated by spaces, or equals signs

[=]

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

Static Limit = 100

is an example of a string used to set an optional parameter. For each optional parameter 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.

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 implied data type (character, integer or real) of each value to set must match that expected by the corresponding optional parameter.

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

Begin * Example options file
   Static Limit = 500 
End

Optional parameter settings are preserved following a call to E05JBF 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 E05JBF.

A complete list of optional parameters, their symbolic names and default values is given in Section 12 in E05JBF.

4 References

None.

5 Parameters

1: $IOPTS$ – INTEGERInput: On entry: the unit number of the option file to be read.

Constraint: IOPTS is a valid unit open for reading.
2: $COMM (LCOMM)$ – REAL (KIND=nag_wp) arrayCommunication Array: On exit: COMM must not be altered between calls to any of the routines E05JBF, E05JCF, E05JDF, E05JEF, E05JFF, E05JGF, E05JHF, E05JJF, E05JKF and E05JLF.
3: $LCOMM$ – INTEGERInput: On entry: the dimension of the array COMM as declared in the (sub)program from which E05JCF is called.

Constraint: $LCOMM \geq 100$ .
4: $IFAIL$ – INTEGERInput/Output: On entry: IFAIL must be set to $0$ , $- 1 or 1$ . If you are unfamiliar with this parameter you should refer to Section 3.3 in the Essential Introduction for details.
For environments where it might be inappropriate to halt program execution when an error is detected, the value $- 1 or 1$ is recommended. If the output of error messages is undesirable, then the value $1$ is recommended. Otherwise, if you are not familiar with this parameter, the recommended value is $0$ . When the value $- 1 or 1$ is used it is essential to test the value of IFAIL on exit.

On exit: $IFAIL = 0$ unless the routine detects an error or a warning has been flagged (see Section 6).

6 Error Indicators and Warnings

If on entry

IFAIL = 0

- 1

, explanatory error messages are output on the current error message unit (as defined by X04AAF).

Errors or warnings detected by the routine:

$IFAIL = 1$: Initialization routine E05JAF has not been called.

On entry, $LCOMM = 〈value〉$ .
Constraint: $LCOMM \geq 100$ .

$IFAIL = 2$: At least one optional parameter from the options file could not be recognized. All optional parameters that were set from the file before this error was encountered will remain set on exit.

BEGIN found, but end-of-file found before END. All optional parameters that were set from the file before this error was encountered will remain set on exit.

Could not read options file.

End-of-file found before BEGIN.

$IFAIL = 3$: Attempt to assign an illegal value of Local Searches ( $lcsrch$ ): $lcsrch = 〈value〉$ .

Attempt to assign an illegal value of Repeatability ( $repeat$ ): $repeat = 〈value〉$ .

Attempt to assign a non-positive value of Function Evaluations Limit ( $nf$ ): $nf = 〈value〉$ .

Attempt to assign a non-positive value of Local Searches Limit ( $loclim$ ): $loclim = 〈value〉$ .

Attempt to assign a non-positive value of Static Limit ( $stclim$ ): $stclim = 〈value〉$ .

Attempt to assign an out-of-bounds value of Infinite Bound Size ( $infbnd$ ): $infbnd = 〈value〉$ .

Attempt to assign too small a value of Local Searches Tolerance ( $loctol$ ): $loctol = 〈value〉$ .

Attempt to assign too small a value of Target Objective Error ( $objerr$ ): $objerr = 〈value〉$ .

Attempt to assign too small a value of Target Objective Safeguard ( $objsfg$ ): $objsfg = 〈value〉$ .

$IFAIL = 5$: One of the numeric values to be set could not be parsed. Check that all such strings specify valid integer or real values.

$IFAIL = - 99$: An unexpected error has been triggered by this routine. Please contact NAG.
See Section 3.8 in the Essential Introduction for further information.

$IFAIL = - 399$: Your licence key may have expired or may not have been installed correctly.
See Section 3.7 in the Essential Introduction for further information.

$IFAIL = - 999$: Dynamic memory allocation failed.
See Section 3.6 in the Essential Introduction for further information.

7 Accuracy

Not applicable.

8 Parallelism and Performance

E05JCF is threaded by NAG for parallel execution in multithreaded implementations of the NAG Library.

Please consult the X06 Chapter Introduction for information on how to control and interrogate the OpenMP environment used within this routine. Please also consult the Users' Note for your implementation for any additional implementation-specific information.

9 Further Comments

E05JDF, E05JEF, E05JFF or E05JGF may also be used to supply optional parameters to E05JBF.

10 Example

This example finds the global minimum of the ‘peaks’ function in two dimensions

F (x, y) = 3 {(1 - x)}^{2} \exp (- x^{2} - {(y + 1)}^{2}) - 10 (\frac{x}{5} - x^{3} - y^{5}) \exp (- x^{2} - y^{2}) - \frac{1}{3} \exp (- {(x + 1)}^{2} - y^{2})

on the box

[- 3, 3] \times [- 3, 3]

The function

F

has several local minima and one global minimum in the given box. The global minimum is approximately located at

(0.23, - 1.63)

, where the function value is approximately

- 6.55

By specifying an initialization list via LIST, NUMPTS and INITPT we can start E05JBF looking close to one of the local minima and check that it really does move away from that point to one of the global minima.

More precisely, we choose

(- 1, 0)

as our initial point (see Section 10.3), and let the initialization list be

(\begin{matrix} - 3 & - 1 & 3 \\ - 3 & 0 & 3 \end{matrix}) .

This example solves the optimization problem using some of the optional parameters described in Section 12 in E05JBF.

NAG Library Routine DocumentE05JCF