Last modified: 2007 Jan 15
Don't overlook the routines in the Run-Time Library that comes with each copy of GAUSS. QNewton() is a function implementing the BFGS method for optimization of an arbitrary function. And there's also an NLP program called SQPSolve() for optimization of functions with general linear/nonlinear, equality/inequality constraints. (Schoenberg recommends using SQPSolve() even when the problem is unconstrained.) These two both have an 'mt' version: Qnewtonmt, and Sqpsolvemt. Qnewtonmt and Sqpsolvemt use the relatively new PV and DS structures. The PV structure handles parameters, and the DS structure handles data. Aptech provides details.