Test that running on output doesn't do anything
For whatever input, the output of general
, expanded
, and general_symbolic
is a block-diagonal matrix.
This means that using the output of these algorithms as an input, should find U=identity
and not do anything.
We could implement this as a consistency test.