Algorithm 738: a software package for unconstrained optimization using tensor methods

This paper describes a software package for finding the unconstrained minimizer of a nonlinear function of <italic>n</italic> variables. The package is intended for problems where <italic>n</italic> is not too large—say, <italic>n</italic> 100—so that the cost of storing one <italic>n</italic> × <italic>n</italic> matrix, and factoring it at each iteration, is acceptable. The software allows the user to choose between a recently developed “tensor method” for unconstrained optimization and an analogous standard method based on a quadratic model. The tensor method bases each iteration on a specially constructed fourth-order model of the objective function not significantly more expensive to form, store, or solve than the standard quadratic model. In our experience, the tensor method requires significantly fewer iterations and function evaluations to solve most unconstrained optimization problems than standard methods based on quadratic models, and also solves a somewhat wider range of problems. For these reasons, it may be a useful addition to numerical software libraries.