Homotopic residual correction processes

From MaRDI portal





The authors discuss a residual correction algorithm for the computation of the inverse or generalized inverse (in the sense of Moore-Penrose) of a matrix. The Newton iteration and some other residual correction processes for matrix inversion are initialized by means of the homotopic (continuation) method. The residual correction processes which are strongly numerically stable, allow effective parallel implementation and converge to the inverse, at least quadratically, if the initial approximation is close to the inverse. For a structured input matrix (e.g. Toeplitz, Hankel, Vandermonde, Cauchy or Pick matrices) the requirement of having a close initial approximation becomes very important. For Hermitian (or real symmetric) positive definite and indefinite matrices, the homotopic residual correction algorithm proposed here guarantees rapid convergence to the inverses, even where a not good initial approximation is available. Some additional practical convergence acceleration can be achieved by combining the homotopic algorithm with a linearly convergent iteration for the linear system. The homotopic residual correction processes do not require any pre-estimation of the smallest singular value of the input matrix. The homotopic algorithm is extended to the numerical computation of the Moore-Penrose generalized inverse and of the inverses of other type of matrices. Numerical experiments performed with general and structured matrices confirm the theory and provide the appropriate values of parameters which define the homotopic residual correction processes. The theoretical estimates for the number of homotopic residual correction steps grow in proportion to the logarithm of the condition number of the input matrix.











This page was built for publication: Homotopic residual correction processes

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q5713233)