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.

Key words： 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, https://doi.org/10.1007/s10589-017-9938-1

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