v0.9.0

v0.9.0 — wide-parameter fitting and the complete output surface

New

• Wide-parameter orbit determination, exposed end to end. Beyond the
  six-element state, a fit can solve for the Marsden non-gravitational
  coefficients, the cometary outgassing time delay DT, and the
  solar-radiation-pressure area-to-mass ratio (AMRAT) — individually or
  together, each partial supplied analytically by the hyperdual
  integrator. Every solved parameter carries its own posterior
  uncertainty in a tagged covariance that names its slot.
• Solar radiation pressure is a first-class force and fit parameter.
  Orbits carry an SRP block that round-trips through propagation,
  ephemerides, impact analysis, and OD without being dropped.
• Post-fit photometry: H and a phase-function slope read back from the
  arc's magnitudes with honest 1σ, climbing the H-only → HG₁₂ → HG₁G₂
  ladder the phase coverage supports. The band→V table covers the
  modern two-character ADES codes; magnitudes in bands with no adopted
  conversion are excluded, counted, and their codes listed — never
  silently used, astrometry unaffected.
• The complete output surface: every value the engine computes now
  crosses the C ABI (version 2; append-only growth, recompile against
  the new header).
  • Ephemeris rows carry the 6×6 sky-plane covariance and the aberrated
    barycentric ICRF state at the photon-emission epoch with its own
    covariance; generate calls return their non-fatal warnings.
  • Residual rows carry radar delay/Doppler residuals, the
    D-optimality information loss on removal, and the full
    along/cross-track 2×2; OD results carry an event-aware
    covariance-trust verdict — absence of a verdict is not trust.
  • Impact records carry the geodetic impact point, the Monte-Carlo
    95% binomial confidence half-width, second-order mean/σ/skewness
    of the miss distance, the distance gradient and Hessian, and the
    mixture component count.
  • Plan evaluation reports radar SNR, range, mode, and link-budget
    provenance plus the predicted optical ephemeris; mixture
    components are basis-tagged with frame and origin (C ABI).
• transform_coordinates carries the full state-plus-non-grav-plus-SRP
  covariance through frame and coordinate changes.

Fixes

• query_sbdb populates AMRAT from JPL's Small-Body Database, so
  SRP-sensitive objects seed with their catalog area-to-mass ratio
  instead of silently missing it.
• A Marsden fit seeded from zero coefficients opens its solve column
  from the covariance prior instead of being read as "no force model".
• Impact probability is evaluated in periapsis space against the
  body's physical radius, so gravitational focusing is not
  double-counted.

Binds empyrean-core v0.9.2 (villeneuve v1.20.2, scott v1.15.0).