A popular and classical method for finding the best rank one approximation of a real tensor is the higher order power method (HOPM). It is known in the literature that the iterative sequence generated by HOPM converges globally, while the convergence rate can be superlinear, linear or sublinear. In this paper, we examine the convergence rate of HOPM in solving the best rank one approximation problem of real tensors. We first show that the iterative sequence of HOPM always converges globally and provide an explicit eventual sublinear convergence rate. The sublinear convergence rate estimate is in terms of the dimension and the order of the underlying tensor space. Then, we examine the concept of nondegenerate singular vector tuples and show that, if the sequence of HOPM converges to a nondegenerate singular vector tuple, then the global convergence rate is R-linear. We show that, for almost all tensors (in the sense of Lebesgue measure), all the singular vector tuples are nondegenerate, and so, the HOPM “typically” exhibits global R-linear convergence rate. Moreover, without any regularity assumption, we establish that the sequence generated by HOPM always converges globally and R-linearly for orthogonally decomposable tensors with order at least 3. We achieved this by showing that each nonzero singular vector tuple of an orthogonally decomposable tensor with order at least 3 is nondegenerate.

In this paper we derive new fractional error bounds for polynomial systems with exponents explicitly determined by the dimension of the underlying space and the number/degree of the involved polynomials. Our major result extends the existing error bounds from the system involving only a single polynomial to a general polynomial system and do not require any regularity assumptions. In this way we resolve, in particular, some open questions posed in the literature. The developed techniques are largely based on variational analysis and generalized differentiation, which allow us to establish, e.g., a nonsmooth extension of the seminal Łojasiewicz’s gradient inequality to maxima of polynomials with explicitly determined exponents. Our major applications concern quantitative Hölderian stability of solution maps for parameterized polynomial optimization problems and nonlinear complementarity systems with polynomial data as well as high-order semismooth properties of the eigenvalues of symmetric tensors.

Image segmentation is an important task in the domain of computer vision and medical imaging. In natural and medical images, intensity inhomogeneity, i.e. the varying image intensity, occurs often and it poses considerable challenges for image segmentation. In this paper, we propose an efficient variational method for segmenting images with intensity inhomogeneity. The method is inspired by previous works on two-stage segmentation and variational Retinex. Our method consists of two stages. In the first stage, we decouple the image into reflection and illumination parts by solving a convex energy minimization model with either total variation or tight-frame regularisation. In the second stage, we segment the original image by thresholding on the reflection part, and the inhomogeneous intensity is estimated by the smoothly varying illumination part. We adopt a primal dual algorithm to solve the convex model in the first stage, and the convergence is guaranteed. Numerical experiments clearly show that our method is robust and efficient to segment both natural and medical images.

In this paper, we present L2 and negative-order norm estimates for the local discontinuous Galerkin (LDG) method solving variable coefficient Schrödinger equations. For these special solutions the LDG method exhibits ‘‘hidden accuracy", and we are able to extract it through the use of a convolution kernel that is composed of a linear combination of B-splines. This technical was initially introduced by Cockburn, Luskin, Shu, and Süli for linear hyperbolic equations and extended by Ryan et al. as a smoothness-increasing accuracy-conserving (SIAC) filter. We demonstrate that it is possible to extend the SIAC filter on Schrödinger equations. When polynomials of degree k are used, we can prove theoretically the LDG method solutions are of order k + 1, whereas the post-processed solutions that convolution with the SIAC filter are of order at least 2k. Additionally, we present numerical results to confirm that the accuracy of LDG solutions can be improved from O(hk+1) to at least O(h2k+1) for Schrödinger equations by using alternating numerical fluxes.

In this paper, we develop the Hamiltonian conservative and L2 conservative local discontinuous Galerkin (LDG) schemes for the Korteweg-de Vries (KdV) type equations with the minimal stencil. For the time discretization, we adopt the semi- implicit spectral deferred correction (SDC) method to achieve the high order accuracy and efficiency. Also we compare the schemes with the dissipative LDG scheme. Stabil- ity of the fully discrete schemes is provided by Fourier analysis for the linearized KdV equation. Numerical examples are shown to illustrate the capability of these schemes. Compared with the dissipative LDG scheme, the numerical simulations also indicate that the conservative LDG scheme with high order time discretization can reduce the long time phase error validly.