C++ radial velocity routines