Search · four archives
Search · four archives
15 papers · ranked by Valyu relevance
Raphael Kruse, Yue Wu
This paper contains an error analysis of two randomized explicit Runge-Kutta schemes for ordinary differential equations (ODEs) with timeirregular coefficient functions. In particular, the methods are applicable to ODEs of Carath´eodory type, whose coefficient functions are only integrable with respect to the time…
Inês A. Ferreira, Juan A. Acebrón, José Monteiro
The Kaczmarz algorithm is an iterative method that solves linear systems of equations. It stands out among iterative algorithms when dealing with large systems for two reasons. First, at each iteration, the Kaczmarz algorithm uses a single equation, resulting in minimal computational work per iteration. Second, solving…
Ruoqi Shen, Yin Tat Lee
Sampling from log-concave distributions is a well researched problem that has many applications in statistics and machine learning. We study the distributions of the form p ∗ ∝ exp(−f(x)), where f : R d → R has an L-Lipschitz gradient and is m-strongly convex. In our paper, we propose a Markov chain Monte Carlo (MCMC)…
Art B. Owen
Randomized quasi-Monte Carlo (RQMC) sampling can bring orders of magnitude reduction in variance compared to plain Monte Carlo (MC) sampling. The extent of the efficiency gain varies from problem to problem and can be hard to predict. This article presents an R function rhalton that produces scrambled versions of…
Chuan-gang Kang, Heng Zhou
The Kaczmarz method is an iterative projection scheme for solving consistent system Ax = b. It is later extended to the inconsistent and ill-posed linear problems. But the classical Kaczmarz method is sensitive to the correlation of the adjacent equations. In order to reduce the impact of correlation on the convergence…
Haochen Jiang, Dongdong Liu, Xianping Wu, Yang Xu
Motivated by the randomized sketch to solve a variety of problems in scientific computation, we improve both the maximal weighted residual Kaczmarz method and the randomized block average Kaczmarz method using two new randomized sketch techniques. Besides, convergence analyses of the proposed methods are provided.…
Jean-Jacques Forneron
Quasi-Monte Carlo (qMC) methods are a powerful alternative to classical Monte-Carlo (MC) integration. Under certain conditions, they can approximate the desired integral at a faster rate than the usual Central Limit Theorem, resulting in more accurate estimates. This paper explores these methods in a simulation-based…
Qiang Ji
Monte Carlo simulations are widely used in many areas including particle accelerators. In this lecture, after a short introduction and reviewing of some statistical backgrounds, we will discuss methods such as direct inversion, rejection method, and Markov chain Monte Carlo to sample a probability distribution…
Guanjie Wang, Qifeng Liao
We present a reduced basis stochastic Galerkin method for partial differential equations with random inputs. In this method, the reduced basis methodology is integrated into the stochastic Galerkin method, such that the cost of solvers for the Galerkin system is significantly reduced. To reduce the main cost of…
Paul Kabaila
Cranley and Patterson put forward the following randomization as the basis for the estimation of the error of a lattice rule for an integral of a one-periodic function over the unit cube in s dimensions. The lattice rule is randomized using independent random shifts in each coordinate direction that are uniformly…
Pelin Çi̇loğlu, Hamdullah Yücel
We investigate numerical behaviour of a convection diffusion equation with random coefficients by approximating statistical moments of the solution. Stochastic Galerkin approach, turning the original stochastic problem to a system of deterministic convection diffusion equations, is used to handle the stochastic domain…
Migran N. Gevorkyan, Anastasia V. Demidova, Anna V. Korolkova, Dmitry S. Kulyabov
'Dmitry S. Kulyabov'] This paper discusses stochastic numerical methods of Runge-Kutta type with weak and strong convergences for systems of stochastic differential equations in Itˆo form. At the beginning we give a brief overview of the stochastic numerical methods and information from the theory of stochastic…
David Arnas, Carl Leake, Daniele Mortari
This work introduces two new techniques for random number generation with any prescribed nonlinear distribution based on the k-vector methodology. The first approach is based on inverse transform sampling using the optimal k-vector to generate the samples by inverting the cumulative distribution. The second approach…
Ramakrishna Tipireddy, Eric Phipps, Roger Ghanem
In this work, solvers and preconditioners based on Gauss-Seidel and Jacobi algorithms are explored for stochastic Galerkin discretization of partial differential equations (PDEs) with random input data. Gauss-Seidel and Jacobi algorithms are formulated such that the existing software is leveraged in the computational…
Jehanzeb H. Chaudhry, Nathanial Burch, Donald Estep
We study the problem of uncertainty quantification for the numerical solution of elliptic partial differential equation boundary value problems posed on domains with stochastically varying boundaries. We also use the uncertainty quantification results to tackle the efficient solution of such problems. We introduce…