Sampling on the sphere from $f(x) \propto x^TAx$

A method for drawing random samples of unit vectors $x$ in $R^p$ with density proportional to $x^TAx$ where $A$ is a symmetric, positive definite matrix. Includes an R function which implements the method.