On computing minimal H-eigenvalue of sign-structured tensors

Abstract

Finding the minimal H-eigenvalue of tensors is an important topic in tensor computation and numerical multilinear algebra. This paper is devoted to a sum-of-squares (SOS) algorithm for computing the minimal H-eigenvalues of tensors with some sign structures called extended essentially nonnegative tensors (EEN-tensors), which includes nonnegative tensors as a subclass. In the even-order symmetric case, we first discuss the positive semi-definiteness of EEN-tensors, and show that a positive semi-definite EEN-tensor is a nonnegative tensor or an M-tensor or the sum of a nonnegative tensor and an M-tensor, then we establish a checkable sufficient condition for the SOS decomposition of EEN-tensors. Finally, we present an efficient algorithm to compute the minimal H-eigenvalues of even-order symmetric EEN-tensors based on the SOS decomposition. Numerical experiments are given to show the efficiency of the proposed algorithm.

Keywords

Extended essentially nonnegative tensor (EEN-tensor) / positive semi-definiteness / H-eigenvalue / sum-of-squares (SOS) polynomial

References

[1] BarmpoutisA, VemuriB C, HowlandD, ForderJ R. Regularized positive-definite fourth order tensor field estimation from DW-MRI.Neuroimage, 2009, 45: 153–162

[2] BasserP J, JonesD K. Diffusion-tensor MRI: theory, experimental design and data analysis a technical review.NMR Biomed, 2002, 15: 456–467

[3] ChenH, HuangZ, QiL. Copositivity detection of tensors: theory and algorithm.J Optim Theory Appl, 2017, 174: 746–761

[4] ChenH, HuangZ, QiL. Copositive tensor detection and its applications in physics and hypergraphs.Comput Optim Appl,2017,

[5] ChenH, LiG, QiL. SOS tensor decomposition: theory and applications.Commun Math Sci, 2016, 14(8): 2073–2100

[6] ChenY, QiL, WangQ. Computing extreme eigenvalues of large scale Hankel tensors.J Sci Comput, 2016, 68(2): 716–738

[7] CooperJ, DutleA. Spectra of uniform hypergraphs.Linear Algebra Appl,2012, 436: 3268–3292

[8] CuiC, DaiY, NieJ. All real eigenvalues of symmetric tensors.SIAM J Matrix Anal Appl, 2014, 35(4): 1582–1601

[9] DingW, QiL, WeiY.M-Tensors and nonsingular M-tensors.Linear Algebra Appl, 2013, 439: 3264–3278

[10] DingW, QiL, WeiY. Fast Hankel tensor-vector product and its application to exponential data fitting.Numer Linear Algebra Appl,2015, 22: 814–832

[11] HuS, LiG, QiL. A tensor analogy of Yuan’s alternative theorem and polynomial optimization with sign structure.J Optim Theory Appl, 2016, 168: 446–474

[12] HuS, LiG, QiL, SongY. Finding the maximum eigenvalue of essentially nonnegative symmetric tensors via sum of squares programming.J Optim Theory Appl,2013, 158: 717–738

[13] HuS, QiL, XieJ. The largest Laplacian and signless Laplacian H-eigenvalues of a uniform hypergraph.Linear Algebra Appl,2015, 469: 1–27

[14] LiuY, ZhouG, IbrahimN F. An always convergent algorithm for the largest eigenvalue of an irreducible nonnegative tensor.J Comput Appl Math, 2010, 235(1): 286–292

[15] LöfbergJ. YALMIP: a toolbox for modeling and optimization in MATLAB.In: IEEE International Symposium on Computer Aided Control Systems Design, Taipei, September24, 2004. 2004, 284–289

[16] LöfbergJ. Pre- and post-processing sums-of-squares programs in practice.IEEE Trans Automat Control, 2009, 54: 1007–1011

[17] NgM, QiL, ZhouG. Finding the largest eigenvalue of a non-negative tensor.SIAM J Matrix Anal Appl, 2009, 31: 1090–1099

[18] QiL. Eigenvalue of a real supersymmetric tensor.J Sci Comput, 2005, 40: 1302–1324

[19] QiL. H+-eigenvalues of Laplacian and signless Laplacian tensors.Commun Math Sci,2014, 12: 1045–1064

[20] QiL, WangF, WangY. Z-eigenvalue methods for a global polynomial optimization problem.Math Program, 2009, 118: 301–316

[21] QiL, YuG,WuE X. Higher order positive semidefinite diffusion tensor imaging.SIAM J Imaging Sci, 2010,3: 416–C 433

[22] SturmJ F. SeDuMi 1.02: a MATLAB toolbox for optimization over symmetric cones.Optim Methods Softw, 1999, 11: 625–653

[23] WangY, CaccettaL, ZhouG. Convergence analysis of a block improvement method for polynomial optimization over unit spheres.Numer Linear Algebra Appl,2015, 22: 1059–1076

[24] WangY, LiuW, CaccettaL, ZhouG. Parameter selection for nonnegative l₁ matrix/tensor sparse decomposition.Oper Res Lett, 2015, 43: 423–426

[25] WangY, QiL, ZhangX. A practical method for computing the largest M-eigenvalue of a fourth-order partially symmetric tensor.Numer Linear Algebra Appl,2009, 16: 589–601

[26] WangY,ZhangK, SunH. Criteria for strong H-tensors.Front Math China, 2016, 11: 577–592

[27] WangY, ZhouG, CaccettaL. Nonsingular H-tensor and its criteria.J Ind Manag Optim, 2016, 12: 1173–1186

[28] ZhangK, WangY. An H-tensor based iterative scheme for identifying the positive definiteness of multivariate homogeneous forms.J Comput Appl Math, 2016, 305: 1–10

[29] ZhangL, QiL, ZhouG. M-tensors and some applications.SIAM J Matrix Anal Appl, 2014, 35: 437–452

RIGHTS & PERMISSIONS

Higher Education Press and Springer-Verlag GmbH Germany