Menu
Home
People
Places
Arts
History
Plants & Animals
Science
Life & Culture
Technology
Reference.org
Mehrotra predictor–corrector method
open-in-new
<p>Mehrotra's predictor–corrector method in <a href="/facts/Optimization_(mathematics)/oRn8Iv5I">optimization</a> is a specific <a href="/facts/Interior_point_method/prIIrceN">interior point method</a> for <a href="/facts/Linear_programming/GduXFQxT">linear programming</a>. It was proposed in 1989 by Sanjay Mehrotra. </p><p>The method is based on the fact that at each <a href="/facts/Iteration/AMc1p7BL">iteration</a> of an interior point algorithm it is necessary to compute the <a href="/facts/Cholesky_decomposition/Q3IjkuTs">Cholesky decomposition</a> (factorization) of a large matrix to find the search direction. The factorization step is the most computationally expensive step in the algorithm. Therefore, it makes sense to use the same decomposition more than once before recomputing it. </p><p>At each iteration of the algorithm, Mehrotra's predictor–corrector method uses the same Cholesky decomposition to find two different directions: a predictor and a corrector. </p><p>The idea is to first compute an optimizing search direction based on a first order term (predictor). The step size that can be taken in this direction is used to evaluate how much centrality correction is needed. Then, a corrector term is computed: this contains both a centrality term and a second order term. </p><p>The complete search direction is the sum of the predictor direction and the corrector direction. </p><p>Although there is no theoretical complexity bound on it yet, Mehrotra's predictor–corrector method is widely used in practice. Its corrector step uses the same <a href="/facts/Cholesky_decomposition/Q3IjkuTs">Cholesky decomposition</a> found during the predictor step in an effective way, and thus it is only marginally more expensive than a standard interior point algorithm. However, the additional overhead per iteration is usually paid off by a reduction in the number of iterations needed to reach an optimal solution. It also appears to converge very fast when close to the optimum. </p>