Loading [MathJax]/jax/output/SVG/jax.js
Research article Special Issues

New error bounds for the tensor complementarity problem


  • This paper discusses new error bounds for the tensor complementarity problem using a P-tensor. A new lower error bound and a global error bound are presented for such a problem. It is proved that the norm of the exact solution of the tensor complementarity problem with a P-tensor has a lower bound and an upper bound. When the order of a tensor is 2, all the results for the tensor complementarity problem obtained reduce to those for the linear complementarity problem.

    Citation: Xin Liu, Guang-Xin Huang. New error bounds for the tensor complementarity problem[J]. Electronic Research Archive, 2022, 30(6): 2196-2204. doi: 10.3934/era.2022111

    Related Papers:

    [1] Yan Li, Yaqiang Wang . Some new results for $ B_1 $-matrices. Electronic Research Archive, 2023, 31(8): 4773-4787. doi: 10.3934/era.2023244
    [2] Haiyan Song, Fei Sun . A numerical method for parabolic complementarity problem. Electronic Research Archive, 2023, 31(2): 1048-1064. doi: 10.3934/era.2023052
    [3] Dongmei Yu, Yifei Yuan, Yiming Zhang . A preconditioned new modulus-based matrix splitting method for solving linear complementarity problem of $ H_+ $-matrices. Electronic Research Archive, 2023, 31(1): 123-146. doi: 10.3934/era.2023007
    [4] Lot-Kei Chou, Siu-Long Lei . High dimensional Riesz space distributed-order advection-dispersion equations with ADI scheme in compression format. Electronic Research Archive, 2022, 30(4): 1463-1476. doi: 10.3934/era.2022077
    [5] Zhengmao Chen . A priori bounds and existence of smooth solutions to a $ L_p $ Aleksandrov problem for Codazzi tensor with log-convex measure. Electronic Research Archive, 2023, 31(2): 840-859. doi: 10.3934/era.2023042
    [6] Peng Zhou, Teng Wang . The parameter-Newton iteration for the second-order cone linear complementarity problem. Electronic Research Archive, 2022, 30(4): 1454-1462. doi: 10.3934/era.2022076
    [7] Jin Wang . Least squares solutions of matrix equation $ AXB = C $ under semi-tensor product. Electronic Research Archive, 2024, 32(5): 2976-2993. doi: 10.3934/era.2024136
    [8] Yiyuan Qian, Kai Zhang, Jingzhi Li, Xiaoshen Wang . Adaptive neural network surrogate model for solving the implied volatility of time-dependent American option via Bayesian inference. Electronic Research Archive, 2022, 30(6): 2335-2355. doi: 10.3934/era.2022119
    [9] Shunjie Bai, Caili Sang, Jianxing Zhao . Localization and calculation for C-eigenvalues of a piezoelectric-type tensor. Electronic Research Archive, 2022, 30(4): 1419-1441. doi: 10.3934/era.2022074
    [10] Yao Sun, Lijuan He, Bo Chen . Application of neural networks to inverse elastic scattering problems with near-field measurements. Electronic Research Archive, 2023, 31(11): 7000-7020. doi: 10.3934/era.2023355
  • This paper discusses new error bounds for the tensor complementarity problem using a P-tensor. A new lower error bound and a global error bound are presented for such a problem. It is proved that the norm of the exact solution of the tensor complementarity problem with a P-tensor has a lower bound and an upper bound. When the order of a tensor is 2, all the results for the tensor complementarity problem obtained reduce to those for the linear complementarity problem.



    The set of all real rth-order n-dimensional tensors is denoted by Tr,n, where r3 and n2 are positive integers. A tensor ATr,n is called a P-tensor [1], if for each vector yRn{0}, there exists an index jJn such that

    yj(Ayr1)j0,

    where Jn defines the index set {1,2,...,n}. For a given vector qRn and a tensor A=(aj1...jr)Tr,n with a multi-array of real entries aj1...jr, where jiJn for i{1,...,r}, the tensor complementarity problem denoted by the TCP(A,q), is to find a real vector yRn such that

    y0,q+Ayr10,andyT(q+Ayr1)=0, (1.1)

    where Ayr1Rn and the jth element of Ayr1 is determined by

    (Ayr1)j:=nj2,...,jr=1ajj2...jryj2yjr

    for jJn. When the tensor A is a matrix, the TCP(A,q) reduces to the linear complementarity problem [2], denoted by the LCP(M,q), which is to find a real vector yRn such that

    y0,q+My0,andyT(q+My)=0,

    where MRn×n and qRn. Error bounds of the LCP(M,q) have been extensively studied in recent years. For example, Mathias and Pang in [3] introduced fundamental quantities associated with an arbitrary P-matrix and derived global upper and lower error bounds for the approximate solution of the LCP(M,q) with a P-matrix. Global upper and lower error bounds to the LCP(M,q) with a P-matrix were given in [2] by Cottle et al. We refer to [4,5,6,7] for more new results on the estimation of error bounds to the LCP(M,q) with a P-matrix. We denote the LCP(M,q) with a P-matrix by the LCP(M,q;p) and the TCP(A,q) with a P-tensor by the TCP(A,q;p).

    In recent years, the properties of several types of structured tensors in [8,9] were used to investigate some properties about the TCP(A,q) solution set. The nonempty of the solution set of the TCP(A,q) and the compactness of the solution set of the TCP(A,q) were explored in [10,11] by Che et al. Huang et al. [12,13] proved the existence of the solution of the TCP(A,q). Yu et al. [14] evaluated the stability features of the TCP(A,q) solution set. Zheng et al. [15] extended the results in [3] to the tensors and presented two TCP(A,q;p) global error bounds. In addition, some applications of the TCP(A,q) were given in [16]. In this paper, we mainly expand the error bounds of the approximate solution of the LCP(M,q;p) in [2,7] to the cases of the TCP(A,q;p). We also extend the bound of the norm of the exact solution of the LCP(M,q;p) in [7] to the case of the TCP(A,q;p).

    The rest of this paper is arranged as follows. Section 2 presents some fundamental concepts and summarizes some related results that will be used later. Some new error bounds of the TCP(A,q;p) are given in Section 3. Section 4 drawn some conclusions.

    In this section, we summarize some results and introduce some denotations that will be used later. The following results are true for a P-tensor.

    Lemma 2.1. ([17]) For any P-tensor ATr,n and any qRn, the solution set of the TCP(A,q) is nonempty and compact.

    Lemma 2.2. ([15,18]) There does not exist an odd order P-tensor.

    For any tensor ATr,n, the infinite norm of A is defined as follows:

    A:=maxjJnnj2,...,jr=1ajj2jr.

    Let T:RnRn be an operator, T is said to be positively homogeneous iff T(ty)=tT(y) holds, for any positive t, y Rn. Song and Qi in [1] introduced two quantities of the form

    α(TA):=miny=1maxjJnyj(TAy)j (2.1)

    for an arbitrary positive even integer r, and

    α(FA):=miny=1maxjJnyj(FAy)j. (2.2)

    Here TA:RnRn and FA:RnRn are two operators that are positively homogenous and defined by

    TAy:={y2r2Ayr1,y0,0,y=0, (2.3)

    and

    FAy:=(Ayr1)[1r1], (2.4)

    respectively. Here y[1r1]=(y1r11,y1r12,...,y1r1n)T for y=(y1,y2,...,yn)T. It is true that

    α(FA)y2maxjJnyj(FAy)j=maxjJnyj(Ayr1)1r1j. (2.5)

    The following result in [1] gives two criterions of a P-tensor based on α(FA) and α(TA).

    Lemma 2.3. ([1]) Let ATr,n, A is a P-tensor iff α(TA)>0. When r is even, it is especially true that A is a P-tensor iff α(FA)>0.

    In this section, new error bounds of the approximate solution of (1.1) are proposed. Denote the solution set of (1.1) by the set of the form

    S:={yRny0,q+Ayr10,yT(q+Ayr1)=0}, (3.1)

    then according to Lemma 2.1, S is nonempty and compact when ATr,n is a P-tensor. Thus for any vRn, there is a vector ˜yS such that

    v˜y=minySvy. (3.2)

    We also need the results as follows.

    Lemma 3.1. For an arbitrary positive integer r, it is true that

    Ayr1≤∥yr1A. (3.3)

    Proof. In fact

    Ayr1=maxjJn(Ayr1)j=maxjJnnj2,...,jr=1ajj2jryj2yjrmaxjJnnj2,...,jr=1ajj2jr∣∥yr1=∥yr1A.

    This completes the proof.

    Theorem 3.1. Suppose qRn, ATr,n is a P-tensor. Let ˜y be the unique solution of the TCP(A,q). We have

    ˜v˜ymax{1,A1r1}v˜y,vRn, (3.4)

    where ˜v˜y is defined by

    ˜v˜y:=min{v,(A(v˜y)r1)[1r1]+(q+A˜yr1)[1r1]}. (3.5)

    Proof. First show that for any v,wRn, the following inequality holds,

    ˜v˜y˜w˜ymax{1,A1r1}(v˜y+w˜y), (3.6)

    where ˜v˜y is defined in (3.5) and ˜w˜y is determined by (3.5) with v being replaced by w. Denote y=˜v˜y and x=˜w˜y, suppose that yjxj∣=yjxj,jJn, then we have

    y=˜v˜y=min{v,(A(v˜y)r1)[1r1]+(q+A˜yr1)[1r1]},
    x=˜w˜y=min{w,(A(w˜y)r1)[1r1]+(q+A˜yr1)[1r1]}.

    Thus,

    yjvjandyj(A(v˜y)r1)1r1j+(q+A˜yr1)1r1j,
    xj=wjorxj=(A(w˜y)r1)1r1j+(q+A˜yr1)1r1j.

    In the following, we prove this result in two cases. If xj=wj, then

    yjxj∣=yjxjvjwjvw=∥(v˜y)(w˜y)v˜y+w˜y.

    Therefore

    yxv˜y+w˜y. (3.7)

    If xj=(A(w˜y)r1)1r1j+(q+A˜yr1)1r1j, then

    yjxj∣=yjxj(A(v˜y)r1)1r1j(A(w˜y)r1)1r1jA(v˜y)r11r1+A(w˜y)r11r1.

    It follows from Lemma 3.1, that

    A(v˜y)r1v˜yr1A,A(w˜y)r1w˜yr1A.

    Thus it holds that

    yxA1r1(v˜y+w˜y). (3.8)

    Combining (3.7) and (3.8) results in (3.6). Let w=˜y in (3.6), then

    ˜v˜ymin{˜y,(q+A˜yr1)[1r1]}max{1,A1r1}v˜y. (3.9)

    Given that ˜y solves the TCP(A,q), then we have

    min{˜y,(q+A˜yr1)[1r1]}=0.

    It follows from (3.9) that

    min{v,(A(v˜y)r1)[1r1]+(q+A˜yr1)[1r1]}max{1,A1r1}v˜y,

    which implies (3.4). This completes the proof.

    We remark that a lower error bound (Theorem 3.1) for the TCP(A,q) with a P-tensor in this paper is sharper than the result (Theorem 3.2) in [15]. Because

    ˜v˜y1+A1r1˜v˜ymax{1,A1r1},˜v˜y:=min{v,(A(v˜y)r1)[1r1]+(q+A˜yr1)[1r1]}.

    The following result gives a global error bound on the approximate solution of (1.1).

    Theorem 3.2. Suppose ATr,n is a P-tensor. Let ˜y be the unique solution of the TCP(A,q). We have

    (1+A1r1)˜v˜y12α(FA)≤∥v˜y(1+A1r1)˜v˜y+12α(FA),vRn, (3.10)

    where

    1=(1+A1r1)2˜v˜y24α(FA)(˜v˜y)2k10

    with k1 satisfying (˜v˜y)k10,

    (A(v˜y)r1)1r1k1(v˜y)k1=maxjJn{(A(v˜y)r1)1r1j(v˜y)j},

    and α(FA) is defined in (2.2).

    Proof. We first show that the following inequality holds

    α(FA)v˜y2(1+A1r1)˜v˜yv˜y+(˜v˜y)2k10. (3.11)

    Suppose w=(q+A˜yr1)[1r1]. Since ˜y solves (1.1), we can prove that

    ˜y0,w0,˜yTw=0.

    Let x=v˜v˜y and z=(A(v˜y)r1)[1r1]+(q+A˜yr1)[1r1]˜v˜y, where ˜v˜y is defined as that in (3.5), then we have that

    x0,z0,xTz=0.

    Therefore, for any jJn,

    0xjzjxjwj˜yjzj+˜yjwj=(zw)j(x˜y)j=[(A(v˜y)r1)[1r1]˜v˜y]j(v˜v˜y˜y)j=(A(v˜y)r1)1r1j(v˜y)j(˜v˜y)j(v˜y)j(A(v˜y)r1)1r1j(˜v˜y)j+(˜v˜y)2j.

    Thus

    (˜v˜y)j(v˜y)j+(A(v˜y)r1)1r1j(˜v˜y)j(˜v˜y)2j(A(v˜y)r1)1r1j(v˜y)j.

    Let

    (A(v˜y)r1)1r1k1(v˜y)k1=maxjJn{(A(v˜y)r1)1r1j(v˜y)j},

    then

    maxjJn{(A(v˜y)r1)1r1j(v˜y)j}(˜v˜y)k1(v˜y)k1+(A(v˜y)r1)1r1k1(˜v˜y)k1(˜v˜y)2k1˜v˜yv˜y+(A(v˜y)r1)1r1˜v˜y(˜v˜y)2k1. (3.12)

    From Lemma 3.1, we have that

    (A(v˜y)r1)1r1v˜yA1r1. (3.13)

    Furthermore, we can deduce from (2.5) that

    α(FA)v˜y2maxjJn{(A(v˜y)r1)1r1j(v˜y)j}. (3.14)

    Therefore, combining (3.12)–(3.14) results in

    α(FA)v˜y2(v˜y)k1(˜v˜y)k1+(A(v˜y)r1)1r1k1(˜v˜y)k1(˜v˜y)2k1v˜y˜v˜y+(A(v˜y)r1)1r1˜v˜y(˜v˜y)2k1=∥v˜y(1+A1r1)˜v˜y(˜v˜y)2k1.

    If (˜v˜y)k1=0, then v=˜y; otherwise, (3.11) holds. According to Lemma 2.3, we have α(FA)>0, here A is a P-tensor. We can derive (3.10) by solving (3.11). The proof is completed.

    If we set v=0 in Theorem 3.2, we obtain a new bound for the norm of the exact solution of (1.1).

    Corollary 1. Let v=0, then under the conditions of Theorem 3.2, we have that

    ((q)+)k2(1+A1r1)22α(FA)≤∥˜y((q)+)k2(1+A1r1)+22α(FA), (3.15)

    where

    (q)+=min{q,0},
    2=((q)+)2k2(1+A1r1)24α(FA)((q)+)2k20

    with k2 satisfying

    ˜yk2(A˜yr1)1r1k2=maxjJn{˜yj(A˜yr1)1r1j},and((q)+)k20,

    and α(FA) is defined in (2.2).

    Proof. The proof is similar to that of Theorem 3.2 and is omitted.

    We have obtained a new lower error bound in Theorem 3.1 and a global error bound in Theorem 3.2 to the TCP(A,q;p).

    We remark that when r=2, Theorem 3.2 and Corollary 1 reduce to Theorem 3.1 in [7] and Corollary 3.1 in [7], respectively.

    When r=2, the tensor ATr,n reduces to a matrix M. Thus, FAy:=(Ayr1)[1r1]=My, we have

    α(FA):=miny=1maxjJnyj(FAy)j=miny=1maxjJnyj(My)j=α(M).

    Furthermore, we can deduce from (3.5) that

    ˜v˜y:=min{v,(A(v˜y)r1)[1r1]+(q+A˜yr1)[1r1]}=min{v,M(v˜y)+q+M˜y}=min{v,q+Mv}=u,

    which implies that

    1=(1+A1r1)2˜v˜y24α(FA)(˜v˜y)2k1=(1+M)2min{v,q+Mv}24α(M)(min{v,q+Mv})2k1=(1+M)2u24α(M)u2k1=△

    with k1 satisfying uk10 and

    (v˜y)k1(M(v˜y))k1=maxjJn{(v˜y)j(M(v˜y))j},

    and that

    2=((q)+)2k2(1+A1r1)24α(FA)((q)+)2k2=((q)+)2k2(1+M)24α(M)((q)+)2k2=

    with k2 satisfying ((q)+)k20 and

    ˜yk2(M˜y)k2=maxjJn{˜yj(M˜y)j}.

    Therefore, Theorem 3.1 reduces to Theorem 5.10.6 in [2].

    Similarly, when r=2, Theorem 3.2 is an extension of Theorem 3.1 in [7], and Corollary 1 is an extension of Corollary 3.1 in [7].

    In this paper, we present a new lower error bound and a result of the global error bound on the approximate solution to the TCP(A,q) with a P-tensor. The new bound on the norm of the exact solution to the TCP(A,q) with a P-tensor is derived.

    As we all know, error bounds have important applications in iterative methods for solving related optimization problems. For example, we use the obtained error bounds to analyze the convergence of the iterative method for solving the TCP(A,q). We can futher study the numerical algorithm to compute the error bounds of the TCP(A,q) with a P-tensor (see [19]).

    Research by G.X. Huang was supported in part by Key Laboratory of bridge nondestructive testing and engineering calculation Open fund projects (grant 2020QZJ03).

    The authors declare there is no conflicts of interest.



    [1] Y. S. Song, L. Qi, Properties of some classes of structured tensors, J. Optim. Theory Appl., 165 (2015), 854–873. https://doi.org/10.1007/s10957-014-0616-5 doi: 10.1007/s10957-014-0616-5
    [2] R. W. Cottle, J. S. Pang, R. E. Stone, The Linear Complementarity Problem, Academic Press, Boston, MA, 1992.
    [3] R. Mathias, J. S. Pang, Error bounds for the linear complementarity problem with a P-matrix, Linear Algebra Appl., 132 (1990), 123–136. https://doi.org/10.1016/0024-3795(90)90058-K doi: 10.1016/0024-3795(90)90058-K
    [4] X. Chen, S. Xiang, Perturbation bounds of P-matrix linear complementarity problems, SIAM J. Optim., 18 (2008), 1250–1265. https://doi.org/10.1137/060653019 doi: 10.1137/060653019
    [5] X. Chen, S. Xiang, Computation of error bounds for P-matrix linear complementarity problems, Math. Program., 106 (2006), 513–525. https://doi.org/10.1007/s10107-005-0645-9 doi: 10.1007/s10107-005-0645-9
    [6] Z. Q. Luo, O. L. Mangasarian, J. Ren, M. V. Solodov, New error bounds for the linear complementarity problem, Math. Oper. Res., 19 (1994), 880–892. https://doi.org/10.1287/moor.19.4.880 doi: 10.1287/moor.19.4.880
    [7] X. M. Fang, Z. J. Qiao, Improved error bounds based on α(M) for the linear complementarity problem, Linear Algebra Appl., 589 (2020), 186–200. https://doi.org/10.1016/j.laa.2019.12.009 doi: 10.1016/j.laa.2019.12.009
    [8] S. Q. Du, L. Y. Cui, Y. Y. Chen, Y. M. Wei, Stochastic tensor complementarity problem with discrete distribution, J. Optim. Theory Appl., 192 (2022), 912–929. https://doi.org/10.1007/s10957-021-01997-7 doi: 10.1007/s10957-021-01997-7
    [9] L. Qi, Z. Y. Luo, Tensor Analysis: Spectral Theory and Special Tensors, SIAM, Philadelphia, PA, 2017. Available from: https://epubs.siam.org/doi/pdf/10.1137/1.9781611974751.bm.
    [10] M. L. Che, L. Q. Qi, Y. M. Wei, Positive definite tensors to nonlinear complementarity problems, J. Optim. Theory Appl., 168 (2016), 475–487. https://doi.org/10.1007/s10957-015-0773-1 doi: 10.1007/s10957-015-0773-1
    [11] Z. H. Huang, L. Q. Qi, Tensor complementarity problems–part Ⅰ: Basic theory, J. Optim. Theory Appl., 183 (2019), 1–23. https://doi.org/10.1007/s10957-019-01566-z doi: 10.1007/s10957-019-01566-z
    [12] L. B. Cui, Y. D. Fan, Y. S. Song, S. L. Wu, The existence and uniqueness of solution for tensor complementarity problem and related systems, J. Optim. Theory Appl., 192 (2022), 312–334. https://doi.org/10.1007/s10957-021-01972-2 doi: 10.1007/s10957-021-01972-2
    [13] Z. H. Huang, Y. Y. Suo, J. Wang, On Q-tensors, preprint, arXiv: 1509.03088.
    [14] W. Yu, C. Ling, H. He, On the properties of tensor complementarity problems, preprient, arXiv: 1608.01735.
    [15] M. M. Zheng, Y. Zhang, Z. H. Huang, Global error bounds for the tensor complementarity problem with a P-tensor, J. Ind. Manage. Optim., 15 (2019), 933–946. https://doi.org/10.3934/jimo.2018078 doi: 10.3934/jimo.2018078
    [16] Z. H. Huang, L. Q. Qi, Tensor complementarity problems–part Ⅲ: Applications, J. Optim. Theory Appl., 183 (2019), 771–791. https://doi.org/10.1007/s10957-019-01573-0 doi: 10.1007/s10957-019-01573-0
    [17] X. L. Bai, Z. H. Huang, Y. Wang, Global uniqueness and solvability for tensor complementarity problems, J. Optim. Theory Appl., 170 (2016), 72–84. https://doi.org/10.1007/s10957-016-0903-4 doi: 10.1007/s10957-016-0903-4
    [18] P. Z. Yuan, L. H. You, Some remarks on P, P0, B and B0 tensors, Linear Algebra Appl., 459 (2014), 511–521. https://doi.org/10.1016/j.laa.2014.07.043 doi: 10.1016/j.laa.2014.07.043
    [19] L. Q. Qi, Z. H. Huang, Tensor complementarity problems–part Ⅱ: Solution methods, J. Optim. Theory Appl., 183 (2019), 365–385. https://doi.org/10.1007/s10957-019-01568-x doi: 10.1007/s10957-019-01568-x
  • This article has been cited by:

    1. Xue-liu Li, Tong-tong Shang, Guo-ji Tang, Lower bounds of the solution set of the polynomial complementarity problem, 2024, 18, 1862-4472, 497, 10.1007/s11590-023-02004-w
    2. Xue-liu Li, Guo-ji Tang, Strict feasibility for the polynomial complementarity problem, 2024, 89, 0925-5001, 57, 10.1007/s10898-023-01339-z
    3. Xue-liu Li, Guo-ji Tang, The Bounds of Solutions to Polynomial Complementarity Problems, 2024, 0022-3239, 10.1007/s10957-024-02511-5
    4. Yi-gong Huang, Tong-tong Shang, Guo-ji Tang, Existence results of solutions to a generalized vertical polynomial complementarity problem in terms of vertical block tensor tuples, 2024, 451, 03770427, 116113, 10.1016/j.cam.2024.116113
    5. Hai-Ying Wang, Zu-Feng Fu, Shi-Liang Wu, On the Bound of the Solution Set for the Vertical Tensor Complementarity Problem, 2025, 204, 0022-3239, 10.1007/s10957-024-02559-3
  • Reader Comments
  • © 2022 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(1601) PDF downloads(59) Cited by(5)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog