16.01.2015 Views

GAMS — The Solver Manuals - Available Software

GAMS — The Solver Manuals - Available Software

GAMS — The Solver Manuals - Available Software

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

MILES<br />

Thomas F. Rutherford, University of Colorado<br />

Contents<br />

1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 313<br />

2 <strong>The</strong> Newton Algorithm . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 314<br />

3 Lemke’s Method with Implicit Bounds . . . . . . . . . . . . . . . . . . . . . . . . . . 316<br />

4 <strong>The</strong> Options File . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 319<br />

5 Log File Output . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 321<br />

6 Status File Output . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 323<br />

7 Termination Messages . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 329<br />

Abstract<br />

MILES is a solver for nonlinear complementarity problems and nonlinear systems of equations. This solver can be<br />

accessed indirectly through <strong>GAMS</strong>/MPSGE or <strong>GAMS</strong>/MCP. This paper documents the solution algorithm, user<br />

options, and program output. <strong>The</strong> purpose of the paper is to provide users of <strong>GAMS</strong>/MPSGE and <strong>GAMS</strong>/MCP<br />

an overview of how the MCP solver works so that they can use the program effectively.<br />

1 Introduction<br />

MILES is a Fortran program for solving nonlinear complementarity problems and nonlinear systems of equations.<br />

<strong>The</strong> solution procedure is a generalized Newton method with a backtracking line search. This code is based<br />

on an algorithm investigated by Mathiesen (1985) who proposed a modeling format and sequential method<br />

for solving economic equilibrium models. <strong>The</strong> method is closely related to algorithms proposed by Robinson<br />

(1975), Hogan (1977), Eaves (1978) and Josephy (1979). In this implementation, subproblems are solved as<br />

linear complementarity problems (LCPs), using an extension of Lemke’s almost-complementary pivoting scheme<br />

in which upper and lower bounds are represented implicitly. <strong>The</strong> linear solver employs the basis factorization<br />

package LUSOL, developed by Gill et al. (1991).<br />

<strong>The</strong> class of problems for which MILES may be applied are referred to as ”generalized” or ”mixed” complementarity<br />

problems, which is defined as follows:<br />

Given: F : R n → R n , l, u ∈ R n<br />

Find:<br />

such that<br />

z, w, v ∈ R n<br />

F (z) = w − v<br />

l ≤ z ≤ u, w ≥ 0, v ≥ 0<br />

w T (z − l) = 0 , v T (u − z) = 0.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!