This calculates the Euler angles of matrix rot in the y-convention. See Elk’s roteuler routine. This will be made redundent when TRIQS uses scipy version 1.4+