Single-pass Nystr\"{o}m approximation in mixed precision

From MaRDI portal
Publication:6400255




Abstract: Low rank matrix approximations appear in a number of scientific computing applications. We consider the Nystr"{o}m method for approximating a positive semidefinite matrix A. The computational cost of its single-pass version can be decreased by running it in mixed precision, where the expensive products with A are computed in a precision lower than the working precision. We bound the extra finite precision error which is compared to the error of the Nystr"{o}m approximation in exact arithmetic and develop a heuristic to identify when the approximation quality is not affected by the low precision computation. Further, the mixed precision Nystr"{o}m method can be used to inexpensively construct a limited memory preconditioner for the conjugate gradient method. We bound the condition number of the resulting preconditioned coefficient matrix, and experimentally show that such a preconditioner can be effective.














This page was built for publication: Single-pass Nystr\"{o}m approximation in mixed precision

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6400255)