aiqbench analysis
This research presents a scalable stochastic algorithm for recovering latent low-rank structures in sparse nonnegative matrices, addressing computational bottlenecks in large-scale matrix completion. By leveraging block-wise operations and CUDA acceleration, the method enables efficient analysis of massive datasets, such as the *Drosophila* connectome, without requiring the storage of dense intermediate products. The work advances machine learning infrastructure for handling high-dimensional sparse data where traditional alternating minimization approaches become computationally prohibitive.
Subzero matrix completion for sparse data analysis: large-scale learning of latent low-rank structure
arXiv:2608.21607v1 Announce Type: new Abstract: We investigate when a sparse nonnegative matrix can be recovered from a real-valued matrix of much lower rank by zeroing out its negative elements. The potential for such decompositions suggests a mathematical connection between sparsity and rank; we analyze a number of sparse matrices with this latent low-rank structure and use them to illustrate the geometric origins of this connection. Previous algorithms have discovered these decompositions via an alternating minimization over the factors of a low-rank matrix, but to do so, they have also needed to compute and store another matrix, neither sparse nor low-rank, that is the size of their product. We develop a stochastic, alternating least-squares algorithm that operates on smaller blocks of this dense matrix and scales as a result to much larger problems. We also show how to further accelerate this algorithm with sparse optimizations and customized CUDA kernels. As one example, we use the algorithm to analyze the sparse matrix of synaptic weights for the recently published $\textit{Drosphilia}$ connectome. The nonzero elements of this matrix, with 139,255 rows and columns, record the number of synapses between cells in the nervous system of a female fruit fly. Despite a slowly decaying spectrum of singular values, this matrix exhibits a latent low-rank structure that is predictive of cell categories across multiple levels of specificity.
Abstract:We investigate when a sparse nonnegative matrix can be recovered from a real-valued matrix of much lower rank by zeroing out its negative elements. The potential for such decompositions suggests a mathematical connection between sparsity and rank; we analyze a number of sparse matrices with this latent low-rank structure and use them to illustrate the geometric origins of this connection. Previous algorithms have discovered these decompositions via an alternating minimization over the factors of a low-rank matrix, but to do so, they have also needed to compute and store another matrix, neither sparse nor low-rank, that is the size of their product. We develop a stochastic, alternating least-squares algorithm that operates on smaller blocks of this dense matrix and scales as a result to much larger problems. We also show how to further accelerate this algorithm with sparse optimizations and customized CUDA kernels. As one example, we use the algorithm to analyze the sparse matrix of synaptic weights for the recently published $\textit{Drosphilia}$ connectome. The nonzero elements of this matrix, with 139,255 rows and columns, record the number of synapses between cells in the nervous system of a female fruit fly. Despite a slowly decaying spectrum of singular values, this matrix exhibits a latent low-rank structure that is predictive of cell categories across multiple levels of specificity.
Machine Learning (cs.LG); Machine Learning (stat.ML)
arXiv:2608.21607 [cs.LG]
(or arXiv:2608.21607v1 [cs.LG] for this version)
https://doi.org/10.48550/arXiv.2608.21607 Focus to learn more arXiv-issued DOI via DataCite (pending registration)
arXivLabs is a framework that allows collaborators to develop and share new arXiv features directly on our website.
Both individuals and organizations that work with arXivLabs have embraced and accepted our values of openness, community, excellence, and user data privacy. arXiv is committed to these values and only works with partners that adhere to them.
Have an idea for a project that will add value for arXiv's community? Learn more about arXivLabs.
Full text under the publisher's open licence · 314 words
Related coverage

Generalization and infinite width
LessWrong · Aug 1, 2026
ROCm 10.1: Breaking the Data-Movement Bottleneck
AMD · Oct 5, 2026

CUDA Toolkit 13.4 Adds Windows on Arm Support and Greater Control over Shared GPUs
NVIDIA · Sep 9, 2026

Co-Designing AI Model Attention for Fast, Interactive Long-Context Inference
NVIDIA · Jul 31, 2026