pub fn random_unitary<A>(n: usize) -> Array2<A>
where
    A: Scalar + Lapack,
{
    let a: Array2<A> = random((n, n));
    let (q, _r) = a.qr_into().unwrap();
    q
}
/// Generate random unitary matrix using QR decomposition
///
/// Be sure that this it **NOT** a uniform distribution. Use it only for test purpose.
Generate random unitary matrix using QR decomposition
Be sure that this it **NOT** a uniform distribution. Use it only for test purpose.
Generate random unitary matrix using QR decomposition Be sure that this it ** NOT ** a uniform distribution . Use it only for test purpose .
pub fn random_unitary<A>(n: usize) -> Array2<A>
where
    A: Scalar + Lapack,
{
    let a: Array2<A> = random((n, n));
    let (q, _r) = a.qr_into().unwrap();
    q
}
pub fn random_unitary < A > ( n : usize ) -> Array2 < A > where A : Scalar + Lapack , { let a : Array2 < A > = random ( ( n , n ) ) ; let ( q , _r ) = a . qr_into ( ) . unwrap ( ) ; q }
Generate random unitary matrix using QR decomposition
Be sure that this it **NOT** a uniform distribution.
Generate random unitary matrix using QR decomposition Be sure that this it ** NOT ** a uniform distribution .
n
usize
n
usize