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 ...

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...

This Licentiate Thesis contains contributions in two different subfields of Computing Science: parallel ScaLAPACK-style algorithms for Sylvester-type matrix equations and periodic eigenvalue reordering in a cyclic product of matrices. Sylvester-type matrix equations, like the continuous-time Sylvester equation AX −XB = C, where A of size m×m, B of size n×n and C of size m×n are general matrices...

Some complex quaternionic equations in the type AX - XB = C are investigated. For convenience, these equations were called generalized Sylvester-quaternion equations, which include the Sylvester equation as special cases. By the real matrix representations of complex quaternions, the necessary and sufficient conditions for the solvability and the general expressions of the solutions are obtained.

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 ...

Recently, recursive blocked algorithms for solving triangular one-sided and two-sided Sylvester-type equations were introduced by Jonsson and K̊agström. This elegant yet simple technique enables an automatic variable blocking that has the potential of matching the memory hierarchies of today’s HPC systems. The main parts of the computations are performed as level 3 general matrix multiply and ad...

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 ...