It follows that there are three levels of iterations, the outer matrix split method, the PRESB and the innermost iterations to solve the basic matrix systems, \(T+W_i\) , \(i=1,2\) . Further, it is most efficient to use a flexible form of outer iteration, see [1]}, [2]}. However, in the numerical tests a direct solver is used for the innermost systems.
