GPU-acceleration of waveform relaxation methods for large differential systems

It is the purpose of this paper to provide an acceleration of waveform relaxation (WR) methods for the numerical solution of large systems of ordinary differential equations. The introduced technique is based on the employ of graphics processing units (GPUs) in order to speed-up the numerical integration process. A CUDA solver based on WR-Picard, WR-Jacobi and red-black WR-Gauss-Seidel iterations is presented and some numerical experiments realized on a multi-GPU machine are provided