In constrained
optimization
Mathematical optimization (alternatively spelled ''optimisation'') or mathematical programming is the selection of a best element, with regard to some criteria, from some set of available alternatives. It is generally divided into two subfiel ...
, a field of
mathematics
Mathematics is a field of study that discovers and organizes methods, Mathematical theory, theories and theorems that are developed and Mathematical proof, proved for the needs of empirical sciences and mathematics itself. There are many ar ...
, a barrier function is a
continuous function
In mathematics, a continuous function is a function such that a small variation of the argument induces a small variation of the value of the function. This implies there are no abrupt changes in value, known as '' discontinuities''. More preci ...
whose value increases to infinity as its argument approaches the boundary of the
feasible region
In mathematical optimization and computer science, a feasible region, feasible set, or solution space is the set of all possible points (sets of values of the choice variables) of an optimization problem that satisfy the problem's constraints, ...
of an optimization problem. Such functions are used to replace inequality
constraints by a penalizing term in the objective function that is easier to handle. A barrier function is also called an interior penalty function, as it is a penalty function that forces the solution to remain within the interior of the feasible region.
The two most common types of barrier functions are
inverse barrier function
Inverse or invert may refer to:
Science and mathematics
* Inverse (logic), a type of conditional sentence which is an immediate inference made from another conditional sentence
* Additive inverse, the inverse of a number that, when added to the ...
s and logarithmic barrier functions. Resumption of interest in logarithmic barrier functions was motivated by their connection with primal-dual
interior point method
Interior-point methods (also referred to as barrier methods or IPMs) are algorithms for solving Linear programming, linear and nonlinear programming, non-linear convex optimization problems. IPMs combine two advantages of previously-known algorit ...
s.
Motivation
Consider the following constrained optimization problem:
:minimize
:subject to
where is some constant. If one wishes to remove the inequality constraint, the problem can be reformulated as
:minimize ,
:where if , and zero otherwise.
This problem is equivalent to the first. It gets rid of the inequality, but introduces the issue that the penalty function , and therefore the objective function , is
discontinuous, preventing the use of
calculus
Calculus is the mathematics, mathematical study of continuous change, in the same way that geometry is the study of shape, and algebra is the study of generalizations of arithmetic operations.
Originally called infinitesimal calculus or "the ...
to solve it.
A barrier function, now, is a continuous approximation to that tends to infinity as approaches from above. Using such a function, a new optimization problem is formulated, viz.
:minimize
where is a free parameter. This problem is not equivalent to the original, but as approaches zero, it becomes an ever-better approximation.
Logarithmic barrier function
For logarithmic barrier functions,
is defined as
when
and
otherwise (in one dimension; see below for a definition in higher dimensions). This essentially relies on the fact that
tends to negative infinity as
tends to 0.
This introduces a gradient to the function being optimized which favors less extreme values of
(in this case values lower than
), while having relatively low impact on the function away from these extremes.
Logarithmic barrier functions may be favored over less computationally expensive
inverse barrier functions
Inverse or invert may refer to:
Science and mathematics
* Inverse (logic), a type of conditional sentence which is an immediate inference made from another conditional sentence
* Additive inverse, the inverse of a number that, when added to the ...
depending on the function being optimized.
Higher dimensions
Extending to higher dimensions is simple, provided each dimension is independent. For each variable
which should be limited to be strictly lower than
, add
.
Formal definition
Minimize
subject to
Assume strictly feasible:
Define logarithmic barrier