We present recursive blocked algorithms for solving triangular Sylvester-type matrix equations. Recursion leads to automatic blocking that is variable and \squarish". The main part of the computations are performed as level 3 general matrix multiply and add (GEMM) operations. We also present new highly optimized superscalar kernels for solving small-sized matrix equations stored in level 1 cach...

We introduce the concept of Stieltjes integral of an operator-valued function with respect to the spectral measure associated with a normal operator. We give sufficient conditions for the existence of this integral and find bounds on its norm. The results obtained are applied to the Sylvester and Riccati operator equations. Assuming that the entry C is a normal operator, the spectrum of the ent...

We present an iterative method based on the Hermitian and skew-Hermitian splitting for solving the continuous Sylvester equation. By using the Hermitian and skew-Hermitian splitting of the coefficient matrices A and B, we establish a method which is practically inner/outer iterations, by employing a CGNR-like method as inner iteration to approximate each outer iterate, while each outer iteratio...

We introduce the notion of Stieltjes integral with respect to the spectral measure corresponding to a normal operator. Sufficient conditions for the existence of this integral are given, and estimates for its norm are established. The results are applied to operator Sylvester and Riccati equations. Assuming that the spectrum of a closed densely defined operator A does not have common points wit...

This paper proposes a fast multi-band image fusion algorithm, which combines a high-spatial low-spectral resolution image and a low-spatial high-spectral resolution image. The well admitted forward model is explored to form the likelihoods of the observations. Maximizing the likelihoods leads to solving a Sylvester equation. By exploiting the properties of the circulant and downsampling matrice...

We consider the linear matrix inequality (LMI) problem of H∞ output feedback control for a generalized plant whose input, measured output, disturbance and controlled are scalar. provide an explicit form optimal value. This is unification some results in literature performance limitation analysis. To obtain value, we focus on non-uniqueness perpendicular matrices, which appear LMI problem. use n...

We present a nested splitting conjugate gradient iteration method for solving large sparse continuous Sylvester equation, in which both coefficient matrices are (non-Hermitian) positive semi-definite, and at least one of them is positive definite. This method is actually inner/outer iterations, which employs the Sylvester conjugate gradient method as inner iteration to approximate each outer it...

Ion mobility spectrometry (IMS) has been widely used for the on-site detection of trace chemicals, but continue to suffer from a low duty cycle ion injection. The Hadamard transform (HT-IMS) technique was employed address problem with increased signal-to-noise ratio (SNR). However, in this work, through simulation, certain deviation between mathematical principle and actual data collection proc...

In various applications, for instance in the detection of a Hopf bifurcation or in solving separable boundary value problems using the two-parameter eigenvalue problem, one has to solve a generalized eigenvalue problem of the form (B1 ⊗A2 −A1 ⊗B2)z = μ(B1 ⊗ C2 − C1 ⊗B2)z, where matrices are 2 × 2 operator determinants. We present efficient methods that can be used to compute a small subset of t...