Function random_unitary 

pub fn random_unitary<A>(n: usize) -> ArrayBase<OwnedRepr<A>, Dim<[usize; 2]>>
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.
• This function uses rand::thread_rng. See random_unitary_using for using another RNG.