Solving stable Sylvester equations via rational iterative schemes

Journal of Scientific Computing

Volumn 28, Issue 1, 2006, Pages 51-83

  Benner, Peter ^a   Quintana Ortí, Enrique S ^b   Quintana Ortí, Gregorio ^b  

^a CHEMNITZ UNIVERSITY OF TECHNOLOGY   (Germany)

^b UNIVERSITAT JAUME I   (Spain)

Author keywords

Halley's method; Image restoration; Matrix sign function; Model reduction; Newton iteration; Newton Schulz iteration; Parallel algorithms; Sylvester equation

Indexed keywords

CONVERGENCE OF NUMERICAL METHODS; IMAGE RECONSTRUCTION; ITERATIVE METHODS; MATHEMATICAL MODELS; NUMERICAL ANALYSIS; PARALLEL ALGORITHMS; PROGRAM PROCESSORS;

HALLEY'S METHOD; MATRIX SIGN FUNCTION; NEWTON ITERATION; NEWTON-SCHULZ ITERATION; SYLVESTER EQUATION;

MATRIX ALGEBRA;

EID: 33746621687     PISSN: 08857474     EISSN: None     Source Type: Journal    
DOI: 10.1007/s10915-005-9007-2     Document Type: Article

References (47)

1. Abels J., and Benner, P. (1999). CAREX - a collection of benchmark examples for continuous-time algebraic Riccati equations (version 2.0). SLICOT Working Note 1999-14. Available from http://www.win.tue.nl/niconet/NIC2/reports. html.
2. Aldhaheri, R. (1991). Model order reduction via real Schur-form decomposition. Int. J. Control, 53(3), 709-716.
3. Anderson, E., Bai, Z., Bischof, C., Demmel, J., Dongarra, J., Du Croz, J., Greenbaum, A., Hammarling, S., McKenney, A., and Sorensen, D. (1999). LAPACK Users' Guide, 3rd ed., SIAM, Philadelphia, PA.
4. Antoulas, A. Approximation of Large-Scale Dynamical Systems. SIAM Publications, Philadelphia, PA, 2005.
5. Bai, Z., and Demmel, J. (1993). Design of a parallel nonsymmetric eigenroutine toolbox, Part I. In R. S. et al., (eds.), Proceedings of the Sixth SIAM Conference on Parallel Processing for Scientific Computing, SIAM, Philadelphia, PA, 391-398, See also: Tech. Report CSD-92-718, Computer Science Division, University of California, Berkeley, CA.
6. Bartels, R., and Stewart, G. (1972). Solution of the matrix equation AX + XB = C: Algorithm 432. Comm. ACM, 15, 820-826.
7. Beavers, A. N., and Denman, E. D. (1975). A new solution method for the Lyapunov matrix equations. SIAM J. Appl. Math. 29, 416-421.
8. Benner, P. (2004). Factorized solution of Sylvester equations with applications in control. In Proc. Intl. Symp. Math. Theory Networks and Syst. MTNS 2004, http://www.mtns2004.be.
9. Benner, P., Claver, v, and Quintana-Orti, E. (1999). Parallel distributed solvers for large stable generalized Lyapunov equations. Parallel Processing Lett., 9(1), 147-158.
10. Benner, P., and Quintana-Orti, E. (1999). Solving stable generalized Lyapunov equations with the matrix sign function. Numer. Algorithms 20(1), 75-100.
11. Bemier, P., Quintana-Ortí, E., and Quintana-Ortí, G. (2002). Numerical solution of discrete stable linear matrix equations on multicomputers. Parallel Algorithms and Appl. 17(1), 127-146.
12. Blackford, L., Choi, J., Cleary, A., D'Azevedo, E., Demmel, J., Dhillon, I., Dongarra, J., Hammarling, S., Henry, G., Petitet, A., Stanley, K., Walker, D., and Whaley, R. (1997). ScaLAPACK Users' Guide. SIAM, Philadelphia, PA.
13. Byers, R. (1987). Solving the algebraic Riccati equation with the matrix sign function. Linear Algebra Appl. 85, 267-279.
14. Calvetti, D., and Reichel, L. (1996). Application of ADI iterative methods to the restoration of noisy images. SIAM J. Matrix Anal. Appl. 17, 165-186.
15. Chan, T. (1987). Rank revealing QR factorizations. Linear Algebra Appl. 88/89, 67-82.
16. Choi, C., and Laub, A. (1990). Efficient matrix-valued algorithms for solving stiff Riccati differential equations. IEEE Trans. Automat. Control, 35, 770-776.
17. Datta, B. (2003). Numerical Methods for Linear Control Systems Design and Analysis. Elsevier Press, Amsterdam.
18. Dieci, L., Osborne, M., and Russell, R., (1988). A. Riccati, transformation method for solving linear BVPs. I: Theoretical aspects. SIAM J. Numer. Anal. 25(5), 1055-1073.
19. Dongarra, J., Croz, J. D., Duff, I., and Hammarling, S. (1990). A set of level 3 Basic linear algebra subprograms. ACM Trans. Math. Soft., 16:1-17.
20. Dongarra, J. Duff, I. S., Sorensen, D. C., and van der Vorst, H. A. (1991). Solving Linear Systems on Vector and Shared Memory Computers. SIAM, Philadelphia, PA.
21. Enright, W (1978). Improving the efficiency of matrix operations in the numerical solution of stiff ordinary differential equations. ACM Trans. Math. Soft. 4, 127-136.
22. Epton, M. (1980). Methods for the solution of AXD - BXC = E and its application in the numerical solution of implicit ordinary differential equations. BIT, 20, 341-345.
23. Fernando, K., and Nicholson, H. (1984). On a fundamental property of the cross-Gramian matrix. IEEE Trans. Circuits Syst. CAS-31(5), 504-505.
24. Gajić, Z., and Qureshi, M. (1995). Lyapunov Matrix Equation in System Stability and Control. Math, in Science and Engineering. Academic Press, San Diego, CA.
25. Gardiner, J., and Laub, A. (1991). Parallel algorithms for algebraic Riccati equations. Int. J. Control, 54(6), 1317-1333.
26. Gardiner, J., Laub, A., Amato, J., and Moler, C. (1992). Solution of the Sylvester matrix equation AXB + CXD = E. ACM Trans. Math. Soft, 18, 223-231.
27. Golub, G., and Van Loan, C. (1996). Matrix Computations, 3rd ed., Johns Hopkins University Press, Baltimore.
28. Golub, G. H., Nash, S., and Van Loan, C. F. (1979). A Hessenberg-Schur method for the problem AX + XB = C. IEEE Trans. Automat. Control. AC-24, 909-913.
29. Goto, K., and van de Geijn, R. (2002). On reducing TLB misses in matrix multiplication. FLAME Working Note 9, Department of Computer Sciences, The University of Texas at Austin, http://www.cs.utexas.edu/users/flame.
30. Grasedyck, L. (2004). Existence of a low rank or H-matrix approximant to the solution of a Sylvester equation. Numer. Lin. Alg. Appl. 11, 371-389.
31. Grasedyck, L., and Hackbusch, W. (2004). A multigrid method to solve large scale Sylvester equations. Preprint 48, Max-Planck Institut für Mathematik in den Naturwissenschaften, Leipzig, Germany.
32. Gropp, W., Lusk, E. and Skjellum, A. (1994). Using MPI: Portable Parallel Programming with the Message-Passing Interface. MIT Press, Cambridge, MA.
33. Halley, E. (1694). Methodus nova, accurata & facilis inveniendi radices aequationum quarumcumque generaliter, sine praevia reductione. Philos. Trans. R Soc. Lond. 18, 136-148.
34. Henry, G., and van de Geijn, R. (1997). Parallelizing the QR algorithm for the unsymmetric algebraic eigenvalue problem: myths and reality. SIAM J. Sci. Comput. 17, 870-883.
35. Higham, N. (1986). Computing the polar decomposition-with applications. SIAM J. Sci. Statist. Comput. 7, 1160-1174.
36. Hoskins, W., Meek, D., and Walton, D. (1977). The numerical solution of the matrix equation XA + AY = F. BIT 17, 184-190.
37. Hu, D., and Reichel, L. (1992). Krylov-subspace methods for the Sylvester equation. Linear Algebra Appl. 172, 283-313.
38. Kagström, B., and Poromaa, P. (1996). Lapack-style algorithms and software for solving the generalized Sylvester equation and estimating the separation between regular matrix pairs. ACM Trans. Math. Soft. 22(1), 78-103.
39. Kenney, C., and Laub, A. (1991). Rational iterative methods for the matrix sign function. SIAM J. Matrix Anal. Appl. 12, 273-291.
40. Kenney, C., and Laub, A. (1995). The matrix sign function. IEEE Trans. Automat. Control 40(8), 1330-1348.
41. Quintana-Orti, E., Quintana-Orti, G., Sun, X., and van de Geijn, R. (2001). A note on parallel matrix inversion. SIAM J. Sci. Comput. 22, 1762-1771.
42. Roberts, J. (1980). Linear model reduction and solution of the algebraic Riccati equation by use of the sign function. Int. J. Control, 32, 677-687, (Reprint of Technical Report No. TR-13, CUED/B-Control, Cambridge University, Engineering Department, 1971).
43. Schulz, G. (1993), Iterative berechnung der reziproken matrix. Z. Angew. Math. Mech. 13, 57-59. In German.
44. Sima, V. (1996). Algorithms for Linear-Quadratic Optimization, Pure and Applied Mathematics, vol. 200, Marcel Dekker Inc., New York, NY
45. Slowik, M., Benner, P., and Sima, V. (2004). Evaluation of the linear matrix equation solvers in SLICOT. SLICOT Working Note 2004-1 Available from http://www.win.tue.nl/niconet/NIC2/reports.html.
46. van de Geijn, R., (1997). Using FLAPACK: Parallel Linear Algebra Package. MIT Press, Cambridge, MA.
47. Wachspress, E. (1988). Iterative solution of the Lyapunov matrix equation. Appl. Math. Lett. 107, 87-90.