Let $\Omega\subset\mathbb{R}^n$ be a bounded strictly convex domain and $\nu:\partial\Omega\rightarrow\mathbb{R}$ be a lower semi-continuous functions. It is well known that the function $\overline{\nu}:\overline{\Omega}\rightarrow\mathbb{R}$ defined by
$$
\overline{\nu}(x):=\sup\Big\{\alpha(x)\,\Big|\,\alpha:\mathbb{R}^n\rightarrow\mathbb{R} \mbox{ is an affine function and }\alpha|\_{\partial\Omega}\leq \nu\Big\}
$$
is not $C^1$ in general but is the generalised convex solution to the Dirichlet problem of Monge-Ampère equation
$$
\det D^2u=0,\quad u|\_{\partial\Omega}=\nu.
$$
For a problem I'm working on, it would be nice if there is a smooth alternative of $\overline{\nu}$: 


**Question.** Is it true that for any $\Omega$, $\nu$ and any $\epsilon>0$ there is a strictly convex smooth function $u$ on $\Omega$ satisfying
$$
\det D^2u<\epsilon,\quad u|\_{\partial\Omega}=\nu\, ?
$$


Any comments or hints for reference are welcome.


**Remark on boundary value.** Here, for a convex function $u$ on $\Omega$, the boundary value $u|\_{\partial\Omega}$ is defined in such a way that $u|\_{\partial\Omega}(x\_0)$ ($x\_0\in\partial\Omega$) is the limit of $u(x)$ when $x$ tends to $x\_0$ along a line segment contained in $\Omega$. This does not depend on the choice of the segment. In general, $u|\_{\partial\Omega}$ is only lower semi-continuous.



Finding a smooth convex function with prescribed boundary value and small Monge-Ampère measure
In the case $n \geq 3$ it is false. One way to see this is to use the Pogorelov example
$$w(x', \,x\_n) = |x'|^{2-2/n}f(x\_n),$$ 
which for the appropriate choice of $f$ solves $\det D^2w = 1$ in $\{|x\_n| < \rho\}$. This example is $C^{1,1-2/n}$ and is $0$ on the $x\_n$-axis. (This example arises from the affine invariance of the Monge-Ampere equation; $w$ is invariant under the scaling $w(x',\,x\_n) = \frac{1}{\lambda^{2-2/n}}w(\lambda x', x\_n)$, which preserves the Monge-Ampere measure).


Take $\nu = w$ on $\partial B\_{\rho}$. By the comparison principle, $u \geq w$ in $B\_{\rho}$, and by convexity $u = 0$ on the $x\_n$ axis. Thus, $u$ is not strictly convex, and its regularity is at best $C^{1,1-2/n}$. (One can use analogues of the Pogorelov example of the form $|x'| + |x'|^{n/2}g(x\_n)$ in the same way to generate Lipschitz singularities for $u$ along a line segment when the boundary data are Lipschitz).


This example is in a sense optimal; [Caffarelli](https://www.jstor.org/stable/1971509?seq=1#page_scan_tab_contents) showed that if $\nu$ is $C^{1,\beta}$ for $\beta > 1-2/n$ then solutions to $\det D^2u = \epsilon, \quad u|\_{\partial \Omega} = \nu$ are strictly convex and smooth in $\Omega$.


In the case $n = 2$, solutions to $\det D^2u = \epsilon$ are locally strictly convex by a classical result of Alexandrov. By taking smooth approximations $\nu\_k$ of the boundary data and solving $\det D^2u\_k = \epsilon, \quad u\_k|\partial \Omega = \nu\_k$ one can hope to obtain the desired approximation in the limit.


