Search

Theory of quasi-orthogonal polynomials is significantly related to constructing vector Padé approximations. The present paper introduces an efficient procedure to compute adjacent families of quasi-orthogonal polynomials as defined in Sadaka (Appl. Numer. Math. 21:57–70, 1996) and Sadaka (Appl. Numer. Math. 24:483–499, 1997). The derived algorithm uses short recursive relations whose coefficients are written in terms of lower triangular block determinants. The strategy of computing such determinants is given and analyzed. [ABSTRACT FROM AUTHOR]

We study the issue of numerically solving the nonnegative inverse eigenvalue problem (NIEP). At first, we reformulate the NIEP as a convex composite optimization problem on Riemannian manifolds. Then we develop a scheme of the Riemannian linearized proximal algorithm (R-LPA) to solve the NIEP. Under some mild conditions, the local and global convergence results of the R-LPA for the NIEP are established, respectively. Moreover, numerical experiments are presented. Compared with the Riemannian Newton-CG method in Z. Zhao et al. (Numer. Math. 140:827–855, 2018), this R-LPA owns better numerical performances for large scale problems and sparse matrix cases, which is due to the smaller dimension of the Riemannian manifold derived from the problem formulation of the NIEP as a convex composite optimization problem. [ABSTRACT FROM AUTHOR]