e04rxf is a part of the NAG optimization modelling suite. It allows a piece of real information to be read from or written to the problem stored in the handle. For example, it may be used to extract the current approximation of the solution during a monitoring step.

2 Specification

Fortran Interface

Subroutine e04rxf (

handle, cmdstr, ioflag, lrarr, rarr, ifail)

Integer, Intent (In)	::	ioflag
Integer, Intent (Inout)	::	lrarr, ifail
Real (Kind=nag_wp), Intent (Inout)	::	rarr(lrarr)
Character (*), Intent (In)	::	cmdstr
Type (c_ptr), Intent (In)	::	handle

C Header Interface

#include <nag.h>

void	e04rxf_ (void *handle, const char cmdstr, const Integer ioflag, Integer lrarr, double rarr[], Integer *ifail, const Charlen length_cmdstr)

The routine may be called by the names e04rxf or nagf_opt_handle_set_get_real.

3 Description

e04rxf adds an additional means of communication to routines within the NAG optimization modelling suite. It allows you to either read or write a piece of information in the handle in the form of a real array. The item is identified by cmdstr and the direction of the communication is set by ioflag.

The following cmdstr are available:

Primal Variables or X: The current value of the primal variables.
Dual Variables or U: The current value of the dual variables (Lagrangian multipliers for bounds, linear, quadratic and nonlinear constraints).
Dual Cone Variables or UC: The current value of the dual cone variables (Lagrangian multipliers for quadratic cone constraints).
Dual Matrix Variables or UA: The current value of the dual matrix variables (Lagrangian multipliers for matrix inequality constraints).
Residual Weights or RW: Array of residual weights.
Covariance Matrix or COV: Covariance matrix of the primal variables.
Hessian Matrix or HES: The current value of the Hessian matrix.
Variance or VAR: Variance for the primal variables.
MCS List or MCSL: Used to provide a custom initialization list for the solver e05kbf.
Slack or SLK: Slack variables.

The functionality is supported by all solvers within the NAG optimization modelling suite for the retrieval of their final solutions. e04rxf, moreover, can also be used by a selection of the solvers in the NAG Library, such as, e04mtf and e04ptf for the retrieval of the approximate solutions within the monitoring step. See the documentation of the individual solvers. In addition, you can retrieve solutions from a specific solver by adding the solver name after the available cmdstr mentioned above, e.g., keywords X e04mt or X handle_solve_lp_ipm return the value of primal variables only from the linear programming solver e04mtf. This is useful when a problem is solved by several different solvers. See Section 3.1 in the E04 Chapter Introduction for more details about the NAG optimization modelling suite.

4 References

None.

5 Arguments

1: $handle$ – Type (c_ptr) Input

On entry: the handle to the problem. It needs to be initialized by e04raf and must not be changed between calls to the NAG optimization modelling suite.

2: $cmdstr$ – Character(*) Input

On entry: a string which identifies the item within the handle to be read or written. The string is case insensitive and space tolerant.

Constraint:

cmdstr =' Primal Variables',' Dual Variables',' Dual Cone Variables',' Dual Matrix Variables',' Residual Weights',' Covariance Matrix',' Hessian Matrix',' Variance',' Slack',' X',' U',' UC',' UA',' RW',' COV',' HES',' VAR', or ​' SLK'

3: $ioflag$ – Integer Input

On entry: indicates the direction of the communication.

$ioflag \neq 0$: e04rxf will extract the requested information from the handle to rarr.
$ioflag = 0$: The writing mode will apply and the content of rarr will be copied to the handle.

4: $lrarr$ – Integer Input/Output

On entry: the dimension of the array rarr.

On exit: the correct expected dimension of rarr if lrarr does not match the item identified by cmdstr (in this case e04rxf returns

ifail = 2

5: $rarr (lrarr)$ – Real (Kind=nag_wp) array Input/Output

On entry: if

ioflag = 0

(write mode), rarr must contain the information to be written to the handle; otherwise, it does not need to be set.

On exit: if

ioflag \neq 0

(read mode), rarr contains the information requested by cmdstr; otherwise, rarr is unchanged.

6: $ifail$ – Integer Input/Output

On entry: ifail must be set to

0

−1

1

to set behaviour on detection of an error; these values have no effect when no error is detected.

A value of

0

causes the printing of an error message and program execution will be halted; otherwise program execution continues. A value of

−1

means that an error message is printed while a value of

1

means that it is not.

If halting is not appropriate, the value

−1

1

is recommended. If message printing is undesirable, then the value

1

is recommended. Otherwise, the value

−1

is recommended since useful values can be provided in some output arguments even when

ifail \neq 0

on exit. 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:

Note: in some cases e04rxf may return useful information.

$ifail = 1$: The supplied handle does not define a valid handle to the data structure for the NAG optimization modelling suite. It has not been properly initialized or it has been corrupted.

$ifail = 2$: On entry, $lrarr = ⟨ value ⟩$ , expected $value = ⟨ value ⟩$ .
Constraint: lrarr must match the size of the data identified in cmdstr.

$ifail = 3$: The provided cmdstr is not recognised.

$ifail = 4$: Reading mode is not supported for the given cmdstr.

$ifail = 5$: Writing mode is not supported for the given cmdstr.

$ifail = 6$: The request cannot be processed at this phase.

The requested information is not available.

$ifail = - 99$: An unexpected error has been triggered by this routine. Please contact NAG.
See Section 7 in the Introduction to the NAG Library FL Interface for further information.

$ifail = - 399$: Your licence key may have expired or may not have been installed correctly.
See Section 8 in the Introduction to the NAG Library FL Interface for further information.

$ifail = - 999$: Dynamic memory allocation failed.
See Section 9 in the Introduction to the NAG Library FL Interface for further information.

7 Accuracy

Not applicable.

8 Parallelism and Performance

Background information to multithreading can be found in the Multithreading documentation.

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

e04rxf makes calls to BLAS and/or LAPACK routines, which may be threaded within the vendor library used by this implementation. Consult the documentation for the vendor library for further information.

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

None.

10 Example

This example demonstrates how to use e04rxf to extract the current approximation of the solution when the monitoring routine monit is called during the solve by e04mtf.

We solve the following linear programming problem:

- 0.02 x_{1} - 0.2 x_{2} - 0.2 x_{3} - 0.2 x_{4} - 0.2 x_{5} + 0.04 x_{6} + 0.04 x_{7}

subject to the bounds

\begin{matrix} - 0.01 \leq x_{1} \leq 0.01 \\ - 0.1 \leq x_{2} \leq 0.15 \\ - 0.01 \leq x_{3} \leq 0.03 \\ - 0.04 \leq x_{4} \leq 0.02 \\ - 0.1 \leq x_{5} \leq 0.05 \\ - 0.01 \leq x_{6} \\ - 0.01 \leq x_{7} \end{matrix}

and the general constraints

\begin{matrix} x_{1} & + & x_{2} & + & x_{3} & + & x_{4} & + & x_{5} & + & x_{6} & + & x_{7} & = & - 0.13 \\ 0.15 x_{1} & + & 0.04 x_{2} & + & 0.02 x_{3} & + & 0.04 x_{4} & + & 0.02 x_{5} & + & 0.01 x_{6} & + & 0.03 x_{7} & \leq & - 0.0049 \\ 0.03 x_{1} & + & 0.05 x_{2} & + & 0.08 x_{3} & + & 0.02 x_{4} & + & 0.06 x_{5} & + & 0.01 x_{6} & \leq & - 0.0064 \\ 0.02 x_{1} & + & 0.04 x_{2} & + & 0.01 x_{3} & + & 0.02 x_{4} & + & 0.02 x_{5} & \leq & - 0.0037 \\ 0.02 x_{1} & + & 0.03 x_{2} & + & 0.01 x_{5} & \leq & - 0.0012 \\ - 0.0992 & \leq & 0.70 x_{1} & + & 0.75 x_{2} & + & 0.80 x_{3} & + & 0.75 x_{4} & + & 0.80 x_{5} & + & 0.97 x_{6} \\ - 0.003 & \leq & 0.02 x_{1} & + & 0.06 x_{2} & + & 0.08 x_{3} & + & 0.12 x_{4} & + & 0.02 x_{5} & + & 0.01 x_{6} & + & 0.97 x_{7} & \leq & 0.002 . \end{matrix}

During the monitoring step of e04mtf, if the three convergence measures are below an acceptable threshold, the approximate solution is extracted with e04rxf and printed on the standard output.

e04rx: FL CL CPP AD PY MB

NAG FL Interfacee04rxf (handle_​set_​get_​real)

▸▿ Contents

1 Purpose

2 Specification

3 Description

4 References

5 Arguments

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

NAG FL Interface
e04rxf (handle_set_get_real)