ePrints@IIScePrints@IISc Home | About | Browse | Latest Additions | Advanced Search | Contact | Help

Concise row-pruning algorithm to invert a matrix

Lakshmikantham, V and Sen, SK and Sivasundaram, S (1994) Concise row-pruning algorithm to invert a matrix. In: Applied mathematics and computation, 61 (1). pp. 17-24.

[img] PDF
concise.pdf - Published Version
Restricted to Registered users only

Download (467Kb) | Request a copy
Official URL: http://dx.doi.org/10.1016/0096-3003(94)90143-0

Abstract

Presented here, in a vector formulation, is an O(mn2) direct concise algorithm that prunes/identifies the linearly dependent (ld) rows of an arbitrary m X n matrix A and computes its reflexive type minimum norm inverse A(mr)-, which will be the true inverse A-1 if A is nonsingular and the Moore-Penrose inverse A+ if A is full row-rank. The algorithm, without any additional computation, produces the projection operator P = (I - A(mr)- A) that provides a means to compute any of the solutions of the consistent linear equation Ax = b since the general solution may be expressed as x = A(mr)+b + Pz, where z is an arbitrary vector. The rank r of A will also be produced in the process. Some of the salient features of this algorithm are that (i) the algorithm is concise, (ii) the minimum norm least squares solution for consistent/inconsistent equations is readily computable when A is full row-rank (else, a minimum norm solution for consistent equations is obtainable), (iii) the algorithm identifies ld rows, if any, and reduces concerned computation and improves accuracy of the result, (iv) error-bounds for the inverse as well as the solution x for Ax = b are readily computable, (v) error-free computation of the inverse, solution vector, rank, and projection operator and its inherent parallel implementation are straightforward, (vi) it is suitable for vector (pipeline) machines, and (vii) the inverse produced by the algorithm can be used to solve under-/overdetermined linear systems.

Item Type: Journal Article
Additional Information: Copyright of this article belongs to Elsevier science.
Department/Centre: Division of Information Sciences > Supercomputer Education & Research Centre
Others
Date Deposited: 25 Mar 2011 09:43
Last Modified: 25 Mar 2011 09:43
URI: http://eprints.iisc.ernet.in/id/eprint/36299

Actions (login required)

View Item View Item