×

Fast solution of a certain Riccati equation through Cauchy-like matrices. (English) Zbl 1189.65075

Summary: We consider a special instance of the algebraic Riccati equation \(XCX-XE-AX+B=0\) encountered in transport theory, where the \(n \times n\) matrix coefficients \(A, B, C, E\) are rank structured matrices. The equation is reduced to unilateral form \(A_1X^2 + A_0X + A_{-1}=0\) and solved by means of Cyclic Reduction (CR). It is shown that the matrices generated by CR are Cauchy-like with respect to a suitable singular operator and their displacement structure is explicitly determined. The application of the Gohberg-Kailath-Olshevsky algorithm provides a method for solving this Riccati equation in \(O(n^2)\) arithmetic operations (ops) with quadratic convergence. The structured doubling algorithm is analyzed in the same framework and accelerated to \(O(n^2)\) ops as well. In critical cases where convergence turns to linear, we present an adaptation of the shift technique which allows us to get rid of the singularity. Numerical experiments and comparisons which confirm the effectiveness of the new approach are reported.

MSC:

65F30 Other matrix algorithms (MSC2010)
65F05 Direct numerical methods for linear systems and matrix inversion
15A24 Matrix equations and identities