In a series of publications ([6, 7, 8, 10, 12]) it has been demonstrated that for such problems a very efficient preconditioner, based on the PREconditioned Square matrix Block (PRESB) method, exists and leads to very tight eigenvalue bounds and, hence, few iterations.