In the article: https://arxiv.org/abs/0906.3217 the authors prove in Lemma 1 a formula which helps compute more easily the integral of the Hessian of a function defined on $\Bbb{S}^2$. More precisely, if $h : \Bbb{S}^2 \to \Bbb{R}$ is a $C^2$ function, $Hess(h)(X,Y) = \langle \nabla_X \nabla h,Y\rangle$ is the Hessian of $H$ and $H(h)$ is the determinant of the Hessian, then $$ \int_{\Bbb{S}^2} H(h) dA = \frac{1}{2} \int_{\Bbb{S}^2} |\nabla h|^2 dA$$ where $dA$ is the area element on $\Bbb{S}^2$. This formula is helpful and simplifies some aspects of my numerical computations. Here and in the following, $\nabla$ and $\Delta$ represent the tangential gradient and Laplace-Beltrami operator on $\Bbb{S}^2$.
The above formula is proved in connection with bodies of constant width, but in the proof in the article they don't seem to use this fact. In my numerical experiments the formula gives the expected result in the general case. However, I work in the case where $H(h)>0$ (I don't know if this is relevant or not...)
I was wondering if it is possible to obtain a similar simplification for the integral $\int_{\Bbb{S}^2} h H(h) dA$? More precisely, is it possible to obtain something of the form $$ \int_{\Bbb{S}^2} h H(h) dA = \int_{\Bbb{S}^2} \mathcal{F}(h,\nabla h,\Delta h)dA $$ where $\mathcal{F}$ is "simple" (polynomial)?