Modulus-Based Matrix Splitting Iteration Methods for a Class of Stochastic Linear Complementarity Problem

For the expected value formulation of stochastic linear complementarity problem, we establish modulus-based matrix splitting iteration methods. The convergence of the new methods is discussed when the coefficient matrix is a positive definite matrix or a positive semi-definite matrix, respectively. The advantages of the new methods are that they can solve the large scale stochastic linear complementarity problem, and spend less computational time. Numerical results show that the new methods are efficient and suitable for solving the large scale problems.


Introduction
The complementarity problems have been widely used in the engineering design, information technology, economic equilibrium, etc. Since some elements may involve uncertain data in practical applications, many problems can be attributed to stochastic variational inequality problems or stochastic linear complementarity problems, and arouse the attention of many researchers. Gurkan et al. [1] proposed the expected value (EV) formulation of stochastic variational inequality by using the sample-path method. Chen and Fukushima [2] proposed the expected residual minimization (ERM) formulation for stochastic linear complementarity problems by quasi-Monte Carlo methods. Lin and Fukushima [3] proposed the stochastic mathematical programs with equilibrium constraints (SMPEC) for the stochastic nonlinear complementarity problems. Zhou and Cacceta [4] transformed the monotone stochastic linear complementarity problem (SLCP) in finite sample space into a constrained minimization problem, and solved it with the Feasible Semi-smooth Newton Method. Mangasarian and Ren [5] given a new error bound for the monotone LCP based on the error bounds.
Chen et al. [6] studied the SLCP involving a random matrix whose expectation matrix is positive semi-definite. Zhang and Chen [7] proposed a smooth projection gradient algorithm to solve the SLCP. However, these methods are only suitable for solving the small-scale SLCP.
In recent years, some scholars have proposed a series of methods for the study of large-scale complementarity problems. Dong and Jiang [8] proposed a class of modified modulus-based method. Bai [9] presented a class of modulus-based matrix splitting iteration methods. Bai and Evans [10] [11] also proposed a class of modulus-based synchronous multi-splitting (MSM) iteration methods. Bai and Zhang [12] further proposed a synchronous two-stage multi-splitting iteration method, which can be applied to solving the large-scale linear complementarity problems. Zhang [13] summarized the latest development and achievements of the modulus-based matrix splitting iteration methods, including the corresponding multi-splitting iteration methods, etc. Zhang [14] improved the convergence theorem of matrix multi-splitting methods for linear complementarity problems. Such methods are easy to be implemented and very efficient in practical applications, and there is no need to project iteration results into space n R + . Li et al. [15] [16] [17] [18] applied a class of modulus-based matrix splitting iteration methods to solving the nonlinear complementarity problem. Numerical results show that the methods are efficient. In the past decade many scholars have made many new achievements in this field, see the literatures [19]- [30].
In this paper, we extend the modulus-based matrix splitting iteration methods to solve the large-scale stochastic linear complementarity problems. We also prove the convergence of these methods when the coefficient matrix is a positive definite matrix or a positive semi-definite matrix. The numerical results show that these methods are efficient.
The outline of the paper is as follows. In Section 2 we present some necessary results and lemmas. In Section 3 we establish the modulus-based matrix splitting iteration methods for solving the SLCP. The convergence of the methods is proved in Section 4. The numerical results are shown in Section 5. Finally, in Section 6, we give some concluding remarks.

Preliminaries
In this section, we briefly introduce some necessary results and lemmas.
and ( ) n q ω ∈  are the rand matrices and vectors for ω ∈ Ω , respectively.
Usually there not exists z for all ω ∈ Ω for Problem (1). In order to get a reasonable solution of (1), in this paper we use the EV formulation proposed by Gurkan et al. [1].
The Expected Value (EV) Formulation [1]:  , and E be the expectation. We consider the following EV formulation: finding a vector We briefly denote it as

Modulus-Based Matrix Splitting Iteration Methods
In this section, we aim at the EV formulation of the stochastic linear complementarity problem (2). We give some corresponding modulus-based matrix splitting iteration methods.
For the strong monotone stochastic linear complementarity problem, the coefficient matrix is positive definite. For this case, we can apply the method proposed by Dong and Jiang [8].
Method 3.1 Step 1: Select an arbitrary initial vector ( ) 0 n x ∈  and set : 0 k = ; Step 2: satisfies the termination rule, then stop; otherwise, set : 1 k k = + and return to Step 2.
Unfortunately, the coefficient matrices of some stochastic linear complementarity problems are positive semi-definite, Method 3.1 is not suitable for solving the problem (2). Cottle et al. [31] presented a regularization method.
Based on this method, we establish a regularized modulus-based matrix splitting iteration method. To simplify the notation, we will denote { } Step 1: Select a positive number 0 R ε ∈ and an arbitrary initial vector x ε ∈  , and set : 0 Step 2: Generate the iteration sequence Step

Convergence
In this section, we analyze the convergence of

The Case of Symmetric Positive Definite Matrix
We first discuss the convergence of Method 3.1 when the coefficient matrix is symmetric positive definite.
Theorem 1 Suppose that the system matrix

The Case of Symmetric Positive Semi-Definite Matrix
We now discuss the convergence of

Numerical Results
In this section, we test some numerical results to show the efficiency of our methods. Let ∈ Ω = =  . The steps to generate test problems can be found in the literature [4]. Numerical experimental results are shown in Tables 1-3. Table 1 shows that Method 3.1 is effective when the coefficient matrix is symmetric positive definite.

Conclusion
In this paper, we study the fast numerical methods for solving the stochastic linear complementarity problems. Firstly, we convert the expected value formulation of stochastic linear complementarity problems into the equivalent fixed point equations, then we establish a class of modulus-based matrix splitting iteration methods, and analyze the convergence of the method. These new methods can be applied to solve the large-scale stochastic linear complementarity problems. The numerical results also show the effectiveness of the new methods.