- About this Journal ·
- Abstracting and Indexing ·
- Advance Access ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Articles in Press ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents

Abstract and Applied Analysis

Volume 2013 (2013), Article ID 836720, 6 pages

http://dx.doi.org/10.1155/2013/836720

## A Projection-Type Method for Multivalued Variational Inequality

Institute of Applied Mathematics, Chongqing University of Posts and Telecommunications, Chongqing 400065, China

Received 18 April 2013; Revised 5 August 2013; Accepted 19 August 2013

Academic Editor: Shawn X. Wang

Copyright © 2013 Changjie Fang et al. 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

We propose a projection-type method for multivalued variational inequality. The iteration sequence generated by the algorithm is proven to be globally convergent to a solution, provided that the multivalued mapping is continuous with nonempty compact convex values. Moreover, we present a necessary and sufficient condition on the nonemptiness of the solution set. Preliminary computational experience is also reported.

#### 1. Introduction

We consider the following multivalued variational inequality, denoted by to find and such that where is a nonempty closed convex set in , is a multivalued mapping from into with nonempty values, and and denote the inner product and norm in , respectively.

Projection-type algorithms have been extensively studied in the literature; see [1–8] and the references therein. Reference [2] proposes a subgradient extragradient algorithm for solving single-valued variational inequality in which the next iterate is a projection onto a halfspace whose bounding hyperplane supports the feasible set at a certain point. Reference [6] proposes a projection method for variational inequality problems in which the hyperplane strictly separates the current iterate from the solution set of (1). Theory and algorithm of multivalued variational inequality have been much studied in the literature [1, 9–18]. Various algorithms for computing the solution of (1) are proposed. The well-known proximal point algorithm [19] requires the multivalued mapping to be monotone. Relaxing the monotonicity assumption, [12] proposed the double projection algorithm for solving (1); also see [5]. Assume that is pseudomonotone; [20] described a combined relaxation method for solving (1); see also [21]. Recently, [22] proposed an extragradient method for generalized variational inequality. In [22], the next iterate is the projection of the current iterate onto the feasible set ; also see [23].

In this paper, we introduce a projection-type method for multivalued variational inequality in which the next iterate is a projection of the initial iterate onto intersection of two halfspaces containing the solution set. We obtain a global convergence theorem, assuming that is pseudomonotone on with respect to the solution set; see (3) in the following. Moreover, we show that the iterative sequence diverges if and only if the solution set is empty. We also present numerical results of the proposed method. Now let us compare our algorithm with algorithms in [5, 12, 22, 23]. First, the next iterate in our method relates to the initial point. In [5, 12], the next iterate is the projection of the current iterate onto the intersection of the hyperplane and the feasible set . Secondly, the next iterate in our method is a projection of the initial point onto the intersection of the two hyperplanes and the feasible set . In addition, our Armijo-type linesearch procedure is also different from those in [12, 22, 23].

The organization of this paper is as follows. In the next section, we present the algorithm details and some lemmas. We prove several preliminary results for convergence analysis in Section 3. Numerical results are reported in the last section.

#### 2. Algorithms

Let us recall the definition of continuous multivalued mapping. is said to be upper semicontinuous at if, for every open set containing , there is an open set containing such that for all . is said to be lower semicontinuous at if, given any sequence converging to and any , there exists a sequence that converges to . is said to be continuous at if, it is both upper semicontinuous and lower semicontinuous at . If is single valued, then both upper semicontinuity and lower semicontinuity reduce to the continuity of .

is called pseudomonotone on in the sense of Karamardian [24], if, for any ,

Let be the solution set of (1), that is, those points satisfying (1). Throughout this paper, we assume that the solution set of the problem (1) is nonempty and is continuous on with nonempty compact convex values satisfying the following property: The property (3) holds if is pseudomonotone on .

Let denote the projector onto , and let be a parameter.

Proposition 1. * and solves the problem (1) if and only if
*

*Algorithm 2. *Choose and three parameters , and . Set .

*Step 1. *If for some , stop; else take arbitrarily .

*Step 2. *For every positive integer , let .

*Step 3. *Let be the smallest nonnegative integer satisfying
Set and .

*Step 4. * Compute , where
Let , and go to Step 1.

*Remark 3. *Let us compare the previous algorithm with Algorithm 3.1 in [11]. First, the parameter is required to be strictly less than 1, and is assumed to be equal to 1 in their Algorithm 3.1. In our Algorithm 2, the parameter can take any positive scalar and . Secondly, since has compact convex values, has closed convex values. Therefore, in Step 2 of our algorithm is uniquely determined by . Hence, it is easy to compute the value of satisfying (5). In Step 1 of their Algorithm 3.1, since is a multivalued mapping, it is very difficult in practice to compute the value of satisfying and at the same time. In addition, we report numerical results concerning our algorithm, while [11] does not present numerical experiments for the proposed algorithm. Finally, we compare the performance of our algorithm and Algorithm 3.1 in [11] (see Table 4).

Lemma 4. *The sequence generated in Step 2 has the following properties:
*

*Proof. *See Lemma 2.1 in [5].

Lemma 5. *For every and ,
*

*Proof. *See Lemma 2.3 in [5].

We show that Algorithm 2 is well defined and implementable.

Lemma 6. *If , there exists satisfying (5).*

*Proof. *In view of Lemma 4, we have . Therefore,
where the first inequality follows from Lemma 5 and the second one follows from and .

Lemma 7. *Let be a closed convex subset of . For any and , the following statements hold:*(i)*,*(ii)*. *

*Proof. *See [25].

Lemma 8. *Let and . If , then the hyperplane
**
strictly separates and the solution set .*

*Proof. *Since ,
where the first inequality follows from (5) and the last one follows from . Since satisfies the property (3), .

Lemma 9. *If , then .*

*Proof. *It follows from Lemma 8 that . Next, it is sufficient to prove that for all . The proof will be given by induction.

Obviously, . Suppose that
Then, .

Let . Since
it follows from Lemma 7 (i) that
Thus, . Therefore, we obtain that for all .

Lemma 10. *Let be a nonempty bounded closed convex set, and let the mapping be lower semicontinuous with nonempty closed convex values; then, the solution set of is nonempty.*

*Proof. *See Lemma 2.9 in [22].

The following lemma says that if the solution set is empty, then is a nonempty set, which implies the feasibility of Algorithm 2.

Lemma 11. *Let be continuous with nonempty compact convex values on , and suppose that ; then, for all .*

*Proof. *On the contrary, suppose that there exists such that . Then, there exists a positive number such that
where
Since is continuous with compact values, Proposition 3.11 in [26] implies that is a bounded set, and so is bounded. Without loss of generality, we assume that
Consider the variational inequality , where
It follows from Lemma 10 that the solution set of , denoted by , is nonempty. We denote the three sequences , and by , and , respectively, when Algorithm 2 is applied to with starting point . We claim that(i)the set has at least elements;(ii) for ;(iii) is not a solution of .

Items (i) and (iii) are obvious. Next we prove the item (ii). It is sufficient to prove that
where .

Since
where the second inequality follows from and Lemma 7 (ii), so , and hence, . Therefore,

Since , it follows from Lemma 9 that . Therefore, for , which contradicts the supposition that .

#### 3. Main Results

Theorem 12. *Let be continuous with nonempty compact convex values on satisfying condition (3). Suppose that Algorithm 2 generates an infinite sequence . If the solution set of is nonempty, then globally converges to a solution of satisfying .*

*Proof. *Since , by Lemma 9 and the definition of projection, it follows that
Therefore, is a bounded sequence.

Since , . Since
it follows that
Thus, it follows from Lemma 7 (ii) that
that is,
Thus, the sequence is nondecreasing and bounded and hence convergent, which implies that
On the other hand, since
and since , we have
where the second inequality follows from (5). Therefore, it follows from (27) that
Since is continuous with compact values, Proposition 3.11 in [26] implies that is a bounded set, and so the sequences and are bounded. Thus, the continuity of implies that is a bounded set. Therefore, is bounded. It follows that
By the boundedness of , there exists a convergent subsequence converging to .

If is a solution of the problem (1), we show next that the whole sequence converges to . Let . Since , by Lemma 9,
Therefore,
Thus,
Letting in (34), we have
where the last inequality follows from Lemma 7 (i) and the fact that and . Therefore,
Thus, the sequence has a unique cluster point , which shows the global convergence of .

Suppose now that is not a solution of the problem (1). We show first that in Algorithm 2 cannot tend to . Since is continuous with compact values, Proposition in [26] implies that is a bounded set, and so the sequence is bounded. Therefore, there exists a subsequence converging to . Since is upper semicontinuous with compact values, Proposition in [26] implies that is closed, and so . By the definition of , we have
If , then . The lower continuity of , in turn, implies the existence of such that converges to . Since , we have and . Therefore and
Letting in (38), we have
with being continuous. It follows from Lemma 5 that
Thus, we obtain the contradiction because . Therefore, is bounded and so is .

By (31) and the boundedness of , we obtain that . Since is continuous and the sequences and are bounded, there exists an accumulation point of such that . This implies that solves the variational inequality (1). Similar to the preceding proof, we obtain that globally converges to .

*Remark 13. *In [11], the mapping is required to be pseudomonotone. Since the pseudomonotonicity implies condition (3), our assumptions of the mapping are more general.

Theorem 14. *Let be continuous with nonempty compact convex values on satisfying condition (3). Suppose that Algorithm 2 generates an infinite sequence . Then the solution set of is empty if and only if the sequence generated by Algorithm 2 diverges.*

*Proof. *In view of Theorem 12, it is sufficient to prove that if the solution set is empty, then the sequence generated by Algorithm 2 diverges. Since inequality (26) also holds in this case, the sequence is still nondecreasing. We claim that
Otherwise, is bounded, and hence it follows from (26) that
A similar discussion as in Theorem 12 would lead to the conclusion that every cluster point of is a solution of , which contradicts the emptiness of the solution set to .

#### 4. Numerical Experiments

In this section, we present some numerical experiments for the proposed algorithm. The MATLAB codes are run on a PC (with CPU Intel P-T2390) under MATLAB Version 7.0.1.24704(R14) Service Pack 1. We compare the performance of our Algorithm 2 and the algorithms in [5, 11, 12, 22]. In Tables 1, 2, 3, and 4, “It.” denotes number of iteration, and “CPU" denotes the CPU time in seconds. The tolerance means that when , the procedure stops.

*Example 15. *Let ,
and be defined by
Then, the set and the mapping satisfy the assumptions of Theorem 12, and is a solution of the multivalued variational inequality. Example 15 is tested in [5, 12, 22]. We choose , and for our algorithm and Algorithm 1 in [5]; , and for Algorithm 1 in [12]; , and for Algorithm 1 in [22]; for Algorithm 3.1 in [11]. We use as the initial point (Tables 1– 4).

#### Acknowledgments

This work is partially supported by Natural Science Foundation Project of CQ CSTC (no. 2010BB9401), Science and Technology Project of Chongqing Municipal Education Committee of China (no. KJ110509), and Foundation of Chongqing University of Posts and Telecommunications for the Scholars with Doctorate (no. A2012-04).

#### References

- T. Q. Bao and P. Q. Khanh, “A projection-type algorithm for pseudomonotone nonlips-chitzian multivalued variational inequalities, in generalized convexity,generalized monotonicity and applications,” in
*Proceedings of the 7th International Symposium on Generalized Convexity and Generalized Monotonicity*, A. Eberhard, N. Hadjisavvas, and D. T. Lus, Eds., pp. 113–129, Springer, 2005. - Y. Censor, A. Gibali, and S. Reich, “Extensions of Korpelevich's extragradient method for the variational inequality problem in Euclidean space,”
*Optimization*, vol. 61, no. 9, pp. 1119–1132, 2012. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - F. Facchinei and J.-S. Pang,
*Finite-Dimensional Variational Inequalities and Complementary Problems*, Springer, New York, NY, USA, 2003. - Y. He, “A new double projection algorithm for variational inequalities,”
*Journal of Computational and Applied Mathematics*, vol. 185, no. 1, pp. 166–173, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - F. Li and Y. He, “An algorithm for generalized variational inequality with pseudomonotone mapping,”
*Journal of Computational and Applied Mathematics*, vol. 228, no. 1, pp. 212–218, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. V. Solodov and B. F. Svaiter, “A new projection method for variational inequality problems,”
*SIAM Journal on Control and Optimization*, vol. 37, no. 3, pp. 765–776, 1999. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - D. Sun, “A class of iterative methods for solving nonlinear projection equations,”
*Journal of Optimization Theory and Applications*, vol. 91, no. 1, pp. 123–140, 1996. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Y. J. Wang, N. H. Xiu, and J. Z. Zhang, “Modified extragradient method for variational inequalities and verification of solution existence,”
*Journal of Optimization Theory and Applications*, vol. 119, no. 1, pp. 167–183, 2003. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - E. Allevi, A. Gnudi, and I. V. Konnov, “The proximal point method for nonmonotone variational inequalities,”
*Mathematical Methods of Operations Research*, vol. 63, no. 3, pp. 553–565, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - A. Auslender and M. Teboulle, “Lagrangian duality and related multiplier methods for variational inequality problems,”
*SIAM Journal on Optimization*, vol. 10, no. 4, pp. 1097–1115, 2000. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - H. Chen, “A new extragradient method for generalized variational inequality in Euclidean space,”
*Fixed Point Theory and Applications*, vol. 2013, article 139, 2013. View at Publisher · View at Google Scholar · View at MathSciNet - C. Fang and Y. He, “A double projection algorithm for multi-valued variational inequalities and a unified framework of the method,”
*Applied Mathematics and Computation*, vol. 217, no. 23, pp. 9543–9551, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - S. C. Fang and E. L. Peterson, “Generalized variational inequalities,”
*Journal of Optimization Theory and Applications*, vol. 38, no. 3, pp. 363–383, 1982. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. Fukushima, “The primal Douglas-Rachford splitting algorithm for a class of monotone mappings with application to the traffic equilibrium problem,”
*Mathematical Programming A*, vol. 72, no. 1, pp. 1–15, 1996. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Y. He, “Stable pseudomonotone variational inequality in reflexive Banach spaces,”
*Journal of Mathematical Analysis and Applications*, vol. 330, no. 1, pp. 352–363, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - R. Saigal, “Extension of the generalized complementarity problem,”
*Mathematics of Operations Research*, vol. 1, no. 3, pp. 260–266, 1976. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - G. Salmon, J.-J. Strodiot, and V. H. Nguyen, “A bundle method for solving variational inequalities,”
*SIAM Journal on Optimization*, vol. 14, no. 3, pp. 869–893, 2003. View at Publisher · View at Google Scholar · View at MathSciNet - M. V. Solodov and B. F. Svaiter, “Forcing strong convergence of proximal point iterations in a Hilbert space,”
*Mathematical Programming A*, vol. 87, no. 1, pp. 189–202, 2000. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - R. T. Rockafellar, “Monotone operators and the proximal point algorithm,”
*SIAM Journal on Control and Optimization*, vol. 14, no. 5, pp. 877–898, 1976. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - I. V. Konnov, “On the rate of convergence of combined relaxation methods,”
*Izvestiya Vysshikh Uchebnykh Zavedeniĭ. Matematika*, vol. 37, no. 12, pp. 89–92, 1993. View at Google Scholar · View at MathSciNet - I. V. Konnov,
*Combined Relaxation Methods for Variational Inequalities*, Springer, Berlin, Germany, 2001. View at Publisher · View at Google Scholar · View at MathSciNet - C. J. Fang and Y. R. He, “An extragradient method for generalized variational inequality,”
*Pacic Journal of Optimization*, vol. 9, no. 1, pp. 47–59, 2013. View at Google Scholar - C. J. Fang, S. L. Chen, and C. D. Yang, “An algorithm for solving multi-valued variational inequality,”
*Journal of Inequalities and Applications*, vol. 2013, article 218, 2013. View at Publisher · View at Google Scholar - S. Karamardian, “Complementarity problems over cones with monotone and pseudomonotone maps,”
*Journal of Optimization Theory and Applications*, vol. 18, no. 4, pp. 445–454, 1976. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - E. H. Zarantonello, “Projections on convex sets in Hilbert space and spectral theory,” in
*Contributions to Nonlinear Functional Analysis*, E. H. Zarantonello, Ed., Academic Press, New York, NY, USA, 1971. View at Google Scholar - J.-P. Aubin and I. Ekeland,
*Applied Nonlinear Analysis*, John Wiley & Sons, New Yor, NY, USA, 1984. View at MathSciNet