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.

Riemannian Preconditioned Coordinate Descent for Low Multilinear Rank Approximation

Riemannian preconditioned coordinate descent for low multi-linear rank approximation

Preconditioned Low-Rank Riemannian Optimization for Symmetric Positive Definite Linear Matrix Equations

Riemannian coordinate descent algorithms on matrix manifolds

Preconditioned low-rank Riemannian optimization for linear systems with tensor product structure

Structure-preserving low multilinear rank approximation of antisymmetric tensors

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

A Preconditioned Riemannian Gradient Descent Algorithm for Low-Rank Matrix Recovery

The Computation of Low Multilinear Rank Approximations of Tensors Via Power Scheme and Random Projection

Projected Gradient Descent Algorithm for Low-Rank Matrix Estimation

Riemannian preconditioned algorithms for tensor completion via tensor ring decomposition

A Quadratically Convergent Algorithm for Structured Low-Rank Approximation

Matrix Nearness Problems with Off-Block-Diagonal Rank Constraints

Decentralized Riemannian natural gradient methods with Kronecker-product approximations

Low-rank approximate inverse for preconditioning tensor-structured linear systems

Learning-Based Low-Rank Approximations

Low-rank plus diagonal approximations for Riccati-like matrix differential equations

Optimization on product manifolds under a preconditioned metric

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

CP decomposition and low-rank approximation of antisymmetric tensors

Jacobi-Davidson method on low-rank matrix manifolds