Wrong computation of relativistic clock correction for one way GNSS measurements
Current implementation is:
dt_rel = -2.0 * (p_station.v_station - p_sat.v_sat) / c²
However, the formula is wrong since there is an inversion between satellite and station coordinates. Furthermore, for one-way range measurements, there is no need to take into account for the station position. Therefore, the formula can just be:
dt_rel = -2.0 * (p_sat.v_sat) / c²