Now, form the m × × m matrix representations of the new block and universe-block operators, form a new block by changing basis with the transformation T , for example: H B = T H B − − l T † † S x B = T S x B − − l T † † At this point the iteration is ended and the algorithm goes back to step 1.