Search papers, labs, and topics across Lattice.
This paper investigates the robustness of randomized sketching techniques, specifically SparseStack, to low-precision rounding on GPUs, which is crucial for optimizing performance in sparse numerical linear algebra. The authors find that the embedding quality remains largely unaffected by the choice of FP16 rounding method, with all variants yielding similar accuracy across various input types. Notably, deterministic rounding offers the best performance-accuracy tradeoff, highlighting the importance of sketch distribution over quantization rules in maintaining subspace distortion and accuracy.
SparseStack's embedding quality remains stable across different FP16 rounding methods, challenging assumptions about low-precision impacts on performance.
Randomized sketching is a core primitive in randomized numerical linear algebra. On modern hardware architectures, in particular on GPUs, the performance of sparse sketches is limited by memory traffic and atomic accumulation rather than floating-point throughput. This makes sketching a natural target for mixed precision, provided that low-precision accumulation does not degrade the embedding quality. We study mixed-precision GPU implementations of sparse oblivious subspace embeddings, focusing on a SparseStack generalization of the GPU CountSketch kernel of Higgins et al. SparseStack improves embedding quality relative to CountSketch on coherent inputs, but its additional nonzeros per column increase atomic-update contention and reduce throughput. We therefore implement FP16 SparseStack variants using deterministic round-to-nearest, exact stochastic rounding, and dithered rounding, and compare them with FP32 SparseStack, CountSketch, mixed-precision CountSketch, and FlashSketch. Our main empirical finding is that, for the tested regimes, SparseStack embedding quality is insensitive to the FP16 rounding rule. Deterministic, stochastic, and dithered rounding FP16 SparseStack produce nearly identical subspace distortion and sketch-and-solve least-squares accuracy across incoherent, coherent, and adversarial test problems. The dominant accuracy factor is the sketch distribution rather than the quantization rule: SparseStack variants substantially improve distortion on coherent inputs, while all methods behave similarly on incoherent inputs. Since deterministic rounding has the lowest overhead, it provides the best performance--accuracy tradeoff among the FP16 SparseStack variants.
