Create two_level_decomposition

So far it only works for square matrices up to a 4x4 matrix. Any ideas automating the scheme for all sizes?