Abstract:This article proposes an efficient numerical method for solving nonlinear partial differential equations (PDEs) based on sparse Gaussian processes (SGPs). Gaussian processes (GPs) have been extensively studied for solving PDEs by formulating the problem of finding a reproducing kernel Hilbert space (RKHS) to approximate a PDE solution. The approximated solution lies in the span of base functions generated by evaluating derivatives of different orders of kernels at sample points. However, the RKHS specified by GPs can result in an expensive computational burden due to the cubic computation order of the matrix inverse. Therefore, we conjecture that a solution exists on a ``condensed" subspace that can achieve similar approximation performance, and we propose a SGP-based method to reformulate the optimization problem in the ``condensed" subspace. This significantly reduces the computation burden while retaining desirable accuracy. The paper rigorously formulates this problem and provides error analysis and numerical experiments to demonstrate the effectiveness of this method. The numerical experiments show that the SGP method uses fewer than half the uniform samples as inducing points and achieves comparable accuracy to the GP method using the same number of uniform samples, resulting in a significant reduction in computational cost. Our contributions include formulating the nonlinear PDE problem as an optimization problem on a ``condensed" subspace of RKHS using SGP, as well as providing an existence proof and rigorous error analysis. Furthermore, our method can be viewed as an extension of the GP method to account for general positive semi-definite kernels.

Localized Method of Approximate Particular Solutions with Polynomial Basis Functions

Localized meshless methods based on polynomial basis functions for solving axisymmetric equations

A Novel Method of Polynomial Approximation for Parametric Problems in Power Systems

The Method of Polynomial Particular Solutions for Solving Ordinary Differential Equations

The Method of Polynomial Particular Solutions for Solving Nonlinear Poisson-Type Equations

A New Spectral Difference Method Using Hierarchical Polynomial Bases for Hyperbolic Conservation Laws

Localized subspace iteration methods for elliptic multiscale problems

Localized Method of Approximate Particular Solutions for Solving Unsteady Navier–Stokes Problem

A Class of Linearization-Based Collocation Methods for Initial Value and Boundary Value Engineering Problems.

Sparse Gaussian processes for solving nonlinear PDEs

Localized method of fundamental solutions for 2D harmonic elastic wave problems

A well conditioned Method of Fundamental Solutions

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

Localized collocation schemes and their applications

Regular polynomial interpolation and approximation of global solutions of linear partial differential equations

Localized Chebyshev and MLS Collocation Methods for Solving 2D Steady State Nonlocal Diffusion and Peridynamic Equations

Solving a linear fractional equation with nonlocal boundary conditions based on multiscale orthonormal bases method in the reproducing kernel space

Accurately Solving the Poisson Equation by Combining Multiscale Radial Basis Functions and Gaussian Quadrature

Solve Partial Differential Equations by Two or More Radial Basis Functions

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

Numerical Analysis for Sturm–Liouville Problems with Nonlocal Generalized Boundary Conditions