Interior-point method


Interior-point methods are a certain class of algorithms that solve linear and nonlinear convex optimization problems.
John von Neumann suggested an interior-point method of linear programming, which was neither a polynomial-time method nor an efficient method in practice. In fact, it turned out to be slower than the commonly used simplex method.
An interior point method, discovered by Soviet mathematician I. I. Dikin in 1967 and reinvented in the U.S. in the mid-1980s.
In 1984, Narendra Karmarkar developed a method for linear programming called Karmarkar's algorithm, which runs in provably polynomial time and is also very efficient in practice. It enabled solutions of linear programming problems that were beyond the capabilities of the simplex method. Contrary to the simplex method, it reaches a best solution by traversing the interior of the feasible region. The method can be generalized to convex programming based on a self-concordant barrier function used to encode the convex set.
Any convex optimization problem can be transformed into minimizing a linear function over a convex set by converting to the epigraph form. The idea of encoding the feasible set using a barrier and designing barrier methods was studied by Anthony V.
Yurii Nesterov and Arkadi Nemirovski came up with a special class of such barriers that can be used to encode any convex set. They guarantee that the number of iterations of the algorithm is bounded by a polynomial in the dimension and accuracy of the solution.
Karmarkar's breakthrough revitalized the study of interior-point methods and barrier problems, showing that it was possible to create an algorithm for linear programming characterized by polynomial complexity and, moreover, that was competitive with the simplex method.
Already Khachiyan's ellipsoid method was a polynomial-time algorithm; however, it was too slow to be of practical interest.
The class of primal-dual path-following interior-point methods is considered the most successful.
Mehrotra's predictor–corrector algorithm provides the basis for most implementations of this class of methods.

Primal-dual interior-point method for nonlinear optimization

The primal-dual method's idea is easy to demonstrate for constrained nonlinear optimization.
For simplicity, consider the all-inequality version of a nonlinear optimization problem:
The logarithmic barrier function associated with is
Here is a small positive scalar, sometimes called the "barrier parameter". As converges to zero the minimum of should converge to a solution of.
The barrier function gradient is
where is the gradient of the original function, and is the gradient of.
In addition to the original variable we introduce a Lagrange multiplier inspired variable
is sometimes called the "perturbed complementarity" condition, for its resemblance to "complementary slackness" in KKT conditions.
We try to find those for which the gradient of the barrier function is zero.
Applying to, we get an equation for the gradient:
where the matrix is the Jacobian of the constraints.
The intuition behind is that the gradient of should lie in the subspace spanned by the constraints' gradients. The "perturbed complementarity" with small can be understood as the condition that the solution should either lie near the boundary, or that the projection of the gradient on the constraint component normal should be almost zero.
Applying Newton's method to and, we get an equation for update :
where is the Hessian matrix of, is a diagonal matrix of, and is a diagonal matrix with.
Because of, the condition
should be enforced at each step. This can be done by choosing appropriate :