scholarly journals Open type quasi-Monte Carlo integration based on Halton sequences in weighted Sobolev spaces

2016 ◽  
Vol 33 ◽  
pp. 169-189 ◽  
Author(s):  
Peter Hellekalek ◽  
Peter Kritzer ◽  
Friedrich Pillichshammer
Keyword(s):  
Monte Carlo ◽  
Sobolev Spaces ◽  
Weighted Sobolev Spaces ◽  
Monte Carlo Integration ◽  
Open Type ◽  
Quasi Monte Carlo ◽  
Halton Sequences

scholarly journals A stochastic analysis of greedy routing in a spatially dependent sensor network

2012 ◽  
Vol 23 (4) ◽  
pp. 485-514 ◽  
Author(s):  
HOLGER P. KEELER
Keyword(s):  
Monte Carlo ◽  
Sensor Network ◽  
Monte Carlo Integration ◽  
Greedy Routing ◽  
Quasi Monte Carlo ◽  
Forwarding Node ◽  
Node Location ◽  
Speed Up ◽  
Halton Sequences ◽  
Random Behaviour

For a sensor network, a tractable spatially dependent node deployment model is presented with the property that the density is inversely proportional to the sink distance. A stochastic model is formulated to examine message advancements under greedy routing in such a sensor network. The aim of this work is to demonstrate that an inhomogeneous Poisson process can be used to model a sensor network with spatially dependent node density. Symmetric elliptic integrals and asymptotic approximations are used to describe the random behaviour of hops. Types of dependence that affect hop advancements are examined. We observe that the dependence between successive jumps in a multi-hop path is captured by including only the previous forwarding node location. We include a simple uncoordinated sleep scheme, and observe that the complexity of the model is reduced when sufficiently many nodes are asleep. All expressions involving multi-dimensional integrals are derived and evaluated with quasi-Monte Carlo integration methods based on Halton sequences and recently developed lattice rules. An importance sampling function is derived to speed up the quasi-Monte Carlo methods. The ensuing results agree extremely well with simulations.

scholarly journals Derandomised lattice rules for high dimensional integration

2019 ◽  
Vol 60 ◽  
pp. C247-C260
Author(s):  
Y. Kazashi ◽  
F. Y. Kuo ◽  
I. H. Sloan
Keyword(s):  
Monte Carlo ◽  
Sobolev Spaces ◽  
Monte Carlo Integration ◽  
High Dimensional ◽  
Lattice Rules ◽  
Worst Case ◽  
Multiple Integration ◽  
Quasi Monte Carlo ◽  
Lattice Methods ◽  
Numer Math

We seek shifted lattice rules that are good for high dimensional integration over the unit cube in the setting of an unanchored weighted Sobolev space of functions with square-integrable mixed first derivatives. Many existing studies rely on random shifting of the lattice, whereas here we work with lattice rules with a deterministic shift. Specifically, we consider 'half-shifted' rules in which each component of the shift is an odd multiple of \(1/(2N)\) where \(N\) is the number of points in the lattice. By applying the principle that there is always at least one choice as good as the average, we show that for a given generating vector there exists a half-shifted rule whose squared worst-case error differs from the shift-averaged squared worst-case error by a term of only order \({1/N^2}\). We carry out numerical experiments where the generating vector is chosen component-by-component (CBC), as for randomly shifted lattices, and where the shift is chosen by a new `CBC for shift' algorithm. The numerical results are encouraging. References J. Dick, F. Y. Kuo, and I. H. Sloan. High-dimensional integration: The quasi-Monte Carlo way. Acta Numer., 22:133–288, 2013. doi:10.1017/S0962492913000044. J. Dick, D. Nuyens, and F. Pillichshammer. Lattice rules for nonperiodic smooth integrands. Numer. Math., 126(2):259–291, 2014. doi:10.1007/s00211-013-0566-0. T. Goda, K. Suzuki, and T. Yoshiki. Lattice rules in non-periodic subspaces of sobolev spaces. Numer. Math., 141(2):399–427, 2019. doi:10.1007/s00211-018-1003-1. F. Y. Kuo. Lattice rule generating vectors. URL http://web.maths.unsw.edu.au/ fkuo/lattice/index.html. D. Nuyens and R. Cools. Fast algorithms for component-by-component construction of rank-1 lattice rules in shift-invariant reproducing kernel Hilbert spaces. Math. Comput., 75:903–920, 2006. doi:10.1090/S0025-5718-06-01785-6. I. H. Sloan and S. Joe. Lattice methods for multiple integration. Oxford Science Publications. Clarendon Press and Oxford University Press, 1994. URL https://global.oup.com/academic/product/lattice-methods-for-multiple-integration-9780198534723. I. H. Sloan and H. Wozniakowski. When are quasi-Monte Carlo algorithms efficient for high dimensional integrals? J. Complex., 14(1):1–33, 1998. doi:10.1006/jcom.1997.0463. I. H. Sloan, F. Y. Kuo, and S. Joe. On the step-by-step construction of quasi-Monte Carlo integration rules that achieve strong tractability error bounds in weighted Sobolev spaces. Math. Comput., 71:1609–1641, 2002. doi:10.1090/S0025-5718-02-01420-5.

Higher order Quasi-Monte Carlo integration for Bayesian PDE Inversion

2019 ◽  
Vol 77 (1) ◽  
pp. 144-172 ◽  
Author(s):  
Josef Dick ◽  
Robert N. Gantner ◽  
Quoc T. Le Gia ◽  
Christoph Schwab
Keyword(s):  
Monte Carlo ◽  
Higher Order ◽  
Monte Carlo Integration ◽  
Quasi Monte Carlo

Quasi-Monte Carlo Integration

1995 ◽  
Vol 122 (2) ◽  
pp. 218-230 ◽  
Author(s):  
William J. Morokoff ◽  
Russel E. Caflisch
Keyword(s):  
Monte Carlo ◽  
Monte Carlo Integration ◽  
Quasi Monte Carlo

Multilevel quasi-Monte Carlo integration with product weights for elliptic PDEs with lognormal coefficients

2019 ◽  
Vol 53 (5) ◽  
pp. 1507-1552 ◽  
Author(s):  
L. Herrmann ◽  
C. Schwab
Keyword(s):  
Monte Carlo ◽  
Random Field ◽  
Gaussian Random Field ◽  
Function Spaces ◽  
Monte Carlo Integration ◽  
Superior Performance ◽  
Elliptic Pdes ◽  
Weighted Function ◽  
Weighted Function Spaces ◽  
Quasi Monte Carlo

We analyze the convergence rate of a multilevel quasi-Monte Carlo (MLQMC) Finite Element Method (FEM) for a scalar diffusion equation with log-Gaussian, isotropic coefficients in a bounded, polytopal domain D ⊂ ℝd. The multilevel algorithm QL* which we analyze here was first proposed, in the case of parametric PDEs with sequences of independent, uniformly distributed parameters in Kuo et al. (Found. Comput. Math. 15 (2015) 411–449). The random coefficient is assumed to admit a representation with locally supported coefficient functions, as arise for example in spline- or multiresolution representations of the input random field. The present analysis builds on and generalizes our single-level analysis in Herrmann and Schwab (Numer. Math. 141 (2019) 63–102). It also extends the MLQMC error analysis in Kuo et al. (Math. Comput. 86 (2017) 2827–2860), to locally supported basis functions in the representation of the Gaussian random field (GRF) in D, and to product weights in QMC integration. In particular, in polytopal domains D ⊂ ℝd, d=2,3, our analysis is based on weighted function spaces to describe solution regularity with respect to the spatial coordinates. These spaces allow GRFs and PDE solutions whose realizations become singular at edges and vertices of D. This allows for non-stationary GRFs whose covariance operators and associated precision operator are fractional powers of elliptic differential operators in D with boundary conditions on ∂D. In the weighted function spaces in D, first order, Lagrangian Finite Elements on regular, locally refined, simplicial triangulations of D yield optimal asymptotic convergence rates. Comparison of the ε-complexity for a class of Matérn-like GRF inputs indicates, for input GRFs with low sample regularity, superior performance of the present MLQMC-FEM with locally supported representation functions over alternative representations, e.g. of Karhunen–Loève type. Our analysis yields general bounds for the ε-complexity of the MLQMC algorithm, uniformly with respect to the dimension of the parameter space.

Sign in / Sign up

Export Citation Format

Share Document