  This research introduce our work on developing Krylov subspace and AMG
solvers on NVIDIA GPUs. As SpMV is a crucial part for these iterative methods,
SpMV algorithms for single GPU and multiple GPUs are implemented. A HEC matrix
format and a communication mechanism are established. And also, a set of
specific algorithms for solving preconditioned systems in parallel environments
are designed, including ILU(k), RAS and parallel triangular solvers. Based on
these work, several Krylov solvers and AMG solvers are developed. According to
numerical experiments, favorable acceleration performance is acquired from our
Krylov solver and AMG solver under various parameter conditions.
