In this paper we consider the numerical solution of large-scale projected generalized continuous-time and discrete-time Sylvester equations with low-rank right-hand sides. First, we present the results on the sufficient conditions for the existence, uniqueness, and analytic formula of the solutions of these equations. Second, we review the low-rank alternating direction implicit method and the ...

We discuss parallel algorithms for solving eight common standard and generalized triangular Sylvester-type matrix equation. Our parallel algorithms are based on explicit blocking, 2D block-cyclic data distribution of the matrices and wavefront-like traversal of the right hand side matrices while solving small-sized matrix equations at different nodes and updating the rest of the right hand side...

An iterative algorithm is constructed to solve the generalized coupled Sylvester matrix equations AXB − CYD,EXF − GYH M,N , which includes Sylvester and Lyapunov matrix equations as special cases, over generalized reflexive matrices X and Y . When the matrix equations are consistent, for any initial generalized reflexive matrix pair X1, Y1 , the generalized reflexive solutions can be obtained b...

An iterative algorithm is presented for solving the extended Sylvester-conjugate matrix equation. By this iterativemethod, the solvability of thematrix equation can be determined automatically. When the matrix equation is consistent, a solution can be obtained within finite iteration steps for any initial values in the absence of round-off errors. The algorithm is also generalized to solve a mo...

We present a new algorithm for solving the Sylvester-Observer Equation: AX ? XH = (0;C). The algorithm embodies two main computational phases: the solution of a series of independent equation systems, and a series of matrix-matrix multiplications. The algorithm is, thus, well suited for parallel and high performance computing. By reducingthe coeecient matrix A to lower Hessenberg form, one can ...

In this paper, Haar wavelet operational matrix method is proposed to solve a class of fractional partial differential equations. We derive the Haar wavelet operational matrix of fractional order integration. Meanwhile, the Haar wavelet operational matrix of fractional order differentiation is obtained. The operational matrix of fractional order differentiation is utilized to reduce the initial ...

The global FOM and GMRES algorithms are among the effective methods to solve Sylvester matrix equations. In this paper, we study these algorithms in the case that the coefficient matrices are real symmetric (real symmetric positive definite) and extract two CG-type algorithms for solving generalized Sylvester matrix equations. The proposed methods are iterative projection methods onto matrix Kr...

Many Krylov subspace methods for shifted linear systems take advantage of the invariance of the Krylov subspace under a shift of the matrix. However, exploiting this fact in the non-Hermitian case introduces restrictions; e.g., initial residuals must be collinear and this collinearity must be maintained at restart. Thus we cannot simultaneously solve shifted systems with unrelated right-hand si...