naginterfaces.library.opt.handle_​set_​linconstr

naginterfaces.library.opt.handle_set_linconstr(handle, bl=None, bu=None, irowb=None, icolb=None, b=None, idlc=0)

handle_set_linconstr is a part of the NAG optimization modelling suite and adds a new block of linear constraints to the problem or modifies an individual linear constraint.

For full information please refer to the NAG Library document for e04rj

https://support.nag.com/numeric/nl/nagdoc_30/flhtml/e04/e04rjf.html

Parameters

handleHandle

The handle to the problem. It needs to be initialized (e.g., by handle_init()) and must not be changed between calls to the NAG optimization modelling suite.

blNone or float, array-like, shape $\left(\text{nclin}\right)$, optional

$bl$ defines the lower bounds of the linear constraints, $l_B$

buNone or float, array-like, shape $\left(\text{nclin}\right)$, optional

$bu$ defines the upper bounds of the linear constraints, $u_B$

irowbNone or int, array-like, shape $\left(\text{nnzb}\right)$, optional

$irowb$ specifies one-based row indices of the $\text{nnzb}$ nonzeros of the sparse matrix $B$ in coordinate storage (CS) format.

icolbNone or int, array-like, shape $\left(\text{nnzb}\right)$, optional

$icolb$ specifies one-based column indices of the $\text{nnzb}$ nonzeros of the sparse matrix $B$ in coordinate storage (CS) format.

bNone or float, array-like, shape $\left(\text{nnzb}\right)$, optional

$b$ stores the $\text{nnzb}$ nonzeros of the sparse matrix $B$ in coordinate storage (CS) format.

idlcint, optional

If $idlc=0$, a new block of linear constraints is added to the model; otherwise, $idlc>0$ refers to the number of an existing linear constraint which will be replaced and $\text{nclin}$ must be set to one.

Returns

idlcint

If $idlc=0$, the number of the last linear constraint added. By definition, it is the number of linear constraints already defined plus $\text{nclin}$. Otherwise, $idlc>0$ stays unchanged.

Raises

NagValueError

(errno $1$)

$handle$ has not been initialized.

(errno $1$)

$handle$ does not belong to the NAG optimization modelling suite, has not been initialized properly or is corrupted.

(errno $1$)

$handle$ has not been initialized properly or is corrupted.

(errno $2$)

The problem cannot be modified right now, the solver is running.

(errno $4$)

On entry, $idlc=⟨value⟩$.

Constraint: $idlc\ge 0$.

(errno $4$)

On entry, $idlc=⟨value⟩$.

The given $idlc$ does not match with any existing linear constraint.

The maximum $idlc$ is $⟨value⟩$.

(errno $6$)

On entry, $\text{nnzb}=⟨value⟩$.

Constraint: $\text{nnzb}\ge 0$.

(errno $6$)

On entry, $\text{nclin}=⟨value⟩$.

Constraint: $\text{nclin}\ge 0$.

(errno $6$)

On entry, $idlc=⟨value⟩$ and $\text{nclin}=⟨value⟩$.

Constraint: If $idlc>0$, $\text{nclin}=1$.

(errno $8$)

On entry, $i=⟨value⟩$, $irowb[\text{i}-1]=⟨value⟩$ and $\text{nclin}=⟨value⟩$.

Constraint: $1\le irowb[\text{i}-1]\le \text{nclin}$.

(errno $8$)

On entry, $i=⟨value⟩$, $icolb[\text{i}-1]=⟨value⟩$ and $n=⟨value⟩$.

Constraint: $1\le icolb[\text{i}-1]\le n$.

(errno $8$)

On entry, more than one element of $b$ has row index $⟨value⟩$ and column index $⟨value⟩$.

Constraint: each element of $b$ must have a unique row and column index.

(errno $10$)

On entry, $j=⟨value⟩$, $bl[j-1]=⟨value⟩$ and $bu[j-1]=⟨value⟩$.

Constraint: $bl[j-1]\le bu[j-1]$.

(errno $10$)

On entry, $j=⟨value⟩$, $bl[j-1]=⟨value⟩$, $\text{bigbnd}=⟨value⟩$.

Constraint: $bl[j-1]<\text{bigbnd}$.

(errno $10$)

On entry, $j=⟨value⟩$, $bu[j-1]=⟨value⟩$, $\text{bigbnd}=⟨value⟩$.

Constraint: $bu[j-1]>-\text{bigbnd}$.

Notes

After the $handle$ has been initialized (e.g., handle_init() has been called), handle_set_linconstr may be used to add to the problem a new block of $m_B$ linear constraints

$$l_B\le Bx\le u_B$$

where $B$ is a general $m_B\times n$ rectangular matrix, $n$ is the current number of decision variables in the model and $l_B$ and $u_B$ are $m_B$-dimensional vectors defining the lower and upper bounds, respectively. The call can be repeated to add multiple blocks to the model.

Note that the bounds are specified for all the constraints of this block. This form allows full generality in specifying various types of constraint. In particular, the $j$th constraint may be defined as an equality by setting $l_j=u_j$. If certain bounds are not present, the associated elements of $l_B$ or $u_B$ may be set to special values that are treated as $-\infty$ or $+\infty$. See the description of the option ‘Infinite Bound Size’ which is common among all solvers in the suite. Its value is denoted as $\text{bigbnd}$ further in this text. Note that the bounds are interpreted based on its value at the time of calling this function and any later alterations to ‘Infinite Bound Size’ will not affect these constraints.

The linear constraints can be edited. To identify the individual constraints, they are numbered starting with $1$, see $idlc$. A single constraint (i.e., a single row of the matrix $B$) can be modified (replaced) by handle_set_linconstr by referring to its $idlc$. An individual coefficient $b_{\text{ij}}$ of the matrix $B$ can be set or modified by handle_set_linconstr_coeff() and bounds of a single constraint can be set or modified by handle_set_bound(). Note that it is also possible to temporarily disable and enable individual constraints in the model by handle_disable() and handle_enable(), respectively.

Linear constraints may be present in many different types of problems, for simplicity of the notation, only one block of linear constraints is presented. For example,

Linear Programming (LP)

$$\begin{array}{c}\begin{array}{ccc}{minimize}_{x\in R^n}& c^{T}x& \text{(a)}\\ \text{subject to}& l_B\le Bx\le u_B\text{,}& \text{(b)}\\ & l_x\le x\le u_x\text{,}& \text{(c)}\end{array}\end{array}$$

Quadratic Programming (QP)

$$\begin{array}{c}\begin{array}{ccc}{minimize}_{x\in R^n}& \frac{1}{2}{x}^{T}Hx+c^{T}x& \text{(a)}\\ \text{subject to}& l_B\le Bx\le u_B\text{,}& \text{(b)}\\ & l_x\le x\le u_x\text{,}& \text{(c)}\end{array}\end{array}$$

Quadratically Constrained Quadratic Programming (QCQP)

$$\begin{array}{c}\begin{array}{ccc}{minimize}_{x\in R^n}& \frac{1}{2}{x}^{T}Hx+c^{T}x& \text{(a)}\\ \text{subject to}& \frac{1}{2}{x}^T{Q}_{k}x+r_k^{T}x+s_k\le 0\text{,}k=1,\dots ,m_Q\text{,}& \text{(b)}\\ & l_B\le Bx\le u_B\text{,}& \text{(c)}\\ & l_x\le x\le u_x\text{,}& \text{(d)}\end{array}\end{array}$$

Nonlinear Programming (NLP)

$$\begin{array}{c}\begin{array}{ccc}{minimize}_{x\in R^n}& f\left(x\right)& \text{(a)}\\ \text{subject to}& l_g\le g\left(x\right)\le u_g\text{,}& \text{(b)}\\ & \frac{1}{2}{x}^T{Q}_{k}x+r_k^{T}x+s_k\le 0\text{,}k=1,\dots ,m_Q\text{,}& \text{(c)}\\ & l_B\le Bx\le u_B\text{,}& \text{(d)}\\ & l_x\le x\le u_x\text{,}& \text{(e)}\end{array}\end{array}$$

or linear Semidefinite Programming (SDP)

$$\begin{array}{c}\begin{array}{ccc}{minimize}_{x\in R^n}& c^{T}x& \text{(a)}\\ \text{subject to}& {\sum }_1^n{x}_i{A}_i^k-A_0^k⪰0\text{,}k=1,\dots ,m_A\text{,}& \text{(b)}\\ & l_B\le Bx\le u_B\text{,}& \text{(c)}\\ & l_x\le x\le u_x\text{.}& \text{(d)}\end{array}\end{array}$$

See the E04 Introduction for more details about the NAG optimization modelling suite.

See also

naginterfaces.library.examples.mip.handle_solve_milp_ex.main()

naginterfaces.library.examples.opt.handle_add_vars_ex.main()

naginterfaces.library.examples.opt.handle_solve_ipopt_ex.main()

naginterfaces.library.examples.opt.handle_solve_lp_ipm_ex.main()

naginterfaces.library.examples.opt.handle_solve_nldf_ex.main()

naginterfaces.library.examples.opt.handle_solve_socp_ipm_ex.main()

naginterfaces.library.examples.opt.handle_solve_ssqp_ex.main()