Resolve "Filter Branch 59"
Closes #65 (closed) . Among the notable changes are the fact that I force the unitary matrices occuring in the eigenvalue decomposition of the operators to be special unitary with det(U) == 1 instead of just a unitary matrix.