Algorithms for numerical integration in high and infinite dimensions: Analysis, applications and implementation

dc.contributor.advisor Kuo, Frances en_US
dc.contributor.advisor Sloan, Ian en_US Gilbert, Alexander en_US 2022-03-22T17:24:33Z 2022-03-22T17:24:33Z 2018 en_US
dc.description.abstract Approximating high and infinite dimensional integrals numerically is in general a very difficult problem. However, it is also one that arises in several applications from statistics, finance and uncertainty quantification, thus motivating a real need for the development and analysis of efficient algorithms. The difficulty lies in the fact that in general high-dimensional problems suffer from the curse of dimensionality where the cost of an approximation rises exponentially with dimension. However, knowing certain properties of the integrands allows one to identify problems that do not suffer from the curse and for which efficient algorithms can be developed. In this thesis we study numerical integration algorithms, specifically Quasi-Monte Carlo (QMC) quadrature rules and the Multivariate Decomposition Method (MDM), when bounds on the first mixed derivatives are known. The focus in this thesis is on analysis and development of algorithms, a new application for QMC methods from the field of uncertainty quantification and efficient strategies for implementing numerical integration algorithms. The main results of this thesis are as follows. First, we present a full error analysis for the application of QMC methods to approximate the expectation of the smallest eigenvalue of an elliptic differential operator with coefficients that are parametrised by infinitely-many stochastic variables. Eigenvalue problems are used to model many physical situations in engineering and the natural sciences, and this problem is motivated by uncertainty quantification of such problems. It also represents a new application for QMC methods. Second, we provide explicit details and numerical results on how to efficiently implement the Multivariate Decomposition Method (MDM) for approximating infinite-dimensional integrals. The third contribution of this thesis is a new method of constructing optimal active sets for use in the MDM. Finally, we present two user-friendly Component-by-Component algorithms for constructing QMC lattice rules, which automatically choose good function space weight parameters. In all cases we present numerical results that display the advantages of the algorithms, and where appropriate substantiate our corresponding theoretical results. en_US
dc.language English
dc.language.iso EN en_US
dc.publisher UNSW, Sydney en_US
dc.rights CC BY-NC-ND 3.0 en_US
dc.rights.uri en_US
dc.subject.other High-dimensional integration en_US
dc.subject.other Numerical integration en_US
dc.subject.other Quasi-Monte Carlo en_US
dc.subject.other Infinite-dimensional integration en_US
dc.title Algorithms for numerical integration in high and infinite dimensions: Analysis, applications and implementation en_US
dc.type Thesis en_US
dcterms.accessRights open access
dcterms.rightsHolder Gilbert, Alexander
dspace.entity.type Publication en_US
unsw.relation.faculty Science
unsw.relation.originalPublicationAffiliation Gilbert, Alexander, Mathematics & Statistics, Faculty of Science, UNSW en_US
unsw.relation.originalPublicationAffiliation Kuo, Frances, Mathematics & Statistics, Faculty of Science, UNSW en_US
unsw.relation.originalPublicationAffiliation Sloan, Ian, Mathematics & Statistics, Faculty of Science, UNSW en_US School of Mathematics & Statistics *
unsw.thesis.degreetype PhD Doctorate en_US
Original bundle
Now showing 1 - 1 of 1
No Thumbnail Available
public version.pdf
1.22 MB
Resource type