Let $\Omega\subset\mathbb{R}^2$ is bounded and convex and $\partial\Omega$ be smooth. Consider the second order elliptic PDE (1) $$ \begin{cases} Lu = f &\text{ on } \Omega\,,\\ u=0 &\text{ on } \partial\Omega\,, \end{cases} $$ where $f\in C^\infty(\Omega)$ satisfies $f(x) > 0$ and $f$ has no isolated critical points on $\Omega$. Then, does the solution $u$ of the above have only one critical point? I am looking for a result similar to Theorem 2 of this paper(ScienceDirect):

Theorem 2: Let $\Omega$ be a bounded, strictly convex domain in $\mathbb{R}^2$ and $u\in C^3(\Omega)\cap C^1(\bar{\Omega})$ a solution to the boundary value problem \begin{align*} \Delta u &= f(u,\nabla u)\quad\text{ in }\Omega\,,\\ u=\text{const}\,,& \quad\nabla u\neq0\,,\quad\text{ on }\partial\Omega\,, \end{align*} where $f\in C^1$, $f_u\geq0$. Then $u$ has exactly one critical point in $\bar{\Omega}$ and there $\det(D^2u)>0$ holds (i.e., a global maximum or minimum).

The paper claims that Theorem 2 is true even if we replace Laplace's operator by another elliptic operator. This means, for example, that if $f$ is constant, then the solution to (1) has a unique critical point.