Parallel Two-Sided Sylvester-Type Matrix Equation Solvers for SMP Systems Using Recursive Blocking

Authors:
Isak Jonsson;Bo Kågström
Affiliations:
-;-
Venue:
PARA '02 Proceedings of the 6th International Conference on Applied Parallel Computing Advanced Scientific Computing
Year:
2002

Citing 13
Cited 0

Distributed and shared memory block algorithms for the triangular Sylvester equation with sep-1 estimators

SIAM Journal on Matrix Analysis and Applications
Solution of the Sylvester matrix equation AXBT + CXDT = E

ACM Transactions on Mathematical Software (TOMS)
Algorithm 705; a FORTRAN-77 software package for solving the Sylvester matrix equation AXBT + CXDT = E

ACM Transactions on Mathematical Software (TOMS)
LAPACK-style algorithms and software for solving the generalized Sylvester equation and estimating the separation between regular matrix pairs

ACM Transactions on Mathematical Software (TOMS)
Recursion leads to automatic variable blocking for dense linear-algebra algorithms

IBM Journal of Research and Development
LAPACK Users' guide (third ed.)

LAPACK Users' guide (third ed.)
Blocked algorithms and software for reduction of a regular matrix pair to generalized Schur form

ACM Transactions on Mathematical Software (TOMS)
Solution of the matrix equation AX + XB = C [F4]

Communications of the ACM
Recursive blocked algorithms for solving triangular systems—Part I: one-sided and coupled Sylvester-type matrix equations

ACM Transactions on Mathematical Software (TOMS)
Recursive blocked algorithms for solving triangular systems—Part II: two-sided and generalized Sylvester and Lyapunov matrix equations

ACM Transactions on Mathematical Software (TOMS)
Recursive Blocked Data Formats and BLAS's for Dense Linear Algebra Algorithms

PARA '98 Proceedings of the 4th International Workshop on Applied Parallel Computing, Large Scale Scientific and Industrial Problems
Parallel Triangular Sylvester-Type Matrix Equation Solvers for SMP Systems Using Recursive Blocking

PARA '00 Proceedings of the 5th International Workshop on Applied Parallel Computing, New Paradigms for HPC in Industry and Academia
Parallel Algorithms for Triangular Sylvester Equations: Design, Scheduling and Saclability Issues

PARA '98 Proceedings of the 4th International Workshop on Applied Parallel Computing, Large Scale Scientific and Industrial Problems

Quantified Score

Hi-index	0.00

Visualization

Abstract

We present recursive blocked algorithms for solving triangular two-sided 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. This is a continuation of the work presented at the PARA2000 conference ([9]), where we presented results for one-sided Sylvester-type matrix equations. The improvements for two-sided Sylvester-type matrix equations are remarkable, and make a substantial impact on solving unreduced matrix equations problems as well. Uniprocessor and SMP parallel performance results are presented and compared with results from existing LAPACK and SLICOT routines for solving this type of matrix equations.