---
title: Min Distance to Riemannian Mean (MDRM)
url: https://www.emergentmind.com/topics/minimum-distance-to-riemannian-mean-mdrm
type: topic
---

# Min Distance to Riemannian Mean (MDRM)

The minimum distance to Riemannian mean (MDRM) quantifies the minimal average squared geodesic distance from empirical data points on a Riemannian manifold to their Fréchet mean. This notion is central to statistical inference, clustering, and geometric data analysis on manifolds. MDRM is rigorously defined via the empirical Fréchet energy and can be computed using both stochastic and deterministic gradient algorithms, with theoretical guarantees under curvature constraints [1111.3120]. Efficient computation of MDRM on complex manifolds is enabled by neural approximations of the metric-constrained Eikonal equation, providing differentiable surrogates of the Riemannian distance [2404.08754].

## 1. Formal Definition and Variational Characterization

Let $(M,g)$ be a complete Riemannian manifold and $\{y_i\}_{i=1}^{n}\subset M$ a finite sample. The empirical measure is $\mu_n = \frac{1}{n}\sum_{i=1}^n \delta_{y_i}$. Define the empirical Fréchet energy:
$$
F_n(x) = \frac{1}{n} \sum_{i=1}^n d^2(x, y_i),
$$
where $d(\cdot,\cdot)$ denotes the geodesic distance. The sample Fréchet mean $\hat\mu_n$ is
$$
\hat\mu_n = \arg\min_{x\in M} F_n(x),
$$
and the minimum distance to mean (MDRM) statistic is
$$
F_n(\hat\mu_n) = \frac{1}{n} \sum_{i=1}^n d^2(\hat\mu_n, y_i).
$$
This functional measures point-cloud spread relative to their intrinsic mean under the manifold geometry [1111.3120].

## 2. Optimality Conditions and Uniqueness

For sufficiently regular distributions and injectivity radius, the gradient of $F$ at $x$ is given by
$$
\nabla F(x) = -2 \int_M \exp_x^{-1}(y)\,dP(y) \in T_xM.
$$
Hence, $x=\mu$ is a Fréchet mean if and only if
$$
\int_M \exp_\mu^{-1}(y)\,dP(y) = 0,
$$
ensuring $\mu$ zeroes the intrinsic mean of log-mapped vectors. On compact Riemannian manifolds and for empirical measures, the Fréchet mean is unique for generic data configurations [1111.3120].

## 3. Curvature Bounds and Convexity Properties

Let $\bar B(x,\rho)$ be a geodesic ball and sectional curvatures $K_\sigma\le\Delta$ on $\bar B(x,\rho)$. When $\rho < \pi/(4\sqrt{\Delta})$ (for $\Delta>0$),
$$
\operatorname{Hess} F(x)[v,v] \ge 2f_{\Delta}(\rho)\|v\|^2,
$$
with $f_\Delta(s) = (\sqrt{\Delta}s)\cot(\sqrt{\Delta}s)$ for $\Delta>0$, $f_\Delta(s)=1$ for $\Delta\le 0$. Thus, $F$ is strongly convex along geodesics, which guarantees uniqueness and stability of the minimizer, thereby affirming the robustness of the MDRM concept under curvature constraints [1111.3120].

## 4. Algorithmic Computation of MDRM

For empirical measures, MDRM can be computed via gradient-based methods:

- **Stochastic Gradient Update:**  
  $X_{k+1} = \exp_{X_k}[ -2 t_{k+1} v_k ]$,  
  with $v_k = \exp^{-1}_{X_k}(P_{k+1})$, $P_{k+1}\sim\mu_n$. Step sizes $t_k$ satisfy $\sum t_k=\infty$, $\sum t_k^2<\infty$ [1111.3120].

- **Deterministic Gradient Descent:**  
  $x_{k+1} = \exp_{x_k}[ -t_k G_k ]$,  
  where $G_k = -2\frac{1}{n}\sum_{i=1}^n \exp_{x_k}^{-1}(y_i)$. Step sizes $t_k$ as above ensure convergence to $\hat\mu_n$ [1111.3120].

- **Practical Recommendations:**  
  Initialize at one of $\{y_i\}$, iterate using chosen scheme until $d(x_{k+1},x_k)<\varepsilon$, then compute MDRM as $F_n(\hat\mu_n)$. Error control is achieved by step-sizes of order $1/k$, yielding $O(1/k)$ convergence in mean-square.

## 5. Neural Eikonal Solvers and Efficient MDRM Evaluation

Deep metric-constrained Eikonal solvers enable efficient approximation of the Riemannian distance function $\phi_\theta$ on manifold $M$ [2404.08754]. The MDRM for a query $x$ is evaluated as
$$
m(x) := d_g(x, \mu^*) \approx \phi_\theta(x, \mu^*),
$$
with computation complexity $O(W)$, where $W$ is the network FLOPs. For geodesic path recovery, backtracking using explicit ODE integration (Euler, RK4, symplectic) from $x$ to $\mu^*$ follows the negative Riemannian gradient field of $\phi_\theta$.

Accuracy guarantees are enforced by training loss $\mathcal{L}$ to maintain $\|\nabla \phi_\theta\|_g \approx 1$ uniformly, yielding relative $\ell^2$ errors of $O(10^{-3} - 10^{-4})$ on benchmarks (Euclid, sphere, Gaussian mixture models). Handling of high-curvature regions involves Ricci-biased sampling; singularities and complex topologies require reparameterisation and reprojection strategies [2404.08754].

## 6. Statistical Consistency and Interpretation

Given probability measure $P$ on $M$, the Fréchet variance is
$$
\operatorname{Var}(P) = \min_{x\in M} \int d^2(x, y)\,dP(y).
$$
Empirically, $F_n(\hat\mu_n)$ strongly consistently estimates $\operatorname{Var}(P)$ as $n\to\infty$. Under uniqueness of the Fréchet mean, almost sure convergence $\hat\mu_n\to\mu$ and $F_n(\hat\mu_n)\to F(\mu)$ is guaranteed [1111.3120]. This affirms MDRM as a fundamental statistic for spread and centrality on Riemannian manifolds.

## 7. Applications and Computational Considerations

MDRM underlies statistical analysis on manifold-valued data in domains such as non-Euclidean machine learning, geometric statistics, and signal processing (e.g., radar target detection using Toeplitz covariance matrices [1111.3120]). The neural Eikonal approach broadens MDRM computation to manifolds with complex metrics or high dimensions, previously intractable with classical solvers. Parameter choices (network depth/width, sampling method) must scale with manifold complexity to ensure accuracy [2404.08754].

---

**References**  
Kelshaw & Magri, “Computing distances and means on manifolds with a metric‐constrained Eikonal approach,” [2404.08754]  
Arnaud et al., “Medians and means in Riemannian geometry: existence, uniqueness and computation,” [1111.3120]

Source: https://www.emergentmind.com/topics/minimum-distance-to-riemannian-mean-mdrm