On matrix factorization and efficient least squares solution.

Astronomy and Astrophysics – Astrophysics

Scientific paper

Rate now

  [ 0.00 ] – not rated yet Voters 0   Comments 0

Details

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

Say what you really think

Search LandOfFree.com for scientists and scientific papers. Rate them and share your experience with other people.

Rating

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.

Rate now

     

Profile ID: LFWR-SCP-O-1826755

  Search
All data on this website is collected from public sources. Our data reflects the most accurate information available at the time of publication.