Introduction to mixed complementarity problems

To be able to reliably solve models with occasionally binding constraints, all equilibrium equations should be represented in RECS as mixed complementarity problems (MCP). Here, we present a short introduction to this kind of problems. For more information, see Rutherford (1995), and Ferris and Pang (1997).

Definition of a mixed complementarity problem

Complementarity problems can be seen as extensions of square systems of nonlinear equations that incorporate a mixture of equations and inequalities. Many economic problems can be expressed as complementarity problems. An MCP is defined as follows (adapted from Munson, 2002):

Definition (Mixed Complementarity Problem) Given a continuously differentiable function $F: \mathbf{R}^n\rightarrow \mathbf{R}^n$, and lower and upper bounds

$$l \in \{\mathbf{R} \cup \{-\infty\}\}^n,$$

$$u \in \{\mathbf{R} \cup \{+\infty\}\}^n.$$

The mixed complementarity problem $l \le x \le u \perp F(x)$ is to find a $x \in \mathbf{R}^n$ such that one of the following holds for each $i \in \{1,\ldots,n\}:$

$F_i(x)=0$ and $l_i\le x_i \le u_i$,

$F_i(x)>0$ and $x_i=l_i$,

$F_i(x)<0$ and $x_i=u_i$.

To summarize, an MCP problem associates each variable, $x_i$, to a lower bound, $l_i$, an upper bound, $u_i$, and an equation, $F_i(x)$. The solution is such that if $x_i$ is between its bounds then $F_i(x)=0$. If $x_i$ is equal to its lower (upper) bound then $F_i(x)$ is positive (negative).

This format encompasses several cases. In particular, it is easy to see that with infinite lower and upper bounds, solving an MCP problem is equivalent to solving a square system of nonlinear equations:

$$-\infty\le x\le +\infty \perp F(x) \Leftrightarrow F(x)=0.$$

A simple example of mixed complementarity

We consider here the traditional consumption/saving problem with borrowing constraint (Deaton, 1991). This problem is solved as a demonstration, see CS1. A consumer with utility $\sum_{t=0}^{\infty} u(C_t)/(1+\delta)^t$ has a stochastic income, $Y_t$ and has to choose each period how much to consume and how much to save. He is prevented from borrowing, but can save at the rate $r$. Without the borrowing constraint, his problem consists of choosing its consumption $C_t$ such that it satisfies the standard Euler equation:

$$u'\left(C_{t}\right)=\frac{1+r}{1+\delta}\mathrm{E}_{t}\left[u'\left(C_{t+1}\right)\right].$$

The borrowing constraint implies the inequality $C_t \le X_t$, where $X_t$ is the cash on hand available in period $t$ and defined as

$$X_{t}=\left(1+r\right)\left(X_{t-1}-C_{t-1}\right)+Y_{t}.$$

When the constraint is binding (i.e., $C_t=X_t$), the Euler equation no longer holds. The consumer would like to borrow but cannot, so its marginal utility of immediate consumption is higher than its discounted marginal utility over future consumption:

$$u'\left(C_{t}\right) \ge \frac{1+r}{1+\delta}\mathrm{E}_{t}\left[u'\left(C_{t+1}\right)\right].$$

Since the Euler equation holds when consumption is not constrained, this behavior can be summarized by the following complementarity equation

$$C_{t}\le X_{t} \quad \perp \quad \frac{1+r}{1+\delta}\mathrm{E}_{t}\left[u'\left(C_{t+1}\right)\right]\le u'\left(C_{t}\right).$$

When do complementarity problems arise?

In addition to encompassing nonlinear systems of equations, complementarity problems appear naturally in economics. Although not exhaustive, we provide here a few situations that lead to complementarity problems:

  • Karush-Kuhn-Tucker conditions of a constrained nonlinear program. The first-order conditions of the following nonlinear programming problem $\min_x f(x)$ subject to $g(x)\le 0$ and $l\le x \le u$ can be written as a system of complementarity equations: $l\le x\le u \perp f'(x)-\lambda g'(x)$ and $\lambda\ge 0 \perp -g(x)\ge 0$, where $\lambda$ is the Lagrange multiplier on the first inequality constraint.
  • Natural representation of regime-switching behaviors. Let us consider two examples. First, a variable $y$ that is determined as a function of another variable $x$ as long as $y$ is superior to a lower bound $a$. Put simply: $y=\max\left[a,\lambda\left(x\right)\right]$. In MCP, we would write this as $y\ge a \perp y\ge \lambda\left(x\right)$. Second, a system of intervention prices in which a public stock is accumulated when prices decrease below the intervention price and sold when they exceed the intervention price (see also STO2) can be represented as $S\ge 0 \perp P-P^I\ge 0$, where $S$, $P$ and $P^I$, respectively, are the public storage, the price and the intervention price.
  • A Walrasian equilibrium can be formulated as a complementarity problem (Mathiesen, 1987) by associating three sets of variables to three sets of inequalities: non-negative levels of activity are associated to zero-profit conditions, non-negative prices are associated to market clearance, and income levels are associated to income balance.

Conventions of notations adopted for representing and solving complementarity problems

To solve complementarity problems, RECS uses several solvers listed in MCP solvers. The convention adopted in most MCP solvers and used by RECS is the one used above in the MCP definition: superior or equal inequalities are associated with the lower bounds and inferior or equal inequalities are associated with the upper bounds. So, when defining your model, take care to respect this convention.

In addition, in Yaml files inequalities should not be written explicitly. As an example, let's consider these 5 equations:

$$ F(x) = 0,$$

$$ y\ge a \quad \perp \quad G\left(y\right)\ge 0,$$

$$ z\le b \quad \perp \quad H\left(z\right)\le 0,$$

$$ a\le k \le b \quad \perp \quad J\left(k\right),$$

$$ l\ge a \quad \perp \quad M\left(l\right)\le 0.$$

They would be written in a Yaml file as

  - F(x)=0 | -inf<=x<=inf
  - G(y)   |    a<=y<=inf
  - H(z)   | -inf<=z<=b
  - J(k)   |    a<=k<=b
  - -M(l)  |    a<=l<=inf

Note that it is necessary to associate lower and upper bounds with every variables, and the "perp" symbol ($\perp$) is substituted by the vertical bar (|). So if there are no finite bounds, one has to associate infinite bounds. The last equation does not respect the convention that associates lower bounds on variables with superior or equal inequality, so, when writing it in the Yaml file, the sign of the equation needs to be reversed.

References

Deaton, A. (1991). Saving and liquidity constraints. Econometrica, 59(5), 1221-1248.

Rutherford, T. F. (1995). Extension of GAMS for complementarity problems arising in applied economic analysis. Journal of Economic Dynamics and Control, 19(8), 1299-1324.

Ferris, M. C. and Pang, J. S. (1997). Engineering and economic applications of complementarity problems. Siam Review, 39(4), 669-713.

Munson, T. (2002). Algorithms and Environments for Complementarity. Unpublished PhD thesis from University of Wisconsin, Madison.

Mathiesen, L. (1987). An algorithm based on a sequence of linear complementarity problems applied to a Walrasian equilibrium model: An example. Mathematical Programming, 37, 1-18.