Sure, I can do that.  I suppose I should also save ``np.cos(theta)`` and ``np.sin(theta)`` since they are calculated many times.