mirror of
https://github.com/cosinekitty/astronomy.git
synced 2025-12-31 03:30:26 -05:00
The reason SSB vector errors were larger than other bodies is because the Sun/Barycenter relationship does not have position and velocity vectors of a "typical" size. The distance of a planet from the SSB is fairly constant, and the speed a planet travels is fairly constant. Therefore, comparing errors by dividing by the magnitude of the correct vector usually makes sense for scaling. But for the barycentric Sun (or the heliocentric SSB), the magnitude of the vectors can become arbitrarily small, nearly zero in fact, resulting in surprisingly large ratios. I compensated for this in all the tests by adding a new rule. When the error thresholds r_thresh and v_thresh are negative, it is a flag that indicates they are absolute, not relative. In other words, r_thresh < 0 indicates that abs(r_thresh) is the maximum number of astronomical units allowed in position errors, and v_thresh < 0 specifies maximum AU/day. This results in more consistent numbers that give confidence the errors are indeed very small and not worth worrying about.