A GPU-based Multi-level Algorithm for Boundary Value Problems

by   J. T. Becerra-Sagredo, et al.

A novel and scalable geometric multi-level algorithm is presented for the numerical solution of elliptic partial differential equations, specially designed to run with high occupancy of streaming processors inside Graphics Processing Units(GPUs). The algorithm consists of iterative, superposed operations on a single grid, and it is composed of two simple full-grid routines: a restriction and a coarsened interpolation-relaxation. The restriction is used to collect sources using recursive coarsened averages, and the interpolation-relaxation simultaneously applies coarsened finite-difference operators and interpolations. The routines are scheduled in a saw-like refining cycle. Convergence to machine precision is achieved repeating the full cycle using accumulated residuals and successively collecting the solution. Its total number of operations scale linearly with the number of nodes. It provides an attractive fast solver for Boundary Value Problems (BVPs), specially for simulations running entirely in the GPU. Applications shown in this work include the deformation of two-dimensional grids, the computation of three-dimensional streamlines for a singular trifoil-knot vortex and the calculation of three-dimensional electric potentials in heterogeneous dielectric media.


page 9

page 11


GPU Methodologies for Numerical Partial Differential Equations

In this thesis we develop techniques to efficiently solve numerical Part...

Multi-level Neural Networks for Accurate Solutions of Boundary-Value Problems

The solution to partial differential equations using deep learning appro...

Semi-analytical calculation of the singular and hypersingular integrals for discrete Helmholtz operators in 2D BEM

Approximate solutions to elliptic partial differential equations with kn...

Convergence analysis of inexact two-grid methods: Multigrid cycles

Multigrid is a popular iterative solver for a large class of linear syst...

Convergent Finite Difference Methods for Fully Nonlinear Elliptic Equations in Three Dimensions

We introduce a generalized finite difference method for solving a large ...

Parallel Skeletonization for Integral Equations in Evolving Multiply-Connected Domains

This paper presents a general method for applying hierarchical matrix sk...

Please sign up or login with your details

Forgot password? Click here to reset