NAG FL Interface
e04mzf (qpconvex1_​sparse_​mps)

Note: this routine is deprecated. Replaced by e04mxf.

Old Code

mpslst = .False.    ! Or = .True.
Call e04mzf(infile,maxn,maxm,maxnnz,xbldef,xbudef,mpslst,n,m,nnz,iobj,   &
            ncolh,a,irowa,iccola,bl,bu,start,pnames,nname,crname,xs,     &
            istate,ifail)

if (ifail == 1 .Or. ifail == 2 .Or. ifail == 3) then
  ! not enough memory, allocate bigger arrays as given in m, n, nnz
  ! and call e04mzf again
else if (ifail >= 4 .And. ifail <= 16) then
  ! mps input file formating error
  ! stop
else if (ifail == 17) then
  ! wrong arguments to e04mzf
  ! stop
else if (ifail == 0) then
  ! data successfully read, call solver
end if

New Code

mpslst = 0      ! or = 1
Call e04mxf(infile,maxn,maxm,maxnnz,maxncolh,maxnnzh,maxlintvar,mpslst,  &
            n,m,nnz,ncolh,nnzh,lintvar,iobj,a,irowa,iccola,bl,bu,pnames, &
            nname,crname,h,irowh,iccolh,minmax,intvar,ifail)

if (ifail == 2) then
  ! not enough memory, allocate bigger arrays as given in m, n, nnz, ncolh,
  ! nnzh, lintvar and call e04mxf again
else if (ifail >= 3 .And. ifail <= 35) then
  ! MPS input file formatting error
  ! stop
else if (ifail == 36) then
  ! wrong input argument
  ! stop
else if (ifail=-999) then
  ! internal memory allocation error
  ! stop
else if (ifail == 0 .Or. ifail == 1) then
  ! data successfully read (with possible warning)
  start = 'C'
  do j=1,n
    xs(j) = min(max(0.0_nag_wp,bl(j)),bu(j))
    istate(j) = 0.0_nag_wp
  end do
  ! call solver
end if
e04mxf has extended the functionality of e04mzf and the interface has changed substantially. If there are Integer variables, a quadratic part of the objective function or OBJSENSE section (see Section 3.1 in e04mxf), e04mxf will read them and return them in the new arguments (lintvar, intvar, ncolh, nnzh, h, irowh, iccolh and minmax), with e04mzf these caused a file formatting error. The new routine e04mxf might also accept a slightly misformatted input file and return a warning ifail=1.
The type of the argument mpslst has changed from logical to integer.
The parameters xbldef and xbudef of e04mzf were removed and fixed in e04mxf to their default values 0 and 1020, respectively. Note that value 1020 within bounds is interpreted in our solvers as + (unconstrained).
Parameters start, xs and istate were present in e04mzf only for the convenience of calling the solver routine e04nkf/​e04nka and have been removed from e04mxf.
Should you need any assistance, please do not hesitate to contact the NAG Technical Support Service.
Settings help

FL Name Style:


FL Specification Language:


1 Purpose

e04mzf reads data for a sparse linear programming or quadratic programming problem from an external file which is in standard or compatible MPSX input format.

2 Specification

Fortran Interface
Subroutine e04mzf ( infile, maxn, maxm, maxnnz, xbldef, xbudef, mpslst, n, m, nnz, iobj, ncolh, a, ha, ka, bl, bu, start, names, nname, crname, xs, istate, ifail)
Integer, Intent (In) :: infile, maxn, maxm, maxnnz
Integer, Intent (Inout) :: ifail
Integer, Intent (Out) :: n, m, nnz, iobj, ncolh, ha(maxnnz), ka(maxn+1), nname, istate(maxn+maxm)
Real (Kind=nag_wp), Intent (In) :: xbldef, xbudef
Real (Kind=nag_wp), Intent (Out) :: a(maxnnz), bl(maxn+maxm), bu(maxn+maxm), xs(maxn+maxm)
Logical, Intent (In) :: mpslst
Character (1), Intent (Out) :: start
Character (8), Intent (InOut) :: names(5)
Character (8), Intent (Out) :: crname(maxn+maxm)
C Header Interface
#include <nag.h>
void  e04mzf_ (const Integer *infile, const Integer *maxn, const Integer *maxm, const Integer *maxnnz, const double *xbldef, const double *xbudef, const logical *mpslst, Integer *n, Integer *m, Integer *nnz, Integer *iobj, Integer *ncolh, double a[], Integer ha[], Integer ka[], double bl[], double bu[], char *start, char names[], Integer *nname, char crname[], double xs[], Integer istate[], Integer *ifail, const Charlen length_start, const Charlen length_names, const Charlen length_crname)
The routine may be called by the names e04mzf or nagf_opt_qpconvex1_sparse_mps.

3 Description

e04mzf reads Linear Programming (LP) or Quadratic Programming (QP) problem data from an external file which is prepared in standard or compatible MPSX (see IBM (1971)) input format and then initializes n (the number of variables), m (the number of general linear constraints), the m×n matrix A, and the vectors l, u and c (stored in row iobj of A) for use with e04nkf, which is designed to solve problems of the form
minimize xRn cTx+12xTHx  subject to  l{ x Ax } u.  
For LP problems, H=0. For QP problems, you must set ncolh>0 (see Section 5) and provide a subroutine to e04nkf to compute Hx for any given vector x. The optional parameter Maximize may be used to specify an alternative problem in which the objective function is maximized (see Section 12.1 in e04nkf/​e04nka).
MPSX input format
The input file of data may only contain two types of lines:
  1. 1.Indicator lines (specifying the type of data which is to follow).
  2. 2.Data lines (specifying the actual data).
The input file must not contain any blank lines. Any characters beyond column 80 are ignored. Indicator lines must not contain leading blank characters (in other words they must begin in column 1). The following displays the order in which the indicator lines must appear in the file:
NAME user-supplied name
ROWS
data line(s)
COLUMNS
data line(s)
RHS
data line(s)
RANGES (optional)
data line(s)
BOUNDS (optional)
data line(s)
ENDATA
The ‘user-supplied name’ specifies a name for the problem and must occupy columns 1522. The name can either be blank or up to a maximum of 8 characters.
A data line follows the same fixed format made up of fields defined below. The contents of the fields may have different significance depending upon the section of data in which they appear.
Field 1 Field 2 Field 3 Field 4 Field 5 Field 6
Columns 23 512 1522 2536 4047 5061
Contents Code Name Name Value Name Value
The names and codes consist of ‘alphanumeric’ characters (i.e., a–z, A–Z, 09, +, -, *, blank ( ), :, $ or full stop (.) only) and the names must not contain leading blank characters. Values are read using Fortran format E12.0. This allows values to be entered in several equivalent forms. For example, 1.2345678, 1.2345678E+0, 123.45678E−2 and 12345678E−07 all represent the same number. It is safest to include an explicit decimal point.
Note that in order to ensure numeric values are interpreted as intended, they should be right-justified in the 12-character field, with no trailing blanks. This is because in some situations trailing blanks may be interpreted as zeros and this can dramatically affect the interpretation of the value. This is relevant if the value contains an exponent, or if it contains neither an exponent nor an explicit decimal point. For example, the fields
%%%%1.23E-2%
%%%%%%%123%%
may be interpreted as 1.23E−20 and 12300 respectively (where % denotes a blank). The actual behaviour is system-dependent.
Comment lines are allowed in the data file. These must have an asterisk (*) in column 1 and any characters in columns 2–80. In any data line, a dollar sign ($) as the first character in Field 3 or 5 indicates that the information from that point through column 80 consists of comments.
Columns outside the six fields must be blank, except for columns 72–80, whose contents are ignored by the routine. These columns may be used to enter a sequence number. A non-blank character outside the predefined six fields and columns 72–80 is considered to be a major error (ifail=13; see Section 6), unless it is part of a comment.
ROWS Data Lines
These lines specify row (constraint) names and their inequality types (i.e., =, or ).
Field 1: defines the constraint type. It may be in column 2 or column 3.
N free row, that is no constraint. It may be used to define the objective row.
G greater than or equal to (i.e., ).
L less than or equal to (i.e., ).
E exactly equal to (i.e., =).
Field 2: defines the row name.
Row type N stands for ‘Not binding’, also known as ‘Free’. It can be used to define the objective row. The objective row is a free row that specifies the vector c in the linear objective term cTx. It is taken to be the first free row, unless some other free row name is specified by the names array (see Section 5). Note that c is assumed to be zero if (for example) the line
%N%%DUMMYROW
(where % denotes a blank) appears in the ROWS section of the MPSX data file, and the row name DUMMYROW is omitted from the COLUMNS section.
COLUMNS Data Lines
These lines specify the names to be assigned to the variables (columns) in the general linear constraint matrix A, and define, in terms of column vectors, the actual values of the corresponding matrix elements.
Field 1: blank (ignored).
Field 2: gives the name of the column associated with the elements specified in the following fields.
Field 3: contains the name of a row.
Field 4: used in conjunction with Field 3 contains the value of the matrix element.
Field 5: is optional (may be used like Field 3).
Field 6: is optional (may be used like Field 4).
Note that only the nonzero elements of A and c need to be specified in the COLUMNS section, as any zero elements of A are removed and any unspecified elements of c are assumed to be zero. In addition, any nonzero elements in the jth column of A must be grouped together before those in the (j+1)th column, for j=1,2,,n-1. Nonzero elements within a column may however appear in any order.
RHS Data Lines
This section specifies the right-hand side values of the general linear constraint matrix A (if any). The lines specify the name to be given to the right-hand side (RHS) vector along with the numerical values of the elements of the vector, which may appear in any order. The data lines have exactly the same format as the COLUMNS data lines, except that the column name is replaced by the RHS name. Only the nonzero elements need be specified. Note that this section may be empty, in which case the RHS vector is assumed to be zero.
RANGES Data Lines (optional)
Ranges are used for constraints of the form lAxu, where both l and u are finite. The range of the constraint is r=u-l. Either l or u must be specified in the RHS section and r must be defined in this section. The data lines have exactly the same format as the COLUMNS data lines, except that the column name is replaced by the RANGES name.
BOUNDS Data Lines (optional)
These lines specify limits on the values of the variables (l and u in lxu). If the variable is not specified in the bound set then it is automatically assumed to lie between default lower and upper bounds (usually 0 and +). Like an RHS column which is given a name, the set of variables in one bound set is also given a name.
Field 1: specifies the type of bound or defines the variable type.
LO lower bound
UP upper bound
FX fixed variable
FR free variable (- to +)
MI lower bound is -
PL upper bound is +. This is the default variable type.
Field 2: identifies a name for the bound set.
Field 3: identifies the column name of the variable belonging to this set.
Field 4: identifies the value of the bound; this has a numerical value only in association with LO, UP, FX in Field 1, otherwise it is blank.
Field 5: is blank and ignored.
Field 6: is blank and ignored.
Note that if RANGES and BOUNDS sections are both present, the RANGES section must appear first.

4 References

IBM (1971) MPSX – Mathematical programming system Program Number 5734 XM4 IBM Trade Corporation, New York

5 Arguments

1: infile Integer Input
On entry: the unit number associated with the MPSX data file.
Constraint: 0infile2147483647.
2: maxn Integer Input
On entry: an upper limit for the number of variables in the problem.
Constraint: maxn1.
3: maxm Integer Input
On entry: an upper limit for the number of constraints (including the objective row) in the problem.
Constraint: maxm1.
4: maxnnz Integer Input
On entry: an upper limit for the number of nonzeros (including the objective row) in the problem.
Constraint: maxnnz1.
5: xbldef Real (Kind=nag_wp) Input
On entry: the default lower bound to be used for the variables in the problem when none is specified in the BOUNDS section of the MPSX data file. For a standard LP or QP problem xbldef would normally be set to zero.
6: xbudef Real (Kind=nag_wp) Input
On entry: the default upper bound to be used for the variables in the problem when none is specified in the BOUNDS section of the MPSX data file. For a standard LP or QP problem xbudef would normally be set to ‘infinity’ (i.e., xbudef1020).
Constraint: xbudefxbldef.
7: mpslst Logical Input
On entry: if mpslst=.TRUE., a listing of the input data is sent to the current advisory message unit (as defined by x04abf). This can be useful for debugging the MPSX data file. If mpslst=.FALSE., no listing is produced.
8: n Integer Output
On exit: n, the actual number of variables in the problem.
9: m Integer Output
On exit: m, the actual number of general linear constraints in the problem (including the objective row).
10: nnz Integer Output
On exit: the actual number of nonzeros in the problem (including the objective row).
11: iobj Integer Output
On exit: if iobj>0, row iobj of A is a free row containing the nonzero coefficients of the vector c.
If iobj=0, the coefficients of c are assumed to be zero.
If iobj=−1, no such row was found and the routine terminates with ifail=4 or 5 (see Section 6).
12: ncolh Integer Output
On exit: ncolh=0. For QP problems, ncolh is the number of leading nonzero columns of the Hessian matrix H and must, therefore, be set >0 before calling e04nkf.
13: a(maxnnz) Real (Kind=nag_wp) array Output
On exit: the nonzero elements of A, ordered by increasing column index.
14: ha(maxnnz) Integer array Output
On exit: the row indices of the nonzero elements stored in a.
15: ka(maxn+1) Integer array Output
On exit: a set of pointers to the beginning of each column of A. More precisely, ka(i) contains the index in a of the start of the ith column, for i=1,2,,n. Note that ka(1)=1 and ka(n+1)=nnz+1.
16: bl(maxn+maxm) Real (Kind=nag_wp) array Output
17: bu(maxn+maxm) Real (Kind=nag_wp) array Output
On exit: bl contains the vector l (the lower bounds) and bu contains the vector u (the upper bounds), for all the variables and constraints in the following order. The first n elements of each array contain the bounds on the variables x and the next m elements contain the bounds for the linear objective term cTx and the general linear constraints Ax (if any). Note that an ‘infinite’ lower bound is indicated by bl(j) = -1.0E+20, an ‘infinite’ upper bound by bu(j)=-1.0E+20 and an equality constraint by bl(j)=bu(j). (The lower bound for cTx, stored in bl(n+iobj), is set to -xbudef. The corresponding upper bound, stored in bu(n+iobj), is set to xbudef.)
Note that e04mzf uses an ‘infinite’ bound size of 1020 in the definition of l and u. In other words, any element of u greater than or equal to 1020 will be regarded as + (and similarly any element of l less than or equal to -1020 will be regarded as -). If this value is deemed to be ‘inappropriate’, you are recommended to reset the value of the optional parameter Infinite Bound Size and make any necessary changes to bl and/or bu before calling e04nkf.
18: start Character(1) Output
On exit: start='C' and an internal Crash procedure will be used by e04nkf to choose an initial basis.
19: names(5) Character(8) array Input/Output
On entry: a set of names associated with the MPSX form of the problem.
names(1)
Must contain either the name of the problem or be blank.
names(2)
Must contain either the name of the objective row or be blank (in which case the first objective free row is used).
names(3)
Must contain either the name of the RHS set to be used or be blank (in which case the first RHS set is used).
names(4)
Must contain either the name of the RANGE set to be used or be blank (in which case the first RANGE set (if any) is used).
names(5)
Must contain either the name of the BOUNDS set to be used or be blank (in which case the first BOUNDS set (if any) is used).
On exit: a set of names associated with the problem as defined in the MPSX data file as follows:
names(1)
Contains the name of the problem (or blank if none).
names(2)
Contains the name of the objective row (or blank if none).
names(3)
Contains the name of the RHS set (or blank if none).
names(4)
Contains the name of the RANGE set (or blank if none).
names(5)
Contains the name of the BOUNDS set (or blank if none).
20: nname Integer Output
On exit: n+m, the total number of variables and constraints in the problem.
21: crname(maxn+maxm) Character(8) array Output
On exit: the MPSX names of all the variables and constraints in the problem in the following order. The first n elements contain the MPSX names for the variables and the next m elements contain the MPSX names for the objective row and general linear constraints (if any). Note that the MPSX name for the objective row is stored in crname(n+iobj).
22: xs(maxn+maxm) Real (Kind=nag_wp) array Output
On exit: a set of initial values for the variables and constraints in the problem. More precisely,
xs(j)=min(max(0.0,bl(j)),bu(j)), for j=1,2,,nname.
23: istate(maxn+maxm) Integer array Output
On exit: a set of initial states for the variables and constraints in the problem. More precisely, istate(j)=1 if xs(j)=bu(j) and 0 otherwise, for j=1,2,,nname.
24: ifail Integer Input/Output
On entry: ifail must be set to 0, −1 or 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 or 1 is recommended. If message printing is undesirable, then the value 1 is recommended. Otherwise, the value 0 is recommended. 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 or −1, explanatory error messages are output on the current error message unit (as defined by x04aaf).
Errors or warnings detected by the routine:
Note: the error exits ifail=4 to ifail=16 are caused by having either a corrupt or a nonstandard MPSX data file. Refer to Section 3 for a detailed description of the MPSX format which can be read by e04mzf. If mpslst=.TRUE., the last line of printed output refers to the line in the MPSX data file which contains the reported error.
ifail=1
Too many rows. Limit is value, but the actual number required is value.
Increase maxm by at least (m-maxm) and rerun e04mzf.
ifail=2
Too many columns. Limit is value, but the actual number required is value.
Increase maxn by at least (n-maxn) and rerun e04mzf.
ifail=3
Too many nonzeros. Limit is value, but the actual number required is value.
Increase maxnnz by at least (nnz-maxnnz) and rerun e04mzf.
ifail=4
No objective row found.
There must be at least one row in the ROWS section with row type N for the objective row.
ifail=5
Objective row name value is not defined in the ROWS section.
ifail=6
No rows specified.
ifail=7
Illegal constraint type value.
An illegal constraint type was detected in the ROWS section. The constraint type must be either N, L, G or E.
ifail=8
Row name with leading blank or non-alphanumeric character: value.
ifail=9
Column name with leading blank or non-alphanumeric character: value.
ifail=10
Illegal bound type value.
An illegal bound type was detected in the BOUNDS section. The bound type must be either LO, UP, FX, FR, MI or PL.
ifail=11
Column name value is not defined in the COLUMNS section.
An unknown column name was detected in the BOUNDS section. All the column names must be specified in the COLUMNS section.
ifail=12
The last line must be the ENDATA indicator line.
ifail=13
Line value is not a comment nor a valid line.
ifail=14
Row name value is not defined in the ROWS section.
An unknown row name was detected in COLUMNS, RHS or RANGES section. All the row names must be specified in the ROWS section.
ifail=15
No columns specified.
ifail=16
BOUNDS name value was not found.
RANGES name value was not found.
RHS name value was not found.
ifail=17
On entry, infile=value.
Constraint: 0infile2147483647.
On entry, maxm=value.
Constraint: maxm1.
On entry, maxn=value.
Constraint: maxn1.
On entry, maxnnz=value.
Constraint: maxnnz1.
On entry, xbldef=value and xbudef=value.
Constraint: xbldefxbudef.
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.
e04mzf is not threaded in any implementation.

9 Further Comments

None.

10 Example

This example solves the quadratic programming problem
minimize cT x + 12 xT H x   subject to   l Ax u, −2 Ax 2,  
where
c= ( -4.0 -1.0 -1.0 -1.0 -1.0 -1.0 -1.0 -0.1 -0.3 ) ,   H= ( 2 1 1 1 1 0 0 0 0 1 2 1 1 1 0 0 0 0 1 1 2 1 1 0 0 0 0 1 1 1 2 1 0 0 0 0 1 1 1 1 2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ) ,   ncolh=5,  
A= ( 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 4.0 1.0 2.0 3.0 4.0 -2.0 1.0 1.0 1.0 1.0 1.0 -1.0 1.0 -1.0 1.0 1.0 1.0 1.0 1.0 ) ,  
l= ( -2.0 -2.0 -2.0 )   and   u= ( 1.5 1.5 4.0 ) .  
The optimal solution (to five figures) is
x*=(2.0,-0.23333,-0.26667,-0.3,-0.1,2.0,2.0,-1.7777,-0.45555)T.  
Three bound constraints and two general linear constraints are active at the solution. Note that, although the Hessian matrix is positive semidefinite, the point x* is unique.
The MPSX representation of the problem is given in Section 10.2.

10.1 Program Text

Program Text (e04mzfe.f90)

10.2 Program Data

Note: the MPSX data which is read by e04mzf begins with the second record of this data file; the first record is a caption which is read by the example program.
Program Options (e04mzfe.opt)

10.3 Program Results

Program Results (e04mzfe.r)