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

Generate random unitary matrix using QR decomposition

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