pub fn solve(
eph: &dyn EphemerisSource,
inputs: &SolveInputs,
with_geodetic: bool,
) -> Result<ReceiverSolution, SppError>Expand description
Run the SPP solve from synthesized/measured pseudoranges.
Uses the core trust-region weighted least-squares solver over the
sqrt(w) * (P_meas - P_hat) residual. The converged position/clock is a
sub-micron solver-agreement result (the linear-algebra step is not
bit-reproducible across BLAS builds), not a 0-ULP claim. The residual /
Jacobian substrate evaluated at recorded states is the 0-ULP target and is
exercised by the trace-replay parity test, not by this entry point.
This is the reference SPP entry point: it runs the legacy
SolverRecipe::NalgebraTrfLegacy trust-region factorization, so its
existing goldens are unchanged. solve_with_solver selects the owned
deterministic kernel.