Thiele-Innes-based solver for speed boost

[sblunt]

[semaphoreP]

Essentially a different basis that makes the orbits linear in 4 parameters (analytical marginalization). Would need new basis class + new chi2 function for analytical marginalization of orbital elements.