“…Since the coefficient matrix is often large and sparse, iterative methods become more attractive than direct methods for the saddle point problem (1). Many efficient methods have been proposed, such as the Uzawa-type schemes ( [10,15,16,19,22]), the Krylov subspace methods ( [25,33,34]), the SOR-like methods ( [9,20,26,30,35]), the Hermitian and skew-Hermitian splitting (HSS) iteration methods ( [2,[4][5][6][7][8]29]) and some other iterative methods [15,28,42], we mention just a few.…”