Article

Recursive blocked algorithms for solving periodic triangular Sylvester-type matrix equations

Authors:

Bo KågströmAuthors Info & Claims

PARA'06: Proceedings of the 8th international conference on Applied parallel computing: state of the art in scientific computing

Pages 531 - 539

Published: 18 June 2006 Publication History

Abstract

Recently, recursive blocked algorithms for solving triangular one-sided and two-sided Sylvester-type equations were introduced by Jonsson and Kågströ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 add (GEMM) operations. We extend and apply the recursive blocking technique to solving periodic Sylvester-type matrix equations. Successive recursive splittings are performed on 3-dimensional arrays, where the third dimension represents the periodicity of a matrix equation.

References

[1]

Bartels, R.H., Stewart, G.W.: Solution of the equation AX + XB = C. Comm. Assoc. Comput. Mach. 15, 820-826 (1972).

Digital Library

[2]

Benner, P., Mehrmann, V., Xu, H.: Perturbation analysis for the eigenvalue problem of a formal product of matrices. BIT 42(1), 1-43 (2002).

[3]

Bittanti, S., Colaneri, P. (eds.): Periodic Control Systems. In: Proceedings volume from the IFAC Workshop, August 27-28, 2001, Elsevier Science & Technology, Cernobbio-Como, Italy (2001).

[4]

Bojanczyk, A.W., Golub, G., Van Dooren, P.: The Periodic Schur Decomposition. Algorithms and Applications. In: Luk, F.T. (ed.) Proceedings SPIE Conference, vol. 1770, pp. 31-42 (1992).

[5]

Elmroth, E., Gustavson, F., Jonsson, I., Kågström, B.: Recursive Blocked Algorithms and Hybrid Data Structures for Dense Matrix Library Software. SIAM Review 46(1), 3-45 (2004).

[6]

Fairweather, G., Gladwell, I.: Algorithms for Almost Block Diagonal Linear Systems. SIAM Review 44(1), 49-58 (2004).

[7]

Granat, R., Kågström, B.: Direct Eigenvalue Reordering in a Product of Matrices in Periodic Schur Form. SIAM J. Matrix Anal. Appl. 28(1), 285-300 (2006).

Digital Library

[8]

Granat, R., Kågström, B., Kressner, D.: Reordering the Eigenvalues of a Periodic Matrix Pair with Applications in Control. In: Proc. of 2006 IEEE Conference on Computer Aided Control Systems Design (CACSD), pp. 25-30 (2006) ISBN: 0- 7803-9797-5.

[9]

Granat, R., Kågström, B., Kressner, D.: Computing Periodic Deflating Subspaces Associated with a Specified Set of Eigenvalues. BIT Numerical Mathematics, December 2006 (submitted).

[10]

Hench, J.J., Laub, A.J.: Numerical solution of the discrete-time periodic Riccati equation. IEEE Trans. Automat. Control 39(6), 1197-1210 (1994).

[11]

Jonsson, I., Kågström, B.: Recursive blocked algorithms for solving triangular systems -- Part I: One-sided and coupled Sylvester-type matrix equations. ACM Trans. Math. Softw. 28(4), 392-415 (2002).

Digital Library

[12]

Jonsson, I., Kågström, B.: Recursive blocked algorithms for solving triangular systems -- Part II: Two-sided and generalized Sylvester and Lyapunov matrix equations. ACM Trans. Math. Softw. 28(4), 416-435 (2002).

Digital Library

[13]

Jonsson, I., Kågström, B.: RECSY -- A High Performance Library for Sylvester-Type Matrix Equations. In: Kosch, H., Böszörményi, L., Hellwagner, H. (eds.) Euro-Par 2003. LNCS, vol. 2790, pp. 810-819. Springer, Heidelberg (2003).

[14]

Kågström, B., Poromaa, P.: LAPACK-Style Algorithms and Software for Solving the Generalized Sylvester Equation and Estimating the Separation between Regular Matrix Pairs. ACM Trans. Math. Software 22, 78-103 (1996).

Digital Library

[15]

Kågström, B.,Westin, L.: Generalized Schur methods with condition estimators for solving the generalized Sylvester equation. IEEE Trans. Automat. Control 34(4), 745-751 (1989).

[16]

Sreedhar, J., Van Dooren, P.: A Schur approach for solving some matrix equations. In: Helmke, U., Menniken, R., Saurer, J. (eds.) Systems and Networks: Mathematical Theory and Applications, Mathematical Research, vol. 77, pp. 339-362 (1994).

[17]

Sun, J.-G.: Perturbation bounds for subspaces associated with periodic eigenproblems. Taiwanese Journal of Mathematics 9(1), 17-38 (2005).

[18]

Varga, A.: Periodic Lyapunov equations: some applications and new algorithms. Internat. J. Control 67(1), 69-87 (1997).

[19]

Varga, A., Van Dooren, P.: Computational methods for periodic systems. In: Prepr. IFAC Workshop on Periodic Control Systems, Como, Italy, pp. 177-182 (2001).

Cited By

Huang BMa C(2019)Finite iterative algorithm for the symmetric periodic least squares solutions of a class of periodic Sylvester matrix equationsNumerical Algorithms10.1007/s11075-018-0553-881:1(377-406)Online publication date: 21-May-2019
https://dl.acm.org/doi/10.1007/s11075-018-0553-8
Hajarian M(2018)A Finite Iterative Method for Solving the General Coupled Discrete-Time Periodic Matrix EquationsCircuits, Systems, and Signal Processing10.1007/s00034-014-9842-134:1(105-125)Online publication date: 27-Dec-2018
https://dl.acm.org/doi/10.1007/s00034-014-9842-1
Hajarian M(2017)Convergence analysis of generalized conjugate direction method to solve general coupled Sylvester discrete-time periodic matrix equationsInternational Journal of Adaptive Control and Signal Processing10.5555/3193602.319360531:7(985-1002)Online publication date: 1-Jul-2017
https://dl.acm.org/doi/10.5555/3193602.3193605
Show More Cited By

Recursive blocked algorithms for solving periodic triangular Sylvester-type matrix equations
1. Computing methodologies
  1. Symbolic and algebraic manipulation
    1. Symbolic and algebraic algorithms

Recommendations

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

Triangular matrix equations appear naturally in estimating the condition numbers of matrix equations and different eigenspace computations, including block-diagonalization of matrices and matrix pairs and computation of functions of matrices. To solve a ...
Recursive blocked algorithms for solving triangular systems—Part II: two-sided and generalized Sylvester and Lyapunov matrix equations

We continue our study of high-performance algorithms for solving triangular matrix equations. They appear naturally in different condition estimation problems for matrix equations and various eigenspace computations, and as reduced systems in standard ...
Parallel algorithms and condition estimators for standard and generalized triangular Sylvester-type matrix equations
PARA'06: Proceedings of the 8th international conference on Applied parallel computing: state of the art in scientific computing

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

Comments

Information & Contributors

Information

Published In

cover image Guide Proceedings

PARA'06: Proceedings of the 8th international conference on Applied parallel computing: state of the art in scientific computing

June 2006

1191 pages

ISBN:3540757546

Editors:
Bo Kågström
Umeå University, Department of Computing Science and High Performance Computing Center North, Umeå, Sweden
,
Erik Elmroth
Umeå University, Department of Computing Science and High Performance Computing Center North, Umeå, Sweden
,
Jack Dongarra
University of Tennessee, Department of Computer Science, Knoxville, TN
,
Jerzy Waśniewski
Technical University of Denmark, Informatics and Mathematical Modelling, Lyngby, Denmark

Sponsors

Umeå University
VR: The Swedish Research Council

Publisher

Springer-Verlag

Berlin, Heidelberg

Publication History

Published: 18 June 2006

Author Tags

Qualifiers

Article

Contributors

Other Metrics

View Article Metrics

Bibliometrics & Citations

Bibliometrics

Article Metrics

4
Total Citations
View Citations
0
Total Downloads

Downloads (Last 12 months)0
Downloads (Last 6 weeks)0

Reflects downloads up to 06 Jan 2025

Other Metrics

View Author Metrics

Citations

Cited By

Huang BMa C(2019)Finite iterative algorithm for the symmetric periodic least squares solutions of a class of periodic Sylvester matrix equationsNumerical Algorithms10.1007/s11075-018-0553-881:1(377-406)Online publication date: 21-May-2019
https://dl.acm.org/doi/10.1007/s11075-018-0553-8
Hajarian M(2018)A Finite Iterative Method for Solving the General Coupled Discrete-Time Periodic Matrix EquationsCircuits, Systems, and Signal Processing10.1007/s00034-014-9842-134:1(105-125)Online publication date: 27-Dec-2018
https://dl.acm.org/doi/10.1007/s00034-014-9842-1
Hajarian M(2017)Convergence analysis of generalized conjugate direction method to solve general coupled Sylvester discrete-time periodic matrix equationsInternational Journal of Adaptive Control and Signal Processing10.5555/3193602.319360531:7(985-1002)Online publication date: 1-Jul-2017
https://dl.acm.org/doi/10.5555/3193602.3193605
Andersson PGranat RJonsson IKågström B(2008)Parallel Algorithms for Triangular Periodic Sylvester-Type Matrix EquationsProceedings of the 14th international Euro-Par conference on Parallel Processing10.1007/978-3-540-85451-7_83(780-789)Online publication date: 26-Aug-2008
https://dl.acm.org/doi/10.1007/978-3-540-85451-7_83

View Options

View options

Media

Figures

Other

Tables

View Table of Contents