# project: Project a vector on the border of the region defined by a set... In lintools: Manipulation of Linear Systems of (in)Equalities

## Description

Compute a vector, closest to x in the Euclidean sense, satisfying a set of linear (in)equality restrictions.

## Usage

 1 2 project(x, A, b, neq = length(b), w = rep(1, length(x)), eps = 0.01, maxiter = 1000L)

## Arguments

 x [numeric] Vector that needs to satisfy the linear restrictions. A [matrix] Coefficient matrix for linear restrictions. b [numeric] Right hand side of linear restrictions. neq [numeric] The first neq rows in A and b are treated as linear equalities. The others as Linear inequalities of the form Ax<=b. w [numeric] Optional weight vector of the same length as x. Must be positive. eps The maximum allowed deviation from the constraints (see details). maxiter maximum number of iterations

## Value

A list with the following entries:

• x: the adjusted vector

• status: Exit status:

• 0: success

• 1: could not allocate enough memory (space for approximately 2(m+n) doubles is necessary).

• 2: divergence detected (set of restrictions may be contradictory)

• 3: maximum number of iterations reached

• eps: The tolerance achieved after optimizing (see Details).

• iterations: The number of iterations performed.

• duration: the time it took to compute the adjusted vector

• objective: The (weighted) Euclidean distance between the initial and the adjusted vector

## Details

The tolerance eps is defined as the maximum absolute value of the difference vector \boldsymbol{Ax}-\boldsymbol{b} for equalities. For inequalities, the difference vector is set to zero when it's value is lesser than zero (i.e. when the restriction is satisfied). The algorithm iterates until either the tolerance is met, the number of allowed iterations is exceeded or divergence is detected.