An algorithm for computing geometric relative velocities through Fermi and observational coordinates

We present a numerical method for computing the \textit{Fermi} and \textit{observational coordinates} of a distant test particle with respect to an observer. We apply this method for computing some previously introduced concepts of relative velocity: \textit{kinematic}, \textit{Fermi}, \textit{spectroscopic} and \textit{astrometric} relative velocities. We also extend these concepts to non-convex normal neighborhoods and we make some convergence tests, studying some fundamental examples in Schwarzschild and Kerr spacetimes. Finally, we show an alternative method for computing the Fermi and astrometric relative velocities.