Abstract:A preconditioning strategy is proposed for the iterative solve of large numbers of linear systems with variable matrix and right-hand side which arise during the computation of solution statistics of stochastic elliptic partial differential equations with random variable coefficients sampled by Monte Carlo. Building on the assumption that a truncated Karhunen-Loève expansion of a known transform of the random variable coefficient is known, we introduce a compact representation of the random coefficient in the form of a Voronoi quantizer. The number of Voronoi cells, each of which is represented by a centroidal variable coefficient, is set to the prescribed number $P$ of preconditioners. Upon sampling the random variable coefficient, the linear system assembled with a given realization of the coefficient is solved with the preconditioner whose centroidal variable coefficient is the closest to the realization. We consider different ways to define and obtain the centroidal variable coefficients, and we investigate the properties of the induced preconditioning strategies in terms of average number of solver iterations for sequential simulations, and of load balancing for parallel simulations. Another approach, which is based on deterministic grids on the system of stochastic coordinates of the truncated representation of the random variable coefficient, is proposed with a stochastic dimension which increases with the number $P$ of preconditioners. This approach allows to bypass the need for preliminary computations in order to determine the optimal stochastic dimension of the truncated approximation of the random variable coefficient for a given number of preconditioners.

A Bernoulli-barycentric rational matrix collocation method with preconditioning for a class of evolutionary PDEs

Barycentric Rational Collocation Methods for a Class of Nonlinear Parabolic Partial Differential Equations

SOLVING PARTIAL DIFFERENTIAL ALGEBRAIC EQUATIONS BY COLLOCATION AND RADIAL BASIS FUNCTIONS

A Multiple-Scale Higher Order Polynomial Collocation Method for 2D and 3D Elliptic Partial Differential Equations with Variable Coefficients

Numerical Solution of a Class of Nonlinear Partial Differential Equations by Using Barycentric Interpolation Collocation Method

A discrete least squares collocation method for two-dimensional nonlinear time-dependent partial differential equations

Multiquadric quasi-interpolation methods for solving partial differential algebraic equations

High Order Numerical Algorithms Based on Biquadratic Spline Collocation for Two-Dimensional Parabolic Partial Differential Equations

A Barycentric Rational Interpolation Collocation Method for Solving the Helmholtz Equation

Reduced Krylov Basis Methods for Parametric Partial Differential Equations

A reduced conjugate gradient basis method for fractional diffusion

Quadratic spline collocation for one-dimensional linear parabolic partial differential equations

Collocation methods for integral fractional Laplacian and fractional PDEs based on radial basis functions

Preconditioners based on Voronoi quantizers of random variable coefficients for stochastic elliptic partial differential equations

Matrix computational collocation approach based on rational Chebyshev functions for nonlinear differential equations

Collocation methods based on barycentric rational interpolation for Volterra integro‐differential equations with weakly singular kernels

A modified domain decomposition spectral collocation method for parabolic partial differential equations

A Multiscale RBF Collocation Method for the Numerical Solution of Partial Differential Equations

A Bernstein Polynomial Collocation Method for the Solution of Elliptic Boundary Value Problems

A Note on Parallel Preconditioning for the All-at-Once Solution of Riesz Fractional Diffusion Equations

Parallel-in-time solution of hyperbolic PDE systems via characteristic-variable block preconditioning