Liqun Qi

Maolin Che, Liqun Qi, Yimin Wei

The main purpose of this note is to investigate some kinds of nonlinear complementarity problems (NCP). For the structured tensors, such as, symmetric positive definite tensors and copositive tensors, we derive the existence theorems on a solution of these kinds of nonlinear complementarity problems. We prove that a unique solution of the NCP exists under the condition of diagonalizable tensors.

Shenglong Hu, Zheng-Hai Huang, Liqun Qi

In this paper, we introduce a new class of nonnegative tensors --- strictly nonnegative tensors. A weakly irreducible nonnegative tensor is a strictly nonnegative tensor but not vice versa. We show that the spectral radius of a strictly nonnegative tensor is always positive. We give some sufficient and necessary conditions for the six well-conditional classes of nonnegative tensors, introduced in the literature, and a full relationship picture about strictly nonnegative tensors with these six classes of nonnegative tensors. We then establish global R-linear convergence of a power method for finding the spectral radius of a nonnegative tensor under the condition of weak irreducibility. We show that for a nonnegative tensor T, there always exists a partition of the index set such that every tensor induced by the partition is weakly irreducible; and the spectral radius of T can be obtained from those spectral radii of the induced tensors. In this way, we develop a convergent algorithm for finding the spectral radius of a general nonnegative tensor without any additional assumption. The preliminary numerical results demonstrate the feasibility and effectiveness of the proposed algorithm.

Shouqiang Du, Liping Zhang, Chiyu Chen et al.

This paper is concerned with solving some structured multi-linear systems, which are called tensor absolute value equations. This kind of absolute value equations is closely related to tensor complementarity problems and is a generalization of the well-known absolute value equations in the matrix case. We prove that tensor absolute value equations are equivalent to some special structured tensor complementary problems. Some sufficient conditions are given to guarantee the existence of solutions for tensor absolute value equations. We also propose a Levenberg-Marquardt-type algorithm for solving some given tensor absolute value equations and preliminary numerical results are reported to indicate the efficiency of the proposed algorithm.

Liping Zhang, Liqun Qi, Ziyan Luo

It is well known that the dominant eigenvalue of a real essentially nonnegative matrix is a convex function of its diagonal entries. This convexity is of practical importance in population biology, graph theory, demography, analytic hierarchy process and so on. In this paper, the concept of essentially nonnegativity is extended from matrices to higher order tensors, and the convexity and log convexity of dominant eigenvalues for such a class of tensors are established. Particularly, for any nonnegative tensor, the spectral radius turns out to be the dominant eigenvalue and hence possesses these convexities. Finally, an algorithm is given to calculate the dominant eigenvalue, and numerical results are reported to show the effectiveness of the proposed algorithm.

Shenglong Hu, Liqun Qi

In 1907, Oskar Perron showed that a positive square matrix has a unique largest positive eigenvalue with a positive eigenvector. This result was extended to irreducible nonnegative matrices by Geog Frobenius in 1912, and to irreducible nonnegative tensors and weakly irreducible nonnegative tensors recently. This result is a fundamental result in matrix theory and has found wide applications in probability theory, internet search engines, spectral graph and hypergraph theory, etc. In this paper, we give a necessary and sufficient condition for the existence of such a positive eigenvector, i.e., a positive Perron vector, for a nonnegative tensor. We show that every nonnegative tensor has a canonical nonnegative partition form, from which we introduce strongly nonnegative tensors. A tensor is called strongly nonnegative, if the spectral radius of each genuine weakly irreducible block is equal to the spectral radius of the tensor, which is strictly larger than the spectral radius of any other block. We prove that a nonnegative tensor has a positive Perron vector if and only if it is strongly nonnegative. The proof is nontrivial. Numerical results for finding a positive Perron vector are reported.

Liping Zhang, Liqun Qi, Guanglu Zhou

We study M-tensors and various properties of M-tensors are given. Specially, we show that the smallest real eigenvalue of M-tensor is positive corresponding to a nonnegative eigenvector. We propose an algorithm to find the smallest positive eigenvalue and then apply the property to study the positive definiteness of a multivariate form.

Z. Luo, L. Qi, Ph. L. Toint

A generalization of the Bernstein matrix concentration inequality to random tensors of general order is proposed. This generalization is based on the use of Einstein products between tensors, from which a strong link can be established between matrices and tensors, in turn allowing exploitation of existing results for the former.

Runyao Yu, Chenhui Gu, Jochen Stiasny et al.

Electricity price forecasting in Europe presents unique challenges due to increasing renewable generation variability, market integration, and the continent's physically interconnected power system. While recent advances in foundation models have led to substantial improvements in general time series forecasting, most existing approaches do not incorporate prior graph knowledge from the transmission topology, which can limit their ability to exploit meaningful cross-region dependencies in interconnected power systems, motivating a domain-specific foundation model. In this paper, we address this gap by first introducing a comprehensive and up-to-date dataset across 24 European countries (38 regions), spanning from 2022-01-01 to 2026-01-01. Building on this groundwork, we propose PriceFM, a probabilistic foundation model pretrained on this large dataset. Specifically, PriceFM maps each region's price and exogenous features, including load, solar, and wind generation forecasts, into a comparable latent embedding via a shared Mixture-of-Experts (MoE) projection layer, then injects prior graph knowledge by constructing a sparse graph mask derived from transmission topology. Across a large-scale European benchmark, PriceFM achieves strong performance and demonstrates superior generalization compared with multiple competitive baselines. The results highlight the value of topology-guided forecasting with increasing renewable generation and strong cross-border interconnections. The methodology is available at: https://runyao-yu.github.io/PriceFM/.

Yannan Chen, Liqun Qi, Qun Wang

In this note, we construct explicit SOS decomposition of A Fourth Order Four Dimensional Hankel Tensor with A Symmetric Generating Vector, at the critical value. This is a supplementary note to Paper [3].

Gaohang Yu, Shaochun Wan, Liqun Qi et al.

Tensor train (TT) factorization and corresponding TT rank, which can well express the low-rankness and mode correlations of higher-order tensors, have attracted much attention in recent years. However, TT factorization based methods are generally not sufficient to characterize low-rankness along each mode of third-order tensor. Inspired by this, we generalize the tensor train factorization to the mode-k tensor train factorization and introduce a corresponding multi-mode tensor train (MTT) rank. Then, we proposed a novel low-MTT-rank tensor completion model via multi-mode TT factorization and spatial-spectral smoothness regularization. To tackle the proposed model, we develop an efficient proximal alternating minimization (PAM) algorithm. Extensive numerical experiment results on visual data demonstrate that the proposed MTTD3R method outperforms compared methods in terms of visual and quantitative measures.

Chenjian Pan, Chen Ling, Hongjin He et al.

Recovering color images and videos from highly undersampled data is a fundamental and challenging task in face recognition and computer vision. By the multi-dimensional nature of color images and videos, in this paper, we propose a novel tensor completion approach, which is able to efficiently explore the sparsity of tensor data under the discrete cosine transform (DCT). Specifically, we introduce two ``sparse + low-rank'' tensor completion models as well as two implementable algorithms for finding their solutions. The first one is a DCT-based sparse plus weighted nuclear norm induced low-rank minimization model. The second one is a DCT-based sparse plus $p$-shrinking mapping induced low-rank optimization model. Moreover, we accordingly propose two implementable augmented Lagrangian-based algorithms for solving the underlying optimization models. A series of numerical experiments including color image inpainting and video data recovery demonstrate that our proposed approach performs better than many existing state-of-the-art tensor completion methods, especially for the case when the ratio of missing data is high.

Chenjian Pan, Chen Ling, Hongjin He et al.

Tensor completion refers to the task of estimating the missing data from an incomplete measurement or observation, which is a core problem frequently arising from the areas of big data analysis, computer vision, and network engineering. Due to the multidimensional nature of high-order tensors, the matrix approaches, e.g., matrix factorization and direct matricization of tensors, are often not ideal for tensor completion and recovery. In this paper, we introduce a unified low-rank and sparse enhanced Tucker decomposition model for tensor completion. Our model possesses a sparse regularization term to promote a sparse core tensor of the Tucker decomposition, which is beneficial for tensor data compression. Moreover, we enforce low-rank regularization terms on factor matrices of the Tucker decomposition for inducing the low-rankness of the tensor with a cheap computational cost. Numerically, we propose a customized ADMM with enough easy subproblems to solve the underlying model. It is remarkable that our model is able to deal with different types of real-world data sets, since it exploits the potential periodicity and inherent correlation properties appeared in tensors. A series of computational experiments on real-world data sets, including internet traffic data sets, color images, and face recognition, demonstrate that our model performs better than many existing state-of-the-art matricization and tensorization approaches in terms of achieving higher recovery accuracy.

Xuezhong Wang, Maolin Che, Liqun Qi et al.

Nonlinear gradient dynamic approach for solving the tensor complementarity problem (TCP) is presented. Theoretical analysis shows that each of the defined dynamical system models ensures the convergence performance. The computer simulation results further substantiate that the considered dynamical system can solve the tensor complementarity problem (TCP).

Zhenghai Huang, Liqun Qi

In this paper, we consider higher order paired symmetric tensors and strongly paired symmetric tensors. Elasticity tensors and higher order elasticity tensors in solid mechanics are strongly paired symmetric tensors. A (strongly) paired symmetric tensor is said to be positive definite if the homogeneous polynomial defined by it is positive definite. Positive definiteness of elasticity and higher order elasticity tensors is strong ellipticity in solid mechanics, which plays an important role in nonlinear elasticity theory. We mainly investigate positive definiteness of fourth order three dimensional and sixth order three dimensional (strongly) paired symmetric tensors. We first show that the concerned (strongly) paired symmetric tensor is positive definite if and only if its smallest $M$-eigenvalue is positive. Second, we propose several necessary and sufficient conditions under which the concerned (strongly) paired symmetric tensor is positive definite. Third, we study the conditions under which the homogeneous polynomial defined by a fourth order three dimensional or sixth order three dimensional (strongly) paired symmetric tensor can be written as a sum of squares of polynomials, and further, propose several necessary and/or sufficient conditions to judge whether the concerned (strongly) paired symmetric tensors are positive definite or not. Fourth, by using semidefinite relaxation we propose a sequential semidefinite programming method to compute the smallest $M$-eigenvalue of a fourth order three dimensional (strongly) paired symmetric tensor, by which we can check positive definiteness of the concerned tensor. The preliminary numerical results demonstrate that our method is effective.

Weiyang Ding, Liqun Qi, Hong Yan

We establish several sufficient conditions for the strong ellipticity of any fourth-order elasticity tensor in this paper. The first presented sufficient condition is an extension of positive definite matrices, which states that the strong ellipticity holds if the unfolding matrix of this fourth-order elasticity tensor can be modified into a positive definite one by preserving the summations of some corresponding entries. An alternating projection algorithm is proposed to verify whether an elasticity tensor satisfies the first condition or not. Conditions for some special cases beyond the first sufficient condition are further investigated, which includes some important cases for the isotropic and some particular anisotropic linearly elastic materials.

Yannan Chen, Liqun Qi, Qun Wang

Large scale tensors, including large scale Hankel tensors, have many applications in science and engineering. In this paper, we propose an inexact curvilinear search optimization method to compute Z- and H-eigenvalues of $m$th order $n$ dimensional Hankel tensors, where $n$ is large. Owing to the fast Fourier transform, the computational cost of each iteration of the new method is about $\mathcal{O}(mn\log(mn))$. Using the Cayley transform, we obtain an effective curvilinear search scheme. Then, we show that every limiting point of iterates generated by the new algorithm is an eigen-pair of Hankel tensors. Without the assumption of a second-order sufficient condition, we analyze the linear convergence rate of iterate sequence by the Kurdyka-Łojasiewicz property. Finally, numerical experiments for Hankel tensors, whose dimension may up to one million, are reported to show the efficiency of the proposed curvilinear search method.