Solution (source code)

= Solution

An affine geodesic follows from
$$
\mathcal L=\frac12\left[r^{-2}(-\dot t^2+\dot r^2)+r^2\dot\theta^2\right].
$$
Its constant norm may be chosen as
$$
2\mathcal L=\varepsilon,
\qquad
\varepsilon=-1,0,+1
$$
for timelike, null and spacelike geodesics respectively.