
arXiv: 1509.06820
The dominant contribution to communication complexity in factorizing a matrix using QR with column pivoting is due to column-norm updates that are required to process pivot decisions. We use randomized sampling to approximate this process which dramatically reduces communication in column selection. We also introduce a sample update formula to reduce the cost of sampling trailing matrices. Using our column selection mechanism we observe results that are comparable in quality to those obtained from the QRCP algorithm, but with performance near unpivoted QR. We also demonstrate strong parallel scalability on shared memory multiple core systems using an implementation in Fortran with OpenMP. This work immediately extends to produce low-rank truncated approximations of large matrices. We propose a truncated QR factorization with column pivoting that avoids trailing matrix updates which are used in current implementations of level-3 BLAS QR and QRCP. Provided the truncation rank is small, avoiding trailing matrix updates reduces approximation time by nearly half. By using these techniques and employing a variation on Stewart's QLP algorithm, we develop an approximate truncated SVD that runs nearly as fast as truncated QR.
low-rank approximation, Numerical solutions to overdetermined systems, pseudoinverses, blocked algorithm, truncated SVD, random sampling, QR factorization, singular value decomposition, column pivoting, Parallel numerical computation, Numerical Analysis (math.NA), Direct numerical methods for linear systems and matrix inversion, FOS: Mathematics, Mathematics - Numerical Analysis, parallel computation, sample update
low-rank approximation, Numerical solutions to overdetermined systems, pseudoinverses, blocked algorithm, truncated SVD, random sampling, QR factorization, singular value decomposition, column pivoting, Parallel numerical computation, Numerical Analysis (math.NA), Direct numerical methods for linear systems and matrix inversion, FOS: Mathematics, Mathematics - Numerical Analysis, parallel computation, sample update
| selected citations These citations are derived from selected sources. This is an alternative to the "Influence" indicator, which also reflects the overall/total impact of an article in the research community at large, based on the underlying citation network (diachronically). | 58 | |
| popularity This indicator reflects the "current" impact/attention (the "hype") of an article in the research community at large, based on the underlying citation network. | Top 10% | |
| influence This indicator reflects the overall/total impact of an article in the research community at large, based on the underlying citation network (diachronically). | Top 10% | |
| impulse This indicator reflects the initial momentum of an article directly after its publication, based on the underlying citation network. | Top 1% |
