arrow
Return

Generalizing Random Butterfly Transforms to Arbitrary Matrix Sizes

delete2024-12-11
delete0
delete
OA
AI
N
Neil Lindquist
P
Piotr Łuszczek
J
Jack Dongarra
DOI:10.1145/3699714delete
deleteOriginal
deleteOriginal request for help
deleteShare
deleteSave
Abstract

Abstract

En 中文
Parker and L & ecirc; introduced random butterfly transforms (RBTs) as a preprocessing technique to replace pivoting in dense LU factorization. Unfortunately, their FFT-like recursive structure restricts the dimensions of the matrix. Furthermore, on multinode systems, efficient management of the communication overheads restricts the matrix's distribution even more. To remove these limitations, we have generalized the RBT to arbitrary matrix sizes by truncating the dimensions of each layer in the transform. We expanded Parker's theoretical analysis to generalized RBT, specifically that in exact arithmetic, Gaussian elimination with no pivoting will succeed with probability 1 after transforming a matrix with full-depth RBTs. Furthermore, we experimentally show that these generalized transforms improve performance over Parker's formulation by up to 62% while retaining the ability to replace pivoting. This generalized RBT is available in the SLATE numerical software library.
Keywords:
Gaussian Elimination
Randomization

Journal

ACM Transactions on Mathematical Software cover
ACM Transactions on Mathematical Software
IF:
3.2
Papers:
34
Citations:
5.1K

Organization

U
Univ Tennessee
Scholars:
1.2K
Papers: 1.8K
Citations: 192