Add test for outer and bug fix

Author: termoshtt

src/vector.rs

@@ -64,11 +64,11 @@ pub fn outer<A, S1, S2>(a: &ArrayBase<S1, Ix1>, b: &ArrayBase<S2, Ix1>) -> Array
 {
     let m = a.len();
     let n = b.len();
-    let mut ab = Array::zeros((m, n));
+    let mut ab = Array::zeros((n, m));
     ImplOuter::outer(m,
                      n,
                      a.as_slice_memory_order().unwrap(),
                      b.as_slice_memory_order().unwrap(),
                      ab.as_slice_memory_order_mut().unwrap());
-    ab
+    ab.reversed_axes()
 }

tests/outer.rs

@@ -0,0 +1,20 @@
+include!("header.rs");
+use ndarray_linalg::vector::outer;
+
+#[test]
+fn outer_() {
+    let dist = RealNormal::<f64>::new(0.0, 1.0);
+    let m = 2;
+    let n = 3;
+    let a = Array::random(m, dist);
+    let b = Array::random(n, dist);
+    println!("a = \n{:?}", &a);
+    println!("b = \n{:?}", &b);
+    let ab = outer(&a, &b);
+    println!("ab = \n{:?}", &ab);
+    for i in 0..m {
+        for j in 0..n {
+            ab[(i, j)].assert_close(a[i] * b[j], 1e-7);
+        }
+    }
+}