Rapid Grassmannian Averaging with Chebyshev Polynomials
Brighton Ancelin, Alex Saad-Falcon, Kason Ancelin, Justin Romberg

TL;DR
This paper introduces fast algorithms for averaging points on Grassmannian manifolds, significantly reducing computational complexity and enabling efficient applications in machine learning, computer vision, and signal processing.
Contribution
The paper presents novel spectral-based algorithms, RGrAv and DRGrAv, for rapid centralized and decentralized Grassmannian averaging, with theoretical guarantees and superior empirical performance.
Findings
Algorithms outperform existing methods in speed and accuracy.
The methods are effective in tasks like K-means clustering on video data.
The algorithms are versatile for various Grassmannian averaging applications.
Abstract
We propose new algorithms to efficiently average a collection of points on a Grassmannian manifold in both the centralized and decentralized settings. Grassmannian points are used ubiquitously in machine learning, computer vision, and signal processing to represent data through (often low-dimensional) subspaces. While averaging these points is crucial to many tasks (especially in the decentralized setting), existing methods unfortunately remain computationally expensive due to the non-Euclidean geometry of the manifold. Our proposed algorithms, Rapid Grassmannian Averaging (RGrAv) and Decentralized Rapid Grassmannian Averaging (DRGrAv), overcome this challenge by leveraging the spectral structure of the problem to rapidly compute an average using only small matrix multiplications and QR factorizations. We provide a theoretical guarantee of optimality and present numerical experiments…
Peer Reviews
Decision·Submitted to ICLR 2025
- The paper addresses an interesting problem and is generally well-written, presenting the proposed approach clearly, though it occasionally lacks detail. - The authors prove that Chebyshev polynomials are optimal for their modified power method. - Their experiments demonstrate that the algorithm outperforms competing methods in synthetic Grassmannian averaging and K-means clustering of subspaces spanned by video sequences for motion clustering. They have also successfully adapted related method
1. Importance of Decentralized Grassmannian Averaging: The paper emphasizes the importance of decentralized averaging, which the authors claim is significant. However, the paper does not provide concrete examples where decentralization is essential and demonstrates it only through a single set of synthetic experiments (emulating a decentralized setup). Could you provide more compelling examples to illustrate the importance of this problem? 2. Comparison and Relation to Existing Literature: How
omitted.
- The contribution at the technical level is incremental, in my opinion. Theorem 1 seems to be a basic variant of the known result for Chebyshev polynomials. The only difference seems to be that the paper adds an extra constraint $f_t(0)=0$. - The assumption on the dual-banded spectrum is also not new. This assumption is used in the AISTATS 2022 paper: Super-Acceleration with Cyclical Step-sizes (https://proceedings.mlr.press/v151/goujaud22a/goujaud22a.pdf); see Eq. 3 there. That paper has a mu
The paper observes that the eigenvalues of projection matrices have a particular structure, which can be exploited to build a more efficient variant of the power iteration to compute eigenvalues/vectors. As far as I can tell, this contribution is novel and is potentially more widely applicable than explored in the submitted paper. The paper is generally well-written and easy to understand. The developed algorithm is shown to be easy to extend to a decentralized implementation. This is a quite
My main concern with the paper is the studied problem. The paper states that points on the Grassmann manifold are "used ubiquitously in machine learning, computer vision, and signal processing". It is true that some years ago (before deep learning took off), some work went into using subspaces (and hence Grassmannians) to represent batches of data. This approach has, however, not caught on and is quite rare today. I do not want to argue that a paper is uninteresting because it explores a (curre
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Taxonomy
TopicsMathematics and Applications · Advanced Topics in Algebra · Polynomial and algebraic computation
Methodsk-Means Clustering
