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 R + n . 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.

In this section, we briefly introduce some necessary results and lemmas.

Let A = ( a i j ) ∈ ℝ n × n , A is said to be positive semi-definite if x T A x ≥ 0 for all x ∈ ℝ n , and positive definite if x T A x > 0 for all x ∈ ℝ n \ { 0 } . A ∈ R n × n is called a P 0 -matrix if all of its principle minors are nonnegative.

Let ( Ω , F , Ρ ) be a probability space, where Ω is a sample subset of ℝ m . Suppose that probability distribution is known, we consider the stochastic linear complementarity problem (SLCP): finding a vector z ∈ ℝ n such that

M ( ω ) z + q ( ω ) ≥ 0 ,     z ≥ 0 ,     z T ( M ( ω ) z + q ( ω ) ) = 0 ,     ω ∈ Ω . (1)

where M ( ω ) ∈ ℝ n × n and q ( ω ) ∈ ℝ n 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] :

Let F ( z , ω ) = M ( ω ) z + q ( ω ) , M ¯ = E [ M ( ω ) ] , q = E [ q ( ω ) ] , and E be the expectation. We consider the following EV formulation: finding a vector z ∈ ℝ n such that

F ¯ ( z ) = E [ F ( z , ω ) ] = M ¯ z + q ≥ 0 ,     z ≥ 0 ,   z T F ¯ ( z ) = 0. (2)

We briefly denote it as LCP ( q , M ¯ ) .

Define

RES ( z ) : = min ( z , M ¯ z + q )

where the min operator denotes the componentwise minimum of two vectors. It is generally known that z * solves the LCP ( q , M ¯ ) if and only if z * solves the equations

R E S ( z ) = 0

The function RES is called the natural residual of the LCP ( q , M ¯ ) and is often used in error analysis.

Lemma 1 (see [8] ) Let α ∈ ( 0 , + ∞ ) be a scalar, then the LCP ( q , M ¯ ) (2) is equivalent to the following fixed-point problem: finding x ∈ ℝ n , satisfying that

( α I + M ¯ ) x = ( α I − M ¯ ) | x | − q (3)

Moreover, if x is the solution of (3), then

r : = α ( | x | − x ) , z : = | x | + x (4)

define a solution pair of Problem (2). On the other hand, if the vector pair z and r solves Problem (2), then x : = 1 / 2 ( z − r / α ) solves the fixed-point problem (3).

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 x ( 0 ) ∈ ℝ n and set k : = 0 ;

Step 2: Calculate x ( k + 1 ) through the iteration scheme

( α I + M ¯ ) x ( k + 1 ) = ( α I − M ¯ ) | x ( k ) | − q

Step 3: Let z ( k + 1 ) = | x ( k + 1 ) | + x ( k + 1 ) , if z ( k + 1 ) satisfies the termination rule, then stop; otherwise, set k : = k + 1 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 { ε } and { x ε k } for { ε k } and { x ε k k } , and denote the regularization problem for LCP ( q , M ¯ + ε I ) .

Method 3.2

Step 1: Select a positive number ε 0 ∈ R and an arbitrary initial vector x ε ( 0 ) ∈ ℝ n , and set k : = 0 ;

Step 2: Generate the iteration sequence x ε ( k + 1 ) through solving the following equations

( α I + M ¯ + ε I ) x ε ( k + 1 ) = ( α I − M ¯ − ε I ) | x ε ( k ) | − q .

Let z ε ( k + 1 ) = | x ε ( k + 1 ) | + x ε ( k + 1 ) .

Step 3: Set ε = α ε , where α ∈ ( 0 , 1 ) is a positive number, k : = k + 1 , and return to Step 2.

In this section, we analyze the convergence of Method 3.1 and Method 3.2 when the coefficient matrix of the LCP ( q , M ¯ ) is a symmetric positive definite matrix and a symmetric positive semi-definite matrix.

In this section, we test some numerical results to show the efficiency of our methods. Let, n be the order of the matrix, IT denote the average iteration steps, and CPU denote the average iteration time.

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.

Table 2 and Table 3 list the numerical experimental results of Method 3.2 when the coefficient matrix is symmetric positive semi-definite. We know that Method 3.2 is effective. (In the following, we briefly denote that Feasible Semi-smooth Newton Method is FSNM.)

Table 4 shows that, Method 3.2 is more effective than FSNM [19] . By Tables 1-4, we know that Method 3.2 improves the computational efficiency and is suitable for solving the large scale problems.

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.

This work was supported by Natural Science Foundation of China (11661027), National Project for Research and Development of Major Scientific Instruments (61627807), and Guangxi Natural Science Foundation (2015 GXNSFAA 139014).

The authors declare no conflicts of interest regarding the publication of this paper.

Lu, Q.Q. and Li, C.L. (2019) Modulus-Based Matrix Splitting Iteration Methods for a Class of Stochastic Linear Complementarity Problem. American Journal of Operations Research, 9, 245-254. https://doi.org/10.4236/ajor.2019.96016