An integral square matrix A is called principally unimodular (PU if every nonsingular principal submatrix is unimodular (that is, has determinant \1). Principal unimodularity was originally studied with regard to skew-symmetric matrices; see [2, 4, 5]; here we consider symmetric matrices. Our main theorem is a generalization of Tutte's excluded minor characterization of totally unimodular matri...

We prove some results on the kernel of the Abel transform on an irreducible Riemannian symmetric space X = G=K with G noncompact and complex, in particular an estimate of this kernel. We also study the behaviour of spherical functions near the walls of Weyl chambers. We show how these harmonic spherical analysis results lead to a new proof of a central limit theorem of Guivarc'h and Raugi in th...

Let $Rin textbf{C}^{mtimes m}$ and $Sin textbf{C}^{ntimes n}$ be nontrivial involution matrices; i.e., $R=R^{-1}neq pm~I$ and $S=S^{-1}neq pm~I$. An $mtimes n$ complex matrix $A$ is said to be an $(R, S)$-symmetric ($(R, S)$-skew symmetric) matrix if $RAS =A$ ($ RAS =-A$). The $(R, S)$-symmetric and $(R, S)$-skew symmetric matrices have a number of special properties and widely used in eng...

We are concerned with the linear-quadratic optimal stochastic control problem where all the coefficients of the control system and the running weighting matrices in the cost functional are allowed to be predictable (but essentially bounded) processes and the terminal state-weighting matrix in the cost functional is allowed to be random. Under suitable conditions, we prove that the value field V...

An efficient numerical method is developed for evaluating φ(A), where A is a symmetric matrix and φ is the function defined by φ(x) = (ex − 1)/x = 1+ x/2 + x2/6+ .... This matrix function is useful in the so-called exponential integrators for differential equations. In particular, it is related to the exact solution of the ODE system dy/dt = Ay + b, where A and b are t-independent. Our method a...

Two canonical forms for skew-symmetric matrix polynomials over arbitrary fields are characterized — the Smith form, and its skew-symmetric variant obtained via unimodular congruences. Applications include the analysis of the eigenvalue and elementary divisor structure of products of two skew-symmetric matrices, the derivation of a Smith-McMillan-like canonical form for skew-symmetric rational m...

We show that any symmetric positive definite homogeneous matrix polynomial M ∈ R[x1, . . . , xn] admits a piecewise semi-certificate, i.e. a collection of identites M(x) = P j fi,j(x)Ui,j(x) T Ui,j(x) where Ui,j(x) is a matrix polynomial and fi,j(x) is a non negative polynomial on a semialgebraic subset Si, where R = ∪ri=1Si. This result generalizes to the setting of biforms. Some examples of c...

We give a matrix factorization for the solution of the linear system Ax = f , when coefficient matrix A is a dense symmetric positive definite matrix. We call this factorization as "WW T factorization". The algorithm for this factorization is given. Existence and backward error analysis of the method are given. The WDWT factorization is also presented. When the coefficient matrix is a symmetric...

Definitions and notation: Let A be a central simple algebra of degree n over a field k of characteristic different from 2. An involution on A is a ring antiautomorphism of order at most 2. An involution σ is of the first kind if σ|k = Idk, and of the second kind if σ|k is a non trivial involution on k, denoted by .̄ In the last case, k is a quadratic extension of the subfield k0 fixed by .̄ So we...