unitary matrix for debugging
close #46