On some Krylov subspace based methods for large-scale nonsymmetric algebraic Riccati problems

被引:8
|
作者
Bentbib, A. [1 ]
Jbilou, K. [2 ]
Sadek, E. M. [1 ,2 ]
机构
[1] Univ Cadi Ayyad, Lab LAMAI, FST, Marrakech, Morocco
[2] ULCO Calais, LMPA, 50 Rue F Buisson, Calais, France
关键词
Extended block Amoldi; Low-rank approximation; Newton method; Nonsymmetric Riccati equation; Transport theory; DOUBLING-ALGORITHM; ITERATIVE SOLUTION; EQUATIONS;
D O I
10.1016/j.camwa.2015.09.025
中图分类号
O29 [应用数学];
学科分类号
070104 ;
摘要
In the present paper, we consider large scale nonsymmetric matrix Riccati equations with low rank right hand sides. These matrix equations appear in many applications such as transport theory, Wiener-Hopf factorization of Markov chains, applied probability and others. We show how to apply directly Krylov methods such as the extended block Arnoldi algorithm to get low rank approximate solutions. We also combine the Newton method and block Krylov subspace methods to get approximations of the desired minimal nonnegative solution. We give some theoretical results and report some numerical experiments for the well known transport equation. (C) 2015 Elsevier Ltd. All rights reserved.
引用
收藏
页码:2555 / 2565
页数:11
相关论文
共 50 条