Added simple verbose output to track convergence
Test that projection operators are idempotent
np.angle ranges from -pi to pi; easier than I thought
Move the W projection to a different class
New W projection class for every component of any given sample to have the same phase

But phase can vary between samples i.e. over time