## Abstract

The G-algorithm was proposed by Bareiss [1] as a method for solving the weighted linear least squares problem. It is a square root free algorithm similar to the fast Givens method except that it triangularizes a rectangular matrix a column at a time instead of one element at a time. In this paper an error analysis of the G-algorithm is presented which shows that it is as stable as any of the standard orthogonal decomposition methods for solving least squares problems. The algorithm is shown to be a competitive method for sparse least squares problems. A pivoting strategy is given for heavily weighted problems similar to that in [14] for the Householder-Golub algorithm. The strategy is prohibitively expensive, but it is not necessary for most of the least squares problems that arise in practice.

Original language | English (US) |
---|---|

Pages (from-to) | 507-520 |

Number of pages | 14 |

Journal | BIT |

Volume | 25 |

Issue number | 3 |

DOIs | |

State | Published - Sep 1 1985 |

## All Science Journal Classification (ASJC) codes

- Software
- Computer Graphics and Computer-Aided Design
- Computational Mathematics
- Applied Mathematics