TY - JOUR
T1 - Convergence of Rump's method for inverting arbitrarily ill-conditioned matrices
AU - Oishi, Shin'ichi
AU - Tanabe, Kunio
AU - Ogita, Takeshi
AU - Rump, Siegfried M.
N1 - Funding Information:
This research was partially supported by Grant-in-Aid for Specially Promoted Research (No. 17002012: Establishment of Verified Numerical Computation) from the Ministry of Education, Science, Sports and Culture of Japan.
PY - 2007/8/1
Y1 - 2007/8/1
N2 - In this paper, the problem of inverting regular matrices with arbitrarily large condition number is treated in double precision defined by IEEE 754 floating point standard. In about 1984, Rump derived a method for inverting arbitrarily ill-conditioned matrices. The method requires the possibility to calculate a dot product in higher precision. Rump's method is of theoretical interest. Rump made it clear that inverting an arbitrarily ill-conditioned matrix in single or double precision does not produce meaningless numbers, but contains a lot of information in it. Rump's method uses such inverses as preconditioners. Numerical experiments exhibit that Rump's method converges rapidly for various matrices with large condition numbers. Why Rump's method is so efficient for inverting arbitrarily ill-conditioned matrices is a little mysterious. Thus, to prove its convergence is an interesting problem in numerical error analysis. In this article, a convergence theorem is presented for a variant of Rump's method.
AB - In this paper, the problem of inverting regular matrices with arbitrarily large condition number is treated in double precision defined by IEEE 754 floating point standard. In about 1984, Rump derived a method for inverting arbitrarily ill-conditioned matrices. The method requires the possibility to calculate a dot product in higher precision. Rump's method is of theoretical interest. Rump made it clear that inverting an arbitrarily ill-conditioned matrix in single or double precision does not produce meaningless numbers, but contains a lot of information in it. Rump's method uses such inverses as preconditioners. Numerical experiments exhibit that Rump's method converges rapidly for various matrices with large condition numbers. Why Rump's method is so efficient for inverting arbitrarily ill-conditioned matrices is a little mysterious. Thus, to prove its convergence is an interesting problem in numerical error analysis. In this article, a convergence theorem is presented for a variant of Rump's method.
KW - Accurate dot product
KW - Ill-conditioned matrix
KW - Matrix inversion
KW - Precondition
UR - http://www.scopus.com/inward/record.url?scp=34247358712&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=34247358712&partnerID=8YFLogxK
U2 - 10.1016/j.cam.2006.05.022
DO - 10.1016/j.cam.2006.05.022
M3 - Article
AN - SCOPUS:34247358712
SN - 0377-0427
VL - 205
SP - 533
EP - 544
JO - Journal of Computational and Applied Mathematics
JF - Journal of Computational and Applied Mathematics
IS - 1
ER -