  The main objective of this work consists in analyzing sub-structuring method
for the parallel solution of sparse linear systems with matrices arising from
the discretization of partial differential equations such as finite element,
finite volume and finite difference. With the success encountered by the
general-purpose processing on graphics processing units (GPGPU), we develop an
hybrid multiGPUs and CPUs sub-structuring algorithm. GPU computing, with CUDA,
is used to accelerate the operations performed on each processor. Numerical
experiments have been performed on a set of matrices arising from engineering
problems. We compare C+MPI implementation on classical CPU cluster with
C+MPI+CUDA on a cluster of GPU. The performance comparison shows a speed-up for
the sub-structuring method up to 19 times in double precision by using CUDA.
