Related papers: Compact Schemes for $A^+B$, $A^+AB$ and $AA^+B$
Let K be an arbitrary field, and a,b,c,d be elements of K such that the polynomials t^2-at-b and t^2-ct-d are split in K[t]. Given a square matrix M with entries in K, we give necessary and sufficient conditions for the existence of two…
A common optimization problem is the minimization of a symmetric positive definite quadratic form $< x,Tx >$ under linear constrains. The solution to this problem may be given using the Moore-Penrose inverse matrix. In this work we extend…
The combined matrix is a very useful concept for many applications. Almost strictly sign regular (ASSR) matrices form an important structured class of matrices with two possible zero patterns, which are either type-I staircase or type-II…
We present scalable distributed-memory algorithms for sparse matrix permutation, extraction, and assignment. Our methods follow an Identify-Exchange-Build (IEB) strategy where each process identifies the local nonzeros to be sent, exchanges…
A formula is given for the propagation of errors during matrix inversion. An explicit calculation for a 2 by 2 matrix using both the formula and a Monte Carlo calculation are compared. A prescription is given to determine when a matrix with…
An algorithm is given for explicitly computing Penrose diagrams for spacetimes of the form $ds^2 = -f(r)\, dt^2 + f(r)^{-1} \, dr^2 + r^2 \, d\Omega^2$. The resulting diagram coordinates are shown to extend the metric continuously and…
We present positive maps and matrix inequalities for variables from the positive cone. These inequalities contain partial transpose and reshuffling operations, and can be understood as positive multilinear maps that are in one-to-one…
The Moore-Penrose pseudo-inverse $X^\dagger$, defined for rectangular matrices, naturally emerges in many areas of mathematics and science. For a pair of rectangular matrices $X, Y$ where the corresponding entries are jointly Gaussian and…
We develop several efficient algorithms for the classical \emph{Matrix Scaling} problem, which is used in many diverse areas, from preconditioning linear systems to approximation of the permanent. On an input $n\times n$ matrix $A$, this…
In this paper, we present some extensions of interpolation between the arithmetic-geometric means inequality. Among other inequalities, it is shown that if $A, B, X$ are $n\times n$ matrices, then \begin{align*}…
A new procedure is presented for computing the matrix cosine and sine simultaneously by means of Taylor polynomial approximations. These are factorized so as to reduce the number of matrix products involved. Two versions are developed to be…
To locate all eigenvalues of a matrix more precisely, we exclude some sets which do not include any eigenvalue of the matrix from the well-known Brauer set to give two new Brauer-type eigenvalue inclusion sets. And it is also shown that the…
In this paper we give different compactifications for the domain and the codomain of an affine rational map $f$ which parametrizes a hypersurface. We show that the closure of the image of this map (with possibly some other extra…
The study of solving the inverse eigenvalue problem for nonnegative matrices has been around for decades. It is clear that an inverse eigenvalue problem is trivial if the desirable matrix is not restricted to a certain structure. Provided…
Let $f$ bea noncommutativepolynomial of degree $m\ge 1$ over an algebraically closed field $F$ of characteristic $0$. If $n\ge m-1$ and $\alpha_1,\alpha_2,\alpha_3$ are nonzero elements from $F$ such that $\alpha_1+\alpha_2+\alpha_3=0$,…
In this paper, we first present an explicit expression for the inverse\emph{} of a type of matrices. As special applications, the inverse of some matrices arising from implicit time integration techniques, such as the well-known implicit…
For an upper bidiagonal matrix $B$ where all the diagonal and the upper subdiagonal entries are positive, two subtraction-free formulae for computation of the traces $J_{M} ( B ) = \textrm{Tr} ( ( B^{\top} B )^{- M} ) = \textrm{Tr} ( ( B…
We propose splitting methods for the computation of the exponential of perturbed matrices which can be written as the sum $A=D+\varepsilon B$ of a sparse and efficiently exponentiable matrix $D$ with sparse exponential $e^D$ and a dense…
A classification of the semisimple subalgebras of the Lie algebra of traceless $3\times 3$ matrices with complex entries, denoted $A_2$, is well-known. We classify its nonsemisimple subalgebras, thus completing the classification of the…
An efficient algorithm for computing eigenvectors of a matrix of integers by exact computation is proposed. The components of calculated eigenvectors are expressed as polynomials in the eigenvalue to which the eigenvector is associated, as…