Sketching with Kerdock's crayons: Fast sparsifying transforms for arbitrary linear maps

Abstract

Given an arbitrary matrix A∈Rn× n, we consider the fundamental problem of computing Ax for any x∈Rn such that Ax is s-sparse. While fast algorithms exist for particular choices of A, such as the discrete Fourier transform, there is currently no o(n2) algorithm that treats the unstructured case. In this paper, we devise a randomized approach to tackle the unstructured case. Our method relies on a representation of A in terms of certain real-valued mutually unbiased bases derived from Kerdock sets. In the preprocessing phase of our algorithm, we compute this representation of A in O(n3 n) operations. Next, given any unit vector x∈Rn such that Ax is s-sparse, our randomized fast transform uses this representation of A to compute the entrywise ε-hard threshold of Ax with high probability in only O(sn + ε-2\|A\|2∞2n n) operations. In addition to a performance guarantee, we provide numerical results that demonstrate the plausibility of real-world implementation of our algorithm.

0

Turn this paper into a lesson

ArcXiv compiles a structured reading guide from this paper's metadata: plain-English importance, contributions, prerequisite concepts, which sections to read first, flashcards, and a quiz. Grounded in the abstract, never invented.

Discussion (0)

Sign in to join the discussion.

Loading comments…