1,655 research outputs found
Computation of Electromagnetic Fields Scattered From Objects With Uncertain Shapes Using Multilevel Monte Carlo Method
Computational tools for characterizing electromagnetic scattering from
objects with uncertain shapes are needed in various applications ranging from
remote sensing at microwave frequencies to Raman spectroscopy at optical
frequencies. Often, such computational tools use the Monte Carlo (MC) method to
sample a parametric space describing geometric uncertainties. For each sample,
which corresponds to a realization of the geometry, a deterministic
electromagnetic solver computes the scattered fields. However, for an accurate
statistical characterization the number of MC samples has to be large. In this
work, to address this challenge, the continuation multilevel Monte Carlo
(CMLMC) method is used together with a surface integral equation solver. The
CMLMC method optimally balances statistical errors due to sampling of the
parametric space, and numerical errors due to the discretization of the
geometry using a hierarchy of discretizations, from coarse to fine. The number
of realizations of finer discretizations can be kept low, with most samples
computed on coarser discretizations to minimize computational cost.
Consequently, the total execution time is significantly reduced, in comparison
to the standard MC scheme.Comment: 25 pages, 10 Figure
A Note on the Importance of Weak Convergence Rates for SPDE Approximations in Multilevel Monte Carlo Schemes
It is a well-known rule of thumb that approximations of stochastic partial
differential equations have essentially twice the order of weak convergence
compared to the corresponding order of strong convergence. This is already
known for many approximations of stochastic (ordinary) differential equations
while it is recent research for stochastic partial differential equations. In
this note it is shown how the availability of weak convergence results
influences the number of samples in multilevel Monte Carlo schemes and
therefore reduces the computational complexity of these schemes for a given
accuracy of the approximations.Comment: 16 pages, 3 figures, updated to version published in the Proceedings
of MCQMC1
Central limit theorems for multilevel Monte Carlo methods
In this work, we show that uniform integrability is not a necessary condition
for central limit theorems (CLT) to hold for normalized multilevel Monte Carlo
(MLMC) estimators and we provide near optimal weaker conditions under which the
CLT is achieved. In particular, if the variance decay rate dominates the
computational cost rate (i.e., ), we prove that the CLT applies
to the standard (variance minimizing) MLMC estimator.
For other settings where the CLT may not apply to the standard MLMC
estimator, we propose an alternative estimator, called the mass-shifted MLMC
estimator, to which the CLT always applies.
This comes at a small efficiency loss: the computational cost of achieving
mean square approximation error is at worst a factor
higher with the mass-shifted estimator than
with the standard one
A multiple replica approach to simulate reactive trajectories
A method to generate reactive trajectories, namely equilibrium trajectories
leaving a metastable state and ending in another one is proposed. The algorithm
is based on simulating in parallel many copies of the system, and selecting the
replicas which have reached the highest values along a chosen one-dimensional
reaction coordinate. This reaction coordinate does not need to precisely
describe all the metastabilities of the system for the method to give reliable
results. An extension of the algorithm to compute transition times from one
metastable state to another one is also presented. We demonstrate the interest
of the method on two simple cases: a one-dimensional two-well potential and a
two-dimensional potential exhibiting two channels to pass from one metastable
state to another one
- …