Factorizations of EP operators

In this paper we characterize EP operators through the existence of different types of factorizations. Our results extend to EP operators existing characterizations for EP matrices and give new characterizations both for EP matrices and EP operators.     

Fast inversion algorithms for diagonal plus semiseparable matrices

Here are considered matrices represented as a sum of diagonal and semiseparable ones. These matrices belong to the class of structured matrices which arises in numerous applications. FastO(N) algorithms for their inversion were developed before under additional restrictions which are a source of instability. Our aim is to eliminate these restrictions and to develop reliable and stable numerical algorithms. In this paper we obtain such algorithms with the only requirement that the considered matrix is invertible and its determinant is not close to zero. The case of semiseparable matrices of order one was considered in detail in an earlier paper of the authors.     

Inherited LU-factorizations of matrices

Various types of LU-factorizations for nonsingular matrices, where L is a lower triangular matrix and U is an upper triangular matrix, are defined and characterized. These types of LU-factorizations are extended to the general m × n case. The more general conditions are considered in the light of the structures of [C.R. Johnson, D.D. Olesky, P. Van den Driessche, Inherited matrix entries: LU factorizations, SIAM J. Matrix Anal. Appl. 10 (1989) 99-104]. Applications to graphs and adjacency matrices are investigated. Conditions for the product of a lower and an upper triangular matrix to be the zero matrix are also obtained.     

Fast state space algorithms for matrix Nehari and Nehari-Takagi interpolation problems

Numerical algorithms with complexityO(n ²) operations are proposed for solving matrix Nehari and Nehari-Takagi problems withn interpolation points. The algorithms are based on explicit formulas for the solutions and on theorems about cascade decomposition of rational matrix function given in a state space form. The method suggests also fast algorithms for LDU factorizations of structured matrices. The numerical behavior of the designed algorithms is studied for a wide set of examples.     

Generalized Leibniz functional matrices and factorizations of some well-known matrices

In this paper, we introduce the generalized Leibniz functional matrices and study some algebraic properties of such matrices. To demonstrate applications of these properties, we derive several novel factorization forms of some well-known matrices, such as the complete symmetric polynomial matrix and the elementary symmetric polynomial matrix. In addition, by applying factorizations of the generalized Leibniz functional matrices, we redevelop the known results of factorizations of Stirling matrices of the first and second kind and the generalized Pascal matrix.     

Polynomial perturbations of bilinear functionals and Hessenberg matrices

This paper deals with symmetric and non-symmetric polynomial perturbations of symmetric quasi-definite bilinear functionals. We establish a relation between the Hessenberg matrices associated with the initial and the perturbed functionals using LU and QR factorizations. Moreover we give an explicit algebraic relation between the sequences of orthogonal polynomials associated with both functionals.     

Multiple LU factorizations of a singular matrix

A singular matrix A may have more than one LU factorizations. In this work the set of all LU factorizations of A is explicitly described when the lower triangular matrix L is nonsingular. To this purpose, a canonical form of A under left multiplication by unit lower triangular matrices is introduced. This canonical form allows us to characterize the matrices that have an LU factorization and to parametrize all possible LU factorizations. Formulae in terms of quotient of minors of A are presented for the entries of this canonical form.     

On the number of P-vertices of some graphs

We provide positive answers to some open questions presented recently by Kim and Shader on a continuity-like property of the P-vertices of nonsingular matrices whose graph is a path. A criterion for matrices associated with more general trees to have at most n − 1 P-vertices is established. The cases of the cycles and stars are also analyzed. Several algorithms for generating matrices with a given number of P-vertices are proposed.     

Matrices related to the Bell polynomials

In this paper, we study the matrices related to the partial exponential Bell polynomials and those related to the Bell polynomials with respect to Ω. As a result, the factorizations of these matrices are obtained, which give unified approaches to the factorizations of many lower triangular matrices. Moreover, some combinatorial identities are also derived from the corresponding matrix representations.     

Numerical stability of fast cosine transforms

In this paper we consider the numerical stability of fast algorithms for discrete cosine transform (DCT) of type III and II, respectively. We show that various fast DCTs can possess a very different behaviour of numerical stability. By matrix factorizations we find that a complex fast DCT which is based mainly on a fast Fouier transform has a better numerical stability than a real fast DCT despite its larger arithmetical complexity. Numerical tests illustrate our theoretical results.     

Functions of a matrix and Krylov matrices

For a given nonderogatory matrix A, formulas are given for functions of A in terms of Krylov matrices of A. Relations between the coefficients of a polynomial of A and the generating vector of a Krylov matrix of A are provided. With the formulas, linear transformations between Krylov matrices and functions of A are introduced, and associated algebraic properties are derived. Hessenberg reduction forms are revisited equipped with appropriate inner products and related properties and matrix factorizations are given.     

Bijections and symmetries for the factorizations of the long cycle

We study the factorizations of the permutation (1,2,…,n)$(1,2,\dots ,n)$ into k factors of given cycle types. Using representation theory, Jackson obtained for each k   an elegant formula for counting these factorizations according to the number of cycles of each factor. In the cases k=2,3$k=2,3$ Schaeffer and Vassilieva gave a combinatorial proof of Jackson?s formula, and Morales and Vassilieva obtained more refined formulas exhibiting a surprising symmetry property. These counting results are indicative of a rich combinatorial theory which has remained elusive to this point, and it is the goal of this article to establish a series of bijections which unveil some of the combinatorial properties of the factorizations of (1,2,…,n)$(1,2,\dots ,n)$ into k factors for all k  . We thereby obtain refinements of Jackson?s formulas which extend the cases k=2,3$k=2,3$ treated by Morales and Vassilieva. Our bijections are described in terms of “constellations”, which are graphs embedded in surfaces encoding the transitive factorizations of permutations.     

Schröder matrix as inverse of Delannoy matrix

Using Riordan arrays, we introduce a generalized Delannoy matrix by weighted Delannoy numbers. It turns out that Delannoy matrix, Pascal matrix, and Fibonacci matrix are all special cases of the generalized Delannoy matrices, meanwhile Schröder matrix and Catalan matrix also arise in involving inverses of the generalized Delannoy matrices. These connections are the focus of our paper. The half of generalized Delannoy matrix is also considered. In addition, we obtain a combinatorial interpretation for the generalized Fibonacci numbers.     

On a class of algebras associated to directed graphs

To any directed graph we associate an algebra with edges of the graph as generators and with relations defined by all pairs of directed paths with the same origin and terminus. Such algebras are related to factorizations of polynomials over noncommutative algebras. We also construct a basis for our algebras associated to layered graphs.     

Accurate computations with Said-Ball-Vandermonde matrices

A generalization of the Vandermonde matrices which arise when the power basis is replaced by the Said-Ball basis is considered. When the nodes are inside the interval (0,1), then those matrices are strictly totally positive. An algorithm for computing the bidiagonal decomposition of those Said-Ball-Vandermonde matrices is presented, which allows us to use known algorithms for totally positive matrices represented by their bidiagonal decomposition. The algorithm is shown to be fast and to guarantee high relative accuracy. Some numerical experiments which illustrate the good behaviour of the algorithm are included.     

Natural density distribution of Hermite normal forms of integer matrices

The Hermite Normal Form (HNF) is a canonical representation of matrices over any principal ideal domain. Over the integers, the distribution of the HNFs of randomly looking matrices is far from uniform. The aim of this article is to present an explicit computation of this distribution together with some applications. More precisely, for integer matrices whose entries are upper bounded in absolute value by a large bound, we compute the asymptotic number of such matrices whose HNF has a prescribed diagonal structure. We apply these results to the analysis of some procedures and algorithms whose dynamics depend on the HNF of randomly looking integer matrices.     

Bases for kernel-based spaces

Since it is well-known (De Marchi and Schaback (2001) [4]) that standard bases of kernel translates are badly conditioned while the interpolation itself is not unstable in function space, this paper surveys the choices of other bases. All data-dependent bases turn out to be defined via a factorization of the kernel matrix defined by these data, and a discussion of various matrix factorizations (e.g. Cholesky, QR, SVD) provides a variety of different bases with different properties. Special attention is given to duality, stability, orthogonality, adaptivity, and computational efficiency. The “Newton” basis arising from a pivoted Cholesky factorization turns out to be stable and computationally cheap while being orthonormal in the “native” Hilbert space of the kernel. Efficient adaptive algorithms for calculating the Newton basis along the lines of orthogonal matching pursuit conclude the paper.     

On unimodular graphs

We study graphs whose adjacency matrices have determinant equal to 1 or −1, and characterize certain subclasses of these graphs. Graphs whose adjacency matrices are totally unimodular are also characterized. For bipartite graphs having a unique perfect matching, we provide a formula for the inverse of the corresponding adjacency matrix, and address the problem of when that inverse is diagonally similar to a nonnegative matrix. Special attention is paid to the case that such a graph is unicyclic.