Abstract

In this paper we apply the parallel hybrid SPIKE algorithm as a preconditioner for a nonstationary iterative method to solve large sparse linear systems. In order to obtain a good preconditioner, we employ several strategies solving combinatorial problems such as matching, reordering, partitioning, and quadratic knapsack. Our SPIKE implementation combines MPI and OpenMP paradigms in a unique multi-core computer. The computational experiments show the influence of each strategy evaluating the number of iterations and CPU time of the iterative solver in a set of large systems from miscellaneous application areas. The experiments suggest that the SPIKE preconditioner is very advantageous when a suitable set of parameters is chosen. The choice of the number of MPI ranks and OpenMP threads is not an easy task, because the SPIKE algorithm can increase the number of iterations when the number of MPI ranks grows. Moreover, the increase in the number of threads does not ensure a better performance.

Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call