Low-rank Sylvester ADI Implementations (v1.0)

This package contains implementations of the the low-rank Sylvester ADI iteration (LR-ADI) to approximately solve large-scale Sylvester equations

 A*X*C + M*X*B = -f*g'

with n x n matrices A,M,  m x m matrices B,C, and f is n x r, g is m x r with  r<<min(n,m). M,C can be identity matrices for standard Sylvester equations.
The spectra (A,M), (B,C) should not intercect and have eigenvalues only in the left or right half plane.. 

Getting Started

Dependencies

MATLAB R2017b and above with the System Control Toolbox. Or Octave 4.0.0 and above with the Control Package. 
The data-files of one demo can be downloaded at mor-wiki.
The M-M.E.S.S. package (https://www.mpi-magdeburg.mpg.de/projects/mess) is not required but can be used to generate some matrices from scratch.

Installation

For using this software, you have to add the folder containing this file to the MATLAB/Octave search path.

Usage

• The main routine is lr_adi_sylv.m which implements the (inexact) low-rank Sylvester ADI iteration.  It has several settings to stear the iterative solution process of the arising linear system (which solver/preconditioner, how to set stopping criteria,....).
• lr_adi_sylv_ex_cpx.m is a much simpler ['exact'] version which only uses direct linear solves with the backslash '\'. These two methods will generate complex low-rank factors, if some of the ADI shift parameters were complex.
• lr_adi_sylv_ex_real.m also uses direct linear solves only, but generates real low-rank factors (provided the sets of shifts are in a suitable ordering). Note the sign-changes in the equation! This version does not yet support inexact linear solves.
• 'demo0.m' shows the basic usage of all the Sylvester ADI versions on an example.
  'demo1.m' illustrates some for finer settings for the inexact solves.
  'demo2.m' to 'demo4.m' run further examples for (generalized) Sylvester equations. 'demo4.m' requires a dataset from https://zenodo.org/records/10017861
• The other '*.m' files contain further routines needed to run the main algorithms. This includes routines for basic ADI shift generation strategies and basic Arnoldi-routines for Ritz value generation.

Please see the comments / help texts for 'lr_adi_sylv.m, lr_adi_sylv_ex_clx.m, lr_adi_sylv_ex_real.m' plus the demos on the usage and possible input parameters.

Authors

Patrick Kürschner

• affiliation: University of Applied Sciences Leipzig (HTWK Leipzig)
• email: patrick.kuerschner@htwk-leipzig.de
• orcid: 0000-0002-6114-8821

License

Copyright (C) 2024 Patrick Kürschner

The software is licensed under the BSD-2-Clause License. See COPYING for a copy of the license.

Citation

DOI

The DOI for version 1.0 is  10.5281/zenodo.10454177

BibTeX

@MISC{Kuersch_Sylvester_ADI,
  author = {K{\"u}rschner, P.},
  title  = {Low-rank Sylvester ADI Implementations (v1.0)},
  year   = 2024,
  doi    = {10.5281/zenodo.10454177}
}

Acknowledgements

Earlier parts of the development of this work (esp. the exact ADI) was supported by Max Planck Society during the time of the author's affiliation with the Max Planck Institute for Dynamics of Complex Technical System, Magdeburg, Germany;