A quadratic-tensor model algorithm for nonlinear least-squares problems with linear constraints A new algorithm is presented for solving nonlinear least squares and nonlinear equation problems. The algorithm is based on approximating the nonlinear functions using the quadratic-tensor model proposed by R. B. Schnabel and P. D. Frank [SIAM J. Numer. Anal. 21, 815-843 (1984; Zbl 0562.65029)]. The problem statement may include simple bounds or more general linear constraints on the unknowns. The algorithm uses a trust-region defined by a box containing the current values of the unknowns. The objective function (Euclidean length of the functions) is allowed to increase at intermediate steps. These increases are allowed as long as our predictor indicates that a new set of best values exists in the trust-region. There is logic provided to retreat to the current best values, should that be required. The computations for the model-problem require a constrained nonlinear least squares solver. This is done using a simpler version of the algorithm. In its present form the algorithm is effective for problems with linear constraints and dense Jacobian matrices. Results on standard test problems are presented in the Appendix. The new algorithm appears to be efficient in terms of function and Jacobian evaluations (Source: http://plato.asu.edu)

This software is also peer reviewed by journal TOMS.