Astronomy and Astrophysics – Astrophysics
Scientific paper
Apr 1995
adsabs.harvard.edu/cgi-bin/nph-data_query?bibcode=1995a%26as..110..405s&link_type=abstract
Astronomy and Astrophysics Supplement, v.110, p.405
Astronomy and Astrophysics
Astrophysics
2
Methods: Numerical, Methods: Data Analysis
Scientific paper
Least squares solution of ill conditioned normal equations by Cholesky-Banachiewicz (ChB) factorization suffers from numerical problems related to near singularity and loss of accuracy. We demonstrate that the near singularity does not arise for correctly posed statistical problems. The accuracy loss is also immaterial since for nonlinear least squares the solution by Newton Raphson iterations yields machine accuracy with no regard for accuracy of an individual iteration (Wilkinson 1963). Since this accuracy may not be achieved using singular value decomposition (SVD) without additional iterations for differential corrections and since SVD is more demanding in terms of number of operations and particularly in terms of required memory, we argue that ChB factorization remains the algorithm of choice for least squares. We present a new, very compact implementation in code of Cholesky (1924) and Banachiewicz (1938b) factorization in an elegant form proposed by Banachiewicz (1942). Source listing of the code is provided. We point out that in the same publication Banachiewicz (1938) discovered LU factorization of square matrices before Crout (1941) and rediscovered factorization of the symmetric matrices after Cholesky (1924). Since the two algorithms became confused, no due credit is given to Banachiewicz in modern literature.
No associations
LandOfFree
On matrix factorization and efficient least squares solution. does not yet have a rating. At this time, there are no reviews or comments for this scientific paper.
If you have personal experience with On matrix factorization and efficient least squares solution., we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and On matrix factorization and efficient least squares solution. will most certainly appreciate the feedback.
Profile ID: LFWR-SCP-O-1826755