I am given real values $p, s, t, u$ and wish to find unknown values $r, v$. As shown in the diagram below, $p$ and $s$ are radii of two given circles, with centers at $(0,-p)$ and $(0,t)$.
At present I'm using Newton's method to solve $a(v)=b(v)$ for $v$, with $a(v) = p - \sqrt{u^2 + v^2}$ [from $(p-a)^2=u^2 + v^2$] and $b(v) = s - \sqrt{u^2 + (t-v)^2}$ [from $(s-a)^2=u^2 + (t-v)^2$], with $a$ and $b$ representing the after-intersecting portions of rays from $(0,-p)$ and $(0,t)$. When $a(v)=b(v)$, we have $a(v)=b(v)=r$. Given $r$ or $v$ it's trivial to find the other and to find tangency point G, etc.
Is there a better method, preferably some closed form, for finding $r$ or $v$ when given $p, s, t, u$?


We can compute the distance from the centers of the two large circles to the center of the small circle in two ways: by the distance formula, and by taking the difference of their radii. This gives us the equations
\begin{align} (s - r)^2 &= u^2 + (v-t)^2 \\ (p-r)^2 &= u^2 + (v+p)^2 \end{align}
These are quadratic in $r$ and $v$, but taking the difference of these two equations and simplifying gives us a linear equation in $r$ and $v$: $$ s^2 + (2p - 2s) r = t^2 - (2p + 2t)v. $$ We can solve for either unknown variable in terms of the other; for example, we can take $$ r = \frac{t^2 - s^2 - (2p+2t)v}{2p-2s}. $$ Now substitute this expression for $r$ into either of the original equations, and you get a quadratic equation for $v$. This can be solved using the quadratic formula; the expression you get is messy, but still gives a closed form for $v$ (and therefore for $r$).