Preprint No. MPIMD/14-21

Title: Low-Rank Newton-ADI methods for Large Nonsymmetric Algebraic Riccati Equations

Author(s): Peter Benner, Patrick K├╝rschner, Jens Saak

Email: kuerschner@mpi-magdeburg.mpg.de

Date: 2014-11-20

Abstract:

The numerical treatment of large-scale, nonsymmetric algebraic Riccati equations (NAREs) by a low-rank variant of Newton\'s method is considered. We discuss a method to compute approximations to the solution of the NARE in a factorized form of low rank. The occurring large-scale Sylvester equations are dealt with using the factored alternating direction implicit iteration (fADI). Several performance enhancing strategies available for the factored ADI as well as the related Newton-ADI for symmetric algebraic Riccati equations are generalized to this combination. This includes the efficient computation of the norm of the residual matrix, adapted shift parameters strategies for fADI, and an acceleration of the Newton\'s scheme by means of a Galerkin projection. Numerical experiments illustrate the capabilities of the proposed method to solve high-dimensional NAREs.

BibTeX:

@TECHREPORT{MPIMD14-21,
author = {Peter Benner and Patrick K├╝rschner and Jens Saak},
title = {Low-Rank Newton-ADI methods for Large Nonsymmetric Algebraic Riccati Equations},
number = {MPIMD/14-21},
month = nov,
year = 2014,
institution = {Max Planck Institute Magdeburg},
type = {Preprint},
note = {Available from \url{http://www.mpi-magdeburg.mpg.de/preprints/}},
}


Download MPIMD/14-21