Journal of Applied Mathematics

Volume 2011 (2011), Article ID 340192, 12 pages

http://dx.doi.org/10.1155/2011/340192

## An Improved Predictor-Corrector Interior-Point Algorithm for Linear Complementarity Problems with -Iteration Complexity

^{1}School of Economics and Management, Wuhan University, Wuhan 430072, China^{2}School of Economics, Wuhan University of Technology, Wuhan 430070, China

Received 18 September 2011; Revised 23 November 2011; Accepted 24 November 2011

Academic Editor: Chong Lin

Copyright © 2011 Debin Fang and Qian Yu. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

This paper proposes an improved predictor-corrector interior-point algorithm for the linear complementarity problem (LCP) based on the Mizuno-Todd-Ye algorithm. The modified corrector steps in our algorithm cannot only draw the iteration point back to a narrower neighborhood of the center path but also reduce the duality gap. It implies that the improved algorithm can converge faster than the MTY algorithm. The iteration complexity of the improved algorithm is proved to obtain which is similar to the classical Mizuno-Todd-Ye algorithm. Finally, the numerical experiments show that our algorithm improved the performance of the classical MTY algorithm.

#### 1. Introduction

Since Karmarkar published the first paper on interior point method [1] in 1984, the interior point methodologies have yielded rich theories and algorithms in the fields of linear programming (LP), quadratic programming (QP), and linear complementarity problems (LCP). Among these interior point methods, predictor-corrector interior-point methods play a special role due to their best polynomial complexity and superlinear convergence.

In 1993, Mizuno et al. [2] proposed the classical representative of predictor-corrector method for linear programming. The Mizuno-Todd-Ye (MTY) algorithm has -iteration complexity which is the best iteration complexity obtained so far for all the interior-point method [3, 4]. Moreover, Ye et al. [5] proved the duality gap of classical MTY algorithm converges to zero quadratically, which dedicated that MTY algorithm has superlinear convergence. The classical MTY algorithm is the first algorithm for LP has both polynomial complexity and superlinear convergence. So the classical MTY algorithm therefore was considered as the most efficient interior point methods for LP.

Afterwards, the MTY algorithm was extended to LCP [6, 7] for its excellent performance. In recent papers, Potra [8, 9] presented several predictor-corrector methods for linear complementarity problems acting in a wide neighborhood of the central path, and Stoer et al. [10] proposed a predictor-corrector algorithm with arbitrarily high order of convergence to degenerate sufficient linear complementarity problems.

Subsequently, Potra proposed a generalization of the MTY algorithm for infeasible starting points for monotone LCP [11]. Based on the so-called “fast step-safe step” strategy, Wright proposed an infeasible interior point method with polynomial and quadratic convergence for nondegenerate LCP [12]. And Ai and Zhang presented an -iteration primal-dual path-following method for monotone LCP based on wide neighborhoods and large updates [13]. These workings improved the MTY algorithm in various aspects. In this paper, the MTY algorithm for LCP will be improved by modifying the corrector of the algorithm to obtain a larger iteration reduce factor than the original, which can guarantee a faster convergence.

The typical iterations of the MTY algorithm operate between two neighborhoods of the central path, and . Before starting the MTY algorithm, an arbitrary initial point will be given. And then, the predictor produces a point by carefully choosing a steplength along the affine-scaling direction at and reduces the primal dual gap by a factor of at least where . In the corrector steps, the corrector produces a point by taking a unit steplength along the centering direction at . It is shown that the corrector put the iteration points back to which is the narrower neighborhood of the central path and maintains the same duality gap. Thus predictor-corrector reduces the duality gap by a factor of at least and holds the iteration points in the neighborhood . It follows that the corresponding iterative procedure has -iteration complexity.

The modified algorithm in this paper has only one corrector step in neighborhood of the central path and one predictor step in a narrow neighborhood same as the MTY algorithm. Moreover, the modified iteration direction of corrector step can make a reduction for duality gap. This improvement attained a larger iteration reduce factor and results in a faster convergence than the classical MTY algorithm.

Section 2 describes the procedure of our predictor-corrector algorithm. Section 3 gives the convergence analysis of our algorithm. It shows that the improved algorithm can generate a sequence in the neighborhood from an arbitrary initial point and converge to optimal solution with -iteration complexity. Finally, the performances of our algorithm are evaluated by numerical experiments in Section 4.

#### 2. Description of the Algorithm

In this paper, LCP is to find a vector pair such that , and , where and is a positive semidefinite matrix.

Denote the set of all feasible points of LCP by and the solution set of LCP by . The relative interior of , is called the set of strictly feasible points, or the set of interior points. We assume in this paper. It is known that if is nonempty, then for any parameter the nonlinear system has a unique positive solution [14]. The set of all such solutions defines the central path of LCP. So the assumption is sufficient to establish the existence of the central path and the existence of a solution to LCP.

The mean value of is denoted as throughout this paper. Then is considered as the neighborhood of central path, where , denotes the vector of ones, and express the Euclidean norm.

The improved predictor-corrector interior-point algorithm (IPCIP) is as follows.

##### 2.1. Main Steps of IPCIP Algorithm

*Step 0. *Choose an initial pair of interior point with , and set the accuracy parameter . Let .

*Step 1. *If , then stop.

*Step 2. *Compute the predictor direction by solving the system (2.1),

*Step 3. *Set where . (We will show later that , where holds for every .)

*Step 4. *Set , and compute the corrector direction by solving the system (2.2),

*Step 5. *Set , update , and return to Step 1.

#### 3. Convergence and Complexity Analysis

Lemma 3.1. *If and , then*(i)*; *(ii)*; *(iii)*.*

*Proof. *The proof of (i) can be seen in Lemma 1 of [2].

From and , we can obtain and . Thus inequalities (ii) and (iii) hold. This completes the proof.

Lemma 3.2. *If the point was generated by the algorithm, then one has , where .*

*Proof. *From Step 2 of our algorithm and is a positive semidefinite matrix, we have . Hence

This completes the proof.

Lemma 3.3. *If , then the point generated by the algorithm satisfies , , and .*

*Proof. *Let us define , where

Note that , then we have from and from .

Furthermore holds since . Therefore we have

So for any .

Since is continuous with respect to and , it is easily seen that , .

Note that , the above argument implies .

It’s well known that the inequality holds for arbitrary vector and arbitrary constant .

So we have from the proof above and Lemma 3.2, thus .

This completes the proof.

*Remark 3.4. *So Lemma 3.3 dedicates that the predictors of our algorithm operate in a wide neighborhood of the central path .

Lemma 3.5. *If then *

*Proof. *

This completes the proof.

Lemma 3.6. *If and the point was generated by the algorithm, then always hold.*

*Proof. *Denote then

From Lemmas 3.1 and 3.5 we have

Note that , we have

It is obvious that when , so .

This completes the proof.

*Remark 3.7. *From Lemmas 3.2 and 3.6, we can obtain immediately that . That means the corrector reduces to . Because the corrector in the classical MTY algorithm just maintains the same duality gap, the improvement of the corrector in our algorithm will make a faster reduction of than MTY algorithm.

Then the polynomial convergence of the algorithm could be established.

Theorem 3.8. *Suppose that the sequence generated by the algorithm, then one has .*

*Proof. *Let us define that , and denote
As discussed in Lemma 3.6, we also have .

From Lemmas 3.1 and 3.5 and , then

So we can write , that is, .

Hence

Since the inequality holds for arbitrary vector and arbitrary constant , thus

Let us define , then .

When , we have , this results in and decreases monotonously, when . Because holds for every thus

Since (otherwise the algorithm will be terminated), it follows that

We have proved that , in Lemma 3.3. So we have

Note that and are continuous with respect to . It implies that and hold for any .

Let , then we have and .

Furthermore, when , thus .

Consider that we can obtain directly from the steps in our algorithm, so holds for every generated by the algorithm.

This completes the proof.

Theorem 3.9. *The iteration complexity of the algorithm is .*

*Proof. *Let then

From Lemma 3.1, we have

From , we have , then ,

Hence , and we can obtain .

So it implies that for each . This means that will decrease at least by a constant factor of at each iteration, which guarantees a -iteration complexity for the algorithm.

This completes the proof.

*Remark 3.10. *Note that the reduce factor in the MTY algorithm is where , but the reduce factor in our algorithm is where , which is larger than . It implies that our algorithm will converge faster than the MTY algorithm, although both have the same iteration complexity.

#### 4. Examples and Numerical Results

Finally, the numerical experiments were carried out to evaluate the performance and practical efficiency of the algorithm,

*Test Problem 1*

Find a vector pair such that , where and is a positive semidefinite matrix,

*Test Problem 2*

Find a vector pair such that , where and M_{2} is a positive semidefinite matrix,

*Test Problem 3*

This example is a general test problem used by Noor et al. [15]. The problem is to find a vector pair such that , where and is a positive semidefinite matrix. To test the efficiency of our algorithm by large-scale problem, we set the dimension of test problem 3 as 100, that is, ,

These test problems are solved by our IPCIP algorithm and the classical MTY algorithm. The experiments are run on a PC (2.6 GHz CPU, 2 G DDR RAM) using MATLAB 7. The results including the corresponding numbers of iterations, , the approximate solutions, and computing times are shown in Tables 1, 2, and 3.

#### 5. Conclusion

This paper modified the MTY algorithm for solving monotone LCPs to strengthen the convergence results. Although the iteration complexity of the improved algorithm was proved to be which is similar to the classical MTY algorithm, but the reduce factor of duality gap was enhanced to 1.168 and results in a faster convergence than the classical MTY algorithm.

The numerical results show that our IPCIP algorithm is more efficient than the MTY algorithm. The number of iterations was decreased, and the computing times could be reduced by nearly 20% to 50%. That indicated IPCIP algorithm has a better performance than the MTY algorithm.

#### Acknowledgments

This paper was supported by the Humanities and Social Science Research Projects of the Ministry of Education of China (no. 11YJCZH221) and the National Natural Science Foundation of China (no. 71103135).

#### References

- N. Karmarkar, “A new polynomial-time algorithm for linear programming,”
*Combinatorica*, vol. 4, no. 4, pp. 373–395, 1984. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - S. Mizuno, M. J. Todd, and Y. Ye, “On adaptive-step primal-dual interior-point algorithms for linear programming,”
*Mathematics of Operations Research*, vol. 18, no. 4, pp. 964–981, 1993. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - T. Illés, M. Nagy, and T. Terlaky, “A polynomial path-following interior point algorithm for general linear complementarity problems,”
*Journal of Global Optimization*, vol. 47, no. 3, pp. 329–342, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - Y. Q. Bai, G. Lesaja, and C. Roos, “A new class of polynomial interior-point algorithms for linear complementary problems,”
*Pacific Journal of Optimization. An International Journal*, vol. 4, no. 1, pp. 19–41, 2008. View at Google Scholar - Y. Ye, O. Güler, R. A. Tapia, and Y. Zhang, “A quadratically convergent $O(\sqrt{nL})$-iteration algorithm for linear programming,”
*Mathematical Programming*, vol. 59, no. 2, pp. 151–162, 1993. View at Publisher · View at Google Scholar - Y. Ye and K. Anstreicher, “On quadratic and $O(\sqrt{nL})$ convergence of a predictor-corrector algorithm for LCP,”
*Mathematical Programming*, vol. 62, no. 3, pp. 537–551, 1993. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - J. Ji, F. A. Potra, and S. Huang, “Predictor-corrector method for linear complementarity problems with polynomial complexity and superlinear convergence,”
*Journal of Optimization Theory and Applications*, vol. 85, no. 1, pp. 187–199, 1995. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - F. A. Potra, “The Mizuno-Todd-Ye algorithm in a larger neighborhood of the central path,”
*European Journal of Operational Research*, vol. 143, no. 2, pp. 257–267, 2002. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - F. A. Potra, “A superlinearly convergent predictor-corrector method for degenerate LCP in a wide neighborhood of the central path with $O(\sqrt{nL})$-iteration complexity,”
*Mathematical Programming A*, vol. 100, no. 2, pp. 317–337, 2004. View at Publisher · View at Google Scholar - J. Stoer, M. Wechs, and S. Mizuno, “High order infeasible-interior-point methods for solving sufficient linear complementarity problems,”
*Mathematics of Operations Research*, vol. 23, no. 4, pp. 832–862, 1998. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - F. A. Potra, “An $O(\sqrt{nL})$ infeasible-interior-point algorithm for LCP with quadratic convergence,”
*Annals of Operations Research*, vol. 62, pp. 81–102, 1996. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - S. J. Wright, “A path-following interior-point algorithm for linear and quadratic problems,”
*Annals of Operations Research*, vol. 62, pp. 103–130, 1996. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - W. Ai and S. Zhang, “An $O(\sqrt{nL})$ iteration primal-dual path-following method, based on wide neighborhoods and large updates, for monotone LCP,”
*SIAM Journal on Optimization*, vol. 16, no. 2, pp. 400–417, 2005. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - M. Kojima, N. Megiddo, T. Noma, and A. Yoshise, “A unified approach to interior point algorithms for linear complementarity problems,”
*Lecture Notes in Computer Science*, vol. 538, 1991. View at Google Scholar - M. A. Noor, Y. J. Wang, and N. Xiu, “Projection iterative schemes for general variational inequalities,”
*Journal of Inequalities in Pure and Applied Mathematics*, vol. 3, no. 3, pp. 1–8, 2002. View at Google Scholar · View at Zentralblatt MATH