Downloads & Free Reading Options - Results

This paper suggests a generalized Gauss-Seidel approach to sparse linear and nonlinear least squares problems. The algorithm, closely related to one given by Elfving, uses the work of Curtis, Powell, and Reid as extended by Coleman and More) to divide the variables into nondisjoint groups of structurally orthogonal columns and then projects the updated residual into each column subspace of the Jacobian in turn. In the linear case, this procedure can be viewed as an alternate ordering of the variables in the Gauss-Seidel method. Preliminary tests indicate that this leads quickly to cheap solutions of limited accuracy for linear problems, and that this approach is promising for an inexact Gauss Newton analog of the inexact Newton approach of Dembo, Eisenstat, and Steihaug.