Gaussian inference on smooth manifolds is central to robotics, but exact marginalization and conditioning are generally non-Gaussian and geometry-dependent. We study tangent-linearized Gaussian inference and derive explicit non-asymptotic W2 stability bounds for projection marginalization and surface-measure conditioning. The bounds separate local second-order geometric distortion from nonlocal tail leakage and, for Gaussian inputs, yield closed-form diagnostics from (μ,Σ) and curvature/reach surrogates. Circle and planar-pushing experiments validate the predicted calibration transition near ∥Σ∥op/R≈1/6 and indicate that normal-direction uncertainty is the dominant failure mode when locality breaks. These diagnostics provide practical triggers for switching from single-chart linearization to multi-chart or sample-based manifold inference. Code and Jupyter notebooks are available at https://github.com/mikigom/StabilityTLGaussian.
Figures & tables
Fig. 1 : W2 stability of marginalizing a Gaussian onto a circular arc. Exact metric projection g vs. the tangent–retraction approximation g^ .
Fig. 2 : W2 stability of conditioning a Gaussian onto a circular arc. Surface-measure conditioning vs. a tangent-plane chart approximation.
Fig. 3 : Circle sweep ( δ=0.2 ) over σ/R : (a) variance ratio ϱ=Varlin/Varexact , (b) realized 95% coverage, and (c) normalized diagnostics on ∥Σ∥op/R : tightened theoretical W2 bound (data-calibrated Ctail with strict upper-envelope safeguard) and empirical W2 proxy. Dashed line marks σ/R=1/6 . Colors indicate R∈{0.5,1,2} .
Fig. 4 : Circle generality stress tests: (a) anisotropic sweep ( η=σn/σt ) versus ∥Σ∥op/R , (b) tail proxy ε versus offset δ/R , and (c) coverage degradation versus δ/R .
tr(Σxy)
λmax(Σxy)
Red.
∥ΔΣxy∥F
k
Unc.
Con.
Unc.
Con.
(%)
1
0.0018
0.0009
0.0009
0.0009
50.5
0.0009
13
0.0234
0.0105
0.0117
0.0105
55.3
0.0118
25
0.0450
0.0190
0.0225
0.0190
57.8
0.0228
37
0.0666
0.0278
0.0333
0.0278
58.2
0.0337
49
0.0882
0.0377
0.0441
0.0377
57.3
0.0446
TABLE I : Quantitative covariance metrics for planar pushing at timesteps
Fig. 5 : Trajectory-wide planar-pushing diagnostics: (a) curvature/reach proxies (κ^k,ρ^k) , (b) spread sk and locality indicator sk/ρ^k , (c) Monte Carlo mismatch (ϱkMC,Δk) .
Fig. 6 : Directional scaling at the terminal planar-pushing pose: isotropic, normal-only, and tangential-only covariance inflation.
Division of Mathematical Sciences, Nanyang Technological University. · School of Statistics and Data Science, Shanghai University of Finance and Economics.
School of Physical and Mathematical Sciences, Nanyang Technological University · School of Statistics and Data Science, Shanghai University of Finance and Economics.