Finding the maximum of a function on $ \Bbb{S}^{7} $. I'm trying to find the maximum of the function
$$2 a^2 h+\sqrt{3} a d f+\sqrt{3} a e g+2 b^2 h-\sqrt{3} b d g+\sqrt{3} b e f\\+2 c^2
   h+\sqrt{3} c d^2+\sqrt{3} c e^2-\sqrt{3} c f^2-\sqrt{3} c g^2\\-d^2 h-e^2 h-f^2
   h-g^2 h-2 h^3$$
on the sphere $$a^2+b^2+c^2+d^2+e^2+f^2+g^2+h^2=1 \text{.}$$
I suspect it is $2$, but I can't come up with a proof. Lagrange multipliers don't work, since the resulting equations don't have an analytic solution.
The value $2$ is attained for $h=-1$.
My current idea is that terms of the form $x^2y$ can't be bigger than $\frac{2}{3\sqrt{3}}$ and terms of the form $xyz$ can't be bigger than $\frac{1}{3\sqrt{3}}$, but I'm not sure this helps.
I would also be satisfied with a slightly higher bound than $2$.
Additional information:
The function is proportional to $\sum_{ijk=1}^8 d_{ijk}p^i p^j p^k$ (the constraint is just $\sum_i (p^i)^2 = 1$), where $ \tau_i\tau_j = -\frac{2}{3}\delta_{ij} +(f_{ijk}-i\, d_{ijk})\tau_k$ holds for the generators $\tau_i$ of $SU(3)$ with $[\tau_i,\tau_j]=f_{ijk} \tau_k$ and the normalization condition $\mathrm{Tr}(\tau_i\tau_j) = -2\delta_{ij}$ is imposed.

Let $$a=p\cos\chi,\, b=p\sin\chi\\
d=m\cos\phi,\,e=m\sin\phi\\
f=k\cos\theta,\,g=k\sin\theta\\
q^2=c^2+q^2$$
Then replace $h^3$ by $h(1-p^2-c^2-m^2-k^2)$
$$\text{Maximise }h\left[4p^2+4c^2+m^2+k^2-2\right]+\sqrt{3}pmk\cos(\chi-\phi+\theta)+\sqrt{3}c(m^2-k^2)$$
Let the angle be zero, $m=n\cos\alpha,k=n\sin\alpha$
$$\text{Maximise } h[4p^2+4c^2+n^2-2]+\sqrt{\frac34}\left[pn^2\sin2\alpha+2cn^2\cos2\alpha\right]$$
Let $p=q\cos\beta,c=q\sin\beta$, 
$$\text{Maximise }h[4q^2+n^2-2]+\sqrt{\frac34}qn^2[\sin2\alpha\cos\beta+2\cos2\alpha\sin\beta]$$
Let $\alpha=0,\beta=\pi/2$
$$\text{Maximize }h[4q^2+n^2-2]+\sqrt{3}qn^2\\
q^2+n^2+h^2=1$$
Use Lagrange Multipliers, I found that either $n=0$, in which case we maximize $h(2-4h^2)$, which is maximized at $h=-1$; or else $h=0$, in which case we maximize $\sqrt{3}(1-q^2)q$, whose maximum is $2/3$.
