Close encounters with major Solar System bodies may bring about a strong amplification of numerical error during inter-planetary orbit propagation. In this work, we reduce global numerical error by integrating regularized equations of motion instead of the classical Newtonian equations in Cartesian coordinates. The integration performance of several sets ofregularized equations is assessed from large-scale numeri-cal propagations of close encounters in the Sun-Earth planar CR3BP. An essential device consists in switching between primary bodies during the propagation, which effectively decomposes a strongly-perturbed heliocentric problem into two weakly-perturbed ones; this propagation approach has been dubbed Online Trajectory Matching (OTM). Through this simple expedient, regularized equations describing the evolution of non-classical orbital elements achieve excellent performances compared to Newtonian equations, even when employing sophisticated adaptive numerical schemes.Further improvements might be expected by carefully selecting the location of the switch of primary bodies during the propagation.