Optimisation Note

The above changes assume a general transformation or rotation is required. As the calculation proceeds, the approximate eigenstates should converge to the true eigenstates and so many of the transformations will approach the identity. It may be possible to exploit this anticipated sparsity to optimise the communication pattern.