Add generate::random(n, m)

Author: termoshtt

src/generate.rs

@@ -19,14 +19,22 @@ pub fn conjugate<A, Si, So>(a: &ArrayBase<Si, Ix2>) -> ArrayBase<So, Ix2>
     a
 }
 
+/// Random matrix
+pub fn random<A, S>(n: usize, m: usize) -> ArrayBase<S, Ix2>
+    where A: RandNormal,
+          S: DataOwned<Elem = A>
+{
+    let mut rng = thread_rng();
+    let v: Vec<A> = (0..n * m).map(|_| A::randn(&mut rng)).collect();
+    ArrayBase::from_shape_vec((n, m), v).unwrap()
+}
+
 /// Random square matrix
 pub fn random_square<A, S>(n: usize) -> ArrayBase<S, Ix2>
     where A: RandNormal,
           S: DataOwned<Elem = A>
 {
-    let mut rng = thread_rng();
-    let v: Vec<A> = (0..n * n).map(|_| A::randn(&mut rng)).collect();
-    ArrayBase::from_shape_vec((n, n), v).unwrap()
+    random(n, n)
 }
 
 /// Random Hermite matrix