Function random_unitary

pub fn random_unitary<A>(n: usize) -> Array2<A>
where
    A: Scalar + Lapack,

Generate random unitary matrix using QR decomposition

Be sure that this it NOT a uniform distribution. Use it only for test purpose.