Hybrid Multi-GPU Solver Based on Schur Complement Method

This paper presents a new hybrid solver based on Schur complement method, in which computations are distributed between multiple CPUs and GPUs. In this solver, the Schur complement is computed either on CPUs for small problem size or on GPUs for large problem sizes. To solve the interface system, we propose a new multi-GPU algorithm that implements conjugate gradient method with explicit preconditioning. Experiments with wrap spring simulation on hybrid multi-core multi-GPU cluster demonstrate efficiency of the proposed method.