Abstract:SIAM Journal on Matrix Analysis and Applications, Volume 45, Issue 2, Page 967-991, June 2024. This paper combines modern numerical computation with theoretical results to improve our understanding of the growth factor problem for Gaussian elimination. On the computational side we obtain lower bounds for the maximum growth for complete pivoting for [math] and [math] using the Julia JuMP optimization package. At [math] we obtain a growth factor bigger than [math]. The numerical evidence suggests that the maximum growth factor is bigger than [math] if and only if [math]. We also present a number of theoretical results. We show that the maximum growth factor over matrices with entries restricted to a subset of the reals is nearly equal to the maximum growth factor over all real matrices. We also show that the growth factors under floating point arithmetic and exact arithmetic are nearly identical. Finally, through numerical search, and stability and extrapolation results, we provide improved lower bounds for the maximum growth factor. Specifically, we find that the largest growth factor is bigger than [math] for [math], and the lim sup of the ratio with [math] is greater than or equal to [math]. In contrast to the old conjecture that growth might never be bigger than [math], it seems likely that the maximum growth divided by [math] goes to infinity as [math].

Gaussian Elimination with Randomized Complete Pivoting

Growth Factors of Orthogonal Matrices and Local Behavior of Gaussian Elimination with Partial and Complete Pivoting

Average-case analysis of the Gaussian Elimination with Partial Pivoting

Distribution of the number of pivots needed using Gaussian elimination with partial pivoting on random matrices

Gaussian elimination corrects pivoting mistakes

A New Upper Bound For the Growth Factor in Gaussian Elimination with Complete Pivoting

Complete pivoting growth of butterfly matrices and butterfly Hadamard matrices

Some New Results on the Maximum Growth Factor in Gaussian Elimination

Randomized Complete Pivoting for Solving Symmetric Indefinite Linear Systems

Adaptive Parallelizable Algorithms for Interpolative Decompositions via Partially Pivoted LU

Formal Derivation of LU Factorization with Pivoting

On a perturbation analysis of Higham squared maximum Gaussian elimination growth matrices

Random Multipliers Numerically Stabilize Gaussian and Block Gaussian Elimination: Proofs and an Extension to Low-rank Approximation

Supporting GENP with Random Multipliers

CholeskyQR with Randomization and Pivoting for Tall Matrices (CQRRPT)

Growth factors of random butterfly matrices and the stability of avoiding pivoting

Rank revealing Gaussian elimination by the maximum volume concept

Iterative methods for solving factorized linear systems

Randomly Pivoted Partial Cholesky: Random How?