Gradient descent now works for the constant amplitude model.  I also wrote print functions for the isotropics and the constants, and re-arranged the array/number conventions in the bessel functions.  This is bothering me, but I got gradient descent to work just by setting the step size to 1/1000.  I think I could do better with an adaptive line search, and I really really want to avoid hard coding anything.