Search · four archives
Search · four archives
15 papers · ranked by Valyu relevance
Riley Murray, James Demmel, Michael W. Mahoney, N. Benjamin Erichson + 9 more
Randomized numerical linear algebra – RandNLA, for short – concerns the use of randomization as a resource to develop improved algorithms for large-scale linear algebra computations. The origins of contemporary RandNLA lay in theoretical computer science, where it blossomed from a simple idea: randomization provides an…
Fabio V. Difonzo, Paweł Przybyłowicz, Yue Wu, Xinheng Xie
In this paper we investigate the existence, uniqueness and approximation of solutions of delay differential equations (DDEs) with the right-hand side functions f = f(t, x, z) that are Lipschitz continuous with respect to x but only H¨older continuous with respect to (t, z). We give a construction of the randomized…
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…
Nawaf Bou‐Rabee, Tore Selland Kleppe
We introduce 5/2- and 7/2-order L 2 -accurate randomized Runge-Kutta-Nyström methods, tailored for approximating Hamiltonian flows within non-reversible Markov chain Monte Carlo samplers, such as unadjusted Hamiltonian Monte Carlo and unadjusted kinetic Langevin Monte Carlo. We establish quantitative 5/2-order L 2…
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…
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…
Nicolas L. Guidotti, Per-Gunnar Martinsson, Juan A. Acebrón, José Monteiro
'José Monteiro'] Abstract. Many scientific applications require the evaluation of the action of the matrix function over a vector and the most common methods for this task are those based on the Krylov subspace. Since the orthogonalization cost and memory requirement can quickly become overwhelming as the basis grows…
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.…
H. P. Wang, Jing Lyu, Zhigang Peng, Min Li
This paper analyzes the SParareal algorithm for stochastic differential equations (SDEs). Compared to the classical Parareal algorithm, the SParareal algorithm accelerates convergence by introducing stochastic perturbations, achieving linear convergence over unbounded time intervals. We first revisit the classical…
Guanjie Wang, Qifeng Liao
equations with random inputs Authors: ['Guanjie Wang' 'Qifeng Liao'] It is known that standard stochastic Galerkin methods face challenges when solving partial differential equations (PDEs) with random inputs. These challenges are typically attributed to the large number of required physical basis functions and…
Julian Clausnitzer, Andreas Kleefeld
We consider the numerical approximation of different ordinary differential equations (ODEs) and partial differential equations (PDEs) with periodic boundary conditions involving a one-dimensional random parameter, comparing the intrusive and non-intrusive polynomial chaos expansion (PCE) method. We demonstrate how to…
Emil Løvbak, Frédéric Blondeel, Adam J. Lee, Lander Vanroye + 2 more
'Andreas Van Barel' 'Giovanni Samaey'] Abstract In PDE-constrained optimization, one aims to find design parameters that minimize some objective, subject to the satisfaction of a partial differential equation. A major challenge is computing gradients of the objective to the design parameters, as applying the chain rule…
Prerna Patil, Hessam Babaee
Low-rank approximation using time-dependent bases (TDBs) has proven effective for reduced-order modeling of stochastic partial differential equations (SPDEs). In these techniques, the random field is decomposed to a set of deterministic TDBs and time-dependent stochastic coefficients. When applied to SPDEs with…
Irina-Beatrice Haas, Michael B. Giles
Multilevel Monte Carlo (MLMC) reduces the total computational cost of financial option pricing by combining SDE approximations with multiple resolutions. This paper explores a further avenue for reducing cost and improving power efficiency through the use of low precision calculations on configurable hardware devices…
Julien Genovese, Francesco Ballarin, Gianluigi Rozza, Claudio Canuto
AbstractIn this manuscript we propose and analyze weighted reduced order methods for stochastic Stokes and Navier-Stokes problems depending on random input data (such as forcing terms, physical or geometrical coefficients, boundary conditions). We will compare weighted methods such as weighted greedy and weighted POD…