An Improved Convergence Condition of the MMS Iteration Method for Horizontal LCP of H + -Matrices

: In this paper, inspired by the previous work in (Appl. Math. Comput., 369 (2020) 124890), we focus on the convergence condition of the modulus-based matrix splitting (MMS) iteration method for solving the horizontal linear complementarity problem (HLCP) with H + -matrices. An improved convergence condition of the MMS iteration method is given to improve the range of its applications, in a way which is better than that in the above published article.


Introduction
As is known, the horizontal linear complementarity problem, for the given matrices A, B ∈ R n×n , is to find that two vectors z, w ∈ R n satisfy Az = Bw + q ≥ 0, z ≥ 0, w ≥ 0 and z T w = 0, (1) where q ∈ R n is given, which is often abbreviated as HLCP.If A = I in (1), the HLCP (1) is no other than the classical linear complementarity problem (LCP) in [1], where I denotes the identity matrix.This implies that the HLCP (1) is a general form of the LCP.
The HLCP (1), used as a useful tool, often arises in a diverse range of fields, including transportation science, telecommunication systems, structural mechanics, mechanical and electrical engineering, and so on, see [2][3][4][5][6][7].In the past several years, some efficient algorithms have been designed to solve the HLCP (1), such as the interior point method [8], the neural network [9], and so on.Particularly, in [10], the modulus-based matrix splitting (MMS) iteration method in [11] was adopted to solve the HLCP (1).In addition, the partial motivation of the present paper is from complex systems with matrix formulation, see [12][13][14] for more details.
Recently, Zheng and Vong [15] further discussed the MMS method, as described below.
The MMS method [10,15].Let Ω be a positive diagonal matrix and r > 0, and let A = M A − N A and B = M B − N B be the splitting of matrices A and B, respectively.Assume that (z (0) , w (0) ) is an arbitrary initial vector.For k = 0, 1, 2, . . .until the iteration sequence (z (k) , w (k) ) converges, compute (z (k+1) , w (k+1) ) by where x (k+1) is obtained by For the later discussion, some preliminaries are gone over.For a square matrix A = (a ij ) ∈ R n×n , |A| = (|a ij |), and A = ( a ij ), where a ii = |a ii | and a ij = −|a ij | for i = j.A matrix A = (a ij ) ∈ R n×n is called a non-singular M-matrix if A −1 ≥ 0 and a ij ≤ 0 for i = j; an H-matrix if its comparison matrix A is a non-singular M-matrix; an H + -matrix if it is an H-matrix with positive diagonals; and a strictly diagonally dominant (s.d.d. For the MMS method with H + -matrix, two new convergence conditions are obtained in [15], which are weaker than the corresponding convergence conditions in [10].One of these is given below.Theorem 1 ([15]).Assume that A, B ∈ R n×n are two H + -matrices and Ω = diag(ω jj ) ∈ R n×n with ω jj > 0, i, 2, . . ., n, Let A = M A − N A be an H-splitting of A, B = M B − N B be an H-compatible splitting of B, and M A + M B Ω be an H + -matrix.Then the MMS method is convergent, provided one of the following conditions holds: with At present, the difficulty in Theorem 1 is to check the condition (4).Besides that, the condition (4) of Theorem 1 is limited by the parameter k.That is to say, if the choice of k is improper, then we cannot use the condition (4) of Theorem 1 to guarantee the convergence of the MMS method.To overcome this drawback, the purpose of this paper is to provide an improved convergence condition of the MMS method, for solving the HLCP of H + -matrices, to improve the range of its applications, in a way which is better than that in Theorem 1 [15].

An Improved Convergence Condition
In fact, by investigating condition (b) of Theorem 1, we know that the left inequality in (4) may have a flaw.Particularly, when the choice of k is improper, we cannot use condition (b) of Theorem 1 to guarantee the convergence of the MMS method.For instance, we consider two matrices To make A and B satisfy the convergence conditions of Theorem 1, we take By the simple computations, For convenience, we take D = D 1 = I, where I denotes the identity matrix.By simple calculations, we have and Further, we have Obviously, when k ≤ 1/3, we naturally do not get that This implies that condition (b) of Theorem 1 may be invalid when we use condition (b) of Theorem 1 to judge the convergence of the MMS method for solving the HLCP.To overcome this disadvantage, we obtain an improved convergence condition for the MMS method, see Theorem 2, whose proof is similar to the proof of Theorem 2.5 in [15].
Theorem 2. Assume that A, B ∈ R n×n are two H + -matrices, and Let A = M A − N A be an H-splitting of A, B = M B − N B be an H-compatible splitting of B, and M A + M B Ω be an H + -matrix.Then the MMS method is convergent, provided one of the following conditions holds: where D is a positive diagonal matrix, such that M A + M B Ω D is an s.d.d.matrix.
Proof.For Case (a), see the proof of Theorem 2.5 in [15].For Case (b), by simple calculations, we have Making use of Equation ( 6), based on the proof of Theorem 2.5 in [15], we have Since M A + M B Ω is an H + -matrix, it follows that Ŝ = M A + M B Ω is a non-singular M-matrix, and the existence of such a matrix D (see [16], p. 137) satisfies From the left inequality in (5), we have Further, based on the inequality (7), we have Thus, based on Lemma 2.3 in [15], we have The proof of Theorem 2 is completed.
Comparing Theorem 2 with Theorem 1, the advantage of the former is that condition (b) of Theorem 2 is not limited by the parameter k of the latter.Besides that, we do not need to find two positive diagonal matrices D and D Incidentally, there exists a simple approach to obtain a positive diagonal matrix D in Theorem 2: first, solving the system Āx = e gives the positive vector x, where Ā = M A + M B Ω ; secondly, we take D = diag( Ā−1 e), which can make M A + M B Ω D an s.d.d.matrix.
In addition, if the H + -matrix M A + M B Ω itself is an s.d.d.matrix, then we can take D = I in Theorem 2. In this case, we can obtain the following corollary.

Corollary 1.
Assume that A, B ∈ R n×n are two H + -matrices, and Let A = M A − N A be an H-splitting of A, B = M B − N B be an H-compatible splitting of B, and the H + -matrix M A + M B Ω be an s.d.d.matrix.Then, the MMS method is convergent, provided one of the following conditions holds:

Numerical Experiments
In this section, we consider a simple example to illustrate our theoretical results in Theorem 2. All the computations are performed in MATLAB R2016B.
In our calculations, we take µ = 4 and ν = 0 for A and B in Example 1, x (0) = (2, 2, . . ., 2) T ∈ R n is used for the initial vector.The modulus-based Jacobi (NMJ) method and Gauss-Seidel (NMGS) method, with r = 2, are adopted.The NMJ and NMGS methods are stopped once the number of iterations is larger than 500 or the norm of residual vectors (RES) is less than 10 −6 , where Here, we consider two cases of Theorem 2. When Ω ≥ D A D −1 B , we take Ω = 2I for the NMJ method and the NMGS method.In this case, Table 1 B , we take D = I, and obtain that I < Ω < 2I and M A + M B Ω D is an s.d.d.matrix.In this case, we take Ω = 1.5I and Ω = 1.2I for the NMJ and NMGS methods, and obtain Tables 2 and 3.The numerical results in Tables 1-3 not only further confirm that the MMS method is feasible and effective, but also show that the convergence condition in Theorem 2 is reasonable.

Conclusions
In this paper, the modulus-based matrix splitting (MMS) iteration method for solving the horizontal linear complementarity problem (HLCP) with H + -matrices, has been further considered.The main aim of this paper is to present an improved convergence condition of the MMS iteration method, to enlarge the range of its applications, in a way which is better than previous work [15].
1 , such that ( M A − |N A |)D and ( M B − |N B |)D 1 are, respectively, s.d.d.matrices, we just find one positive diagonal matrix D, such that M A + M B Ω D is an s.d.d.matrix.