Improving accuracy of parallel SLICOT model reduction routines for stable systems

© 2015 IEEE. Personal use of this material is permitted. Permission from IEEE must be obtained for all other uses, in any current or future media, including reprinting/republishing this material for advertising or promotional purposes, creating new collective works, for resale or redistribution to servers or lists, or reuse of any copyrighted component of this work in other works.
This paper shows part of the work carried out to develop parallel versions of the SLICOT routines for model reduction of stable systems. In particular, the routines that have been parallelised are those based on the solution of Lyapunov equations. The goal is to be able to work with larger unreduced models and also to obtain better performance in the reduction process. New routines have been developed using standard libraries to improve portability and efficiency. A preliminary version was released previously by the authors, which achieved high performance. However, accuracy improvements have been necessary in order to make the new routines similar to the sequential ones in this aspect. Routines presented in this paper preserve good performance obtained by the previous parallel implementation while maintaining high accuracy of sequential SLICOT routines.