2,397 research outputs found
Lectures on Randomized Numerical Linear Algebra
This chapter is based on lectures on Randomized Numerical Linear Algebra from
the 2016 Park City Mathematics Institute summer school on The Mathematics of
Data.Comment: To appear in the edited volume of lectures from the 2016 PCMI summer
schoo
Random Sampling in Computational Algebra: Helly Numbers and Violator Spaces
This paper transfers a randomized algorithm, originally used in geometric
optimization, to computational problems in commutative algebra. We show that
Clarkson's sampling algorithm can be applied to two problems in computational
algebra: solving large-scale polynomial systems and finding small generating
sets of graded ideals. The cornerstone of our work is showing that the theory
of violator spaces of G\"artner et al.\ applies to polynomial ideal problems.
To show this, one utilizes a Helly-type result for algebraic varieties. The
resulting algorithms have expected runtime linear in the number of input
polynomials, making the ideas interesting for handling systems with very large
numbers of polynomials, but whose rank in the vector space of polynomials is
small (e.g., when the number of variables and degree is constant).Comment: Minor edits, added two references; results unchange
Fast randomized iteration: diffusion Monte Carlo through the lens of numerical linear algebra
We review the basic outline of the highly successful diffusion Monte Carlo
technique commonly used in contexts ranging from electronic structure
calculations to rare event simulation and data assimilation, and propose a new
class of randomized iterative algorithms based on similar principles to address
a variety of common tasks in numerical linear algebra. From the point of view
of numerical linear algebra, the main novelty of the Fast Randomized Iteration
schemes described in this article is that they work in either linear or
constant cost per iteration (and in total, under appropriate conditions) and
are rather versatile: we will show how they apply to solution of linear
systems, eigenvalue problems, and matrix exponentiation, in dimensions far
beyond the present limits of numerical linear algebra. While traditional
iterative methods in numerical linear algebra were created in part to deal with
instances where a matrix (of size ) is too big to store, the
algorithms that we propose are effective even in instances where the solution
vector itself (of size ) may be too big to store or manipulate.
In fact, our work is motivated by recent DMC based quantum Monte Carlo schemes
that have been applied to matrices as large as . We
provide basic convergence results, discuss the dependence of these results on
the dimension of the system, and demonstrate dramatic cost savings on a range
of test problems.Comment: 44 pages, 7 figure
- …