I need to maximise a quadratic form in 3D with the constraint that the solution vector, $\mathbf{x}$, lies on the unit sphere. I know that for a simple quadratic form:
$$ \mathbf{x}^T\mathbf{A}\mathbf{x} $$
this would simply be the eigenvector with the highest eigenvalue, and it can be proved simply enough. However, my situation is a more general form
$$ \mathbf{x}^T\mathbf{A}\mathbf{x} + \mathbf{b}\mathbf{x} + c $$
(not even sure if I can call it quadratic any more at this point since there's linear terms mixed into it...), basically corresponding to a translation from the origin. In these conditions, I can't see how to find the constrained maximum in an elegant way. Of course I could just optimise it numerically but it would be nice if there was an analytical method, even just to find the stationary points (and then check them one by one). Thanks!
This is the trust-region subproblem, see e.g. https://www8.cs.umu.se/kurser/5DA001/HT07/lectures/trust-handouts.pdf