Improved Sparse Recovery for Approximate Matrix Multiplication
Abstract
We present a simple randomized algorithm for approximate matrix multiplication (AMM) whose error scales with the *output* norm . Given any matrices and a runtime parameter , the algorithm produces in time, a matrix with total squared error , per-entry variance and bias . Alternatively, the algorithm can compute an *unbiased* estimation with expected total squared error , recovering the state-of-art AMM error obtained by Pagh's TensorSketch algorithm (Pagh, 2013). Our algorithm is a log-factor faster. The key insight in the algorithm is a new variation of pseudo-random rotation of the input matrices (a Fast Hadamard Transform with asymmetric diagonal scaling), which redistributes the Frobenius norm of the *output* uniformly across its entries.
Cite
@article{arxiv.2602.04386,
title = {Improved Sparse Recovery for Approximate Matrix Multiplication},
author = {Yahel Uffenheimer and Omri Weinstein},
journal= {arXiv preprint arXiv:2602.04386},
year = {2026}
}