Solution of large sparse linear systems is the most lime consuming part in many power system simulations. Direct solvers based on LU factorization although robust are known to have limited satiability on parallel platforms. Thus. Krylov subspace based iterative methods (i.e. Conjugate Gradient method Generalized Minimal Residuals (GMRES) method) can be used as alternatives. To achieve competitive performance and robustness however the Krylov subspace methods need a suitable preconditioner. In this ...