\documentclass[reqno]{amsart}
\usepackage{hyperref}

\AtBeginDocument{{\noindent\small
\emph{Electronic Journal of Differential Equations},
Vol. 2010(2010), No. 77, pp. 1--12.\newline
ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu
\newline ftp ejde.math.txstate.edu}
\thanks{\copyright 2010 Texas State University - San Marcos.}
\vspace{9mm}}

\begin{document}
\title[\hfilneg EJDE-2010/77\hfil Infinity Laplace equation]
{Infinity Laplace equation with non-trivial right-hand side}

\author[G. Lu, P. Wang\hfil EJDE-2010/77\hfilneg]
{Guozhen Lu, Peiyong Wang}  % in alphabetical order

\address{Guozhen Lu \newline
 Department of Mathematics\\
 Wayne State University\\
 656 W. Kirby, 1150 FAB\\
 Detroit, MI 48202, USA}
\email{gzlu@math.wayne.edu}

\address{Peiyong Wang \newline
 Department of Mathematics\\
 Wayne State University\\
 656 W. Kirby, 1150 FAB\\
 Detroit, MI 48202, USA}
\email{pywang@math.wayne.edu}


\thanks{Submitted July 2, 2009. Published June 8, 2010.}
\thanks{G. Lu is partially supported by US NSF grant DMS0901761.}
\subjclass[2000]{35J70, 35B35}
\keywords{Infinity Laplace equation; inhomogeneous equation;
\hfill\break\indent  viscosity solutions; least solution;
greatest solution; strict comparison principle;
\hfill\break\indent existence; uniqueness; local Lipschitz continuity}

\begin{abstract}
 We analyze the set of continuous viscosity solutions of
 the infinity Laplace equation $-\Delta^N_{\infty}w(x) = f(x)$,
 with generally sign-changing right-hand side in a bounded domain.
 The existence of a least and a greatest continuous viscosity
 solutions, up to the boundary, is proved through a Perron's
 construction by means of a strict comparison principle. These
 extremal solutions are proved to be absolutely extremal solutions.
\end{abstract}

\maketitle
\numberwithin{equation}{section}
\newtheorem{theorem}{Theorem}[section]
\newtheorem{lemma}[theorem]{Lemma}

\section{Introduction}

In this article, we consider a nonlinear differential operator
known as the normalized infinity Laplacian and symbolically defined as
\begin{equation}\label{infinitylaplacian}
-\Delta^N_{\infty}w(x) = -\frac{1}{|\nabla w(x)|^2}
\sum^n_{i,j=1}\partial_{x_i}w(x)\partial_{x_j}w(x)
\partial^2_{ x_ix_j}w(x),
\end{equation}
which is abbreviated as
\begin{equation}\label{abbreviation}
-\Delta^N_{\infty}w(x) = -\frac{1}{|\nabla w(x)|^2}
\langle D^2w(x) \nabla w(x), \nabla w(x)\rangle.
\end{equation}
Here $\langle\cdot,\cdot\rangle$ denotes the inner product in
the Euclidean space $\Re^n$. The expression $\langle D^2 w(x) \nabla w(x), \nabla
w(x)\rangle $ stands for the un-normalized infinity Laplacian of $w$ at
$x$, sometimes denoted by $\Delta_{\infty}w(x)$.

We assume $\Omega\in\Re^n$ is a bounded open set
and consider the  boundary-value problem
\begin{equation}\label{dirichlet}
 \begin{gathered}
-\Delta^N_{\infty}w(x) = f(x) \quad \text{in }\Omega \\
w(x) = g(x) \quad \text{on }\partial\Omega.
\end{gathered}
\end{equation}
Here we assume that $f\in C(\Omega)$ and $g\in C(\partial\Omega)$.
Next, we make clear the meaning of $\Delta^N_{\infty}w(x)$. For
a twice differentiable function $\varphi$, the normalized infinity
Laplacian is
\begin{equation}\label{normalized}
\Delta^N_{\infty}\varphi(x_0)
= \begin{cases}
\frac{1}{|\nabla \varphi(x_0)|^2}\langle D^2\varphi(x_0)\nabla\varphi
(x_0), \nabla \varphi(x_0)\rangle &\text{if }\ \nabla
\varphi(x_0)\neq 0 \\[3pt]
[\lambda_{\rm min}(D^2\varphi(x_0)), \lambda_{\rm max}(D^2\varphi
(x_0))] &\text{if } \nabla \varphi(x_0) = 0,
\end{cases}
\end{equation}
where $\lambda_{\rm min}(M)$ and $\lambda_{\rm max}(M)$ denote
respectively the least and greatest eigenvalues of a square matrix
$M$. Another pair of symbols are $\Delta^+_{\infty}
\varphi$ and $\Delta^-_{\infty}\varphi$ which are given
equivalently by $\Delta^+_{\infty}\varphi(x) =
\Delta^-_{\infty}\varphi(x) = \Delta^N_{\infty}
\varphi$ if $\nabla\varphi(x)\neq 0$, and $\Delta^+_{\infty}
\varphi(x) = \lambda_{\rm max}(D^2\varphi(x))$ and
$\Delta^-_{\infty}\varphi(x) = \lambda_{\rm min}(D^2(\varphi(x)))$
if $\nabla\varphi(x) = 0$. In case $\nabla \varphi(x) = 0$,
$\Delta^N_{\infty}\varphi(x) \geq f(x)$ means that
$\lambda_{\rm max}(D^2\varphi(x)) \geq f(x)$ and
$\Delta^N_{\infty}\varphi(x) \leq f(x)$ means that
$\lambda_{\rm min}(D^2\varphi(x)) \leq f(x)$. Equivalently,
$\Delta^N_{\infty}\varphi(x) \geq f(x)$ means that
$\Delta^+_{\infty}\varphi(x) \geq f(x)$ and
$\Delta^N_{\infty}\varphi(x) \leq f(x)$ means that
$\Delta^-_{\infty}\varphi(x) \leq f(x)$. For a detailed
explanation of this definition, we refer to \cite{LW2}.

An upper semi-continuous function ,$u\in USC(\Omega)$, is
\textit{a viscosity sub-solution}
of the infinity Laplace equation
\begin{equation}\label{infinitylaplace}
-\Delta^N_{\infty}u(x) = f(x)
\end{equation}
if the condition $u\prec_{x_0}\varphi$ for $x_0\in\Omega$ and
$\varphi\in C^2 (\Omega)$ implies
$-\Delta^N_{\infty}\varphi(x_0)\leq f(x_0)$. Here
$USC(\Omega)$ and $LSC(\Omega)$ denote the sets of upper
semi-continuous and lower semi-continuous functions in $\Omega$
respectively, and $u\prec_{x_0}\varphi$ means $u - \varphi$
attains a local maximum at $x_0$. Similarly, \textit{a viscosity
super-solution} of \eqref{infinitylaplace} is a function $u\in
LSC(\Omega)$ which satisfies the condition that
$\varphi\prec_{x_0}u$ for $x_0\in\Omega$ and $\varphi\in C^2
(\Omega)$ implies $-\Delta^N_{\infty}\varphi(x_0)\geq
f(x_0)$. \textit{A viscosity solution} of \eqref{infinitylaplace}
is both a viscosity sub-solution and super-solution.

The study of the infinity Laplace equation
$\Delta_{\infty} u(x) = 0$ was initiated in the 1960s by
Aronsson in \cite{A1,A2,A3}, where he deduced the
infinity Laplace equation $\Delta_{\infty}u(x) = 0$ as the
Euler-Aronsson equation for smooth absolute minimizers. Partly due
to the lack of a proper notion of solutions of the highly
degenerate nonlinear infinity Laplace equation, the study had been
dormant for quite a while until the introduction of viscosity
solutions by  Evans, Crandall, Ishii,Lions, et al
(see \cite{CIL} and the references therein). The
existence of a solution of the equation was proved by Bhattacharya,
DiBenedetto and Manfredi in \cite{BDM}.
Jensen presented the first proof of the uniqueness of a
viscosity solution of the Dirichlet problem for the homogeneous
infinity Laplace equation in a bounded domain in 1993 in \cite{J},
which revived the study of the infinity Laplacian. Since then, the
Dirichlet problem for the infinity Laplace equation has received
extensive attention. The works
\cite{Ju1,CE,CEG,BB,C,E2,LM1,BJW,ACJ,CW,CGW}
give a partial list of the references in
the literature. Among them, \cite{BB} contains a second proof of
the uniqueness of a viscosity solutions of the Dirichlet problem
for the homogeneous infinity Laplacian in a bounded domain. A
third uniqueness proof is given in \cite{CGW} which works for
unbounded domains. Meanwhile, the study of the eigenvalue problem
for the infinity Laplacian and the evolution problem for the
infinity Laplacian were also taken up
(see \cite{JLM,CW,JK,Ju2}). The authors of the current
paper investigated in
the well-posedness of the inhomogeneous problems
$\Delta_{\infty}u(x) = f(x)$ and
$\Delta^N_{\infty}u(x) = f(x)$, with $f(x) > 0$, in
\cite{LW1,LW2}, where the existence and uniqueness of a viscosity solution
of the Dirichlet problem are proved. Peres-Schramm-Sheffield-Wilson
provided interpretation of the normalized infinity Laplacian from
the point of view of the differential game theory in
\cite{PSSW}. Quoted from \cite{PSSW},
the continuum value of a differential game called the
``tug-of-war" verifies the inhomogeneous infinity Laplace equation
$-\Delta^N_{\infty}u(x) = 2f(x)$, where $f$ is the running
payoff function which satisfies $\inf f > 0$ in the domain. A
counter-example was also provided in \cite{PSSW} to show the
uniqueness of a viscosity solution of the Dirichlet problem for
the inhomogeneous equation fails if $f$ could change sign. It is
unclear what one can say about the multiple viscosity solutions of
the Dirichlet problem \eqref{dirichlet} for a general ``payoff"
function $f$, though. The theme of this paper is to answer at
least partially this question. In fact, we prove that there always
exist continuous viscosity solutions of the Dirichlet problem
\eqref{dirichlet} for the normalized infinity Laplacian and for
any continuous right-hand-side $f$ (Theorem \ref{existence}).
Moreover, the greatest and least viscosity solutions are
constructed (Theorem \ref{existence}) through the Perron's method
combined with a strict comparison theorem (Theorem
\ref{comparison2}).

This article is organized as follows. The second section is devoted
to the derivation of the local Lipschitz continuity of a viscosity
sub-solution (Lemma \ref{lemma1.1}) and a strict comparison
principle (Theorem \ref{comparison2}). The third section contains
the construction of the least and the greatest solutions,
i.\,e.\,the main theorem (Theorem \ref{existence}). The last section
contains closely related problems yet to be solved.

In this article, especially when the inhomogeneous term $f$ is not
continuous in its arguments, the strict differential inequality
$-\Delta^N_{ \infty}w(x) < f(x, \nabla w(x))$ in $\Omega$
in the viscosity sense is understood in \textit{the locally uniform sense}
that for any $x_0\in \Omega$, there exist a neighborhood $N$ of
$x_0$ in $\Omega$ and a $\delta
> 0$ such that $-\Delta^N_{\infty}w(x) \leq f(x, \nabla
w(x)) - \delta$ in the viscosity sense in $N$. The differential
inequality $-\Delta^N_{ \infty}w(x) > f(x, \nabla w(x))$
is similarly understood.

We recall \cite[Lemma 1.10]{LW2}, the proof of which may also be
found therein.

\begin{lemma}\label{lemma0}
Assume $\Omega$ is an open subset of $\Re^n$ and $f\in C(\Omega)$.
$\Lambda$ is an index set.

(a) Suppose $u(x) = \sup_{\lambda\in\Lambda}u_{\lambda}(x) < \infty$,
$x\in\Omega$, where $-\Delta^N_{\infty}u_{\lambda}\leq f$ in
$\Omega$ in the viscosity sense for every $\lambda\in\Lambda$. If
$u\in C(\Omega)$, then $-\Delta^N_{\infty}u\leq f$ in $\Omega$
in the viscosity sense.

(b) Similarly, if $u(x) = \inf_{\lambda\in\Lambda}u_{\lambda}(x) >
-\infty$, $x\in\Omega$, where
$-\Delta^N_{\infty}u_{\lambda}\geq f$ in $\Omega$ in the
viscosity sense for every $\lambda\in\Lambda$. Then $u\in
C(\Omega)$ implies that $-\Delta^N_{\infty}u\geq f$ in
$\Omega$ in the viscosity sense.
\end{lemma}

A similar result holds for the infinity Laplace equation
$-\Delta_{\infty}u = f$, the proof of which is simpler as the
singularity caused by $\nabla u = 0$ does not present in this case.

\section{A Comparison Theorem}

For a nonzero vector $x$, $\hat{x} = x/|x|$ denotes its
normalized vector. The notation $C_b(\Omega)$ denotes the set of
bounded continuous functions defined in $\Omega$. For two sets $V$
and $U$, $V\subset\subset U$ means that $V$ is compactly contained
in $U$.

We start out to prove a lifting lemma stated as follows.

\begin{lemma}\label{liftlemma}
If $u\in USC(\Omega)$ is a viscosity sub-solution of $\Delta_{\infty}u
= k_1$ in $\Omega$, and $v\in C^2(\Sigma)$ verifies $\Delta_{\infty}v \geq k_2$
in $\Sigma$, for constants $k_1$ and $k_2$, then the function $w: (x,y)\mapsto
u(x)+v(y)$ is a viscosity sub-solution of $\Delta_{\infty}w(x,y) =
k_1 + k_2$ in $\Omega\times\Sigma$.
\end{lemma}

\begin{proof}
It suffices to prove $\Delta_{\infty}\varphi(x_0,y_0)\geq k_1 + k_2$ for
any $\varphi\in C^2(\Omega)$ and $(x_0,y_0)\in\Omega$ such that $u(x) + v(y)
\prec_{(x_0, y_0)}\varphi(x,y)$. Without the loss of generality, we may
assume $(x_0,y_0) = (0,0)$, $u(0) = 0$, $v(0) = 0$, $\varphi(0,0) = 0$, and
$\varphi$ is a quadratic polynomial.
Denote
$$\nabla \varphi(0,0) = \begin{pmatrix}
\varphi_{x}(0,0)\\
\varphi_{y}(0,0)\end{pmatrix}
$$
and
$$
D^2\varphi(0,0) = \begin{pmatrix}
\varphi_{xx}(0,0) &\varphi_{xy}(0,0) \\
\varphi_{yx}(0,0) &\varphi_{yy}(0,0)\end{pmatrix}
$$
Then
\begin{equation} \label{(1)}
\begin{aligned}
u(x) + v(y)
&\leq \varphi_x(0,0)\cdot x + \varphi_y(0,0)\cdot y
+ \frac{1}{2}\langle \varphi_{xx}(0,0)x,x\rangle  \\
&\quad + \langle \varphi_{xy}(0,0)x,y\rangle +
\frac{1}{2}\langle \varphi_{yy}(0,0)y,y\rangle.
\end{aligned}
\end{equation}
We write
$$
v(y) = \nabla v(0)\cdot y + \frac{1}{2}\langle D^2v(0)y,y\rangle +
\circ(|y|^2).
$$
 Replacing this in  \eqref{(1)}, we obtain
\begin{align*}
& u(x) + \nabla v(0)\cdot y + \frac{1}{2}\langle D^2v(0)y,y\rangle
+ \circ (|y|^2)\\
&\leq  \varphi_x(0,0)\cdot x + \varphi_y(0,0)\cdot y
+ \frac{1}{2}\langle \varphi_{xx}(0,0)x,x\rangle \\
&\quad +  \langle \varphi_{xy}(0,0)x,y\rangle +
\frac{1}{2}\langle \varphi_{yy}(0,0)y,y\rangle
\end{align*}
or equivalently
\begin{align*}
u(x) &\leq  \varphi_x(0,0)\cdot x + (\varphi_y(0,0) - \nabla v(0))\cdot
y +  \frac{1}{2}\langle \varphi_{xx}(0,0)x, x\rangle \\
&\quad + \langle \varphi_{xy}(0,0)x,y\rangle
+\frac{1}{2}\langle (\varphi_{yy}(0,0) - D^2v(0))y, y\rangle
+ \circ(|y|^2)
\end{align*}
 for any small $x$ and $y$.

It is clear that $\varphi_y(0,0) = \nabla v(0)$ and $\varphi_{yy}(0,0)
- D^2v(0) \geq 0$. Denote $B = \varphi_{yy}(0,0) - D^2v(0)$. Then
\begin{equation} \label{(2)}
\begin{aligned}
u(x)&\leq \varphi_x(0,0)\cdot x + \frac{1}{2}\langle
\varphi_{xx}(0,0)x, x\rangle\\
&\quad + \langle \varphi_{xy}(0,0)x,y\rangle
+ \frac{1}{2}\langle By,y\rangle + \circ(|y|^2).
\end{aligned}
\end{equation}
First, we assume the matrix $B$ is invertible. So $B = A^2$ for some
symmetric invertible matrix $A$.
Then the right-hand-side of
\eqref{(2)} is equal to
\begin{equation}
\begin{aligned}
&\varphi_x(0,0)\cdot x + \frac{1}{2}\langle \varphi_{xx}(0,0)x,x\rangle
+ \langle A^{-1} \varphi_{xy}(0,0)x, Ay\rangle  \\
&+ \frac{1}{2}\langle Ay, Ay\rangle + \circ(|y|^2) \\
&= \varphi_x(0,0)\cdot x + \frac{1}{2}\langle (\varphi_{xx}(0,0)
- \varphi_{xy}(0,0) B^{-1}\varphi_{xy}(0,0))x, x\rangle  \\
&\quad + \frac{1}{2}|A^{-1}\varphi_{xy}(0,0)x + Ay|^2 + \circ(|y|^2),
\end{aligned}
\end{equation}
where $x$ and $y$ are any small vectors. Take
$y = -B^{-1}\varphi_{xy}(0,0)x$
for each small $x$. Then
$$
u(x)\leq \varphi_x(0,0)\cdot x + \frac{1}{2}\langle (\varphi_{xx}(0,0)
- \varphi_{xy}(0,0)
B^{-1}\varphi_{xy}(0,0))x, x\rangle + \circ (|x|^2)
$$
 for all small vector $x$. Therefore,
on account of the fact that $\Delta_{\infty}u \geq k_1$ in $\Omega$, we obtain
\begin{equation} \label{(3)}
\langle (\varphi_{xx}(0,0) - \varphi_{xy}(0,0)B^{-1}\varphi_{xy}(0,0))\varphi_x(0,0),
\varphi_x(0,0)\rangle \geq k_1.
\end{equation}
As a result, the following equalities and inequalities hold at $(0,0)$:
\begin{align*}
&\langle \varphi_{xx}\varphi_x,\varphi_x\rangle
 + 2\langle \varphi_{xy}\varphi_x,\varphi_y\rangle
 + \langle \varphi_{yy}\varphi_y, \varphi_y\rangle  \\
&= \langle \varphi_{xx}\varphi_x,\varphi_x\rangle
 + 2\langle \varphi_{xy}\varphi_x,\varphi_y\rangle
 + \langle B\varphi_y, \varphi_y\rangle
 + \langle D^2v\nabla v, \nabla v\rangle \\
&\geq \langle \varphi_{xx}\varphi_x,\varphi_x\rangle
 + 2\langle \varphi_{xy}\varphi_x,\varphi_y\rangle
 + \langle B\varphi_y, \varphi_y\rangle +  k_2 \\
&= \langle \varphi_{xx}\varphi_x, \varphi_x\rangle
 + 2\langle A^{-1}\varphi_{xy}\varphi_x, A\varphi_y\rangle
 + \langle A\varphi_y, A\varphi_y\rangle + k_2\\
&= \langle (\varphi_{xx}-\varphi_{xy}B^{-1}\varphi_{xy})\varphi_x, \varphi_x\rangle + |A^{-1}
\varphi_{xy}\varphi_x + A\varphi_y|^2 + k_2 \\
&\geq  k_1 + k_2,
\end{align*}
 according to \eqref{(3)}.
In general, when $B$ is not invertible, we define
$B^{\varepsilon} = B + \varepsilon I$
for every small $\varepsilon > 0$. Then $B^{\varepsilon}$ is
invertible and the
inequalities \eqref{(2)} and \eqref{(3)} still hold with $B$
replaced by $B^{
\varepsilon}$. Let $B^{\varepsilon} = A^2$ for a positive definite
matrix $A$.
In the end, we have, at $(0,0)$,
\begin{align*}
&\langle \varphi_{xx}\varphi_x,\varphi_x\rangle
 + 2\langle \varphi_{xy}\varphi_x,\varphi_y\rangle
 + \langle \varphi_{yy}\varphi_y, \varphi_y\rangle
 +  \varepsilon |\varphi_y|^2\\
&= \langle \varphi_{xx}\varphi_x,\varphi_x\rangle
 + 2\langle \varphi_{xy}\varphi_x,\varphi_y\rangle
 + \langle B^{\varepsilon}\varphi_y, \varphi_y\rangle
 + \langle D^2v\nabla v, \nabla v\rangle  \\
&\geq \langle \varphi_{xx}\varphi_x,\varphi_x\rangle
 + 2\langle \varphi_{xy}\varphi_x,\varphi_y\rangle
 + \langle B^{\varepsilon}\varphi_y, \varphi_y\rangle +  k_2\\
&= \langle \varphi_{xx}\varphi_x, \varphi_x\rangle
 + 2<A^{-1}\varphi_{xy}\varphi_x, A\varphi_y\rangle
 + \langle A\varphi_y, A\varphi_y\rangle +  k_2\\
&= \langle (\varphi_{xx}-\varphi_{xy}(B^{\varepsilon})^{-1}
 \varphi_{xy})\varphi_x, \varphi_x\rangle + |A^{-1}
 \varphi_{xy}\varphi_x + A\varphi_y|^2 +  k_2\\
&\geq  k_1 + k_2,
\end{align*}
 for every small $\varepsilon > 0$.
Then $\langle \varphi_{xx}\varphi_x,\varphi_x\rangle
+ \langle 2\varphi_{xy}\varphi_x,\varphi_y\rangle
+ \langle \varphi_{yy}\varphi_y, \varphi_y\rangle
 \geq k_1 + k_2$. The proof is complete.
\end{proof}


\begin{lemma}\label{lemma1.1}
Suppose $f\in C(\Omega)$, and $w\in USC(\Omega)$ is locally
bounded.
\begin{itemize}
\item[(a)] If $-\Delta_{\infty}w(x) \leq f(x)$ in $\Omega$, then
$w: \Omega\to \Re$ is locally Lipschitz continuous.

\item[(b)] If $-\Delta^N_{\infty}w(x) \leq f(x)$ in $\Omega$,
then $w: \Omega\to \Re$ is locally Lipschitz continuous.
\end{itemize}
Furthermore, the Lipschitz constant of $w$ over $\Omega'\subset
\subset\Omega$ may be taken as $C(1+\|w\|_{L^{\infty}(\tilde{
\Omega})})$, where $\Omega'\subset\subset \tilde{\Omega} \subset
\subset \Omega$ and $C$ depends on $\|f\|_{L^{\infty}(\Omega')}$.
\end{lemma}

\begin{proof}
Assume $f(x)\leq K$, $x\in\Omega'\subset\subset\Omega$. Without
loss of generality, we also assume $\|w\|_{L^{\infty}(\Omega)} <
\infty$.

(a) Define $u(x,y) = w(x)
+ Cy^{4/3}$, for $x\in\Omega$, $1\leq y\leq 2$. We notice that
$Cy^{4/3}$ is a $C^2$ solution of the equation
$\Delta_{\infty}w(y) = \frac{64}{81}C^3$ for $y\neq 0$.
Then the preceding lift
lemma \eqref{liftlemma} implies that $u$ is an infinity sub-harmonic
function in $\Omega\times (1, 2)$, if $C$ is sufficiently large.
A well-known fact about
semi-continuous infinity sub-harmonic functions (see, for example,
\cite[Lemma 2.9]{ACJ} for continuous functions, and \cite{CEG},
or \cite{LM2} for semi-continuous functions) states that $u$ is
Lipschitz continuous on $\Omega' \times [1.1, 1.9]$ with some
Lipschitz constant $L = L(\|w\|_{ L^{\infty}(\Omega)})$. As a
result, $w$ is Lipschitz continuous on $\Omega'$ with Lipschitz
constant $L$.

(b) Clearly, that $-\Delta^N_{\infty}w(x)\leq f(x)$ in
$\Omega$ in the viscosity sense implies that $-\Delta_{
\infty}w(x) \leq |\nabla w(x)|^2f(x)$ in the viscosity sense in
$\Omega$ (but not the converse). Without the loss of generality,
we assume that $w > 0$ in $\Omega$. Take $\lambda > 0$ small so
that $\lambda\|w\|_{L^{\infty}(\Omega)} < \frac{1}{2}$. Define $u
= G(w) = w + \frac{ \lambda}{2}w^2$ in $\Omega$. For simplicity,
we assume $w$ is $C^2$. All steps of the following computation can
be made rigorous by means of viscosity solutions. We leave the
details to the reader. Then $G'(w) = 1 + \lambda w$ and $G''(w) =
\lambda$. In particular, $2 > G'(w) > \frac{1}{2}$. Moreover,
$\nabla u = G'(w)\nabla w$, $D^2u = G'(w)D^2w + G''(w) \nabla
w\otimes \nabla w$, and
\begin{align*}
-\Delta_{\infty}u &= -(G'(w))^3\Delta_{\infty} w
- (G'(w))^2G''(w)|\nabla w|^4 \\
&\leq (G'(w))^3\{f(x) - \frac{G''(w)}{G'(w)}|\nabla w|^2\}
|\nabla w|^2 \\
&= \frac{(1+\lambda w)^4}{\lambda}\{f(x) - \frac{\lambda}{1 +
\lambda w}|\nabla w|^2\}\frac{\lambda |\nabla w|^2}{1+\lambda w}\\
&\leq \frac{(1+\lambda w)^4}{4\lambda}(f(x))^2, \quad
\text{(due to the Cauchy-Schwarz inequality)} \\
&< \frac{4K^2}{\lambda}
\end{align*}
By (a), one deduces that $u$ is locally Lipschitz continuous in
$\Omega$. So $w = \frac{2u}{1+\sqrt{1+2\lambda u}}$ is also
locally Lipschitz continuous in $\Omega$.
\end{proof}

The following comparison theorem is a generalization of a strict
comparison principle stated in \cite[Theorem 3.1]{LW2}.

\begin{theorem}\label{comparison}
Assume $f\in C(\Omega\times\Re^n)$, and the modulus of continuity
of the function $x\mapsto f(x,p)$ is independent of $p\in\Re^n$.
Suppose $u_j\in C(\Omega)$, $j = 1, 2$, verify in the viscosity
sense either
$$
\Delta^N_{\infty}u_1(x) < f(x, \nabla u_1)\quad \text{and}\quad
\Delta^N_{\infty}u_2(x) \geq f(x, \nabla u_2)
$$
or
$$\Delta^N_{\infty}u_1(x) \leq f(x, \nabla u_1)\quad \text{and}\quad
 \Delta^N_{\infty}u_2(x) > f(x, \nabla u_2)
$$
in $\Omega$.

If $\limsup_{x\in\Omega\to z}(u_2(x) - u_1(x))\leq 0$ for
any $z\in\partial\Omega$, then $u_2(x) \leq u_1(x)$ in $\Omega$.
\end{theorem}

\begin{proof}
One may follow the proof of \cite[Theorem 3.1]{LW2} which can be
simplified substantially with the application of a fourth order
penalty function $w_{\varepsilon}(x,y) = u_2(x) - u_1(y) -
\frac{1}{4 \varepsilon}|x - y|^4$, $(x,y)\in\Omega\times\Omega$,
used in \cite{CGG} and \cite{JK}. We leave the details to the
reader.
\end{proof}

The preceding comparison theorem and the lemma \eqref{lemma1.1}
imply the following theorem immediately.

\begin{theorem}\label{comparison2}
Assume $f\in C(\Omega)$. Suppose $u_1\in LSC(\Omega)$, $u_2\in
USC(\Omega)$, and they verify either
$$
\Delta^N_{\infty}u_1(x) < f(x)\quad \text{and}\quad
\Delta^N_{\infty}u_2(x) \geq f(x)
$$
or
$$
\Delta^N_{\infty}u_1(x) \leq f(x)\quad \text{and}\quad
 \Delta^N_{\infty}u_2(x) > f(x)
$$
in the viscosity sense in $\Omega$.

If $\limsup_{x\in\Omega\to z}(u_2(x) - u_1(x))\leq 0$ for
any $z\in\partial\Omega$, then $u_2(x) \leq u_1(x)$ in $\Omega$.
\end{theorem}

\section{Continuous Solutions Of The Dirichlet Problem }

There are different approaches to the existence of a viscosity
solution of the boundary value problem \eqref{dirichlet}. The
approach used here is the Perron's method combined with a delicate
albeit elementary analysis which depends essentially on the strict
comparison theorem \eqref{comparison2}.

For $f\in C(\Omega)$ and $g\in C(\partial\Omega)$, we define the set of
strict super-solutions
\begin{equation}
\mathcal{A}^+_{f,g} = \{v\in C(\bar{\Omega}):
-\Delta^N_{\infty} v(x)
> f(x)\text{ in } \Omega, \text{ and } v\geq g \text{ on }\
\partial\Omega\}
\end{equation}
and the set of strict sub-solutions
\begin{equation}
\mathcal{A}^-_{f,g} = \{v\in C(\bar{\Omega}):
-\Delta^N_{\infty} v(x) < f(x) \text{ in } \Omega,
\text{ and } v\leq g \text{ on }\partial\Omega\}.
\end{equation}
Whenever there is little confusion, we will write $\mathcal{A}^+$
and $\mathcal{A}^-$ for $\mathcal{A}^+_{f,g}$ and
$\mathcal{A}^-_{f,g}$, respectively.
Obviously, $\mathcal{A}^+$ and $\mathcal{A}^-$ are both nonempty.

Define $w^+(x) = \inf_{v\in\mathcal{A}^+}v(x)$ and
$w^-(x) = \sup_{v\in\mathcal{ A}^-}v(x)$, $x\in\bar{\Omega}$.
By definition, $w^+$ is upper semi-continuous and $w^-$ is
lower semi-continuous on $\bar{\Omega}$. Obviously, $w^+$
and $w^-$ are both bounded on
$\bar{\Omega}$, and $w^-(x)\leq g(x) \leq w^+(x)$ on
$\partial\Omega$ according to the preceding comparison theorem
\eqref{comparison2}.

Take a super-solution $\phi$ in $\mathcal{A}^+_{f,g}$. For example,
$\phi(x) = - C|x-z|^2 + D$ for suitable $C$ and $D$. Define
\begin{equation}
\mathcal{A}^+_{f,g,\phi} = \{\min(v,\phi): v\in\mathcal{A}^+_{f,g}\}.
\end{equation}
Clearly, $\mathcal{A}^+_{f,g,\phi}\subset\mathcal{A}^+_{f,g}$ and
$w^+(x) = \inf_{v\in\mathcal{A}^+_{f,g,\phi}}v(x)$. For every $v$
in $\mathcal{A}^+_{f,g,\phi}$, Lemma \eqref{lemma1.1} says $v$ is
locally Lipschitz continuous with Lipschitz constant $\leq C(1 +
\|\phi\|_{L^{\infty}(\Omega)})$, i.\,e.\,$\mathcal{A}^+_{f,g,\phi}$
is locally Lipschitz equi-continuous.

On the other hand, one may pick a sequence $\{v_k\}$ in
$\mathcal{A}^+_{f,g,\phi}$ such that $v_k$ converges to $w^+$
on a countable dense
subset $E$ of $\bar{\Omega}$. Define $\tilde{v}_k = \min\{v_1, v_2, \dots,
v_k\}$, $k = 1, 2, \dots$. Then $\tilde{v}_k\in\mathcal{A}^+_{f,g,\phi}$
and $\tilde{v}_k$ converges to $w^+$ on $E$.
Replacing $v_k$ by $\tilde{v}_k$,
one may assume that $v_k\geq v_{k+1}$ for all $k$.
Consequently, a subsequence
of $\{v_k\}$, which will still be denoted by $\{v_k\}$, converges to
some $v$ locally uniformly on $\Omega$. Then
$v\in C(\Omega)$ and
$v_k\geq v\geq w^+$ on $\Omega$. Clearly, $v = w^+$ on $E\cap\Omega$.
As $w^+$ is
upper semi-continuous on $\bar{\Omega}$, for any $x\in\Omega$,
\begin{equation}
w^+(x)\geq \limsup_{z\in E\to x}w^+(z) = \limsup_{z\in E \to
x}v(z) = v(x).
\end{equation}
So $w^+ = v$ on $\Omega$ and whence $\{v_k\}$ converges to $w^+$
locally uniformly on $\Omega$. Therefore $w^+\in C(\Omega)$ and
$-\Delta^N_{\infty}w^+(x)\geq f(x)$ on $\Omega$ on account of
Lemma \eqref{lemma0}. Similarly,
$w^-\in C(\Omega)$ and $-\Delta^N_{\infty}w^-(x)\leq f(x)$ on
$\Omega$.

Next, we show $w^+ = w^- = g$ on $\partial\Omega$. For any $z\in
\partial\Omega$ and any $\varepsilon > 0$, there exists $r > 0$
such that $g(x) \leq g(z) + \varepsilon$ for all $x\in\partial
\Omega$ with $|x-z| \leq r$. Take $C > \|f\|_{L^{\infty}(\Omega)}$
and $D\geq C diam(\Omega) + 2\|g\|_{L^{\infty}(\partial\Omega)}$.
Define $v\in C(\bar{\Omega})$ by
\begin{equation}
v(x) = g(z) + \varepsilon - C|x-z|^2 + D|x-z|.
\end{equation}
Then $-\Delta^N_{\infty}v(x) = 2C > f(x)$ for
$x\in\Omega$. For $x\in\partial\Omega$ with $|x-z|\leq r$, $v(x)
\geq g(z) + \varepsilon + |x-z|\{D - Cr\} \geq g(z) + \varepsilon
\geq g(x)$, while for $x\in\partial\Omega$ with $|x-z| > r$, $v(x)
\geq g(z) + \varepsilon + r\{D - C \mathop{\rm diam}(\Omega)\} \geq g(z) +
\varepsilon + 2\|g\|_{L^{\infty}(\partial\Omega)} \geq g(x)$. So
$v\in\mathcal{A}^+_{f,g}$. As a result, $w^+(z)\leq v(z) = g(z) +
\varepsilon$, for all $\varepsilon > 0$. So $w^+ = g$ on $\partial
\Omega$. Similarly $w^- = g$ on $\partial\Omega$.

Since $w^+$ is upper semi-continuous and $w^-$ is lower
semi-continuous on $\bar{\Omega}$, for any $z\in\partial\Omega$,
the following inequalities hold
\begin{equation}
g(z) = w^+(z)\geq\limsup_{x\in\Omega\to z}w^+(x) \geq
\liminf_{x\in\Omega\to z}w^+(x) \geq \liminf_{x\in\Omega
\to z}w^-(x) \geq w^-(z) = g(z).
\end{equation}
Consequently, all the above inequalities are indeed equalities. So
\begin{equation}
\lim_{x\in\Omega\to z}w^+(x) = w^+(z) = g(z).
\end{equation}
Similarly, one obtains
\begin{equation}
\lim_{x\in\Omega\to z}w^-(x) = w^-(z).
\end{equation}
Therefore $w^+$ and $w^-$ are in $C(\bar{\Omega})$.

We now show that $-\Delta^N_{\infty}w^+(x) = f(x)$ and
$-\Delta^N_{\infty}w^-(x) = f(x)$ in $\Omega$. We need
only prove that $-\Delta^N_{\infty}w^+(x) \leq f(x)$ in
$\Omega$. Suppose the contrary that there exist a $C^2$ function
$\varphi$ and a point $x_0\in\Omega$ such that
$w^+\prec_{x_0}\varphi$ and $\Delta^+_{\infty}\varphi
(x_0) < f(x_0)$.

For any small $\varepsilon > 0$, we define
\begin{equation}
\varphi_{\varepsilon}(x) = \varphi(x_0) + \nabla \varphi(x_0)\cdot
(x-x_0) + \frac{1}{2}\langle D^2\varphi(x_0)(x-x_0),x-x_0\rangle
+ \varepsilon |x-x_0|^2
\end{equation}
so that $x_0$ is a strict local maximum point of $w^+ -
\varphi_{\varepsilon}$. We claim that $-\Delta^+_{\infty}
\varphi_{\varepsilon}(x) < f(x)$ for all $x$ sufficiently close to
$x_0$ if $\varepsilon$ is small enough.

In fact, if $\nabla \varphi(x_0)\neq 0$, then
$\nabla \varphi(x)\neq 0$ in a neighborhood of $x_0$, and in this
neighborhood,
\begin{equation}
\Delta^+_{\infty}\varphi_{\varepsilon}(x)
= \langle D^2\varphi_{ \varepsilon}(x)\hat{\nabla
\varphi}_{\varepsilon}(x), \hat{\nabla
\varphi}_{\varepsilon}(x)> = \Delta^+_{\infty}\varphi(x_0)
+ O(\varepsilon).
\end{equation}
The claim follows from the continuity of
$\Delta^+_{\infty}\varphi$ and $f$.

If $\nabla \varphi(x_0) = 0$, then $\lambda_{\rm max}(D^2\varphi(x_0)
) = \Delta^+_{\infty}\varphi(x_0) < f(x_0)$.
As $\lambda_{\rm max}(D^2\varphi_{\varepsilon}(x))\leq
\lambda_{\rm max}(D^2\varphi(x)) + C\varepsilon$,
\begin{equation}
\Delta^+_{\infty}\varphi_{\varepsilon}(x)\leq
\lambda_{\rm max}(D^2\varphi_{\varepsilon}(x)) < f(x)
\end{equation}
holds for $x$ sufficiently close to $x_0$.

We take $\delta > 0$ small enough so that the function
$\hat{\varphi}(x) := \varphi_{\varepsilon}(x) - \delta$ satisfies
$\hat{\varphi} < w^+$ in a neighborhood of $x_0$ which is
contained in the set $\{x\in\Omega:
\Delta^+_{\infty}\varphi_{ \varepsilon}(x) < f(x)\}$, and
$\hat{\varphi} \geq w^+$ outside this neighborhood of $x_0$.

We know from the previous part of the proof that there exists a
sequence $\{v_k\}$ in $\mathcal{A}^+_{f,g}$ that converges to
$w^+$ locally uniformly in $\Omega$. Therefore there is an element
$v$ of $\mathcal{A}^+_{f,g}$ such that $\hat{\varphi} < v$ in a
neighborhood $N$ of $x_0$ which is a subset of the set
$\{x\in\Omega: \Delta^+_{\infty}\varphi_{ \varepsilon}(x)
< f(x)\}$, and $\hat{\varphi}\geq v$ outside $N$ and in some
$\Omega'\subset\subset\Omega$, if $\delta$ is taken smaller as
needed. We may without loss of generality modify the values of
$\hat{\varphi}$ near $\partial\Omega$ so that $\hat{\varphi} \geq
v$ in $\Omega\backslash N$.

Take $\hat{v} = \min\{\hat{\varphi}, v\}$. Then $\hat{v} =
\hat{\varphi}$ in the neighborhood $N$ of $x_0$ and $\hat{v} = v$
elsewhere. So $\hat{v}\in \mathcal{A}^+_{f,g}$. But $\hat{v} =
\hat{\varphi} < w^+$ in a neighborhood of $x_0$, which is a
contradiction to the definition of $w^+$. So $-\Delta^N_{
\infty}w^+(x) \leq f(x)$ in $\Omega$. Similarly,
$-\Delta^N_{\infty}w^-(x) \geq f(x)$ in $\Omega$.

Furthermore, the comparison theorem \eqref{comparison2} implies
that for any solution $w\in C(\bar{\Omega})$ of the Dirichlet
problem
\begin{gather*}
-\Delta^N_{\infty}w(x) = f(x) \quad \text{in } \Omega \\
w(x) = g(x) \quad \text{on } \partial\Omega
\end{gather*}
$w^-\leq w\leq w^+$ holds on $\bar{\Omega}$, as it holds in
$\bar{\Omega}$ that $v_2\leq w\leq v_1$ for any $v_1\in
\mathcal{A}^+$ and any $v_2\in\mathcal{A}^-$.
We have proved the following existence theorem.

\begin{theorem}\label{existence}
There exists at least one solution in $C(\bar{\Omega})$ of the boundary
value problem \eqref{dirichlet}. Every continuous solution of \eqref{dirichlet}
is locally Lipschitz continuous in $\Omega$. Among all the continuous
solutions of the boundary value problem \eqref{dirichlet}, there
are one least solution $w^-$ and one greatest solution $w^+$ as constructed
above.
\end{theorem}

Furthermore, we can acquire a clearer picture of the set of
continuous solutions of the Dirichlet problem \eqref{dirichlet} by
inspecting the solutions in the following \textit{absolute} way.
First, the construction of $w^+$ and $w^-$ and the above theorem
imply the following theorem.

\begin{theorem}\label{absolutemaxmin}
For any open set $V\subset\subset\Omega$, if $w\in C(\bar{V})$
satisfies
\begin{equation}
\begin{gathered}
-\Delta^N_{\infty}w(x) = f(x) \quad (x\in V) \\
w(x) = w^+(x) \quad (x\in\partial V),
\end{gathered}
\end{equation}
then $w\leq w^+$ in $\bar{V}$.

Similarly, if $w\in C(\bar{V})$ satisfies
\begin{equation}
\begin{gathered}
-\Delta^N_{\infty}w(x) = f(x) \quad (x\in V) \\
w(x) = w^-(x) \quad (x\in\partial V),
\end{gathered}
\end{equation}
then $w\geq w^-$ in $\bar{V}$.
\end{theorem}

\begin{proof}
According to the preceding Theorem \ref{existence}, we may assume
that $w$ is the greatest solution in the region $V$ with boundary
data $w^+$ on $\partial V$. Then $w^+\leq w$ on $\bar{V}$. We need
to prove the reverse inequality $w\leq w^+$.
Define $\tilde{w}$ on $\bar{\Omega}$ by
$$
\tilde{w}(x) = \begin{cases}
w(x), & x\in V \\
w^+(x), & x\in \bar{\Omega}\backslash V.
\end{cases}
$$
Then $-\Delta^N_{\infty}\tilde{w}(x) \leq f(x)$ in
$\Omega$ in the viscosity sense. In fact, if $\tilde{w}\prec_{x_0}
\varphi$ for a point $x_0\in\Omega$ and a $C^2$ function
$\varphi$, and if $x_0 \not\in \partial V$, then clearly
$-\Delta^N_{\infty}\varphi(x_0) \leq f(x_0)$. If
$x_0\in\partial V$, then $w^+\prec_{x_0}\varphi$ as $w^+\leq w$ in
$V$ and $w^+ = w$ on $\partial V$. As a result,
$-\Delta^N_{\infty}\varphi(x_0) \leq f(x_0)$ holds.

For any $v\in \mathcal{A}^+_{f,g}$, Theorem \ref{comparison2}
implies that $v\geq \tilde{w}$. Consequently, $w^+(x)\geq
\tilde{w}(x)$, $x\in\Omega$, and in particular $w^+(x) \geq w(x)$,
$x\in \bar{V}$.

The proof of the second part is similar.
\end{proof}

Define the set of viscosity solutions of the Dirichlet problem
\eqref{dirichlet} by
\begin{equation}
\mathcal{A}_{f,g} = \{u\in C(\bar{\Omega}):
-\Delta^N_{\infty} u(x) = f(x) \text{ in } \Omega,
\text{ and } u = g \text{ on }\partial\Omega.\}
\end{equation}
According to the preceding theorem \eqref{absolutemaxmin},
$w^+$ and $w^-$ are the extremal solutions in $\mathcal{A}_{f,h}$
in an absolute sense as mentioned above.

We conclude this section with a lemma which will be used in the
next section. The proof of the following partial continuity of the
infinity Laplacian lemma is straightforward if one observes that
if $\nabla \varphi(x_0) = 0$ for a smooth function $\varphi$, then
$\Delta^+_{\infty}\varphi(x_0) = \lambda_{\rm max}(D^2
\varphi(x_0))$.

\begin{lemma}
Suppose $\varphi$ is a $C^2$ function, and $x_k\to x_0$.
\begin{itemize}
\item[(i)] If $\nabla \varphi(x_0) \neq 0$, then $\Delta^N_{
\infty}\varphi(x_k)\to \Delta^N_{\infty}
\varphi(x_0)$.

\item[(ii)] If $\nabla \varphi(x_0) = 0$, then $\Delta^+_{\infty}
\varphi(x_0)\geq \limsup_k\Delta^+_{\infty}\varphi(x_k)$.
\end{itemize}
\end{lemma}

\noindent\textbf{Remark:}
In Lemma 3.1(ii), the inequality
holds obviously. In many cases, the inequality is indeed an
equality. However, in general, the equality is not true. For
example, in 2D, take $\varphi(x,y) = \frac{1}{2}x^2 -
\frac{1}{2}y^2$. Then $\Delta^+_{ \infty}\varphi(0,0) = 1$
but $\Delta^+_{\infty}\varphi(x,y) = \frac{ x^2 -
y^2}{x^2+y^2}$ does not necessarily converge to 1 as
$(x,y)\to (0,0)$.

\section{Unanswered Questions}

Following the proof of the existence of the maximum
and minimum solutions of the Dirichlet problem \eqref{dirichlet}
with non-trivial right-hand-side in this work, some closely related
problems need to be answered.

Naturally, one would ask when the uniqueness
of a viscosity solution of the Dirichlet problem \eqref{dirichlet}
holds even if $f$ changes sign. More precisely, what is the
necessary and sufficiency condition on $f$ (and possibly on $g$ as
well) and on the domain $\Omega$ that ensures the Dirichlet
problem \eqref{dirichlet} has a unique continuous solution? Are there always more than
one viscosity solutions of the Dirichlet problem if $f$ changes
sign? A recent work by Armstrong and Smart, \cite{AS}, answered
part of the questions. Interested reader may read their work for
up-to-date development.

One may also ask at most how many distinct solutions can the
Dirichlet problem \eqref{dirichlet} have for any non-trivial
right-hand-side? Under what condition are there infinitely many solutions?
In case there exist multiple solutions, what is the structure of the set of the
continuous solutions of the Dirichlet problem \eqref{dirichlet}?
Do the extremal solutions $w^+$ and $w^-$ determine all the solutions
of the Dirichlet problem in some
way? Or parallelly, ``What is a criterion for a
continuous function to be an element of $\mathcal{A}_{f,g}$?"

We will be more precise in our notations below and hope that the following
discussion will justify our use of multiple subscripts. Let $\mathcal{A}_{
f,g}(\Omega)$ denote the set of the viscosity solutions of the Dirichlet
problem \eqref{dirichlet} in a bounded open set $\Omega$. $w^+_{f,g,\Omega}$
and $w^-_{f,g,\Omega}$ denote the maximum and minimum solutions in
$\mathcal{A}_{f,g}(\Omega)$. The following theorem is a criterion which is
not quite up to the authors' satisfaction in that it depends on the maximum
and minimum solutions for every open subset and does not give enough
information about the solution $u$ solely in terms of  $w^+_{f,g,\Omega}$
and $w^-_{f,g,\Omega}$.

\begin{theorem}
Suppose $u\in C(\Omega)$. Then $-\Delta^N_{\infty}u(x) = f(x)$ in $\Omega$
if and only if for every open set $V\subset\subset\Omega$,
$$w^-_{f,g,V}(x) \leq u(x) \leq w^+_{f,g,V}(x),\ \ \text{for\ }\ x\in V,$$
where $g = u|_{\partial V}$.
\end{theorem}

\begin{proof}
The necessity follows from the Theorem \eqref{existence}.

To show the sufficiency, we only prove $-\Delta^N_{\infty}
u\leq f$ in $\Omega$, as the proof of $-\Delta^N_{\infty}
u\geq f$ is similar. Suppose $u\prec_{x_0}\varphi$ for some
$x_0\in\Omega$
and some $C^2$ function $\varphi$. For any small $r>0$,
let $V = B_r(x_0)$
and $w^+_r$ be the maximum solution of the Dirichlet problem in $V$. As
$w^+_r\geq u$ in $V$ and $w^+_r = u$ on $\partial V$, it is clear that
$w^+_r\prec_{x_r} \varphi$ for some point $x_r\in V$.
So $-\Delta^N_{\infty}\varphi(x_r)\leq f(x_r)$. Sending $r$ to 0,
one obtains $-\Delta^N_{\infty}\varphi(x_0)\leq f(x_0)$ on account
of the continuity of $f$ and the smoothness of $\varphi$,
noticing the fact
that $-\lambda_{\rm max}(D^2\varphi(x_0))\leq \liminf_{r\downarrow 0}-
\Delta^N_{\infty}\varphi(x_r)$ if $\nabla \varphi(x_0) = 0$.
\end{proof}

Clearly, $u\in C(\bar{\Omega})$ is an element of
$\mathcal{A}_{f,g}(\Omega)$
if and only if $u$ verifies the condition stated in the preceding theorem
and $u = g$ on $\partial\Omega$. On the other hand, it is unknown if the
comparison property $\sup_V(u - w^+_{f,g,\Omega})\leq\max_{\partial V}(u - w^+_{
f,g,\Omega})$ and $\inf_V(u - w^-_{f,g,\Omega})\geq\min_{\partial V}(u - w^-_{f,
g,\Omega})$ for every open subset $V\subset\Omega$ alone implies $u\in\mathcal{A}_{
f,g}(\Omega)$.

In addition, can we anticipate a differential game theory
interpretation of the Dirichlet problem \eqref{dirichlet} with the
nontrivial right-hand-side $f$ as we do with the case
$\sup_{\Omega}f(x) < 0$ (\cite{PSSW}, \cite{BEJ} and
\cite{E1})? This question has been partially answered by Armstrong and
Smart in \cite{AS}. Furthermore, one may still ask the questions such
as ``Are there any connections between the maximal and minimal solutions
and the value functions of the players II and I in the generalized
`tug-of-war' game?"

In the end, one may also consider the inverse problem of
the Dirichlet problem \eqref{dirichlet}, ``For what continuous functions $u$, are there
continuous functions $f$ such that $-\Delta^N_{\infty}u =
f$?" The uniqueness of $f$ was initially considered in \cite{PSSW}
and has recently been proved by Y.\,Yu (\cite{Y}).

\subsection*{Acknowledgments}
We would like to thank the referees for very valuable suggestions
and corrections, especially for their comments on open questions
and the suggestion to revise part of the proofs.
After we submitted this paper,
we received a preprint from  S.N. Armstrong and C.K. Smart \cite{AS}
in which they independently established similar results to ours
among other things through a finite difference method.

\begin{thebibliography}{00}

\bibitem{AS} Armstrong, S.\,N. and Smart, C.\,K.,``A finite difference
approach to the infinity Laplace equation and tug-of-war games", preprint 2009.

\bibitem{A1} Aronsson,G., ``Minimization problems for
the functional $\sup_x F(x, f(x), f'(x))$ ", { \em Ark. Mat.}, 6,
1965, 33-53 (1965).

\bibitem{A2} Aronsson, G., ``Minimization problems for
the functional $\sup_x F(x, f(x), f'(x))$. II." {\em Ark. Mat.},
6, 1966, 409-431 (1966).

\bibitem{A3} Aronsson, G., ``Extension of functions
satisfying Lipschitz conditions", { \em Ark. Mat.}, 6, 1967,
551-561 (1967).

\bibitem{ACJ} Aronsson, G., Crandall, M. and Juutinen, P.,
``A tour of the theory of absolute minimizing functions", { \em
Bull. A.M.S.}, 41(2004), no.4, 439-505.


\bibitem{BB} Barles, G.\, and Busca, J.\,, ``Existence and
comparison results for fully nonlinear degenerate elliptic
equations without zeroth-order term", {\em Comm. Partial Diff.
Equations}, 26(2001), 2323-2337.

\bibitem{BDM} Bhattacharya, T., DiBenedetto, E.\,and Manfredi,
J.``Limits as $p\to\infty$ of $\Delta_pu_p = f$
and related extremal problems, Some topics in nonlinear PDEs
(Turin, 1989)", {\em Rend. Sem. Mat. Univ. Politec. Torino 1989},
Special Issue, 15-68(1991).

\bibitem{BEJ} Barron, E.\,N., Evans, L.\,C.\, and Jensen,
R.\,R.\,, ``The infinity Laplacian, Aronsson's equation and their
generalizations", {\em Trans. Amer. Math. Soc.}, 360(2008), no.1,
77-101(electronic).

\bibitem{BJW} Barron, E.\,N.\,,Jensen, R.\,R.\, and Wang,
C.\,Y.\,, ``The Euler equation and absolute minimizers of
$L^{\infty}$ funcitonals", Arch. Ration. Mech. Anal. 157(2001),
no.4, 255-283.

\bibitem{C} Crandall, M.\,G.\,, ``An efficient derivation of the
Aronsson equation", Arch. Rational Mech. Anal., 167(2003), no.4,
271-279.

\bibitem{CE} Crandall,M.\,G.\, and Evans,L.\,C.\,, ``A remark on
infinity harmonic functions", Proceeding of the USA-Chile Workshop
on Nonlinear Analysis (Vi\~{a} del Mar-Valparaiso, 2000), 123-129
(electronic), Electon. J. Differ. Equ. Conf., 6, Southwest Texas
State Univ., San Marcos, TX, 2001.

\bibitem{CEG}Crandall, M.\,G.\,, Evans, L.\,C.\, and Gariepy,
R.\,F.\,, ``Optimal Lipschitz extensions and the infinity
Laplacian", Calc. Var. Partial Differential Equations 13(2001),
no. 2, 123-139.

\bibitem{CGG} Chen, Y., Giga, Y. and Goto, S., ``Uniqueness
and existence of viscosity solutions of generalized mean curvature
flow equations", {\em J. Differential Geom.}, 33(1991), 749-786.

\bibitem{CGW} Crandall, M.\,G.\,, Gunnarsson, G.\, and Wang,
P.\,, ``Uniqueness of $\infty$-harmonic functions and the eikonal
equation", {\em Comm. Partial Diff. Equations.}, 32(2007),
1587-1615.

\bibitem{CIL} Crandall, M.\,G., Ishii, H. and Lions, P.\,L.,
``User's guide to viscosity solutions of second-order partial
differential equations", { \em Bull. A.M.S.}, 27(1992), 1-67.

\bibitem{CW} Crandall, M.\,G. and Wang, P. ``Another Way To Say
Caloric", {\em J. Evol. Equ.} 3(2004), 653-672.

\bibitem{E1} Evans, L.\,C.\,, ``The 1-Laplacian, the infinity
Laplacian and differential games", {\em Perspectives in nonlinear
partial differential equations, Comtemp. Math. 446},
Amer.\,Math.\,Soc.\,, Providence, RI, (2007), 245-254.

\bibitem{E2} Evans, L.\,C.\,, ``Estimates for smooth absolutely
minimizing Lipschitz extensions," {\em Electron. J. Differential
Equations} (1993), no.03, approx.9pp.

\bibitem{J} Jensen, R.\, ``Uniqueness of Lipschitz
extensions: minimizing the sup norm of the gradient", {\em Arch.
Rational Mech. Anal.} 123(1993), no.1, 51-74.

\bibitem{Ju1} Juutinen, P., ``Minimization problems for Lipschitz
functions via viscosity solutions", Ann. Acad. Sci. Fenn. Math.
Diss. No.115 (1998), 53 pp.

\bibitem{Ju2} Juutinen, P., ``Principal eigenvalue of a very
badly degenerate operator and applications," {\em J.Differential
Equations} 236(2007), no.2, 532-550.

\bibitem{JK} Juutinen, P. and Kawohl, B., ``On the evolution
governed by the infinity Laplacian," {\em Math. Ann.} 335(2006),
819-851.

\bibitem{JLM} Juutinen, P.\,,Lindqvist, P.\,,Manfredi, J.,
``The $\infty$-eigenvalue problem," {\em Arch. Ration. Mech. Anal.}
148(1999), no.2, 89-105.

\bibitem{LM1} Lindqvist, P., and Manfredi, J., ``The Harnack
inequality for $\infty$-harmonic functions", {\em Electron. J.
Differential Equations} 1995, No.04, approx. 5 pp.

\bibitem{LM2} Lindqvist, P., and Manfredi, J., ``Note on
$\infty$-superharmonic functions", {\em Rev. Mat. Univ. Complut.
Madrid} 10(1997), no.2, 471-480.

\bibitem{LW1} Lu, G.\, and Wang, P.,``Inhomogeneous infinity
Laplace equation," {\em Adv. Math.} 217(2008), no.4, 1838-1868.

\bibitem{LW2} Lu, G. and Wang, P.,``A PDE perspective of the
normalized infinity Laplacian", {\em Comm. P.D.E.} 10(2008),
1788-1817.

\bibitem{PSSW}Peres, Y., Schramm, O., Sheffield, S., and
Wilson, D.\,, ``Tug of war and the infinity Laplacian",
{\em J. Amer. Math. Soc.} 22(2009), no.1, 167-210.

\bibitem{Y} Yu, Y., ``Uniqueness of values of Aronsson operators
and running costs in `tug-of-war' games", preprint.

\end{thebibliography}

\end{document}
