Does this P.D.E: $$\nabla\cdot\left( \frac{ \nabla u}{u} \right)+a\, \Delta u+b\,u=0 \hspace{3cm} (*)$$ have a variational structure? Here $a$ and $b$ are constants.
In other words, the question I am asking is:
Does there exist a functional such that the corresponding Euler–Lagrange equation is (*)?
I)
Sketched proof: The functional derivative is
$$ \frac{\delta S[u]}{\delta u} ~=~\frac{\partial {\cal L}(u,\nabla u)}{\partial u} -\nabla \cdot \frac{\partial {\cal L}(u,\nabla u)}{\partial (\nabla u)}$$ $$\tag{3}~=~-\left(a+\frac{1}{u}\right) \left[-\frac{1}{u^2}(\nabla u)^2 +\left(a+\frac{1}{u}\right)\nabla^2u +bu\right],$$
which leads to the Euler-Lagrange equation
$$\tag{4}\nabla\cdot\left(\frac{1}{u}\nabla u\right) +a\nabla^2u +bu~=~0\quad \vee\quad u~=~-\frac{1}{a}.$$
II)
Let $F$ and ${\cal V}$ be antiderivatives (aka. primitive or indefinite integrals) of $f$ and $ge^{2F}$, respectively, i.e.
$$\tag{6}F^{\prime}(u)~=~f(u) \quad\text{and}\quad {\cal V}^{\prime}(u)~=~g(u)e^{2F(u)}.$$
We now multiply the normalized equation (5) with an integrating factor
$$\tag{7}\lambda(u)~:=~ e^{2F(u)}. $$
Sketched proof:
$$ \frac{\delta S[u]}{\delta u} ~=~\frac{\partial {\cal L}(u,\nabla u)}{\partial u} -\nabla \cdot \frac{\partial {\cal L}(u,\nabla u)}{\partial (\nabla u)}$$ $$\tag{9}~=~-e^{2F(u)}\left[ \nabla^2u+ f(u)(\nabla u)^2+g(u) \right].$$