Theorems on Schur complement of block diagonally dominant matrices and their application in reducing the order for the solution of large scale linear systems

We firstly consider the block dominant degree for I-(II-)block strictly diagonally dominant matrix and their Schur complements, showing that the block dominant degree for the Schur complement of an I-(II-)block strictly diagonally dominant matrix is greater than that of the original grand block matrix. Then, as application, we present some disc theorems and some bounds for the eigenvalues of the Schur complement by the elements of the original matrix. Further, by means of matrix partition and the Schur complement of block matrix, based on the derived disc theorems, we give a kind of iteration called the Schur-based iteration, which can solve large scale linear systems though reducing the order by the Schur complement and the numerical example illustrates that the iteration can compute out the results faster.     

Absolute and relative Weyl theorems for generalized eigenvalue problems

Weyl-type eigenvalue perturbation theories are derived for Hermitian definite pencils A-λB, in which B is positive definite. The results provide a one-to-one correspondence between the original and perturbed eigenvalues, and give a uniform perturbation bound. We give both absolute and relative perturbation results, defined in the standard Euclidean metric instead of the chordal metric that is often used.     

On Padé-type model order reduction of J-Hermitian linear dynamical systems

A simple, yet powerful approach to model order reduction of large-scale linear dynamical systems is to employ projection onto block Krylov subspaces. The transfer functions of the resulting reduced-order models of such projection methods can be characterized as Padé-type approximants of the transfer function of the original large-scale system. If the original system exhibits certain symmetries, then the reduced-order models are considerably more accurate than the theory for general systems predicts. In this paper, the framework of J-Hermitian linear dynamical systems is used to establish a general result about this higher accuracy. In particular, it is shown that in the case of J-Hermitian linear dynamical systems, the reduced-order transfer functions match twice as many Taylor coefficients of the original transfer function as in the general case. An application to the SPRIM algorithm for order reduction of general RCL electrical networks is discussed.     

Structured eigenvalue condition numbers and linearizations for matrix polynomials

This work is concerned with eigenvalue problems for structured matrix polynomials, including complex symmetric, Hermitian, even, odd, palindromic, and anti-palindromic matrix polynomials. Most numerical approaches to solving such eigenvalue problems proceed by linearizing the matrix polynomial into a matrix pencil of larger size. Recently, linearizations have been classified for which the pencil reflects the structure of the original polynomial. A question of practical importance is whether this process of linearization significantly increases the eigenvalue sensitivity with respect to structured perturbations. For all structures under consideration, we show that this cannot happen if the matrix polynomial is well scaled: there is always a structured linearization for which the structured eigenvalue condition number does not differ much. This implies, for example, that a structure-preserving algorithm applied to the linearization fully benefits from a potentially low structured eigenvalue condition number of the original matrix polynomial.     

On a matrix decomposition of Hartwig and Spindelböck

The main task of the paper is to demonstrate that Corollary 6 in [R.E. Hartwig, K. Spindelböck, Matrices for which A^∗ and A^† commute, Linear and Multilinear Algebra 14 (1984) 241-256] provides a powerful tool to investigate square matrices with complex entries. This aim is achieved, on the one hand, by obtaining several original results involving square matrices, and, on the other hand, by reestablishing some of the facts already known in the literature, often in extended and/or generalized forms. The particular attention is paid to the usefulness of the aforementioned corollary to characterize various classes of matrices and to explore matrix partial orderings.     

Generalized even and odd totally positive matrices

A generalization of the definition of an oscillatory matrix based on the theory of cones is given in this paper. The positivity and simplicity of all the eigenvalues of a generalized oscillatory matrix are proved. Classes of generalized even and odd oscillatory matrices are introduced. Spectral properties of the obtained matrices are studied. Criteria of generalized even and odd oscillation are given. Examples of generalized even and odd oscillatory matrices are presented.     

Hyponormal and strongly hyponormal matrices in inner product spaces

Complex matrices that are structured with respect to a possibly degenerate indefinite inner product are studied. Based on earlier works on normal matrices, the notions of hyponormal and strongly hyponormal matrices are introduced. A full characterization of such matrices is given and it is shown how those matrices are related to different concepts of normal matrices in degenerate inner product spaces. Finally, the existence of invariant semidefinite subspaces for strongly hyponormal matrices is discussed.     

Arithmetic-Geometric Mean determinantal identity

In this paper, we give a generalization of a determinantal identity posed by Charles R. Johnson about minors of a Toeplitz matrix satisfying a specific matrix identity. These minors are those appear in the Dodgson’s condensation formula.     

Permanents of doubly stochastic matrices with fixed zero pattern

The doubly stochastic matrices with a given zero pattern which are closest in Euclidean norm to J_nn, the matrix with each entry equal to 1/n, are identified. If the permanent is restricted to matrices having a given zero pattern confined to one row or to one column, the permanent achieves a local minimum at those matrices with that zero pattern which are closest to J_nn. This need no longer be true if the zeros lie in more than one row or column.     

Maps on matrices that preserve the spectrum

Let T be a continuous map of the space of complex n×n matrices into itself satisfying T(0)=0 such that the spectrum of T(x)-T(y) is always a subset of the spectrum of x-y. There exists then an invertible n×n matrix u such that either T(a)=uau^-1 for all a or T(a)=ua^tu^-1 for all a. We arrive at the same conclusion by supposing that the spectrum of x-y is always a subset of the spectrum of T(x)-Tt(y), without the continuity assumption on T.     

Max algebraic powers of irreducible matrices in the periodic regime: An application of cyclic classes

In max algebra it is well known that the sequence of max algebraic powers A^k, with A an irreducible square matrix, becomes periodic after a finite transient time T(A), and the ultimate period γ is equal to the cyclicity of the critical graph of A.In this connection, we study computational complexity of the following problems: (1) for a given k, compute a periodic power A^r with and r?T(A), (2) for a given x, find the ultimate period of {A^l⊗x}. We show that both problems can be solved by matrix squaring in O(n³logn) operations. The main idea is to apply an appropriate diagonal similarity scaling A?X^-1AX, called visualization scaling, and to study the role of cyclic classes of the critical graph.     

Some results on symplectic matrices

The principal results are that if A is an integral matrix such that AA^T is symplectic then A = CQ, where Q is a permutation matrix and C is symplectic; and that if A is a hermitian positive definite matrix which is symplectic, and B is the unique hermitian positive definite pth.root of A, where p is a positive integer, then B is also symplectic.     

Determinant and Pfaffian of sum of skew symmetric matrices

We completely describe the determinants of the sum of orbits of two real skew symmetric matrices, under similarity action of orthogonal group and the special orthogonal group respectively. We also study the Pfaffian case and the complex case.     

Some results on symplectic matrices

The principal results are that if A is an integral matrix such that AA^T is symplectic then A = CQ, where Q is a permutation matrix and C is symplectic; and that if A is a hermitian positive definite matrix which is symplectic, and B is the unique hermitian positive definite pth.root of A, where p is a positive integer, then B is also symplectic.     

The Center Map of an Affine Immersion

We study the center map of an equiaffine immersion which is introduced using the equiaffine support function. The center map is a constant map if and only if the hypersurface is an equiaffine sphere. We investigate those immersions for which the center map is affine congruent with the original hypersurface. In terms of centroaffine geometry, we show that such hypersurfaces provide examples of hypersurfaces with vanishing centroaffine Tchebychev operator. We also characterize them in equiaffine differential geometry using a curvature condition involving the covariant derivative of the shape operator. From both approaches, assuming the dimension is 2 and the surface is definite, a complete classification follows. Received: May 24, 2006. Revised: July 26, 2006. Accepted: July 28, 2006.     

A note on matrix equivalence

Let R be a principal ideal ringR_n the ring of n × n it matrices over R. It is shown that if A, B, X, Y are elements of R* such that A = XB, B = YA, then A and B are left equivalent. Some consequences are given.     

The ratio field of values

The ratio field of values, a generalization of the classical field of values to a pair of n-by-n matrices, is defined and studied, primarily from a geometric point of view. Basic functional properties of the ratio field are developed and used. A decomposition of the ratio field into line segments and ellipses along a master curve is given and this allows computation. Primary theoretical results include the following. It is shown (1) for which denominator matrices the ratio field is always convex, (2) certain other cases of convex pairs are given, and (3) that, at least for n=2, the ratio field obeys a near convexity property that the intersection with any line segment has at most n components. Generalizations of the ratio field of values involving more than one matrix in both the numerator and denominator are also investigated. It is shown that generally such extensions need not be convex or even simply connected.     

Computation of eigenpair partial derivatives by Rayleigh–Ritz procedure

Based on Rayleigh–Ritz procedure, a new method is proposed for a few eigenpair partial derivatives of large matrices. This method simultaneously computes the approximate eigenpairs and their partial derivatives. The linear systems of equations that are solved for eigenvector partial derivatives are greatly reduced from the original matrix size. And the left eigenvectors are not required. Moreover, errors of the computed eigenpairs and their partial derivatives are investigated. Hausdorff distance and containment gap are used to measure the accuracy of approximate eigenpair partial derivatives. Error bounds on the computed eigenpairs and their partial derivatives are derived. Finally numerical experiments are reported to show the efficiency of the proposed method.     

Commutative matrix subalgebras and length function

It is proved that if the length of a commutative matrix subalgebra is maximal then this subalgebra is maximal under inclusion. The examples are given showing that the converse does not hold. To establish this result, we prove several fundamental properties of the length function.