Input: A , B , w ( 0 ) = [ z 1 ( 0 ) , z 2 ( 0 ) , α ( 0 ) , β ( 0 ) ] T , k m a x .
for k = 0 , 1 , 2 , ⋯ , until convergence
Choose τ ( k ) .
Solve M ( k ) u ( k ) = B 1 w ( k ) inexactly,
‖ P − 1 M ( k ) u ( k ) − P − 1 B 1 w ( k ) ‖ ≤ τ ( k ) .
end for
Output: u ( k m a x ) .