Complete pivoting growth of butterfly matrices and butterfly Hadamard matrices
read the original abstract
The growth problem in Gaussian elimination (GE) remains a foundational question in numerical analysis and numerical linear algebra. Wilkinson resolved the growth problem in GE with partial pivoting (GEPP) in his initial analysis from the 1960s, while he was only able to establish an upper bound for the GE with complete pivoting (GECP) growth problem. The GECP growth problem has seen a spike in recent interest, culminating in improved lower and upper bounds established by Bisain, Edelman, and Urschel in 2023, but still remains far from being fully resolved. Due to the complex dynamics governing the location of GECP pivots, analysis of GECP growth for particular input matrices often estimates the actual growth rather than computes the growth exactly. We present a class of dense random butterfly matrices for which we can compute the exact GECP growth. We extend previous results that established exact growth computations for butterfly matrices when using GEPP and GE with rook pivoting (GERP) to now also include GECP for structured subclasses of inputs. Moreover, we present a new method to construct random Hadamard matrices using butterfly matrices.
This paper has not been read by Pith yet.
Forward citations
Cited by 1 Pith paper
-
The Horton-Strahler number of butterfly trees
Butterfly trees have Horton-Strahler numbers scaling linearly with depth: exactly pq/(1-pq) per level in the simple model via an 8-state Markov chain, and approximately 0.445 in the general uniform model per simulations.
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.