Abstract:The problem of determining the configuration of points from partial distance information, known as the Euclidean Distance Geometry (EDG) problem, is fundamental to many tasks in the applied sciences. In this paper, we propose two algorithms grounded in the Riemannian optimization framework to address the EDG problem. Our approach formulates the problem as a low-rank matrix completion task over the Gram matrix, using partial measurements represented as expansion coefficients of the Gram matrix in a non-orthogonal basis. For the first algorithm, under a uniform sampling with replacement model for the observed distance entries, we demonstrate that, with high probability, a Riemannian gradient-like algorithm on the manifold of rank-$r$ matrices converges linearly to the true solution, given initialization via a one-step hard thresholding. This holds provided the number of samples, $m$, satisfies $m \geq \mathcal{O}(n^{7/4}r^2 \log(n))$. With a more refined initialization, achieved through resampled Riemannian gradient-like descent, we further improve this bound to $m \geq \mathcal{O}(nr^2 \log(n))$. Our analysis for the first algorithm leverages a non-self-adjoint operator and depends on deriving eigenvalue bounds for an inner product matrix of restricted basis matrices, leveraging sparsity properties for tighter guarantees than previously established. The second algorithm introduces a self-adjoint surrogate for the sampling operator. This algorithm demonstrates strong numerical performance on both synthetic and real data. Furthermore, we show that optimizing over manifolds of higher-than-rank-$r$ matrices yields superior numerical results, consistent with recent literature on overparameterization in the EDG problem.

Euclidean Representation of Low-Rank Matrices and Its Geometric Properties

Metric Transforms and Low Rank Matrices via Representation Theory of the Real Hyperrectangle

Recovery of Low-Rank Matrices Based on the Rank Null Space Properties.

Riemannian Optimization for Non-convex Euclidean Distance Geometry with Global Recovery Guarantees

Permutation-Invariant Log-Euclidean Geometries on Full-Rank Correlation Matrices

On the geometric interpretation of the nonnegative rank

A Quadratically Convergent Algorithm for Structured Low-Rank Approximation

A general error analysis for randomized low-rank approximation with application to data assimilation

Data-Driven Linear Complexity Low-Rank Approximation of General Kernel Matrices: A Geometric Approach

High curvature means low-rank: On the sectional curvature of Grassmann and Stiefel manifolds and the underlying matrix trace inequalities

When big data actually are low-rank, or entrywise approximation of certain function-generated matrices

Multilinear Low-Rank Tensors on Graphs & Applications

Exploiting the Structure Effectively and Efficiently in Low Rank Matrix Recovery.

Using oriented matroids to find low rank structure in presence of nonlinearity

Generalized Matrix Local Low Rank Representation by Random Projection and Submatrix Propagation

Riemannian Optimization on Tensor Products of Grassmann Manifolds: Applications to Generalized Rayleigh-Quotients

Bures-Wasserstein minimizing geodesics between covariance matrices of different ranks

No Spurious Local Minima in Nonconvex Low Rank Problems: A Unified Geometric Analysis

The Global Geometry of Centralized and Distributed Low-rank Matrix Recovery without Regularization

General Low-rank Matrix Optimization: Geometric Analysis and Sharper Bounds

Scalable Algorithms for Locally Low-Rank Matrix Modeling