Abstract:Interpolative and CUR decompositions involve "natural bases" of row and column subsets, or skeletons, of a given matrix that approximately span its row and column spaces. These low-rank decompositions preserve properties such as sparsity or non-negativity, and are easily interpretable in the context of the original data. For large-scale problems, randomized sketching to sample the row or column spaces with a random matrix can serve as an effective initial step in skeleton selection to reduce computational cost. A by now well-established approach has been to extract a randomized sketch, followed by column-pivoted QR CPQR) on the sketch matrix. This manuscript describes an alternative approach where CPQR is replaced by LU with partial pivoting (LUPP). While LUPP by itself is not rank-revealing, it is demonstrated that when used in a randomized setting, LUPP not only reveals the numerical rank, but also allows the estimation of the residual error as the factorization is built. The resulting algorithm is both adaptive and parallelizable, and attains much higher practical speed due to the lower communication requirements of LUPP over CPQR. The method has been implemented for both CPUs and GPUs, and the resulting software has been made publicly available.

LinApart: optimizing the univariate partial fraction decomposition

A Fast Algorithm for Partial Fraction Decompositions

Efficient Recursive Methods for Partial Fraction Expansion of General Rational Functions.

Modified Affine Arithmetic in Tensor Form for Trivariate Polynomial Evaluation and Algebraic Surface Plotting

Partial fraction decompositions and an algorithm for computing the vector partition function

Computation of polynomial and rational approximations in complex domains by the -method

Blade: A package for block-triangular form improved Feynman integrals decomposition

Rational Function Simplification for Integration-by-Parts Reduction and Beyond

On Exact and Approximate Interpolation of Sparse Rational Functions

Parallel Optimization for Accelerating the Generation of Correctly Rounded Elementary Functions

A unified presentation of certain families of non-Fuchsian differential equations via fractional calculus operators

Adaptive Parallelizable Algorithms for Interpolative Decompositions via Partially Pivoted LU

Recursive sparse LU decomposition based on nested dissection and low rank approximations

Approximate factorization of multivariate polynomials via differential equations

Decomposing Solution Sets of Polynomial Systems: A New Parallel Monodromy Breakup Algorithm

On the partial realization problem

Efficient reduction of Feynman integrals on supercomputers

Efficient MPI parallel interpolative decomposition

Rigorous numerical integration of algebraic functions

A new approach to MacMahon's Partition Analysis and associated applications

Symbolic Computation of Analytic Approximate Solutions for Nonlinear Differential Equations with Initial Conditions