Add initial orbit determination

[moeyensj]

An implementation of Gauss' method with Herrick+Gibbs or Gibbs to determine the velocity needs to be added so we can validate clusters on the fly.

This also requires quality of fit calculations using ephemeris generation (residuals, chi2, etc..).