An efficient volume-preserving MBO scheme for data clustering and classification
Fabius Krämer, Tim Laux
TL;DR
The paper advances clustering and classification on graphs by formulating a volume-preserving MBO scheme that enforces exact cluster volumes via a vector-valued $V$-order statistic. It provides an exact, efficient algorithm to compute the volume-constrained thresholding step, with worst-case complexity ${O}(N(\log N + P)P^2)$ and improved, data-driven running times under a gradient-flow big-data regime, achieving ${O}(\sqrt{h}\,N\log N)$ per iteration in favorable settings. A rigorous variational analysis connects the discrete scheme to volume-preserving mean curvature flow, establishing convergence of discrete order statistics to continuous counterparts on manifolds and deriving ${L^2}$-bounds for Lagrange multipliers. The approach is complemented by extensive numerics across multiple diffusion kernels and datasets, showing competitive accuracy with favorable running times, and the authors release public code to facilitate adoption. Overall, the work delivers a principled, scalable framework for volume-aware graph clustering and semi-/unsupervised classification with strong theoretical and practical implications.
Abstract
We propose and study a novel efficient algorithm for clustering and classification tasks based on the famous MBO scheme. On the one hand, inspired by Jacobs et al. [J. Comp. Phys. 2018], we introduce constraints on the size of clusters leading to a linear integer problem. We prove that the solution to this problem is induced by a novel order statistic. This viewpoint allows us to develop exact and highly efficient algorithms to solve such constrained integer problems. On the other hand, we prove an estimate of the computational complexity of our scheme, which is better than any available provable bounds for the state of the art. This rigorous analysis is based on a variational viewpoint that connects this scheme to volume-preserving mean curvature flow in the big data and small time-step limit.
