Formal derivation of algorithms
- 1 June 2003
- journal article
- Published by Association for Computing Machinery (ACM) in ACM Transactions on Mathematical Software
- Vol. 29 (2) , 218-243
- https://doi.org/10.1145/779359.779365
Abstract
In this paper we apply a formal approach for the derivation of dense linear algebra algorithms to the triangular Sylvester equation. The result is a large family of provably correct algorithms. By using a coding style that reflects the algorithms as they are naturally presented, the correctness of the algorithms carries through to the correctness of the implementations. Analytically motivated heuristics are used to subsequently choose members from the family that can be expected to yield high performance. Finally, we report performance on the Intel (R) Pentium (R) III processor that is competitive with that of recursive algorithms reported previously in the literature for this operation.Keywords
This publication has 11 references indexed in Scilit:
- Accuracy and Stability of Numerical AlgorithmsPublished by Society for Industrial & Applied Mathematics (SIAM) ,2002
- FLAMEACM Transactions on Mathematical Software, 2001
- Application of ADI Iterative Methods to the Restoration of Noisy ImagesSIAM Journal on Matrix Analysis and Applications, 1996
- Solution of the Sylvester matrix equation AXB
T
+ CXD
T
= EACM Transactions on Mathematical Software, 1992
- Distributed and Shared Memory Block Algorithms for the Triangular Sylvester Equation with $\operatorname{sep}^{ - 1} $ EstimatorsSIAM Journal on Matrix Analysis and Applications, 1992
- A set of level 3 basic linear algebra subprogramsACM Transactions on Mathematical Software, 1990
- An extended set of FORTRAN basic linear algebra subprogramsACM Transactions on Mathematical Software, 1988
- Numerical Solution of the Stable, Non-negative Definite Lyapunov Equation Lyapunov EquationIMA Journal of Numerical Analysis, 1982
- A Hessenberg-Schur method for the problem AX + XB= CIEEE Transactions on Automatic Control, 1979
- Algorithm 432 [C2]: Solution of the matrix equation AX + XB = C [F4]Communications of the ACM, 1972